num4diff 0.0.22 → 0.3.4

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (7) hide show
  1. checksums.yaml +4 -4
  2. data/ext/num4diff/CNum4Diff.c +81 -0
  3. data/ext/num4diff/CNum4Diff.h +25 -0
  4. data/ext/num4diff/Rakefile +5 -0
  5. data/lib/num4diff.rb +28 -27
  6. metadata +34 -11
  7. data/ext/num4diff/libnum4diff.so +0 -0
checksums.yaml CHANGED
@@ -1,7 +1,7 @@
1
1
  ---
2
2
  SHA256:
3
- metadata.gz: 11beab66abec9bf696ee9b5c236dd40e78c5b98d79eea57839149153c01829b9
4
- data.tar.gz: 8389c5b61a387a0eef3cfa7fe251743a25028ba937d74eabf6bd6ab53d81e7d0
3
+ metadata.gz: b15fb62d70f51c8f9975fa35f530b134aa9aa5f642ec061a5c4d33222c04201d
4
+ data.tar.gz: ff5ca49fe1bf82ea9823f95a64953b07311b540e6b0c1ca4b695d37cbdf4fe9a
5
5
  SHA512:
6
- metadata.gz: 47df3270e03fd56b3ff55ac29ddcd5b0ed8a1a9d48d2ff113c3bf6081966cbf266fb649f0246fabd3b851dd3db5a80adcd5695702e3f2d4cfad3b84314326c1d
7
- data.tar.gz: df43a0e2609a1d6b8d5f545d38316805bdfb9f006cdf29c8e41c9619f58fe9ef3435ea8d28866c334084a639ca7ec67e10a7d846ae5511358c7544825710684d
6
+ metadata.gz: 9a54382e8876e73fd52278bf583f7d2040e2534535bcec46d132da5ece1c506dcd31beecc87ab49016c5a92a3dfabdb2dd2e217143fa86b9b867f8869cb5d157
7
+ data.tar.gz: 5b5ec4e57a03980aa36676751dc88c2b44ec28ad53f723427a29eb311f4ab589a38620ea320b319f97821b1eb9a344369e3d6b64d3d20068dbd465dfbb17cd48
data/ext/num4diff/CNum4Diff.c ADDED
@@ -0,0 +1,81 @@
1
+ #include <stdio.h>
2
+ #include <stdlib.h>
3
+ #include <assert.h>
4
+ #include "CNum4Diff.h"
5
+
6
+ static double CNum4Diff_doEulerMethod(double yi, double xi, double h, Func func);
7
+ static double CNum4Diff_doHeunMethod(double yi, double xi, double h, Func func);
8
+ static double CNum4Diff_doRungeKuttaMethod(double yi, double xi, double h, Func func);
9
+ static CNum4Diff _cNum4Diff = {
10
+ .FP_eulerMethod = CNum4Diff_doEulerMethod,
11
+ .FP_heunMethod = CNum4Diff_doHeunMethod,
12
+ .FP_rungeKuttaMethod = CNum4Diff_doRungeKuttaMethod,
13
+ };
14
+ /**************************************/
15
+ /* InterFface部 */
16
+ /**************************************/
17
+ /**************************************/
18
+ /* Class部 */
19
+ /**************************************/
20
+ double CNum4Diff_eulerMethod(double yi, double xi, double h, Func func)
21
+ {
22
+ assert(func != 0);
23
+
24
+ return _cNum4Diff.FP_eulerMethod(yi, xi, h, func);
25
+ }
26
+ double CNum4Diff_heunMethod(double yi, double xi, double h, Func func)
27
+ {
28
+ assert(func != 0);
29
+
30
+ return _cNum4Diff.FP_heunMethod(yi, xi, h, func);
31
+ }
32
+ double CNum4Diff_rungeKuttaMethod(double yi, double xi, double h, Func func)
33
+ {
34
+ assert(func != 0);
35
+
36
+ return _cNum4Diff.FP_rungeKuttaMethod(yi, xi, h, func);
37
+ }
38
+ /**************************************/
39
+ /* 処理実行部 */
40
+ /**************************************/
41
+ /*
42
+ * オイラー法
43
+ */
44
+ static double CNum4Diff_doEulerMethod(double yi, double xi, double h, Func func)
45
+ {
46
+ double y = func(xi);
47
+
48
+ return yi + h * y;
49
+ }
50
+ /*
51
+ * ホイン法
52
+ */
53
+ static double CNum4Diff_doHeunMethod(double yi, double xi, double h, Func func)
54
+ {
55
+ double y = func(xi);
56
+ double k1 = 0.0;
57
+ double k2 = 0.0;
58
+
59
+ k1 = h * y; // k1 = h * f(xi, y)
60
+ k2 = h * (yi + k1); // k2 = h * f(xi_h, h * f(xi + h, yi + k1)
61
+ return yi + (k1 + k2) / 2.0;
62
+ }
63
+ /*
64
+ * ルンゲクッタ法
65
+ */
66
+ static double CNum4Diff_doRungeKuttaMethod(double yi, double xi, double h, Func func)
67
+ {
68
+ double y = 0.0;
69
+ double k1 = 0.0;
70
+ double k2 = 0.0;
71
+ double k3 = 0.0;
72
+ double k4 = 0.0;
73
+
74
+ y = func(xi);
75
+ k1 = h * y; // k1 = h * f(xi, y)
76
+ k2 = h * (yi + k1 / 2.0); // k2 = h * f(xi + h / 2, yi + k1 / 2)
77
+ k3 = h * (yi + k2 / 2.0); // k3 = h * f(xi + h / 2, y1 + k2 / 2)
78
+ k4 = h * (yi + k3); // k4 = h * f(xi + h, yi + k3)
79
+ return yi + (k1 + 2 * k2 + 2 * k3 + k4) / 6.0;
80
+ }
81
+
data/ext/num4diff/CNum4Diff.h ADDED
@@ -0,0 +1,25 @@
1
+ #ifndef _CNum4Diff_H_
2
+ #define _CNum4Diff_H_
3
+
4
+ /**************************************/
5
+ /* 構造体宣言 */
6
+ /**************************************/
7
+ typedef struct _CNum4Diff CNum4Diff;
8
+ typedef double (*Func)(double x);
9
+
10
+ struct _CNum4Diff
11
+ {
12
+ double (*FP_eulerMethod)(double yi, double xi, double h, Func func);
13
+ double (*FP_heunMethod)(double yi, double xi, double h, Func func);
14
+ double (*FP_rungeKuttaMethod)(double yi, double xi, double h, Func func);
15
+ };
16
+ /**************************************/
17
+ /* define宣言 */
18
+ /**************************************/
19
+ /**************************************/
20
+ /* プロトタイプ宣言 */
21
+ /**************************************/
22
+ double CNum4Diff_eulerMethod(double yi, double xi, double h, Func func);
23
+ double CNum4Diff_heunMethod(double yi, double xi, double h, Func func);
24
+ double CNum4Diff_rungeKuttaMethod(double yi, double xi, double h, Func func);
25
+ #endif
data/ext/num4diff/Rakefile ADDED
@@ -0,0 +1,5 @@
1
+ require 'ffi-compiler/compile_task'
2
+
3
+ FFI::Compiler::CompileTask.new('num4diff') do |c|
4
+ c.have_header?('num4diff.h', '.')
5
+ end
data/lib/num4diff.rb CHANGED
@@ -1,47 +1,48 @@
1
1
  require 'ffi'
