num4equ 0.0.7

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 846b376182aec92ffed0f060ef0d57d1d096d12697d4538e2197202b87afc3a5
+  data.tar.gz: bb2ac4a150591ad1e536195a76632a8fce5c986a39e5568f0d5ad3c79b809b59
+SHA512:
+  metadata.gz: 30628042787f597a38cc35524a0e242dee82e8c02d9dbeb802281c02d346be92074b2226244c9faa574df021a53ccbfbb108ff9eba9657b8af8a8225d227781a
+  data.tar.gz: b4ee2e0cdd81b78544a2750598f9c10c987b9e7f11031b4d02ad6b5ee1555de74d8633ac53c516a543772abe9db361041f89ac145847cfe9d237f7b73a6495a6

data/Gemfile

@@ -0,0 +1,2 @@
+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/num4equ/CNum4Equ.c

@@ -0,0 +1,84 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <assert.h>
+#include <math.h>
+#include "CNum4Equ.h"
+
+static double CNum4Equ_doNewtonMethodFFI(double a, Func func, DFunc dFunc, int *ok);
+static double CNum4Equ_doBisectionMethodFFI(double a, double b, Func func, int *ok);
+
+static CNum4Equ _cNum4Equ = {
+    .FP_newtonMethodFFI = CNum4Equ_doNewtonMethodFFI,
+    .FP_bisectionMethodFFI = CNum4Equ_doBisectionMethodFFI,
+};
+
+/**************************************/
+/* InterFface部                       */
+/**************************************/
+/**************************************/
+/* Class部                            */
+/**************************************/
+double CNum4Equ_newtonMethodFFI(double a, Func func, DFunc dFunc, int *ok)
+{
+    assert(func != 0);
+    assert(dFunc != 0);
+
+    return _cNum4Equ.FP_newtonMethodFFI(a, func, dFunc, ok);
+}
+double CNum4Equ_bisectionMethodFFI(double a, double b, Func func, int *ok)
+{
+    assert(func != 0);
+    assert(a < b);
+
+    return _cNum4Equ.FP_bisectionMethodFFI(a, b, func, ok);
+}
+/**************************************/
+/* 処理実行部                         */
+/**************************************/
+static double CNum4Equ_doNewtonMethodFFI(double a, Func func, DFunc dFunc, int *ok)
+{
+    double xn;
+    double x0 = a;
+    long idx = 0;
+
+    *ok = 0;
+    xn = x0;
+    do {
+        double f;
+        double df;
+
+        x0 = xn;
+        f = func(x0);
+        df = dFunc(x0);
+        xn = -1 * f / df + x0;
+        idx++;
+        if (100000 < idx) {
+            *ok = -1;
+            break;
+        } 
+    } while(fabs(xn - x0) > EPS);
+    return xn;
+}
+static double CNum4Equ_doBisectionMethodFFI(double a, double b, Func func, int *ok)
+{
+    double fa = func(a);
+    double fb = func(b);
+    double fxc;
+    double xc;
+
+    *ok = (fa * fb) < 0 ? 0 : -1; 
+    if (0 != *ok) { return 0; }
+    do {
+        xc = (a + b) / 2.0;
+        fxc = func(xc);
+        if (0.0 >= fa * fxc) {
+            b = xc;
+            fb = fxc;
+        } else {
+            a = xc;
+            fa = fxc;
+        }
+    } while(fxc != 0 && fabs(b - a) > EPS);
+    return xc; 
+}
+

data/ext/num4equ/CNum4Equ.h

@@ -0,0 +1,25 @@
+#ifndef _CNum4Equ_H_
+#define _CNum4Equ_H_
+
+/**********************************************/
+/* 構造体宣言                                 */
+/**********************************************/
+typedef struct _CNum4Equ CNum4Equ;
+typedef double (*Func)(double x);
+typedef double (*DFunc)(double x);
+
+struct _CNum4Equ
+{
+    double (*FP_newtonMethodFFI)(double a, Func func, DFunc dFunc, int *ok);
+    double (*FP_bisectionMethodFFI)(double a, double b, Func func, int *ok);
+};
+/**********************************************/
+/* definen宣言                                */
+/**********************************************/
+#define EPS (1.0e-6)
+/**********************************************/
+/* プロトタイプ宣言                           */
+/**********************************************/
+double CNum4Equ_newtonMethodFFI(double a, Func func, DFunc dFunc, int *ok);
+double CNum4Equ_bisectionMethodFFI(double a, double b, Func func, int *ok);
+#endif

data/ext/num4equ/Rakefile

@@ -0,0 +1,5 @@
+require 'ffi-compiler/compile_task'
+
+FFI::Compiler::CompileTask.new('num4equ') do |c|
+    c.ldflags << "-lm"
+end

data/lib/num4equ.rb

@@ -0,0 +1,56 @@
+require 'ffi'
+require 'ffi-compiler/loader'
+
+#
+# 数値計算による方程式の解法ライブラリ
+#
+module Num4EquLib
+    extend FFI::Library
+
+    ffi_lib FFI::Compiler::Loader.find('num4equ') 
+    callback   :f, [:double], :double
+    callback   :df, [:double], :double
+
+    attach_function :newtonMethodFFI,
+        :CNum4Equ_newtonMethodFFI, [:double, :f, :df, :buffer_out], :double
+    attach_function :bisectionMethodFFI,
+        :CNum4Equ_bisectionMethodFFI, [:double, :double, :f, :buffer_out], :double
+    class << self
+        # @overload newtonMethod(a, func, dfunc)
+        #   ニュートン法による解法
+        #   @param [double] a aの値
+        #   @param [callback] func aに対する値を計算
+        #   @param [callback] dfunc fに対する微分関数
+        #   @return [double] xの値
+        #   @raise RangeError
+        def newtonMethod(a, func, dfunc)
+            ok_ptr = FFI::MemoryPointer.new :int 
+            x = newtonMethodFFI(a, func, dfunc, ok_ptr)
+            ok = ok_ptr.read_int
+            ok_ptr.free()
+            if ok < 0 then
+                raise RangeError.new("a:" + a.to_s)
+            end
+            return x
+        end
+        # @overload bisection(a, b, func)
+        #   二分法による解法
+        #   @param [double] a aの値
+        #   @param [double] b bの値
+        #   @param [callback] func aに対する値を計算
+        #   @return [double] xの値
+        #   @raise RangeError
+        def bisectionMethod(a, b, func)
+            ok_ptr = FFI::MemoryPointer.new :int 
+            x = bisectionMethodFFI(a, b, func, ok_ptr)
+            ok = ok_ptr.read_int
+            ok_ptr.free()
+            if ok < 0 then
+                raise RangeError.new("a:" + a.to_s + " " + "b:" + b.to_s)
+            end
+            return x
+        end
+
+        private :newtonMethodFFI, :bisectionMethodFFI
+    end
+end

metadata

@@ -0,0 +1,89 @@
+--- !ruby/object:Gem::Specification
+name: num4equ
+version: !ruby/object:Gem::Version
+  version: 0.0.7
+platform: ruby
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2023-03-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 equation
+email: siranovel@gmail.com
+executables: []
+extensions:
+- ext/num4equ/Rakefile
+extra_rdoc_files: []
+files:
+- Gemfile
+- LICENSE
+- ext/num4equ/CNum4Equ.c
+- ext/num4equ/CNum4Equ.h
+- ext/num4equ/Rakefile
+- lib/num4equ.rb
+homepage: http://github.com/siranovel/num4equation
+licenses:
+- MIT
+metadata: {}
+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 equation
+test_files: []