bigdecimal-math_r 0.2.1 → 0.2.2

data/ext/bigdecimal/math_r/internal/algo/bigmath/sinh/ser.h

@@ -1,11 +1,12 @@
 VALUE
 sinh_ser(VALUE x, VALUE prec)
 {
-	VALUE y = Qundef;
+	VALUE y = Qundef, n;
 
 	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
 
-	sinhcosh_ser(x, prec, &y, NULL);
+	sinhcosh_ser(x, n, &y, NULL);
 
-	return y;
+	return rb_num_round(y, prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/sinhcosh/ser.h

@@ -10,9 +10,9 @@ sinhcosh_inline(VALUE x, VALUE prec, VALUE *sinh, VALUE *cosh)
 	w = BIG_ONE;
 	a = BIG_ONE;
 	y = BIG_ZERO;
-	
 	if (sinh != NULL)  *sinh = BIG_ZERO;
 	if (cosh != NULL)  *cosh = BIG_ZERO;
+
 	while (rb_numdiff_condition_p(y, a, n, &m))
 	{
 		a = rb_funcall(w, div, 2, f, m);
@@ -27,10 +27,12 @@ sinhcosh_inline(VALUE x, VALUE prec, VALUE *sinh, VALUE *cosh)
 		w = rb_funcall1(w, '*', x);
 		f = rb_funcall1(f, '*', LONG2NUM(++i));
 	}
+
 	RB_GC_GUARD(f);
 	RB_GC_GUARD(w);
 	RB_GC_GUARD(a);
 	RB_GC_GUARD(y);
+
 	if (sinh != NULL)  RB_GC_GUARD(*sinh);
 	if (cosh != NULL)  RB_GC_GUARD(*cosh);
 }
@@ -38,8 +40,13 @@ sinhcosh_inline(VALUE x, VALUE prec, VALUE *sinh, VALUE *cosh)
 void
 sinhcosh_ser(VALUE x, VALUE prec, VALUE *sinh, VALUE *cosh)
 {
+	VALUE n;
 	int sign;
-	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
+
+	rb_check_precise(prec);
+
+	n = rb_numdiff_make_n(prec);
+	x = rb_num_canonicalize(x, n, ARG_REAL, ARG_RAWVALUE);
 	
 	if (sinh == NULL && cosh == NULL)  return;
 	
@@ -63,7 +70,7 @@ sinhcosh_ser(VALUE x, VALUE prec, VALUE *sinh, VALUE *cosh)
 	}
 	else
 	{
-		sinhcosh_inline(x, prec, sinh, cosh);
+		sinhcosh_inline(x, n, sinh, cosh);
 		if (sinh != NULL)  *sinh = rb_num_round(*sinh, prec);
 		if (cosh != NULL)  *cosh = rb_num_round(*cosh, prec);
 	}

data/ext/bigdecimal/math_r/internal/algo/bigmath/sqrt/builtin.h

@@ -5,8 +5,8 @@ sqrt_builtin(VALUE x, VALUE prec)
 
 	rb_check_precise(prec);
 
-	if (rb_num_negative_p(x) || rb_num_notequal_p(x, x))
+	if (rb_num_negative_p(x) || rb_num_nan_p(x))
 		return BIG_NAN;
 
-	return rb_funcall(x, sqrt, 1, prec);
+	return rb_num_round(rb_funcall(x, sqrt, 1, prec), prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/sqrt/edom.h

@@ -5,5 +5,6 @@ sqrt_edom(VALUE x, VALUE prec)
 	RUBY_ASSERT(rb_num_negative_p(x));
 	x = rb_num_uminus(x);
 	imag = rb_bigmath_sqrt(x, prec);
+
 	return rb_num_round(rb_Complex(BIG_ZERO, imag), prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/tan/eulerf.h

@@ -1,11 +1,13 @@
 VALUE
 tan_eulerf(VALUE z, VALUE prec)
 {
-	VALUE w;
+	VALUE w, n;
 
 	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
 
-	w = trig_f_euler(mf_tan, rb_ImaginaryZ(z, SIGN_PLUS), prec);
+	w = trig_f_euler(mf_tan, rb_ImaginaryZ(z, SIGN_PLUS), n);
+	w = rb_num_real(w);
 
-	return rb_num_real(w);
+	return rb_num_round(w, prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/tan/ser.h

@@ -2,17 +2,12 @@ VALUE
 tan_ser(VALUE x, VALUE prec)
 {
 	const ID div = rb_intern("div");
-	VALUE sin, cos;
+	VALUE sin, cos, n;
 
 	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
 
-	if (!rb_num_finite_p(x))
-		return BIG_NAN;
-	rb_check_precise(prec);
-	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
-	rb_bigmath_sincos(x, prec, &sin, &cos);
-	sin = rb_num_round(sin, prec);
-	cos = rb_num_round(cos, prec);
+	rb_bigmath_sincos(x, n, &sin, &cos);
 
 	return rb_funcall(sin, div, 2, cos, prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/tanh/eulerf.h

@@ -1,11 +1,13 @@
 VALUE
 tanh_eulerf(VALUE x, VALUE prec)
 {
-	VALUE y;
+	VALUE y = Qundef, n;
 
 	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
 
-	y = rb_ImaginaryZ(trig_f_euler(mf_tan, x, prec), SIGN_PLUS);
+	y = rb_ImaginaryZ(trig_f_euler(mf_tan, x, n), SIGN_PLUS);
+	y = rb_num_real(y);
 
-	return rb_num_real(y);
+	return rb_num_round(y, prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/tanh/formula.h

@@ -4,11 +4,15 @@ tanh_formula(VALUE x, VALUE prec)
 	const ID pred = rb_intern("pred");
 	const ID succ = rb_intern("succ");
 	const ID div = rb_intern("div");
-	VALUE n = rb_numdiff_make_n(prec);
-	VALUE two = rb_BigDecimal1(INT2FIX(2));
-        VALUE t = rb_bigmath_exp(rb_funcall1(x, '*', two), n), y = Qundef;
+	VALUE n, two, t, y = Qundef;
 	int sign;
 
+	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
+
+	two = rb_BigDecimal1(INT2FIX(2));
+        t = rb_bigmath_exp(rb_funcall1(x, '*', two), n);
+
 	if (TYPE(t) == T_FIXNUM || TYPE(t) == T_BIGNUM)
 	{
 		y = rb_Rational(rb_funcall(t, pred, 0), rb_funcall(t, succ, 0));

data/ext/bigdecimal/math_r/internal/algo/bigmath/tanh/ser.h

@@ -1,13 +1,13 @@
 VALUE
 tanh_ser(VALUE x, VALUE prec)
 {
-	VALUE y = Qundef, sinh = Qundef, cosh = Qundef;
+	VALUE y = Qundef, sinh = Qundef, cosh = Qundef, n;
 	const ID div = rb_intern("div");
 
 	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
 
-	sinhcosh_ser(x, prec, &sinh, &cosh);
-
+	sinhcosh_ser(x, n, &sinh, &cosh);
 	y = rb_funcall(sinh, div, 2, cosh, prec);
 
 	return y;

data/ext/bigdecimal/math_r/internal/algo/bigmath/trig/f_euler.h

@@ -263,18 +263,24 @@ f_euler_esc(const ID func, VALUE z, VALUE prec)
 VALUE
 trig_f_euler(const ID func, VALUE theta, VALUE prec)
 {
-	VALUE h = Qundef, sc = Qundef;
+	VALUE h = Qundef, sc = Qundef, n;
+
+	rb_check_precise(prec);
+	n = rb_numdiff_make_n(prec);
+
 	if (rb_num_nan_p(theta))
 		return rb_Complex(BIG_NAN, BIG_NAN);
 	if (rb_num_zero_p(rb_num_imag(theta)))
 	{
-		h = f_euler_h(func, theta, prec);
+		h = f_euler_h(func, theta, n);
+		h = rb_num_round(h, prec);
 		if (rb_num_nan_p(h))
 			return h;
 	}
 	if (rb_num_zero_p(rb_num_real(theta)))
 	{
-		sc = f_euler_sc(func, theta, prec);
+		sc = f_euler_sc(func, theta, n);
+		sc = rb_num_round(sc, prec);
 		if (rb_num_nan_p(sc))
 			return sc;
 	}
@@ -284,7 +290,7 @@ trig_f_euler(const ID func, VALUE theta, VALUE prec)
 		return sc;
 	else
 	{
-		VALUE y = f_euler_esc(func, theta, prec);
-		return y;
+		VALUE y = f_euler_esc(func, theta, n);
+		return rb_num_round(y, prec);
 	}
 }

data/ext/bigdecimal/math_r/log.c

@@ -201,7 +201,7 @@ __impl_clog_formula(VALUE unused_obj, VALUE z, VALUE prec)
 static VALUE
 __impl_clog2_formula(VALUE unused_obj, VALUE z, VALUE prec)
 {
-	return clog_branch(z, prec, clog2_formula);
+	return clog2_branch(z, prec, clog2_formula);
 }
 
 /**
@@ -220,7 +220,7 @@ __impl_clog2_formula(VALUE unused_obj, VALUE z, VALUE prec)
 static VALUE
 __impl_clog10_formula(VALUE unused_obj, VALUE z, VALUE prec)
 {
-	return clog_branch(z, prec, clog10_formula);
+	return clog10_branch(z, prec, clog10_formula);
 }
 
 

data/ext/bigdecimal/math_r/math_r/bigdecimal.h

@@ -1,8 +1,12 @@
 #ifndef INTERNAL_API_BIGDECIMAL
 #define INTERNAL_API_BIGDECIMAL
 
-/* API for BigDecimal */
+/**
+ *  Query the Parameter `prec` that is in valid range $p \lt 0$.
+ *  if invalid range, error occurrd.
+ */
 void rb_check_precise(VALUE prec);
+#define CheckPrecise(prec) rb_check_precise(prec) 
 
 VALUE rb_big_zero(void);
 #define BIG_ZERO rb_big_zero()

data/ext/bigdecimal/math_r/math_r/bigmath.h

@@ -5,6 +5,7 @@
 #define DEF_PREC 1000
 #define BIGMATH_PI     "$M_PI"
 #define BIGMATH_E      "$M_E"
+#define BIGMATH_EulerGamma "$M_EulerGamma"
 #define BIGMATH_LOG2   "$M_LOG2"
 #define BIGMATH_LOG10  "$M_LOG10"
 #define BIGMATH_LOG_PI "$M_LOG_PI"
@@ -21,6 +22,7 @@ VALUE rb_bigmath_const(const char *gv, VALUE (*func)(VALUE prec), VALUE prec);
 VALUE rb_bigmath_const_pi(VALUE prec);
 VALUE rb_bigmath_const_log_pi(VALUE prec);
 VALUE rb_bigmath_const_e(VALUE prec);
+VALUE rb_bigmath_const_euler_gamma(VALUE prec);
 VALUE rb_bigmath_const_log2(VALUE prec);
 VALUE rb_bigmath_const_log10(VALUE prec);
 VALUE rb_bigmath_const_sqrt2(VALUE prec);

data/ext/bigdecimal/math_r/math_r/globals.h

@@ -23,7 +23,7 @@ RUBY_EXT_EXTERN VALUE rb_mComplexPlane;
 RUBY_EXT_EXTERN VALUE rb_mEulerF;
 RUBY_EXT_EXTERN VALUE rb_mLogRep;
 RUBY_EXT_EXTERN VALUE rb_mPower;
-RUBY_EXT_EXTERN VALUE rb_mGammaR;
+RUBY_EXT_EXTERN VALUE rb_mGammaBeta;
 RUBY_EXT_EXTERN VALUE rb_mErf;
 RUBY_EXT_EXTERN VALUE rb_mSolver;
 

data/ext/bigdecimal/math_r/math_r.c

@@ -21,7 +21,7 @@ void InitVM_EulerF(void);
 void InitVM_LogRep(void);
 void InitVM_Power(void);
 void InitVM_Erf(void);
-void InitVM_GammaR(void);
+void InitVM_GammaBeta(void);
 void InitVM_Solver(void);
 
 /**
@@ -36,8 +36,10 @@ void InitVM_Solver(void);
 void
 Init_math_r(void)
 {
+
 	rb_gv_set(BIGMATH_PI, Qnil);
 	rb_gv_set(BIGMATH_E, Qnil);
+	rb_gv_set(BIGMATH_EulerGamma, Qnil);
 	rb_gv_set(BIGMATH_LOG2, Qnil);
 	rb_gv_set(BIGMATH_LOG10, Qnil);
 	rb_gv_set(BIGMATH_LOG_PI, Qnil);
@@ -130,8 +132,8 @@ Init_math_r(void)
 	rb_mErf = rb_define_module_under(rb_mBigMathR, "Erf");
 	InitVM(Erf);
 
-	rb_mGammaR = rb_define_module_under(rb_mBigMathR, "GammaR");
-	InitVM(GammaR);
+	rb_mGammaBeta = rb_define_module_under(rb_mBigMathR, "GammaBeta");
+	InitVM(GammaBeta);
 
 	rb_mSolver = rb_define_module_under(rb_mBigMathR, "Solver");
 	InitVM(Solver);

data/lib/bigdecimal/math_r/const/E_euler.rb

@@ -0,0 +1,46 @@
+# frozen_string_literal: true
+
+module BigMathR
+  module Const
+    module E
+      module_function
+      # Implement by continued fraction.
+      # 
+      # @example
+      #  BigMathR::Const::E.euler(20)
+      #  #=> 0.271828182845904523536e1
+      # @param prec [Integer] Arbitrary precision
+      # @return [BigDecimal] Constant e.
+      # @raise [TypeError] not an Integer
+      # @raise [RangeError] Zero or negative precision
+      def euler(prec)
+        raise TypeError, "precision must be in Integer" unless prec.class == Integer
+        raise RangeError, "Zero or negative precision" if prec <= 0
+        n = BigDecimal.double_fig + prec
+
+        zero = BigDecimal(0)
+        one = BigDecimal(1)
+        two = BigDecimal(2)
+
+        c = one
+        b = one
+        p1 = one;  q1 = zero;  p2 = two;  q2 = one;
+        i = one
+        loop do
+          prev = p2
+          t = p1 * c + p2 * b;  p1 = p2;  p2 = t;
+          t = q1 * c + q2 * b;  q1 = q2;  q2 = t;
+          if q2.nonzero?
+            p1 = p1.div(q2, n);  q1 = q1.div(q2, n);  p2 = p2.div(q2, n);  q2 = one;
+          end
+          i += one
+          c = one
+          b = (i + one) % 3 == 0 ? ((i.to_i.succ) / 3) * two : one
+          break unless q2.nonzero? && prev != p2
+        end
+        p2.round(prec)
+      end
+    end
+  end
+end
+

data/lib/bigdecimal/math_r/version.rb

@@ -2,6 +2,6 @@
 
 module BigMathR
   # Module version.
-  VERSION = "0.2.1"
+  VERSION = "0.2.2"
 end
 

data/lib/bigdecimal/math_r.rb

@@ -4,6 +4,7 @@ require 'bigdecimal/math'
 require_relative './math_r.so'
 
 require_relative 'math_r/version'
+require_relative 'math_r/const/E_euler'
 require_relative 'math_r/const/PI_chudnovsky'
 require_relative 'math_r/const/PI_euler'
 require_relative 'math_r/const/PI_ramanujan1'

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: bigdecimal-math_r
 version: !ruby/object:Gem::Version
-  version: 0.2.1
+  version: 0.2.2
 platform: ruby
 authors:
 - tribusonz-2
@@ -27,16 +27,16 @@ dependencies:
   name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
-        version: '0'
+        version: '13.0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
-        version: '0'
+        version: '13.0'
 - !ruby/object:Gem::Dependency
   name: rake-compiler
   requirement: !ruby/object:Gem::Requirement
@@ -59,7 +59,6 @@ extensions:
 - ext/bigdecimal/math_r/extconf.rb
 extra_rdoc_files: []
 files:
-- ".ruby-version"
 - CODE_OF_CONDUCT.md
 - LICENSE
 - README.md
@@ -87,6 +86,7 @@ files:
 - ext/bigdecimal/math_r/api/bigmath/clog.h
 - ext/bigdecimal/math_r/api/bigmath/const.h
 - ext/bigdecimal/math_r/api/bigmath/const_e.h
+- ext/bigdecimal/math_r/api/bigmath/const_euler_gamma.h
 - ext/bigdecimal/math_r/api/bigmath/const_log10.h
 - ext/bigdecimal/math_r/api/bigmath/const_log2.h
 - ext/bigdecimal/math_r/api/bigmath/const_log_pi.h
@@ -140,9 +140,10 @@ files:
 - ext/bigdecimal/math_r/euler_f.c
 - ext/bigdecimal/math_r/exp.c
 - ext/bigdecimal/math_r/extconf.rb
-- ext/bigdecimal/math_r/gamma_r.c
+- ext/bigdecimal/math_r/gamma_beta.c
 - ext/bigdecimal/math_r/hyperb.c
-- ext/bigdecimal/math_r/internal/algo/bigmath/E/ser.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/E/napier.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/EulerGamma/borwein_bailey.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/LOG10/ser.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/LOG2/BBP2002.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/LOG_PI/builtin.h
@@ -213,6 +214,7 @@ files:
 - ext/bigdecimal/math_r/internal/algo/bigmath/catan/logrep.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/catanh/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/catanh/logrep.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/ccbrt/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ccbrt/formula.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ccos/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ccos/eulerf.h
@@ -230,10 +232,14 @@ files:
 - ext/bigdecimal/math_r/internal/algo/bigmath/cexp/eulerf.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/cexp2/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/cexp2/eulerf.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/cexpm1/branch.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/cexpm1/identity.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/clog/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/clog/formula.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/clog10/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/clog10/formula.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/clog1p/branch.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/clog2/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/clog2/formula.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/cos/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/cos/eulerf.h
@@ -265,11 +271,13 @@ files:
 - ext/bigdecimal/math_r/internal/algo/bigmath/csin/eulerf.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/csinh/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/csinh/eulerf.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/csqrt/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/csqrt/formula.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ctan/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ctan/eulerf.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ctanh/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/ctanh/eulerf.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/cuberoot/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/cuberoot/newton.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/erf/algo911.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/erf/branch.h
@@ -285,6 +293,8 @@ files:
 - ext/bigdecimal/math_r/internal/algo/bigmath/exp/edf.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/exp2/branch.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/exp2/edf.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/expm1/branch.h
+- ext/bigdecimal/math_r/internal/algo/bigmath/expm1/identity.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/expxt/edf.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/f_euler/formula.h
 - ext/bigdecimal/math_r/internal/algo/bigmath/hypot/branch.h
@@ -356,7 +366,7 @@ files:
 - ext/bigdecimal/math_r/solver.c
 - ext/bigdecimal/math_r/trig.c
 - lib/bigdecimal/math_r.rb
-- lib/bigdecimal/math_r.so
+- lib/bigdecimal/math_r/const/E_euler.rb
 - lib/bigdecimal/math_r/const/EulerGamma_engel.rb
 - lib/bigdecimal/math_r/const/LOG2_bbp2007.rb
 - lib/bigdecimal/math_r/const/PI_chudnovsky.rb
@@ -387,7 +397,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.6.7
+rubygems_version: 3.7.2
 specification_version: 4
 summary: Refined Math for BigDecimal
 test_files: []

data/.ruby-version

@@ -1 +0,0 @@
-3.4.3