num4simdiff 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 84111ebd98c9c8f7f000733a780973aa030185287c716b15d4a6d7ef5c62424d
+  data.tar.gz: 86d00d4725e4af07cfb5aff14647be96645f8c8fa17080f310da93817e0ca4cd
+SHA512:
+  metadata.gz: 93a012aff1ff1fd8f8614c483f12836e3ff5e8c0d39adbf1f8c277503c9669d36f4665b7a9d9397cd5ce5d572ceefd8b8a49e13b4c073074e3176fa9be0d0ad4
+  data.tar.gz: c6c1d8bb3a3f35f82becdf0a584d7c9971593d3435296e39c9cb349a8b82095193c979a6172f2423b719780430a5afc1edf7f20e2c984f1e6d14a72094207cba

data/CHANGELOG.md

@@ -0,0 +1,9 @@
+# Change Log
+
+## Unreleased
+
+## [0.0.1] - 2023-04-29
+
+### Fixed
+- fix first fixed.
+

data/Gemfile

@@ -0,0 +1,3 @@
+source "https://rubygems.org"
+
+gemspec

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2023 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/num4simdiff/CNum4SimDiff.c

@@ -0,0 +1,110 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <assert.h>
+#include <stdlib.h>
+#include "CNum4SimDiff.h"
+
+static double* CNum4SimDiff_doEulerMethod(double *yi, double xi, double h, Func func1, Func func2);
+static double* CNum4SimDiff_doHeunMethod(double *yi, double xi, double h, Func func1, Func func2);
+static double* CNum4SimDiff_doRungeKuttaMethod(double *yi, double xi, double h, Func func1, Func func2);
+static CNum4SimDiff _cNum4SiDiff = {
+    .FP_eulerMethod      = CNum4SimDiff_doEulerMethod,
+    .FP_heunMethod       = CNum4SimDiff_doHeunMethod,
+    .FP_rungeKuttaMethod = CNum4SimDiff_doRungeKuttaMethod,
+};
+/**************************************/
+/* InterFface部                       */
+/**************************************/
+/**************************************/
+/* Class部                            */
+/**************************************/
+double* CNum4SimDiff_eulerMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    assert(func1 != 0);
+    assert(func2 != 0);
+
+    return _cNum4SiDiff.FP_eulerMethod(yi, xi, h, func1, func2);
+}
+double* CNum4SimDiff_heunMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    assert(func1 != 0);
+    assert(func2 != 0);
+
+    return _cNum4SiDiff.FP_heunMethod(yi, xi, h, func1, func2);
+}
+double* CNum4SimDiff_rungeKuttaMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    assert(func1 != 0);
+    assert(func2 != 0);
+
+    return _cNum4SiDiff.FP_rungeKuttaMethod(yi, xi, h, func1, func2);
+}
+/**************************************/
+/* 処理実行部                         */
+/**************************************/
+/* 
+ * オイラー法
+ */
+static double* CNum4SimDiff_doEulerMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    int i;
+    double *f = malloc(sizeof(double) * N);
+    double *yi_1 = malloc(sizeof(double) * N);
+
+    f[0] = func1(xi, yi[1]);
+    f[1] = func2(xi, yi[1]);
+    // yi_1 = yi + h * f(xi, y)
+    for (i = 0; i < N; i++) {
+        yi_1[i] = yi[i] + h * f[i];
+    }
+    return yi_1;
+}
+/*
+ * ホイン法
+ */
+static double* CNum4SimDiff_doHeunMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    int i;
+    double *f = malloc(sizeof(double) * N);
+    double *yi_1 = malloc(sizeof(double) * N);
+    double *k1 = malloc(sizeof(double) * N);
+    double *k2 = malloc(sizeof(double) * N);
+
+    f[0] = func1(xi, yi[1]);
+    f[1] = func2(xi, yi[1]);
+    // k1 = h * f(xi, y)
+    // k2 = h * f(xi + h, h * f(xi + h, yi + k1)
+    // yi_1 = yi + (k1 + k2) / 2.0;
+    for (i = 0; i < N; i++) {
+        k1[i] = h * f[i];
+        k2[i] = h * (yi[i] + k1[i]);
+        yi_1[i] = yi[i] + (k1[i] + k2[i]) / 2.0;
+    }
+    return yi_1;
+}
+static double* CNum4SimDiff_doRungeKuttaMethod(double *yi, double xi, double h, Func func1, Func func2)
+{
+    int i;
+    double *f = malloc(sizeof(double) * N);
+    double *yi_1 = malloc(sizeof(double) * N);
+    double *k1 = malloc(sizeof(double) * N);
+    double *k2 = malloc(sizeof(double) * N);
+    double *k3 = malloc(sizeof(double) * N);
+    double *k4 = malloc(sizeof(double) * N);
+
+    f[0] = func1(xi, yi[1]);
+    f[1] = func2(xi, yi[1]);
+    // k1 = h * f(xi, y)
+    // k2 = h * f(xi + h / 2, yi + k1 / 2)
+    // k3 = h * f(xi + h / 2, y1 + k2 / 2)
+    // k4 = h * f(xi + h, yi + k3)
+    for (i = 0; i < N; i++) {
+        k1[i] = h * f[i];
+        k2[i] = h * (yi[i] + k1[i] / 2.0);
+        k3[i] = h * (yi[i] + k2[i] / 2.0); 
+        k4[i] = h * (yi[i] + k3[i]);
+        yi_1[i] = yi[i] + (k1[i] + 2 * k2[i] + 2 * k3[i] + k4[i]) / 6.0;
+    }
+    return yi_1;
+}
+

