bigdecimal-math_r 0.2.0 → 0.2.2

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);
 }
 
 
@@ -243,9 +243,9 @@ InitVM_Log(void)
 	rb_define_module_function(rb_mLog, "log_builtin", __impl_log_builtin, 2);
 	rb_define_module_function(rb_mLog, "log_ser_okumura", __impl_log_ser_okumura, 2);
 
-	rb_define_module_function(rb_mEDF, "log_edf", __impl_log_edf, 2);
-	rb_define_module_function(rb_mEDF, "log2_edf", __impl_log2_edf, 2);
-	rb_define_module_function(rb_mEDF, "log10_edf", __impl_log10_edf, 2);
+	rb_define_module_function(rb_mLog, "log_edf", __impl_log_edf, 2);
+	rb_define_module_function(rb_mLog, "log2_edf", __impl_log2_edf, 2);
+	rb_define_module_function(rb_mLog, "log10_edf", __impl_log10_edf, 2);
 
 	rb_define_module_function(rb_mLog, "log1p_ser_mercator", __impl_log1p_ser_mercator, 2);
 	rb_define_module_function(rb_mLog, "log1p_p_adic", __impl_log1p_p_adic, 2);

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/numeric.h

@@ -26,6 +26,7 @@ bool rb_num_nan_p(VALUE z);
 
 bool rb_num_numeric_p(VALUE self);
 
