num4inte 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 391c7af590564ae676a706c71d1c3d67ba1d8f3c3d01cece8bc5b298d13de1c9
+  data.tar.gz: fae05e502bb9acc4a4456fb285a9cae64d43d895e8822e21b68a7b7b31be7bdd
+SHA512:
+  metadata.gz: 19add3204882ae70fe3173567f13e1f268da893df1a46f73dd58833191f5ec194d9200cc4769527a4a74995f5a0257a794acf4c1b27692e31763c5254ff4240d
+  data.tar.gz: 9d81760ab9526a4ddc621ba374f1e30ae35d339903f9c301e3aee15bf64b857768044d74f09893011f11e81c38b526fd6f9fe75356475372bcd10f843abd0799

data/CHANGELOG.md

@@ -0,0 +1,9 @@
+# Change Log
+
+## Unreleased
+
+## [0.0.1] - 2023-04-25
+
+### 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/num4inte/CNum4Inte.c

@@ -0,0 +1,86 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <assert.h>
+#include "CNum4Inte.h"
+
+static double CNum4Inte_doLeftReimannSumMetod(double a, double b, double h, Func func);
+static double CNum4Inte_doRigtReimannSumMetod(double a, double b, double h, Func func);
+static double CNum4Inte_doNewtonMethod(double a, double b, double h, Func func);
+static CNum4Inte _cNum4Inte = {
+    .FP_leftReimannSumMetod = CNum4Inte_doLeftReimannSumMetod,
+    .FP_rigtReimannSumMetod = CNum4Inte_doRigtReimannSumMetod,
+    .FP_newtonMethod        = CNum4Inte_doNewtonMethod,
+};
+/**************************************/
+/* InterFface部                       */
+/**************************************/
+/**************************************/
+/* Class部                            */
+/**************************************/
+double CNum4Inte_leftReimannSumMetod(double a, double b, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Inte.FP_leftReimannSumMetod(a, b, h, func);
+}
+double CNum4Inte_rigtReimannSumMetod(double a, double b, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Inte.FP_rigtReimannSumMetod(a, b, h, func);
+}
+double CNum4Inte_newtonMethod(double a, double b, double h, Func func)
+{
+    assert(func != 0);
+
+    return _cNum4Inte.FP_newtonMethod(a, b, h, func);
+}
+/**************************************/
+/* 処理実行部                          */
+/**************************************/
+/*
+ * 左リーマン和法
+ */
+static double CNum4Inte_doLeftReimannSumMetod(double a, double b, double h, Func func)
+{
+    double rimann = 0.0;
+    double x = 0.0;
+
+    for (x = a; x < b; x += h) {
+        rimann += func(x);
+    }
+    return rimann * h;
+}
+/*
+ * 右リーマン和法
+ */
+static double CNum4Inte_doRigtReimannSumMetod(double a, double b, double h, Func func)
+{
+    double rimann = 0.0;
+    double x = 0.0;
+
+    for (x = a; x < b; x += h) {
+        rimann += func(x + h);
+    }
+    return rimann * h;
+}
+/*
+ * ニュートン・シンプソン法
+ */
+static double CNum4Inte_doNewtonMethod(double a, double b, double h, Func func)
+{
+    double newton = 0.0;
+    double x = 0.0;
+    double fx = func(a);
+    double fxh = 0.0;
+
+    for (x = a; x < b; x += h) {
+        fxh = func(x + h);
+        newton += (fx + fxh) /2;
+        fx = fxh;
+    }    
+    return newton * h;
+}
+
+
+

data/ext/num4inte/CNum4Inte.h

@@ -0,0 +1,25 @@
+#ifndef _CNum4Inte_H_
+#define _CNum4Inte_H_
+
+/**************************************/
+/* 構造体宣言                         */
+/**************************************/
+typedef struct _CNum4Inte CNum4Inte;
+typedef double (*Func)(double x);
+
+struct _CNum4Inte
+{
+    double (*FP_leftReimannSumMetod)(double a, double b, double h, Func func);
+    double (*FP_rigtReimannSumMetod)(double a, double b, double h, Func func);
+    double (*FP_newtonMethod)(double a, double b, double h, Func func);
+};
+/**************************************/
+/* define宣言                         */
+/**************************************/
+/**************************************/
+/* プロトタイプ宣言                   */
+/**************************************/
+double CNum4Inte_leftReimannSumMetod(double a, double b, double h, Func func);
+double CNum4Inte_rigtReimannSumMetod(double a, double b, double h, Func func);
+double CNum4Inte_newtonMethod(double a, double b, double h, Func func);
+#endif

data/ext/num4inte/Rakefile

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

data/lib/num4inte.rb

@@ -0,0 +1,55 @@
+require 'ffi'
+require 'ffi-compiler/loader'
+
+#
+# 数値計算による数値積分の解法するライブラリ
+#
+module Num4InteLib
+    extend FFI::Library
+
+    ffi_lib FFI::Compiler::Loader.find('num4inte') 
+    # @overload f(xi)
+    #   dy = f(xi)
+    #   @yield [x] dy = f(x)
+    #   @yieldparam [double] xi xiの値
+    #   @return [double] xiに対するyの値
+    callback   :f, [:double], :double
+
+    #
+    # 左リーマン和法
+    # @overload leftReimannSumMetod(a, b, h, func)
+    #   y = leftReimannSumMetod(a, b, h, func)
+    #   @param [double] a aの値
+    #   @param [double] b bの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func xiに対する傾きを計算する関数
+    #   @return [double] [a,b]の間の積分値
+    #
+    attach_function :leftReimannSumMetod,
+        :CNum4Inte_leftReimannSumMetod, [:double, :double, :double, :f], :double
+    #
+    # 右リーマン和法
+    # @overload rigtReimannSumMetod(a, b, h, func)
+    #   y = rigtReimannSumMetod(a, b, h, func)
+    #   @param [double] a aの値
+    #   @param [double] b bの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func xiに対する傾きを計算する関数
+    #   @return [double] [a,b]の間の積分値
+    #
+    attach_function :rigtReimannSumMetod,
+        :CNum4Inte_rigtReimannSumMetod, [:double, :double, :double, :f], :double
+    #
+    # ニュートン・シンプソン法
+    # @overload newtonMethod(a, b, h, func)
+    #   y = newtonMethod(a, b, h, func)
+    #   @param [double] a aの値
+    #   @param [double] b bの値
+    #   @param [double] h  刻み幅
+    #   @param [callback] func xiに対する傾きを計算する関数
+    #   @return [double] [a,b]の間の積分値
+    #
+    attach_function :newtonMethod,
+        :CNum4Inte_newtonMethod, [:double, :double, :double, :f], :double
+end
+

metadata

@@ -0,0 +1,93 @@
+--- !ruby/object:Gem::Specification
+name: num4inte
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2023-04-25 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 integral
+email: siranovel@gmail.com
+executables: []
+extensions:
+- ext/num4inte/Rakefile
+extra_rdoc_files: []
+files:
+- CHANGELOG.md
+- Gemfile
+- LICENSE
+- ext/num4inte/CNum4Inte.c
+- ext/num4inte/CNum4Inte.h
+- ext/num4inte/Rakefile
+- lib/num4inte.rb
+homepage: http://github.com/siranovel/num4integral
+licenses:
+- MIT
+metadata:
+  changelog_uri: http://github.com/siranovel/num4integral/blob/main/CHANGELOG.md
+  documentation_uri: https://rubydoc.info/gems/num4inte/0.0.1
+  homepage_uri: http://github.com/siranovel/num4integral
+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 integral
+test_files: []