num4diff 0.5.1 → 0.5.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 236c79b20fa3b5aff7de7ac41181c458a3d42cf2c846e8a1d9e01265729fd762
-  data.tar.gz: d0ace9350577ba3cadeaa37f05061a406274194357b8912df2890efc1c0f4e38
+  metadata.gz: 29001c8ab562e1167f790ad60e1aad72bb7911a1fdd93f6124437818bf1d6ace
+  data.tar.gz: 252f28a67131a94080d1226462ec8759c3362a19abfa31c3d315069cadc74dbe
 SHA512:
-  metadata.gz: d125cbc9bfdc8cc65ff57c9057c579ce1030c167f60af69c5723dc3c08d71fcb2c36822df3671dab5b0ca944b62de885e21e795943739e31b22d8505464ed3c9
-  data.tar.gz: ee2914de197552cd9be54a17e135da56cfb2a0ef36076fb14d66bedf379b3acf31d07f4a1db95a14d9d2aa0e39fdf38b8326979f200f59e4cb8b6fed44306b67
+  metadata.gz: ac10c3cbcee6872cfeacf33830aea5515cf9a3756dcc8388cb2c5af1da463118d9d9181bbbfb715ff53cd649075f013cff6fbad244ae9e1ce1cacf4c9c3825ac
+  data.tar.gz: 6707b0b1feac33281df4d61d25c8b89c475fd208f62ad2d14e45af0b0cccb443beac22b59e89071c0126876dcd843778ee5bf84fe2d54862d2cf467a4cab8176

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.5.2] - 2023-07-02
+
+### Fixed
+- change bug for  adamsBashforthMethod and adamsMoultonMethod.
+
+## [0.5.2] - 2023-06-08
+
+### Fixed
+- Fix refactor to adamsBashforthMethod and adamsMoultonMethod.
+
 ## [0.5.1] - 2023-05-18
 
 ### Fixed

data/ext/num4diff/CNum4Diff.c

@@ -141,19 +141,21 @@ static double CNum4Diff_doAdamsBashforthMethod(int k, double a, double b, double
     double *f = malloc(sizeof(double) * k);
     int i;
     double bk;
+    AdamsTbl *ptBash = &bash[k - 2];
 
     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);
+        xi = xi + h;
     }
-    for (xi = xi + h; xi < b; xi += h) {
+    for (xi = xi; xi < b; xi += h) {
         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];
+            bk += ptBash->bv[i] * f[i];
         }
-        y = f[0] + h * bk / bash[k - 2].s;
+        y = f[0] + h * bk / ptBash->s;
         // f値をずらす
         for (i = 0; i < k - 1; i++) {
             f[k - (i + 1)] = f[k - (i + 2)];
@@ -173,27 +175,30 @@ static double CNum4Diff_doAdamsMoultonMethod(int k, double a, double b, double y
     double *f2 = malloc(sizeof(double) * (k + 1));
     int i;
     double bk;
+    AdamsTbl *ptBash = &bash[k - 2];
+    AdamsTbl *ptMoulton = &moulton[k - 2];
 
     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);
+        xi = xi + h;
     }
-    for (xi = xi + h; xi < b; xi += h) {
+    for (xi = xi; xi < b; xi += h) {
         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];
+            bk += ptBash->bv[i] * f[i];
             f2[i + 1] = f[i];
         }
-        y_pred = f[0] + h * bk / bash[k - 2].s;
+        y_pred = f[0] + h * bk / ptBash->s;
         f2[0] = y_pred;
         // 修正子
         bk = 0.0;
         for (i = 0; i < k; i++) {
-            bk += moulton[k - 2].bv[i] * f2[i];
+            bk += ptMoulton->bv[i] * f2[i];
         }
-        y = f[0] + h * bk / moulton[k - 2].s;
+        y = f[0] + h * bk / ptMoulton->s;
         // f値をずらす
         for (i = 0; i < k - 1; i++) {
             f[k - (i + 1)] = f[k - (i + 2)];

metadata

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