num4partialdiff 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 16d60c15e1c6da491e4b8485ecfa5024fcb4a5985ddb05f8b8e65999332ff889
+  data.tar.gz: a50c7000c50c87f3471a2e8dbd3a31dadfd45b471fd682b4b13831b48b9073fa
+SHA512:
+  metadata.gz: a1498e064591a3865a5500f4fcb69421b72bbcd46f1cba7ba27bdd55082b6da1fed44eb1b9a8cf98659e83bfd222f87b797d42e2733e5b6610dcfda12579b99b
+  data.tar.gz: 5658eebcbc846cb902a2a32e1b048b1ae32ed0031df36d0c802cf665894898b8c59f71b204a55a61d144b076a4f61676ab981633e2ceaf08f15c3716c2cae6a1

data/.yardopts

@@ -0,0 +1 @@
+--no-private

data/CHANGELOG.md

@@ -0,0 +1,9 @@
+# Change Log
+
+## Unreleased
+
+## [0.0.1] - 2023-08-12
+
+### 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) 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/num4partialdiff/CNum4PartialDiff.c

@@ -0,0 +1,42 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <assert.h>
+#include "CNum4PartialDiff.h"
+
+static double* CNum4PartialDiff_doFiniteDiffMethod(int n, double *yi, double h, double *f);
+static CNum4PartialDiff _cNum4PDE = {
+    .FP_finiteDiffMethod = CNum4PartialDiff_doFiniteDiffMethod,
+};
+/**************************************/
+/* InterFface部                       */
+/**************************************/
+/**************************************/
+/* Class部                            */
+/**************************************/
+double* CNum4PartialDiff_finiteDiffMethod(int n, double *yi, double h, double *f)
+{
+    assert(f != 0);
+    assert(yi != 0);
+    assert( n > 0);
+
+    return _cNum4PDE.FP_finiteDiffMethod(n, yi, h, f);
+}
+/**************************************/
+/* 処理実行部                         */
+/**************************************/
+/*
+ * 前進差分法
+ */
+static double* CNum4PartialDiff_doFiniteDiffMethod(int n, double *yi, double h, double *f)
+{
+    int i;
+    double *yi_1 = malloc(sizeof(double) * n);
+
+    assert(yi_1 != 0);
+    // yi_1 = yi + h * f(xi, y)
+    for (i = 0; i < n; i++) {
+        yi_1[i] = yi[i] + h * f[i];
+    }
+    return yi_1;
+}
+

data/ext/num4partialdiff/CNum4PartialDiff.h

@@ -0,0 +1,20 @@
+#ifndef _CNum4PartialDiff_H_
+#define _CNum4PartialDiff_H_
+
+/**************************************/
+/* 構造体宣言                         */
+/**************************************/
+typedef struct _CNum4PartialDiff CNum4PartialDiff;
+
+struct _CNum4PartialDiff
+{
+    double* (*FP_finiteDiffMethod)(int n, double *yi, double h, double *f);
+};
+/**************************************/
+/* define宣言                         */
+/**************************************/
+/**************************************/
+/* プロトタイプ宣言                   */
+/**************************************/
+double* CNum4PartialDiff_finiteDiffMethod(int n, double *yi, double h, double *f);
+#endif

data/ext/num4partialdiff/Rakefile

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

data/lib/num4partialdiff.rb

@@ -0,0 +1,91 @@
+require 'ffi'
+require 'ffi-compiler/loader'
+require 'fiddle'
+
+
+#
+# 数値計算によ一次偏微分方程式の解法するライブラリ
+module Num4PartialDiffLib
+    extend FFI::Library
+
+    ffi_lib FFI::Compiler::Loader.find('num4partialdiff')
+    attach_function :finiteDiffMethodFFI,
+        :CNum4PartialDiff_finiteDiffMethod,[:int, :buffer_in, :double, :buffer_in], :pointer
+    class << self
+        #
+        # @private
+        class Counter
+            @idx = 0
+            def self.getIdx
+                return @idx
+            end
+            def self.countup
+               @idx = @idx + 1
+            end
+        end
+        # 有限差分法による数値計算
+        #
+        # @overload finiteDiffMethod(yi, x0, h, func)
+        #   @param [double[]] yi 各変数の値
+        #   @param [double[]] x0 各変数の初期値
+        #   @param [double] h  刻み幅
+        #   @param [callback] func xiに対する傾きを計算する関数
+        #   @return [double[]] x0+nhに対するfi_1の値(配列)
+        # @example
+        #   yi = [0.0, 0.0, 0.0]
+        #   x0 = [0.0, 0.0, 0.0]
+        #   h = 0.001
+        #   func = Proc.new do | n, xi |
+        #     f = 2 * xi[0] + 3 * xi[1] + 4 * xi[2] + 5
+        #     next f
+        #   end
+        #   yi_1 = Num4PartialDiffLib.finiteDiffMethod(yi, x0, h, func)
+        # @note 各変数の刻み幅は、一定です。
+        # 
+        def finiteDiffMethod(yi, x0, h, func)
+            n = yi.size
+            fi_ptr = calcVectFi(n, x0, h, func)
+            yi_ptr = cnvRbAry2pt(n, yi)
+            yi_1_ptr = finiteDiffMethodFFI(n, yi_ptr, h, fi_ptr)
+            yi_1 = cnvPt2RbAry(n, yi_1_ptr)
+            return yi_1
+        end
+        #
+        # @private
+        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
+        #
+        # @private
+        def cnvPt2RbAry(n, pt)
+          rbAry = n.times.map { |i|
+             pt.get_double(i * Fiddle::SIZEOF_DOUBLE)
+          }
+          return rbAry  
+        end
+        #
+        # @private
+        def calcVectFi(n, x0, h, func)
+            xi = []
+            fi_ptr = FFI::MemoryPointer.new(:double, n)
+            idx = Counter.getIdx()
+            n.times.map { |i|
+                xi = x0
+                xi[i] = x0[i] + idx * h
+                fi_ptr.put_double(i * Fiddle::SIZEOF_DOUBLE, func.call(n, xi))
+            }
+            Counter.countup
+            return fi_ptr
+        end
+
+        private :cnvRbAry2pt
+        private :cnvPt2RbAry
+        private :calcVectFi
+        private :finiteDiffMethodFFI
+    end
+end
+

metadata

@@ -0,0 +1,95 @@
+--- !ruby/object:Gem::Specification
+name: num4partialdiff
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2023-08-13 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 partial differential equations
+email: siranovel@gmail.com
+executables: []
+extensions:
+- ext/num4partialdiff/Rakefile
+extra_rdoc_files: []
+files:
+- ".yardopts"
+- CHANGELOG.md
+- Gemfile
+- LICENSE
+- ext/num4partialdiff/CNum4PartialDiff.c
+- ext/num4partialdiff/CNum4PartialDiff.h
+- ext/num4partialdiff/Rakefile
+- lib/num4partialdiff.rb
+homepage: http://github.com/siranovel/num4partialdifferent
+licenses:
+- MIT
+metadata:
+  changelog_uri: http://github.com/siranovel/num4partialdifferent/blob/main/CHANGELOG.md
+  documentation_uri: https://rubydoc.info/gems/num4partialdiff/0.0.1
+  homepage_uri: http://github.com/siranovel/num4partialdifferent
+  wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4partidifferent
+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 partial different
+test_files: []