bigdecimal-math_r 0.2.5 → 0.2.7

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: d14324805ebc5939b65cb1da2f7c5ffb351d7093883f6329633ad6813adfb42d
-  data.tar.gz: e435ef471f9880fafb5cd05866b11869d6c627ace4fa2f41f1f4059a7641eacf
+  metadata.gz: 9afea6caf50f7f5fb373b89cbb3014c57cfd14377db98462b681dbfe27e288d6
+  data.tar.gz: 27f3ebfbdf90fbce30895a10be0b34318be3f0e2eae23b66aaf14851e7fc6973
 SHA512:
-  metadata.gz: 70d16c4a4e4fe6d4e2f647e3cbf782477327c8e8cda71fd50b08676212a54cddb1af3fd093e8cc75968bae124a519fe602cd87169334f532244eacc592c2d85a
-  data.tar.gz: 770a02f513ae60c965bde0a34cba8a6c74a68a1b5301006d9652da6897ccaf58fecafd938a8960135d439588c307d36fa32bfcca617a36cc99d37caeff99da43
+  metadata.gz: 1b0492824d423b2673ff0c6faa70bb443ea92086b88b63a6791e02412e52e1cb6f6075815025724ad029d13da86ca85509bbeb8596694ae5c6c46093cc6f309f
+  data.tar.gz: 3916b4ad816aa0736561e22f9d56c522d07d97b5ebc089021254fd7dd40c8edc71a55f9a19784cf78be3d719ad8f3512d2546d79752da3100fb53287a5d1f059

data/changelog.md

@@ -7,6 +7,39 @@ The format is based on [Keep a Changelog](http://keepachangelog.com/en/1.0.0/) a
 
 ## [Unreleased]
 
+## [0.2.7] - 2025-10-XX
+
+### Fix
+- rcm2_edf(): Ignore Increase/Declease for exponent
+
+### Added
+- BigMathR#EulerGamma
+- BigMathR::Const::EulerGamma#borwein_bailey
+
+## [0.2.6] - 2025-10-XX
+
+### Fix
+- A some function's branch-cut set open-ended range
+- BigMathR#asech invalid soltion
+- `solver_triginv()` and `solver_hyperbinv` an error occurs when the value is a Complex with an imaginary part of 0.
+- rcm2_edf(): type conversion error
+
+## [0.2.5] - 2025-10-XX
+
+### Fix
+- Change require `./math_r.so` to `math_r`
+
+## [0.2.4] - 2025-10-XX
+
+### Fix
+- Delete `BigMathR.EulerGamma`
+- Delete `BigMathR::Const::EulerGamma#borwein_bailey`
+
+## [0.2.3] - 2025-10-XX
+
+### Fix
+- unload `math_r.so`.
+
 ## [0.2.2] - 2025-10-XX
 
 ### Changed

data/ext/bigdecimal/math_r/const.c

@@ -54,7 +54,6 @@ __impl_const_e(VALUE unused_obj, VALUE prec)
 	return rb_bigmath_const_e(prec);
 }
 
