num4diff 0.4.5 → 0.5.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: fc6954bb9f3eb7e7b35429c7f4429aa60f5a998ce1934dbe79b0117ae1013a51
-  data.tar.gz: 3caa0e4e7b083a435b374c0739134ac583ab149f30306327783cf239896c263b
+  metadata.gz: 236c79b20fa3b5aff7de7ac41181c458a3d42cf2c846e8a1d9e01265729fd762
+  data.tar.gz: d0ace9350577ba3cadeaa37f05061a406274194357b8912df2890efc1c0f4e38
 SHA512:
-  metadata.gz: b9e5780c5bbfecb0dd284e561775f78a8cff1530194817f4ac58f22db6f0103e851558a548e3cf04d54e4f4d009a2cad2508d1f5bb6a2e2cf33f8aa59a50174e
-  data.tar.gz: 2d4736748fa20f35f7ba11fa0e614f165d9b30a73038d4942262c30f3951d33db45bf0ffb46970197b45e3d98e51d639365a3845e69a8ebc52c4db309be9569b
+  metadata.gz: d125cbc9bfdc8cc65ff57c9057c579ce1030c167f60af69c5723dc3c08d71fcb2c36822df3671dab5b0ca944b62de885e21e795943739e31b22d8505464ed3c9
+  data.tar.gz: ee2914de197552cd9be54a17e135da56cfb2a0ef36076fb14d66bedf379b3acf31d07f4a1db95a14d9d2aa0e39fdf38b8326979f200f59e4cb8b6fed44306b67

data/CHANGELOG.md

@@ -2,6 +2,19 @@
 
 ## Unreleased
 
+## [0.5.1] - 2023-05-18
+
+### Fixed
+- Fix correspondence of k step.
+
+### Changed
+- chg wiki-uri
+
+## [0.4.6] - 2023-04-24
+
+### Fixed
+- Fix Split into single-step and multi-step solutions.
+
 ## [0.4.5] - 2023-04-14
 
 ### Fixed

data/ext/num4diff/CNum4Diff.c

@@ -6,14 +6,47 @@
 static double CNum4Diff_doEulerMethod(double yi, double xi, double h, Func func);
 static double CNum4Diff_doHeunMethod(double yi, double xi, double h, Func func);
 static double CNum4Diff_doRungeKuttaMethod(double yi, double xi, double h, Func func);
-static double CNum4Diff_doAdamsBashforthMethod(double a, double b, double y0, double h, Func func);
-static double CNum4Diff_doAdamsMoultonMethod(double a, double b, double y0, double h, Func func);
+static double CNum4Diff_doAdamsBashforthMethod(int k, double a, double b, double y0, double h, Func func);
+static double CNum4Diff_doAdamsMoultonMethod(int k, double a, double b, double y0, double h, Func func);
 static CNum4Diff _cNum4Diff = {
+    .tierMethod = {
 	.FP_eulerMethod          = CNum4Diff_doEulerMethod,
 	.FP_heunMethod           = CNum4Diff_doHeunMethod,
 	.FP_rungeKuttaMethod     = CNum4Diff_doRungeKuttaMethod,
+    },
+    .multistageMethod = {
 	.FP_adamsBashforthMethod = CNum4Diff_doAdamsBashforthMethod,
 	.FP_adamsMoultonMethod   = CNum4Diff_doAdamsMoultonMethod,
+    },
+};
+/* アダムステーブル */
+static AdamsTbl bash[] = {
+    [0] = {.s = 2,                     // k = 2
+           .bv = {[0] = 3,   [1] = -1}
+          },          
+    [1] = {.s = 12,                    // k = 3
+           .bv = {[0] = 23,  [1] = -16,  [2] = 5}
+          },
+    [2] = {.s = 24,                    // k = 4
+           .bv = {[0] = 55,  [1] = -59,  [2] = 37,  [3] = -9}
+          }, 
+    [3] = {.s = 720,                   // k = 5
+           .bv = {[0] = 1901,[1] = -2774,[2] = 2616,[3] = -1274, [4] = 251}
+          }, 
+};
+static AdamsTbl moulton[] = {
+    [0] = {.s = 2,                     // k = 2
+           .bv = {[0] = 1,  [1] = 1}
+          },           
+    [1] = {.s = 12,                    // k = 3
+           .bv = {[0] = 5,  [1] = 8,  [2] = -1}
+          },
+    [2] = {.s = 24,                    // k = 4
+           .bv = {[0] = 9,  [1] = 19, [2] = -5,  [3] = 1}
+          }, 
+    [3] = {.s = 720,                   // k = 5
+          .bv = {[0] = 251,[1] = 646,[2] = -264,[3] = 106,[4] = -19}
+          }, 
 };
 /**************************************/
 /* InterFface部                       */