2
+ require 'ffi-compiler/loader'
2
3
 
3
4
  #
4
5
  # 数値計算による常微分方程式の解法するライブラリ
5
6
  #
6
7
  module Num4DiffLib
7
8
  extend FFI::Library
8
- spec = Gem.loaded_specs["num4diff"]
9
- gem_root = spec.gem_dir
10
- ffi_lib ["#{gem_root}/ext/num4diff/libnum4diff.so"]
11
- # y = f(xi)
12
- # @param [double] xi xiの値
13
- # @return [double] y xiに対するyの値
9
+
10
+ ffi_lib FFI::Compiler::Loader.find('num4diff')
11
+ # @overload f(xi)
12
+ # @param [double] xi xiの値
13
+ # @return [double] xiに対するyの値
14
14
  callback :f, [:double], :double
15
15
 
16
16
  #
17
- # オイラー法による数値計算
18
- # yi_1 = eulerMethod(yi, xi, h, func)
19
- # @param [double] yi xiに対するyiの値
20
- # @param [double] xi xiの値
21
- # @param [double] h 刻み幅
22
- # @param [callback] func xiに対する微分を計算する関数
23
- # @return [double] yi_1 xi+hに対するyiの値
17
+ # @overload eulerMethod(yi, xi, h, func)
18
+ # オイラー法による数値計算
19
+ # @param [double] yi xiに対するyiの値
20
+ # @param [double] xi xiの値
21
+ # @param [double] h 刻み幅
22
+ # @param [callback] func xiに対する傾きを計算する関数
23
+ # @return [double] xi+hに対するyiの値
24
24
  #
25
25
  attach_function :eulerMethod,
26
26
  :CNum4Diff_eulerMethod, [:double, :double, :double, :f], :double