data/ext/num4simdiff/CNum4SimDiff.h

@@ -0,0 +1,26 @@
+#ifndef _CNum4SimDiff_H_
+#define _CNum4SimDiff_H_
+
+/**************************************/
+/* 構造体宣言                          */
+/**************************************/
+typedef struct _CNum4SimDiff CNum4SimDiff;
+typedef double (*Func)(double x, double dx);
+
+struct _CNum4SimDiff
+{
+    double* (*FP_eulerMethod)(double *yi, double xi, double h, Func func1, Func func2);
+    double* (*FP_heunMethod)(double *yi, double xi, double h, Func func1, Func func2);
+    double* (*FP_rungeKuttaMethod)(double *yi, double xi, double h, Func func1, Func func2);
+};
+/**************************************/
+/* define宣言                         */
+/**************************************/
+#define N 2
+/**************************************/
+/* プロトタイプ宣言                    */
+/**************************************/
+double* CNum4SimDiff_eulerMethod(double *yi, double xi, double h, Func func1, Func func2);
+double* CNum4SimDiff_heunMethod(double *yi, double xi, double h, Func func1, Func func2);
+double* CNum4SimDiff_rungeKuttaMethod(double *yi, double xi, double h, Func func1, Func func2);
+#endif

data/ext/num4simdiff/Rakefile

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

data/lib/num4simdiff.rb