+bool rb_num_domain_p(VALUE range, VALUE x);
 
 /**
  * Rounds the argument +x+ to an arbitrary precision +prec+. <br>

data/ext/bigdecimal/math_r/math_r.c

@@ -7,6 +7,7 @@
 #include <ruby.h>
 #define DEFINE_SYMBOL
 #include "math_r/globals.h"
+#include "math_r/bigmath.h"
 
 void InitVM_Overrides(void);
 void InitVM_Const(void);
@@ -20,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);
 
 /**
@@ -35,6 +36,16 @@ 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);
+	rb_gv_set(BIGMATH_SQRT2, Qnil);
+	rb_gv_set(BIGMATH_SQRT3, Qnil);
+
 	mf_exp = rb_intern_const("exp");
 	mf_exp2 = rb_intern_const("exp2");
 	mf_expm1 = rb_intern_const("expm1");
@@ -121,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/ext/bigdecimal/math_r/solver.c

@@ -74,15 +74,15 @@ bigmath_func1 cb_casech;
 bigmath_func1 cb_acoth;
 bigmath_func1 cb_cacoth;
 
+bigmath_func1 cb_sqrt;
+bigmath_func1 cb_csqrt;
+bigmath_func1 cb_cbrt;
+bigmath_func1 cb_ccbrt;
+
+bigmath_func1 cb_erf;
+bigmath_func1 cb_erfc;
 
 
-static bool
-domain_p(VALUE range, VALUE x)
-{
-	const ID cover_p = rb_intern("cover?");
-	VALUE ans = rb_funcall(range, cover_p, 1, x);
-	return ans == Qtrue ? true : false;
-}
 
 static void
 rb_id_includes(int n, const ID *funcs, ID func)
@@ -867,7 +867,7 @@ solver_triginv(ID func, VALUE z, VALUE prec)
 			if (func == mf_asin)
 			{
 				VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), 0);
-				if (domain_p(domain, z))
+				if (rb_num_domain_p(domain, z))
 					w = asin_branch(z, prec, cb_asin);
 				else
 					w = asin_edom(z, prec);
@@ -875,7 +875,7 @@ solver_triginv(ID func, VALUE z, VALUE prec)
 			else if (func == mf_acos)
 			{
 				VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), 0);
-				if (domain_p(domain, z))
+				if (rb_num_domain_p(domain, z))
 					w = acos_branch(z, prec, cb_acos);
 				else
 					w = acos_edom(z, prec);
@@ -886,7 +886,7 @@ solver_triginv(ID func, VALUE z, VALUE prec)
 			{
 				VALUE m_domain = rb_range_new(BIG_MINUS_INF, BIG_MINUS_ONE, 0);
 				VALUE p_domain = rb_range_new(BIG_ONE, BIG_INF, 0);
-				if (domain_p(m_domain, z) || domain_p(p_domain, z))
+				if (rb_num_domain_p(m_domain, z) || rb_num_domain_p(p_domain, z))
 					w = acsc_branch(z, prec, cb_acsc);
 				else
 					w = acsc_edom(z, prec);
@@ -895,7 +895,7 @@ solver_triginv(ID func, VALUE z, VALUE prec)
 			{
 				VALUE m_domain = rb_range_new(BIG_MINUS_INF, BIG_MINUS_ONE, 0);
 				VALUE p_domain = rb_range_new(BIG_ONE, BIG_INF, 0);
-				if (domain_p(m_domain, z) || domain_p(p_domain, z))
+				if (rb_num_domain_p(m_domain, z) || rb_num_domain_p(p_domain, z))
 					w = asec_branch(z, prec, cb_asec);
 				else
 					w = asec_edom(z, prec);
@@ -1121,9 +1121,9 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 			{
 				VALUE m_domain = rb_range_new(BIG_MINUS_INF, BIG_MINUS_ONE, 0);
 				VALUE p_domain = rb_range_new(BIG_ONE, BIG_INF, 0);
-				if (domain_p(p_domain, z))
+				if (rb_num_domain_p(p_domain, z))
 					w = acosh_branch(z, prec, cb_acosh);
-				else if (domain_p(m_domain, z))
+				else if (rb_num_domain_p(m_domain, z))
 				{
 					z = rb_num_abs(z);
 					w = acosh_branch(z, prec, cb_acosh);
@@ -1135,7 +1135,7 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 			else if (func == mf_atanh)
 			{
 				VALUE domain = rb_range_new(INT2FIX(-1), INT2FIX(1), 0);
-				if (domain_p(domain, z))
+				if (rb_num_domain_p(domain, z))
 					w = atanh_branch(z, prec, cb_atanh);
 				else
 					w = atanh_edom(z, prec);
@@ -1148,9 +1148,9 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 			{
 				VALUE m_domain = rb_range_new(BIG_MINUS_ONE, BIG_ZERO, 0);
 				VALUE p_domain = rb_range_new(BIG_ZERO, BIG_ONE, 0);
-				if (domain_p(p_domain, z))
+				if (rb_num_domain_p(p_domain, z))
 					w = asech_branch(z, prec, cb_asech);
-				else if (domain_p(m_domain, z))
+				else if (rb_num_domain_p(m_domain, z))
 				{
 					w = asech_branch(z, prec, cb_asech);
 					w = rb_Complex(w, rb_bigmath_const_pi(prec));
@@ -1162,7 +1162,7 @@ solver_hyperbinv(ID func, VALUE z, VALUE prec)
 			{
 				VALUE m_domain = rb_range_new(BIG_MINUS_INF, BIG_MINUS_ONE, 0);
 				VALUE p_domain = rb_range_new(BIG_ONE, BIG_INF, 0);
-				if (domain_p(m_domain, z) || domain_p(p_domain, z))
+				if (rb_num_domain_p(m_domain, z) || rb_num_domain_p(p_domain, z))
 					w = acoth_branch(z, prec, cb_acoth);
 				else
 					w = acoth_edom(z, prec);
@@ -1337,12 +1337,12 @@ __impl_bigmath_sqrt(VALUE unused_obj, VALUE z, VALUE prec)
 		if (rb_num_negative_p(z))
 			w = sqrt_edom(z, prec);
 		else
-			w = sqrt_newton(z, prec);
+			w = cb_sqrt(z, prec);
 	}
 	else
 	{
 		z = rb_num_canonicalize(z, prec, ARG_COMPLEX, ARG_RAWVALUE);
-		w = csqrt_formula(z, prec);
+		w = cb_csqrt(z, prec);
 	}
 	return w;
 }
@@ -1374,16 +1374,58 @@ __impl_bigmath_cbrt(VALUE unused_obj, VALUE z, VALUE prec)
 	if (rb_num_real_p(z))
 	{
 		z = rb_num_canonicalize(z, prec, ARG_REAL, ARG_RAWVALUE);
-		w = cuberoot_newton(z, prec);
+		w = cb_cbrt(z, prec);
 	}
 	else
 	{
 		z = rb_num_canonicalize(z, prec, ARG_COMPLEX, ARG_RAWVALUE);
-		w = ccbrt_formula(z, prec);
+		w = cb_ccbrt(z, prec);
 	}
 	return w;
 }
 
+/**
+ *  Computes error function of +z+.
+ *  
+ *  @param z [Numeric] Numerical argument
+ *  @param prec [Integer] Arbitrary precision
+ *  @return [BigDecimal] Real solution
+ *  @raise [ArgumentError] Occurs when +prec+ is not a positive integer.
+ *  @raise [TypeError] Occurs when +z+ is not a numeric class.
+ *  @since 0.2.1
+ *  @example
+ *   BigMathR.erf(5, 20)
+ *   #=> 0.99999999999846254021e0
+ *   BigMathR.erf(5, 40)
+ *   #=> 0.9999999999984625402055719651498116565146e0
+ */
+static VALUE
+__impl_bigmath_erf(VALUE unused_obj, VALUE z, VALUE prec)
+{
+	return erf_branch(z, prec, cb_erf);
+}
+
+/**
+ *  Computes complementary error function of +z+.
+ *  
+ *  @param z [Numeric] Numerical argument
+ *  @param prec [Integer] Arbitrary precision
+ *  @return [BigDecimal] Real solution
+ *  @raise [ArgumentError] Occurs when +prec+ is not a positive integer.
+ *  @raise [TypeError] Occurs when +z+ is not a numeric class.
+ *  @since 0.2.1
+ *  @example
+ *   BigMathR.erfc(5, 20)
+ *   #=> 0.153745979e-11
+ *   BigMathR.erfc(5, 40)
+ *   #=> 0.15374597944280348501883434854e-11
+ */
+static VALUE
+__impl_bigmath_erfc(VALUE unused_obj, VALUE z, VALUE prec)
+{
+	return erfc_branch(z, prec, cb_erfc);
+}
+
 /**
  *  Document-module: BigMathR::Solver
  *  
@@ -1463,6 +1505,15 @@ InitVM_Solver(void)
 	cb_casech = casech_logrep;
 	cb_cacoth = cacoth_logrep;
 
+	cb_sqrt = sqrt_newton;
+	cb_csqrt = csqrt_formula;
+	cb_cbrt = cuberoot_newton;
+	cb_ccbrt = ccbrt_formula;
+
+	cb_erf = erf_algo911;
+	cb_erfc = erfc_algo911;
+
+
 	rb_define_singleton_method(rb_mSolver, "exp", __impl_solver_exp, 3);
 
 	rb_define_singleton_method(rb_mBigMathR, "exp", __impl_bigmath_exp, 2);
@@ -1513,4 +1564,7 @@ InitVM_Solver(void)
 
 	rb_define_singleton_method(rb_mBigMathR, "sqrt", __impl_bigmath_sqrt, 2);
 	rb_define_singleton_method(rb_mBigMathR, "cbrt", __impl_bigmath_cbrt, 2);
+
+	rb_define_singleton_method(rb_mBigMathR, "erf", __impl_bigmath_erf, 2);
+	rb_define_singleton_method(rb_mBigMathR, "erfc", __impl_bigmath_erfc, 2);
 }

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.0"
+  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.0
+  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,10 +59,8 @@ extensions:
 - ext/bigdecimal/math_r/extconf.rb
 extra_rdoc_files: []
 files:
-- ".ruby-version"
 - CODE_OF_CONDUCT.md
 - LICENSE
-- LICENSE.txt
 - README.md
 - Rakefile
 - changelog.md
@@ -88,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
@@ -118,6 +117,7 @@ files:
 - ext/bigdecimal/math_r/api/numeric/abs.h
 - ext/bigdecimal/math_r/api/numeric/canonicalize.h
 - ext/bigdecimal/math_r/api/numeric/cmpeql.h
+- ext/bigdecimal/math_r/api/numeric/domain_p.h
 - ext/bigdecimal/math_r/api/numeric/equal_p.h
 - ext/bigdecimal/math_r/api/numeric/finite_p.h
 - ext/bigdecimal/math_r/api/numeric/imag.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

data/LICENSE.txt

@@ -1,21 +0,0 @@
-The MIT License (MIT)
-
-Copyright (c) 2025 tribusonz-2
-
-Permission is hereby granted, free of charge, to any person obtaining a copy
-of this software and associated documentation files (the "Software"), to deal
-in the Software without restriction, including without limitation the rights
-to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
-copies of the Software, and to permit persons to whom the Software is
-furnished to do so, subject to the following conditions:
-
-The above copyright notice and this permission notice shall be included in
-all copies or substantial portions of the Software.
-
-THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
-IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
-FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
-AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
-LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
-OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
-THE SOFTWARE.