num4diff 0.0.17 → 0.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f446c837fe8380d46e1467a9a2c1365fd6b6e5b5571f44be19d94e3f102d566e
-  data.tar.gz: f5618c68783af3470c5f8f3304c94d4b9307a9eabcba53f470507a72f4d32d56
+  metadata.gz: 0d3136cc91be5e19847b49e534e7fb8292c8e8d6e2e6ea7b64c19a8e4f0f5e87
+  data.tar.gz: 5687b8ae545c9df2da7b39fa2ebe6ae6017276edd8d1e5be750c583d7809c5ff
 SHA512:
-  metadata.gz: 6ccacc59e40dcab11980802247733c4efce0bfa553a93a9d204a972e574a49e422c36fbf51a46509bf8c6dd931f85b7c9ddfd89c28bcd5fe33841dde735e54e6
-  data.tar.gz: 8c51b60e4d77cabdac3c4b49c2f41a2b7787f8b6ecf9b7779a98f427dedb92c49a2376ed759fba43eb615cc8e05410e3474f3c996b3e0a7c582b1a449c45a783
+  metadata.gz: 1707e0ebc7e844bc39c64baff378a57ed2cffe80ad7ff3b1ebe51b9d6a2095b31b6018996cb8179ce298b4d7c7c04da5b758dcb896279db8b6703e4c5a3e189e
+  data.tar.gz: b9a91f1ebeccff6370a47204ec5de6b2a3e2b42e08934e609f512b313aeff78e53658d5b511585bb87ac82d6ecaced36e463446b70fdbf36efa92122607b437a

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2022 siranovel
+
+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.

data/ext/num4diff/CNum4Diff.c

@@ -0,0 +1,81 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <assert.h>
+#include "CNum4Diff.h"
+
+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 CNum4Diff _cNum4Diff = {
+	.FP_eulerMethod = CNum4Diff_doEulerMethod,
+	.FP_heunMethod = CNum4Diff_doHeunMethod,
+	.FP_rungeKuttaMethod = CNum4Diff_doRungeKuttaMethod,
+};
+/**************************************/
+/* InterFface部                       */
+/**************************************/
+/**************************************/
+/* Class部                            */
+/**************************************/
+double CNum4Diff_eulerMethod(double yi, double xi, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Diff.FP_eulerMethod(yi, xi, h, func);
+}
+double CNum4Diff_heunMethod(double yi, double xi, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Diff.FP_heunMethod(yi, xi, h, func);
+}
+double CNum4Diff_rungeKuttaMethod(double yi, double xi, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Diff.FP_rungeKuttaMethod(yi, xi, h, func);
+}
+/**************************************/
+/* 処理実行部                         */
+/**************************************/
+/* 
+ * オイラー法
+ */
+static double CNum4Diff_doEulerMethod(double yi, double xi, double h, Func func)
+{
+    double y = func(xi);
+
+    return yi + h * y;
+}
+/*
+ * ホイン法
+ */
+static double CNum4Diff_doHeunMethod(double yi, double xi, double h, Func func)
+{
+    double y = func(xi);
+    double k1 = 0.0;
+    double k2 = 0.0;
+
+    k1 = h * y;                        // k1 = h * f(xi, y)
+    k2 = h * (yi + k1);                // k2 = h * f(xi_h, h * f(xi + h, yi + k1)
+    return yi + (k1 + k2) / 2.0;
+}
+/*
+ * ルンゲクッタ法
+ */
+static double CNum4Diff_doRungeKuttaMethod(double yi, double xi, double h, Func func)
+{
+    double y = 0.0;
+    double k1 = 0.0;
+    double k2 = 0.0;
+    double k3 = 0.0;
+    double k4 = 0.0;
+
+    y = func(xi);
+    k1 = h * y;                        // k1 = h * f(xi, y)
+    k2 = h * (yi + k1 / 2.0);          // k2 = h * f(xi + h / 2, yi + k1 / 2) 
+    k3 = h * (yi + k2 / 2.0);          // k3 = h * f(xi + h / 2, y1 + k2 / 2)
+    k4 = h * (yi + k3);                // k4 = h * f(xi + h, yi + k3)
+    return yi + (k1 + 2 * k2 + 2 * k3 + k4) / 6.0; 
+}
+