@@ -0,0 +1,97 @@
+require 'ffi'
+require 'ffi-compiler/loader'
+require 'fiddle'
+
+module Num4SimDiffLib
+    extend FFI::Library
+
+    ffi_lib FFI::Compiler::Loader.find('num4simdiff')
+
+    # @overload f(x, dx)
+    #   dy = f(x, dx)
+    #   @yield [x] dy = f(x, dx)
+    #   @yieldparam [double] x xiの値
+    #   @yieldparam [double] dx xiの値
+    #   @return [double] xiに対するyの値
+    callback   :f, [:double, :double], :double
+
+    #
+    # オイラー法による数値計算
+    # @overload eulerMethod(yi, xi, h, func1, func2)
+    #   yi_1 = eulerMethod(yi, xi, h, func1, func2)
+    #   @param [double] yi xiに対するyiの値(配列)
+    #   @param [double] xi xiの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func1 xiに対する傾きを計算する関数
+    #   @param [callback] func2 xiに対する傾きを計算する関数
+    #   @return [double] xi+hに対するyi_1の値(配列)
+    #
+    attach_function :eulerMethodFFI,
+        :CNum4SimDiff_eulerMethod, [:buffer_in, :double, :double, :f, :f], :pointer
+    #
+    # ホイン法による数値計算
+    # @overload heunMethod(yi, xi, h, func1, func2)
+    #   yi_1 = heunMethod(yi, xi, h, func1, func2)
+    #   @param [double] yi xiに対するyiの値(配列)
+    #   @param [double] xi xiの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func1 xiに対する傾きを計算する関数
+    #   @param [callback] func2 xiに対する傾きを計算する関数
+    #   @return [double] xi+hに対するyi_1の値(配列)
+    #
+    attach_function :heunMethodFFI,
+        :CNum4SimDiff_heunMethod, [:buffer_in, :double, :double, :f, :f], :pointer
+    #
+    # 4次のルンゲ＝クッタ法による数値計算
+    # @overload  rungeKuttaMethod(yi, xi, h, func1, func2)
+    #   yi_1 = rungeKuttaMethod(yi, xi, h, func1, func2)
+    #   @param [double] yi xiに対するyiの値(配列)
+    #   @param [double] xi xiの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func1 xiに対する傾きを計算する関数
+    #   @param [callback] func2 xiに対する傾きを計算する関数
+    #   @return [double] xi+hに対するyi_1の値(配列)
+    #
+    attach_function :rungeKuttaMethodFFI,
+        :CNum4SimDiff_rungeKuttaMethod, [:buffer_in, :double, :double, :f, :f], :pointer
+    class << self
+        def eulerMethod(yi, x, h, func1, func2)
+          yi_ptr = cnvRbAry2pt(2, yi)
+          yi_1_ptr = eulerMethodFFI(yi_ptr, x, h, func1, func2)
+          yi_1 = cnvPt2RbAry(2, yi_1_ptr)
+          return 1 + yi_1[0] + 0.5 * yi_1[1]
+        end
+        def heunMethod(yi, x, h, func1, func2)
+          yi_ptr = cnvRbAry2pt(2, yi)
+          yi_1_ptr = heunMethodFFI(yi_ptr, x, h, func1, func2)
+          yi_1 = cnvPt2RbAry(2, yi_1_ptr)
+          return 1 + yi_1[0] + 0.5 * yi_1[1] 
+        end
+        def rungeKuttaMethod(yi, x, h, func1, func2)
+          yi_ptr = cnvRbAry2pt(2, yi)
+          yi_1_ptr = rungeKuttaMethodFFI(yi_ptr, x, h, func1, func2)
+          yi_1 = cnvPt2RbAry(2, yi_1_ptr)
+          return 1 + yi_1[0] + 0.5 * yi_1[1]
+        end
+
+
+        def cnvRbAry2pt(n, ary)
+            yi_ptr = FFI::MemoryPointer.new(:double, n)
+            n.times.map { |i|
+                yi_ptr.put_double(i * Fiddle::SIZEOF_DOUBLE, ary[i].to_f)
+            }
+            return yi_ptr   
+        end
+        def cnvPt2RbAry(n, pt)
+          rbAry = n.times.map { |i|
+             pt.get_double(i * Fiddle::SIZEOF_DOUBLE)
+          }
+          return rbAry  
+        end
+        private :eulerMethodFFI
+        private :heunMethodFFI
+        private :rungeKuttaMethodFFI
+        private :cnvRbAry2pt
+        private :cnvPt2RbAry
+    end
+end

metadata

@@ -0,0 +1,93 @@
+--- !ruby/object:Gem::Specification
+name: num4simdiff
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2023-04-17 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 simultaneous ordinaray differential equations
+email: siranovel@gmail.com
+executables: []
+extensions:
+- ext/num4simdiff/Rakefile
+extra_rdoc_files: []
+files:
+- CHANGELOG.md
+- Gemfile
+- LICENSE
+- ext/num4simdiff/CNum4SimDiff.c
+- ext/num4simdiff/CNum4SimDiff.h
+- ext/num4simdiff/Rakefile
+- lib/num4simdiff.rb
+homepage: http://github.com/siranovel/num4simdiffrent
+licenses:
+- MIT
+metadata:
+  changelog_uri: http://github.com/siranovel/num4simdiffrent/blob/main/CHANGELOG.md
+  documentation_uri: https://rubydoc.info/gems/num4simdiff/0.0.1
+  homepage_uri: http://github.com/siranovel/num4simdiffrent
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.3.7
+signing_key: 
+specification_version: 4
+summary: num for simultaneous different
+test_files: []