-#if 0
 /**
  * The Euler-Mascheroni constant gamma.
  * @example
@@ -69,7 +68,6 @@ __impl_const_euler_gamma(VALUE unused_obj, VALUE prec)
 {
 	return rb_bigmath_const_euler_gamma(prec);
 }
-#endif
 
 /**
  * The natural logarithm of 2.
@@ -165,7 +163,6 @@ __impl_const_e_napier(VALUE unused_obj, VALUE prec)
 	return E_napier(prec);
 }
 
-#if 0
 /**
  * Implement by Borwein-Bailey's formula (Brent-McMillan type formula)
  * 
@@ -182,7 +179,6 @@ __impl_const_euler_gamma_borwein_bailey(VALUE unused_obj, VALUE prec)
 {
 	return EulerGamma_borwein_bailey(prec);
 }
-#endif
 
 /**
  * Implement by BBP's formula (Borwein and Bailey 2002)
@@ -208,9 +204,7 @@ InitVM_Const(void)
 {
 	rb_define_module_function(rb_mBigMathR, "PI", __impl_const_pi, 1);
 	rb_define_module_function(rb_mBigMathR, "E", __impl_const_e, 1);
-#if 0
 	rb_define_module_function(rb_mBigMathR, "EulerGamma", __impl_const_euler_gamma, 1);
-#endif
 	rb_define_module_function(rb_mBigMathR, "LOG2", __impl_const_log2, 1);
 	rb_define_module_function(rb_mBigMathR, "LOG_PI", __impl_const_log_pi, 1);
 	rb_define_module_function(rb_mBigMathR, "LOG10", __impl_const_log10, 1);
@@ -220,12 +214,10 @@ InitVM_Const(void)
 	rb_define_module_function(rb_mConstPI, "machin", __impl_const_pi_machin, 1);
 	rb_define_module_function(rb_mConstE, "napier", __impl_const_e_napier, 1);
 	rb_define_module_function(rb_mConstE, "euler_number", __impl_const_e_napier, 1);
-#if 0
 	rb_define_module_function(rb_mConstEulerGamma, "borwein_bailey",
 		__impl_const_euler_gamma_borwein_bailey, 1);
 	rb_define_module_function(rb_mConstEulerGamma, "brent_mcmillan",
 		__impl_const_euler_gamma_borwein_bailey, 1);
-#endif
 	rb_define_module_function(rb_mConstLOG2, "bbp2002", __impl_const_log2_bbp2002, 1);
 }
 

data/ext/bigdecimal/math_r/internal/algo/bigmath/acos/branch.h

@@ -4,7 +4,7 @@ acos_branch(VALUE x, VALUE prec, bigmath_func1 acos_cb)
 	const ID leq = rb_intern("<=");
 	VALUE y = Qundef;
 	VALUE p001 = rb_BigDecimal1(rb_str_new_cstr("0.001"));
-	VALUE domain = rb_range_new(INT2FIX(0), INT2FIX(1), true);
+	VALUE domain = rb_range_new(INT2FIX(0), INT2FIX(1), false);
 
 	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
 

data/ext/bigdecimal/math_r/internal/algo/bigmath/acsc/branch.h

@@ -3,7 +3,7 @@ acsc_branch(VALUE x, VALUE prec, bigmath_func1 acsc_cb)
 {
 	VALUE y = Qundef;
 	VALUE p_domain = rb_range_new(INT2FIX(1), DBL2NUM(HUGE_VAL), false);
-	VALUE m_domain = rb_range_new(DBL2NUM(-HUGE_VAL), INT2FIX(-1), true);
+	VALUE m_domain = rb_range_new(DBL2NUM(-HUGE_VAL), INT2FIX(-1), false);
 
 	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
 

data/ext/bigdecimal/math_r/internal/algo/bigmath/asech/branch.h

@@ -2,7 +2,7 @@ VALUE
 asech_branch(VALUE x, VALUE prec, bigmath_func1 asech_cb)
 {
 	VALUE y;
-	VALUE domain = rb_range_new(INT2FIX(0), INT2FIX(1), true);
+	VALUE domain = rb_range_new(INT2FIX(0), INT2FIX(1), false);
 
 	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
 

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

@@ -11,13 +11,12 @@ asech_edom(VALUE x, VALUE prec)
 	sqr = rb_funcall1(
 		rb_funcall1(BIG_ONE, '-', absx), '/', 
 		rb_funcall1(absx, '+', BIG_ONE));
-	w = rb_Complex(BIG_ZERO, rb_bigmath_sqrt(sqr, n));
+	w = rb_Complex(BIG_ZERO, rb_bigmath_sqrt(rb_num_uminus(sqr), n));
 	w = rb_funcall1(w, '*', rb_funcall1(recip_x, '+', BIG_ONE));
 	w = rb_funcall1(w, '+', recip_x);
 	y = rb_bigmath_carg(w, n);
 	if (rb_num_negative_p(x))
 		y =rb_funcall1(rb_bigmath_const_pi(n), '-', y);
 	y = rb_Complex(BIG_ZERO, y);
-
 	return rb_num_round(y, prec);
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/asech/logrep.h

@@ -12,6 +12,5 @@ asech_logrep(VALUE x, VALUE prec)
 	y = rb_funcall1(y, '+', BIG_ONE);
 	y = rb_funcall1(y, '/', x);
 	y = rb_bigmath_log(y, prec);
-
 	return y;
 }

data/ext/bigdecimal/math_r/internal/algo/bigmath/asin/branch.h

@@ -4,7 +4,7 @@ asin_branch(VALUE x, VALUE prec, bigmath_func1 asin_cb)
 	const ID leq = rb_intern("<=");
 	VALUE y = Qundef;
 	VALUE p001 = rb_BigDecimal1(rb_str_new_cstr("0.001"));
-	VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), true);
+	VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), false);
 
 	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
 

data/ext/bigdecimal/math_r/internal/algo/bigmath/atanh/branch.h

@@ -4,7 +4,7 @@ atanh_branch(VALUE x, VALUE prec, bigmath_func1 atanh_cb)
 	const ID leq = rb_intern("<=");
 	VALUE y = Qundef;
 	VALUE p0999 = rb_BigDecimal1(rb_str_new_cstr("0.999"));
-	VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), true);
+	VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), false);
 
 	x = rb_num_canonicalize(x, prec, ARG_REAL, ARG_RAWVALUE);
 

data/ext/bigdecimal/math_r/internal/algo/bigmath/rcm2/edf.h

@@ -19,7 +19,7 @@ rcm2_edf(VALUE x, VALUE *reso)
 			while (RTEST(rb_num_coerce_cmp(INT2FIX(2), fra, rb_intern("<="))))
 			{
 				fra = rb_funcall1(fra, '/', rat_two);
-				TYPE(resov) == T_FIXNUM ? 
+				resov = TYPE(resov) == T_FIXNUM ? 
 					LONG2FIX(FIX2LONG(resov) + 1) : 
 					rb_funcall(resov, succ, 0);
 			}
@@ -29,14 +29,14 @@ rcm2_edf(VALUE x, VALUE *reso)
 			while (RTEST(rb_num_coerce_cmp(INT2FIX(1), fra, '>')))
 			{
 				fra = rb_funcall1(fra, '*', rat_two);
-				TYPE(resov) == T_FIXNUM ? 
+				resov = TYPE(resov) == T_FIXNUM ? 
 					LONG2FIX(FIX2LONG(resov) - 1) : 
 					rb_funcall(resov, pred, 0);
 			}
 		}
 		if (has_sign)
 			fra = rb_num_uminus(fra);
-		*reso = LONG2FIX(resov);
+		*reso = resov;
 	}
 	return fra;
 }

data/ext/bigdecimal/math_r/solver.c

@@ -864,6 +864,8 @@ solver_triginv(ID func, VALUE z, VALUE prec)
 	{
 		if (rb_num_real_p(z) || rb_num_zero_p(rb_num_imag(z)))
 		{
+			if (RB_TYPE_P(z, T_COMPLEX))
+				z = rb_num_real(z);
 			if (func == mf_asin)
 			{
 				VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), 0);
@@ -1113,6 +1115,8 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 	{
 		if (rb_num_real_p(z) || rb_num_zero_p(rb_num_imag(z)))
 		{
+			if (RB_TYPE_P(z, T_COMPLEX))
+				z = rb_num_real(z);
 			if (func == mf_asinh)
 			{
 				w = asinh_branch(z, prec, cb_asinh);
@@ -1152,7 +1156,7 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 					w = asech_branch(z, prec, cb_asech);
 				else if (rb_num_domain_p(m_domain, z))
 				{
-					w = asech_branch(z, prec, cb_asech);
+					w = asech_branch(rb_num_uminus(z), prec, cb_asech);
 					w = rb_Complex(w, rb_bigmath_const_pi(prec));
 				}
 				else

data/lib/bigdecimal/math_r/version.rb

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

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: bigdecimal-math_r
 version: !ruby/object:Gem::Version
-  version: 0.2.5
+  version: 0.2.7
 platform: ruby
 authors:
 - tribusonz-2
@@ -366,7 +366,6 @@ 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