27
- # ホイン法による数値計算
28
- # yi_1 = heunMethod(yi, xi, h, func)
29
- # @param [double] yi xiに対するyiの値
30
- # @param [double] xi xiの値
31
- # @param [double] h 刻み幅
32
- # @param [callback] func xiに対する微分を計算する関数
33
- # @return [double] yi_1 xi+hに対するyiの値
27
+ # @overload heunMethod(yi, xi, h, func)
28
+ # ホイン法による数値計算
29
+ # @param [double] yi xiに対するyiの値
30
+ # @param [double] xi xiの値
31
+ # @param [double] h 刻み幅
32
+ # @param [callback] func xiに対する傾きを計算する関数
33
+ # @return [double] xi+hに対するyiの値
34
34
  #
35
35
  attach_function :heunMethod,
36
36
  :CNum4Diff_heunMethod, [:double, :double, :double, :f], :double
37
37
  #
38
- # 4次のルンゲ＝クッタ法による数値計算
39
- # yi_1 = rungeKuttaMethod(yi, xi, h, func)
40
- # @param [double] yi xiに対するyiの値
41
- # @param [double] xi xiの値
42
- # @param [double] h 刻み幅
43
- # @param [callback] func xiに対する微分を計算する関数
44
- # @return [double] yi_1 xi+hに対するyiの値
38
+ # @overload rungeKuttaMethod(yi, xi, h, func)
39
+ # 4次のルンゲ＝クッタ法による数値計算
40
+ # @param [double] yi xiに対するyiの値
41
+ # @param [double] xi xiの値
42
+ # @param [double] h 刻み幅
43
+ # @param [callback] func xiに対する傾きを計算する関数
44
+ # @return [double] xi+hに対するyiの値
45
+ #
45
46
  attach_function :rungeKuttaMethod,
46
47
  :CNum4Diff_rungeKuttaMethod, [:double, :double, :double, :f], :double
47
48
  end
metadata CHANGED
@@ -1,43 +1,66 @@
1
1
  --- !ruby/object:Gem::Specification
2
2
  name: num4diff
3
3
  version: !ruby/object:Gem::Version
4
- version: 0.0.22
4
+ version: 0.3.4
5
5
  platform: ruby
6
6
  authors:
7
7
  - siranovel
8
8
  autorequire:
9
9
  bindir: bin
10
10
  cert_chain: []
11
- date: 2022-04-04 00:00:00.000000000 Z
11
+ date: 2023-02-19 00:00:00.000000000 Z
12
12
  dependencies:
13
13
  - !ruby/object:Gem::Dependency
14
- name: ffi
14
+ name: ffi-compiler
15
15
  requirement: !ruby/object:Gem::Requirement
16
16
  requirements:
17
- - - ">="
18
- - !ruby/object:Gem::Version
19
- version: 1.15.0
20
17
  - - "~>"
21
18
  - !ruby/object:Gem::Version
22
- version: '1.15'
19
+ version: '1.0'
20
+ - - ">="
21
+ - !ruby/object:Gem::Version
22
+ version: 1.0.1
23
23
  type: :runtime
24
24
  prerelease: false
25
25
  version_requirements: !ruby/object:Gem::Requirement
26
26
  requirements:
27
+ - - "~>"
28
+ - !ruby/object:Gem::Version
29
+ version: '1.0'
27
30
  - - ">="
28
31
  - !ruby/object:Gem::Version
29
- version: 1.15.0
32
+ version: 1.0.1
33
+ - !ruby/object:Gem::Dependency
34
+ name: rake
35
+ requirement: !ruby/object:Gem::Requirement
36
+ requirements:
30
37
  - - "~>"
31
38
  - !ruby/object:Gem::Version
32
- version: '1.15'
39
+ version: '12.3'
40
+ - - ">="
41
+ - !ruby/object:Gem::Version
42
+ version: 12.3.3
43
+ type: :runtime
44
+ prerelease: false
45
+ version_requirements: !ruby/object:Gem::Requirement
46
+ requirements:
47
+ - - "~>"
48
+ - !ruby/object:Gem::Version
49
+ version: '12.3'
50
+ - - ">="
51
+ - !ruby/object:Gem::Version
52
+ version: 12.3.3
33
53
  description: numerical solution for ordinaray differential equations
34
54
  email: siranovel@gmail.com
35
55
  executables: []
36
- extensions: []
56
+ extensions:
57
+ - ext/num4diff/Rakefile
37
58
  extra_rdoc_files: []
38
59
  files:
39
60
  - LICENSE
40
- - ext/num4diff/libnum4diff.so
61
+ - ext/num4diff/CNum4Diff.c
62
+ - ext/num4diff/CNum4Diff.h
63
+ - ext/num4diff/Rakefile
41
64
  - lib/num4diff.rb
42
65
  homepage: http://github.com/siranovel/num4different
43
66
  licenses:
data/ext/num4diff/libnum4diff.so DELETED
Binary file