data/ext/num4diff/CNum4Diff.h

@@ -0,0 +1,25 @@
+#ifndef _CNum4Diff_H_
+#define _CNum4Diff_H_
+
+/**************************************/
+/* 構造体宣言                         */
+/**************************************/
+typedef struct _CNum4Diff CNum4Diff;
+typedef double (*Func)(double x);
+
+struct _CNum4Diff
+{
+	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);
+};
+/**************************************/
+/* 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);
+#endif

data/ext/num4diff/Rakefile

@@ -0,0 +1,5 @@
+require 'ffi-compiler/compile_task'
+
+FFI::Compiler::CompileTask.new('num4diff') do |c|
+  c.have_header?('num4diff.h', '.')
+end

data/lib/num4diff.rb

@@ -1,19 +1,48 @@
 require 'ffi'
+require 'ffi-compiler/loader'
 
+#
+# 数値計算による常微分方程式の解法するライブラリ
+#
 module Num4DiffLib
     extend FFI::Library
-    spec = Gem.loaded_specs["num4diff"]
-    gem_root = spec.gem_dir
-    ffi_lib ["#{gem_root}/ext/num4diff/libnum4diff.so"]
-    
+
+    ffi_lib FFI::Compiler::Loader.find('num4diff') 
+    # dy = f(xi)
+    # @param [double] xi xiの値
+    # @return [double] xiに対するyの値
     callback   :f, [:double], :double
-    # yi_1 = eulerMethod(yi, x, h, func)
+
+    #
+    # オイラー法による数値計算
+    #   yi_1 = eulerMethod(yi, xi, h, func)
+    # @param [double] yi xiに対するyiの値
+    # @param [double] xi xiの値
+    # @param [double] h  刻み幅
+    # @param [callback] func xiに対する傾きを計算する関数
+    # @return [double] xi+hに対するyiの値
+    #
     attach_function :eulerMethod,
         :CNum4Diff_eulerMethod, [:double, :double, :double, :f], :double
-    # yi_1 = heunMethod(yi, x, h, func)
+    # ホイン法による数値計算
+    #   yi_1 = heunMethod(yi, xi, h, func)
+    # @param [double] yi xiに対するyiの値
+    # @param [double] xi xiの値
+    # @param [double] h  刻み幅
+    # @param [callback] func xiに対する傾きを計算する関数
+    # @return [double] xi+hに対するyiの値
+    #
     attach_function :heunMethod,
         :CNum4Diff_heunMethod, [:double, :double, :double, :f], :double
-    # yi_1 = rungeKuttaMethod(yi, x, h, func)
+    #
+    # 4次のルンゲ＝クッタ法による数値計算
+    #   yi_1 = rungeKuttaMethod(yi, xi, h, func)
+    # @param [double] yi xiに対するyiの値
+    # @param [double] xi xiの値
+    # @param [double] h  刻み幅
+    # @param [callback] func xiに対する傾きを計算する関数
+    # @return [double] xi+hに対するyiの値
+    #
     attach_function :rungeKuttaMethod,
         :CNum4Diff_rungeKuttaMethod, [:double, :double, :double, :f], :double
 end

metadata

@@ -1,22 +1,66 @@
 --- !ruby/object:Gem::Specification
 name: num4diff
 version: !ruby/object:Gem::Version
-  version: 0.0.17
+  version: 0.1.3
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2022-03-17 00:00:00.000000000 Z
-dependencies: []
+date: 2022-04-14 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: ffi-compiler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.0.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.0.1
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
 description: numerical solution for ordinaray differential equations
 email: siranovel@gmail.com
 executables: []
-extensions: []
+extensions:
+- ext/num4diff/Rakefile
 extra_rdoc_files: []
 files:
-- ext/num4diff/libnum4diff.so
+- LICENSE
+- ext/num4diff/CNum4Diff.c
+- ext/num4diff/CNum4Diff.h
+- ext/num4diff/Rakefile
 - lib/num4diff.rb
 homepage: http://github.com/siranovel/num4different
 licenses:
@@ -40,5 +84,5 @@ requirements: []
 rubygems_version: 3.0.3.1
 signing_key: 
 specification_version: 4
-summary: num for diffrent!
+summary: num for different!
 test_files: []