@@ -21,37 +54,39 @@ static CNum4Diff _cNum4Diff = {
 /**************************************/
 /* Class部                            */
 /**************************************/
-double CNum4Diff_eulerMethod(double yi, double xi, double h, Func func)
+double CNum4Diff_Tier_eulerMethod(double yi, double xi, double h, Func func)
 {
     assert(func != 0);
 
-    return _cNum4Diff.FP_eulerMethod(yi, xi, h, func);
+    return _cNum4Diff.tierMethod.FP_eulerMethod(yi, xi, h, func);
 }
-double CNum4Diff_heunMethod(double yi, double xi, double h, Func func)
+double CNum4Diff_Tier_heunMethod(double yi, double xi, double h, Func func)
 {
     assert(func != 0);
 
-    return _cNum4Diff.FP_heunMethod(yi, xi, h, func);
+    return _cNum4Diff.tierMethod.FP_heunMethod(yi, xi, h, func);
 }
-double CNum4Diff_rungeKuttaMethod(double yi, double xi, double h, Func func)
+double CNum4Diff_Tier_rungeKuttaMethod(double yi, double xi, double h, Func func)
 {
     assert(func != 0);
 
-    return _cNum4Diff.FP_rungeKuttaMethod(yi, xi, h, func);
+    return _cNum4Diff.tierMethod.FP_rungeKuttaMethod(yi, xi, h, func);
 }
-double CNum4Diff_adamsBashforthMethod(double a, double b, double y0, double h, Func func)
+double CNum4Diff_Multistage_adamsBashforthMethod(int k, double a, double b, double y0, double h, Func func)
 {
     assert(func != 0);
     assert(a < b);
+    assert((1 < k) && (k < 6));
 
-    return _cNum4Diff.FP_adamsBashforthMethod(a, b, y0, h, func);
+    return _cNum4Diff.multistageMethod.FP_adamsBashforthMethod(k, a, b, y0, h, func);
 }
-double CNum4Diff_adamsMoultonMethod(double a, double b, double y0, double h, Func func)
+double CNum4Diff_Multistage_adamsMoultonMethod(int k, double a, double b, double y0, double h, Func func)
 {
     assert(func != 0);
     assert(a < b);
+    assert((1 < k) && (k < 6));
 
-    return _cNum4Diff.FP_adamsMoultonMethod(a, b, y0, h, func);
+    return _cNum4Diff.multistageMethod.FP_adamsMoultonMethod(k, a, b, y0, h, func);
 }
 /**************************************/
 /* 処理実行部                         */
@@ -97,50 +132,72 @@ static double CNum4Diff_doRungeKuttaMethod(double yi, double xi, double h, Func
     return yi + (k1 + 2 * k2 + 2 * k3 + k4) / 6.0; 
 }
 /*
- * アダムス・バッシュフォース法(3段)
+ * アダムス・バッシュフォース法(k段)
  */
-static double CNum4Diff_doAdamsBashforthMethod(double a, double b, double y0, double h, Func func)
+static double CNum4Diff_doAdamsBashforthMethod(int k, double a, double b, double y0, double h, Func func)
 {
     double xi = a;
-    double fi;
-    double fi1 = 0;
-    double fi2 = 0;
     double y = 0.0;
+    double *f = malloc(sizeof(double) * k);
+    int i;
+    double bk;
 
-    fi2  = y0;
-    xi = xi + h;
-    fi1  = CNum4Diff_rungeKuttaMethod(fi2, xi, h, func);
-    y = fi1;
+    f[k - 1] = y0;
+    for (i = 0; i < k - 2; i++) {
+        f[k - (i + 2)] =  CNum4Diff_Tier_rungeKuttaMethod(f[k - (i + 1)], xi, h, func);
+    }
     for (xi = xi + h; xi < b; xi += h) {
-        fi  = CNum4Diff_rungeKuttaMethod(fi1, xi, h, func);
-        y = y + h * (23 * fi - 16 * fi1 + 5 * fi2) / 12.0;
-        fi2 = fi1;
-        fi1 = fi;
+        f[0]  = CNum4Diff_Tier_rungeKuttaMethod(f[1], xi, h, func);
+        // 予測子
+        bk = 0.0;
+        for (i = 0; i < k; i++) {
+            bk += bash[k - 2].bv[i] * f[i];
+        }
+        y = f[0] + h * bk / bash[k - 2].s;
+        // f値をずらす
+        for (i = 0; i < k - 1; i++) {
+            f[k - (i + 1)] = f[k - (i + 2)];
+        }
     }
     return y;
 }
 /*
- * アダムス・ムルトン法(3段)
+ * アダムス・ムルトン法(k段)
  */
-static double CNum4Diff_doAdamsMoultonMethod(double a, double b, double y0, double h, Func func)
+static double CNum4Diff_doAdamsMoultonMethod(int k, double a, double b, double y0, double h, Func func)
 {
     double xi = a;
-    double fi;
-    double fi1 = 0;
-    double fi2 = 0;
     double y_pred = 0.0;
     double y = 0.0;
+    double *f = malloc(sizeof(double) * k);
+    double *f2 = malloc(sizeof(double) * (k + 1));
+    int i;
+    double bk;
 
-    fi2  = y0;
-    xi = xi + h;
-    fi1  = CNum4Diff_rungeKuttaMethod(fi2, xi, h, func);
-    y = fi1;
+    f[k - 1] = y0;
+    for (i = 0; i < k - 2; i++) {
+        f[k - (i + 2)] =  CNum4Diff_Tier_rungeKuttaMethod(f[k - (i + 1)], xi, h, func);
+    }
     for (xi = xi + h; xi < b; xi += h) {
-        fi  = CNum4Diff_rungeKuttaMethod(fi1, xi, h, func);
-        y_pred = y + h * (23 * fi - 16 * fi1 + 5 * fi2) / 12.0;
-        y = y + h * (5 * y_pred  + 8 * fi -1 * fi1) / 12.0;
-        fi2 = fi1;
-        fi1 = fi;
+        f[0]  = CNum4Diff_Tier_rungeKuttaMethod(f[1], xi, h, func);
+        // 予測子
+        bk = 0.0;
+        for (i = 0; i < k; i++) {
+            bk += bash[k - 2].bv[i] * f[i];
+            f2[i + 1] = f[i];
+        }
+        y_pred = f[0] + h * bk / bash[k - 2].s;
+        f2[0] = y_pred;
+        // 修正子
+        bk = 0.0;
+        for (i = 0; i < k; i++) {
+            bk += moulton[k - 2].bv[i] * f2[i];
+        }
+        y = f[0] + h * bk / moulton[k - 2].s;
+        // f値をずらす
+        for (i = 0; i < k - 1; i++) {
+            f[k - (i + 1)] = f[k - (i + 2)];
+        }
     }
     return y;
 }

data/ext/num4diff/CNum4Diff.h

@@ -4,27 +4,43 @@
 /**************************************/
 /* 構造体宣言                         */
 /**************************************/
-typedef struct _CNum4Diff CNum4Diff;
+typedef struct _CNum4Diff           CNum4Diff;
+typedef struct _CNum4DiffTier       CNum4DiffTier;
+typedef struct _CNum4DiffMultistage CNum4DiffMultistage;
 typedef double (*Func)(double x);
+typedef struct _AdamsTbl            AdamsTbl;
 
-struct _CNum4Diff
+struct _AdamsTbl
+{
+    int s;
+    int bv[5];
+};
+struct _CNum4DiffTier
+{
+    double (*FP_eulerMethod)(double yi, double xi, double h, Func func);
+    double (*FP_heunMethod)(double yi, double xi, double h, Func func);
+    double (*FP_rungeKuttaMethod)(double yi, double xi, double h, Func func);
+};
+struct _CNum4DiffMultistage
 {
-	double (*FP_eulerMethod)(double yi, double xi, double h, Func func);
-	double (*FP_heunMethod)(double yi, double xi, double h, Func func);
-	double (*FP_rungeKuttaMethod)(double yi, double xi, double h, Func func);
-	double (*FP_adamsBashforthMethod)(double a, double b, double y0, double h, Func func);
-	double (*FP_adamsMoultonMethod)(double a, double b, double y0, double h, Func func);
+    double (*FP_adamsBashforthMethod)(int k, double a, double b, double y0, double h, Func func);
+    double (*FP_adamsMoultonMethod)(int k, double a, double b, double y0, double h, Func func);
 ;
 };
+struct _CNum4Diff
+{
+    CNum4DiffTier tierMethod;
+    CNum4DiffMultistage multistageMethod;
+};
 /**************************************/
 /* define宣言                         */
 /**************************************/
 /**************************************/
 /* プロトタイプ宣言                   */
 /**************************************/
-double CNum4Diff_eulerMethod(double yi, double xi, double h, Func func);
-double CNum4Diff_heunMethod(double yi, double xi, double h, Func func);
-double CNum4Diff_rungeKuttaMethod(double yi, double xi, double h, Func func);
-double CNum4Diff_adamsBashforthMethod(double a, double b, double y0, double h, Func func);
-double CNum4Diff_adamsMoultonMethod(double a, double b, double y0, double h, Func func);
+double CNum4Diff_Tier_eulerMethod(double yi, double xi, double h, Func func);
+double CNum4Diff_Tier_heunMethod(double yi, double xi, double h, Func func);
+double CNum4Diff_Tier_rungeKuttaMethod(double yi, double xi, double h, Func func);
+double CNum4Diff_Multistage_adamsBashforthMethod(int k, double a, double b, double y0, double h, Func func);
+double CNum4Diff_Multistage_adamsMoultonMethod(int k, double a, double b, double y0, double h, Func func);
 #endif

data/lib/num4diff.rb

@@ -26,7 +26,7 @@ module Num4DiffLib
     #   @return [double] xi+hに対するyi_1の値
     #
     attach_function :eulerMethod,
-        :CNum4Diff_eulerMethod, [:double, :double, :double, :f], :double
+        :CNum4Diff_Tier_eulerMethod, [:double, :double, :double, :f], :double
     #
     # ホイン法による数値計算
     # @overload heunMethod(yi, xi, h, func)
@@ -38,7 +38,7 @@ module Num4DiffLib
     #   @return [double] xi+hに対するyi_1の値
     #
     attach_function :heunMethod,
-        :CNum4Diff_heunMethod, [:double, :double, :double, :f], :double
+        :CNum4Diff_Tier_heunMethod, [:double, :double, :double, :f], :double
     #
     # 4次のルンゲ＝クッタ法による数値計算
     # @overload  rungeKuttaMethod(yi, xi, h, func)
@@ -50,10 +50,11 @@ module Num4DiffLib
     #   @return [double] xi+hに対するyi_1の値
     #
     attach_function :rungeKuttaMethod,
-        :CNum4Diff_rungeKuttaMethod, [:double, :double, :double, :f], :double
+        :CNum4Diff_Tier_rungeKuttaMethod, [:double, :double, :double, :f], :double
     #
-    # アダムス・バッシュフォース法(3段)による数値計算
-    # @overload adamsBashforthMethod(a, b, y0, h, func)
+    # アダムス・バッシュフォース法(k段)による数値計算
+    # @overload adamsBashforthMethod(k, a, b, y0, h, func)
+    #   @param [int]    k k段アダムス法
     #   @param [double] a 下限値
     #   @param [double] b 上限値
     #   @param [double] y0 y(x0=a)の値
@@ -62,10 +63,11 @@ module Num4DiffLib
     #   @return [double] [a,b]の積分値
     #
     attach_function :adamsBashforthMethod,
-        :CNum4Diff_adamsBashforthMethod, [:double, :double, :double, :double, :f], :double
+        :CNum4Diff_Multistage_adamsBashforthMethod, [:int, :double, :double, :double, :double, :f], :double
     #
-    # アダムス・ムルトン法(3段)による数値計算
-    # @overload adamsMoultonMethod(a, b, y0, xi, h, func)
+    # アダムス・ムルトン法(k段)による数値計算
+    # @overload adamsMoultonMethod(k, a, b, y0, xi, h, func)
+    #   @param [int]    k k段アダムス法
     #   @param [double] a 下限値
     #   @param [double] b 上限値
     #   @param [double] y0 y0 y(x0=a)の値
@@ -74,5 +76,5 @@ module Num4DiffLib
     #   @return [double] [a,b]の積分値
     #
     attach_function :adamsMoultonMethod,
-        :CNum4Diff_adamsMoultonMethod, [:double, :double, :double, :double, :f], :double
+        :CNum4Diff_Multistage_adamsMoultonMethod, [:int, :double, :double, :double, :double, :f], :double
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4diff
 version: !ruby/object:Gem::Version
-  version: 0.4.5
+  version: 0.5.1
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-04-13 00:00:00.000000000 Z
+date: 2023-05-18 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: ffi-compiler
@@ -69,7 +69,9 @@ licenses:
 - MIT
 metadata:
   changelog_uri: http://github.com/siranovel/num4different/blob/main/CHANGELOG.md
+  documentation_uri: https://rubydoc.info/gems/num4diff/0.5.1
   homepage_uri: http://github.com/siranovel/num4different
+  wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4different
 post_install_message: 
 rdoc_options: []
 require_paths: