num4diff 0.5.1 → 0.5.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 236c79b20fa3b5aff7de7ac41181c458a3d42cf2c846e8a1d9e01265729fd762
-  data.tar.gz: d0ace9350577ba3cadeaa37f05061a406274194357b8912df2890efc1c0f4e38
+  metadata.gz: 4032b47cfd83b41e99e62fb6d2f7085af8c17692de3a7f40abc3480073213d28
+  data.tar.gz: 725a5365017a298863de4504382f47ef2340e0b52976eb925164269509809991
 SHA512:
-  metadata.gz: d125cbc9bfdc8cc65ff57c9057c579ce1030c167f60af69c5723dc3c08d71fcb2c36822df3671dab5b0ca944b62de885e21e795943739e31b22d8505464ed3c9
-  data.tar.gz: ee2914de197552cd9be54a17e135da56cfb2a0ef36076fb14d66bedf379b3acf31d07f4a1db95a14d9d2aa0e39fdf38b8326979f200f59e4cb8b6fed44306b67
+  metadata.gz: 810201bbeef874ac1c946cf454815fb0eb7dcc72cadd5b4ffd5f1150c3bcb5942e342127ef1b2440794a44e797e40922ab5aae55d8ba7be42fb63ab69ec7bcef
+  data.tar.gz: e330e6e32a55a08798682bdcd83e47b4e9b591d2489491884f07a43b86b530ae95628c2b08578ef112adbd13526a8ddc96645afa437ecfb078436597e3ffd2c5

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [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,6 +141,7 @@ 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++) {
@@ -151,9 +152,9 @@ static double CNum4Diff_doAdamsBashforthMethod(int k, double a, double b, double
         // 予測子
         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,6 +174,8 @@ 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++) {
@@ -183,17 +186,17 @@ static double CNum4Diff_doAdamsMoultonMethod(int k, double a, double b, double y
         // 予測子
         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.2
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-05-18 00:00:00.000000000 Z
+date: 2023-06-08 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.2
   homepage_uri: http://github.com/siranovel/num4different
   wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4different
 post_install_message: