pspline 5.0.4 → 5.0.5

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: dd676129ee15c3a24f0aba8101aa3c4e2f0a6ac1
-  data.tar.gz: 28632b037254bcef8dee32470ee0ab513864e21f
+  metadata.gz: 8dacedf68daad7cf3e2b1e3f8aa1e59ccce2408d
+  data.tar.gz: 948d586d9d991821a7e0d640ad9ef63fd98fc074
 SHA512:
-  metadata.gz: b029fa56366d6557ae1d97fc2f9f4124b393b9c35afe0ccdbca1c8da5e51aacb64da8f00ca3f7446bf857eee0c61a068f9bce55d2d8ea56ba961c56d6ecea19d
-  data.tar.gz: f5f54d6eabcfb247a31f3d91bf415b84cc4834dccb6c151e1dd90bbd41846639e1adf4279a9e79fd4176d710f407acf1e2e03ef1bff846253dc1b4892d002986
+  metadata.gz: c61c3af1bfaed26c57d94e144ad9e6546a5e7a21e2b194d3a703f90640b2beafc87450e4752b1399cc452ecba6cf41e154dee27359912ef5977ba396946dd5d5
+  data.tar.gz: 386924860c51388e0ddbd74fc6e95dca67dd771556c13b2d58500e3b833bb476003cd40b539bd496862012f11371e9ecd4e8a1999cb05a4bea2baaba376244bb

data/Gemfile

@@ -1,5 +1,5 @@
-source 'https://rubygems.org'
-
-# Specify your gem's dependencies in pspline.gemspec
-gemspec
-gem "rake-compiler"
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in pspline.gemspec
+gemspec
+gem "rake-compiler"

data/README.md

@@ -29,6 +29,7 @@ See pspline.wiki at https://github.com/vsp2old/pspline .
 5.0.2	Upload for release.   
 5.0.3	Add class PSPLINE::Rfft, class PSPLINE::Cfft.   
 5.0.4	Add method line_integral.   
+5.0.5	Add Multiple Bspline.   
 
 ## Development
 

data/Rakefile

@@ -1,6 +1,6 @@
-require "bundler/gem_tasks"
-require "rake/extensiontask"
-
-Rake::ExtensionTask.new "pspline" do |ext|
-  ext.lib_dir = "lib/pspline"
-end
+require "bundler/gem_tasks"
+require "rake/extensiontask"
+
+Rake::ExtensionTask.new "pspline" do |ext|
+  ext.lib_dir = "lib/pspline"
+end

data/bin/console

@@ -1,14 +1,14 @@
-#!/usr/bin/env ruby
-
-require "bundler/setup"
-require "pspline"
-
-# You can add fixtures and/or initialization code here to make experimenting
-# with your gem easier. You can also use a different console, if you like.
-
-# (If you use this, don't forget to add pry to your Gemfile!)
-# require "pry"
-# Pry.start
-
-require "irb"
-IRB.start
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "pspline"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start

data/bin/setup

@@ -1,8 +1,8 @@
-#!/usr/bin/env bash
-set -euo pipefail
-IFS=$'\n\t'
-set -vx
-
-bundle install
-
-# Do any other automated setup that you need to do here
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/ext/pspline/example/exbspline.rb

@@ -1,36 +1,57 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-Ad = [[-4.0, 1.34095e-3],[-3.6, 2.98189e-3],[-3.2, 6.62420e-3],[-2.8, 1.46827e-2],
-	  [-2.4, 3.23838e-2],[-2.0, 7.06508e-2],[-1.6, 1.50527e-1],[-1.2, 3.05020e-1],
-	  [-0.8, 5.59055e-1],[-0.4, 8.55639e-1],[ 0.4, 8.55639e-1],[ 0.8, 5.59055e-1],
-	  [ 1.2, 3.05020e-1],[ 1.6, 1.50527e-1],[ 2.0, 7.06508e-2],[ 2.4, 3.23838e-2],
-	  [ 2.8, 1.46827e-2],[ 3.2, 6.62420e-3],[ 3.6, 2.98189e-3],[ 4.0, 1.34095e-3]]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Bs = Bspline.new(Ad, 20, 5)
-
-puts "# Interpolation of the sec2(x) function"
-
-vv = []
-Ad.each do |p|
-	printf "% .2f % f\n", p[0], p[1]
-	vv.push p[0];
-end
-printf "# value of interpolation points, Dp = %d", Dp
-if Jbn == 0
-	print "\n"
-elsif Jbn > 0
-	printf ", Jbn = %d\n", Jbn
-elsif Jbn < 0
-	printf ", Jsk = %d\n", -Jbn
-end
-
-s = Bs.plot(vv, Dp, Jbn) do |a,b|
-	printf "% .2f % f\n", a, b[0]
-end
-# STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+Ad = [[-4.0, 1.34095e-3],[-3.6, 2.98189e-3],[-3.2, 6.62420e-3],[-2.8, 1.46827e-2],
+	  [-2.4, 3.23838e-2],[-2.0, 7.06508e-2],[-1.6, 1.50527e-1],[-1.2, 3.05020e-1],
+	  [-0.8, 5.59055e-1],[-0.4, 8.55639e-1],[ 0.4, 8.55639e-1],[ 0.8, 5.59055e-1],
+	  [ 1.2, 3.05020e-1],[ 1.6, 1.50527e-1],[ 2.0, 7.06508e-2],[ 2.4, 3.23838e-2],
+	  [ 2.8, 1.46827e-2],[ 3.2, 6.62420e-3],[ 3.6, 2.98189e-3],[ 4.0, 1.34095e-3]]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Bs = Bspline.new(Ad, 20, 5)
+
+puts "# Interpolation of the sec2(x) function"
+
+vv = []
+Ad.each do |p|
+	printf "% .2f % f\n", p[0], p[1]
+	vv.push p[0];
+end
+printf "# value of interpolation points, Dp = %d", Dp
+if Jbn == 0
+	print "\n"
+elsif Jbn > 0
+	printf ", Jbn = %d\n", Jbn
+elsif Jbn < 0
+	printf ", Jsk = %d\n", -Jbn
+end
+
+s = Bs.plot(vv, Dp, Jbn) do |a,b|
+	printf "% .2f % f\n", a, b[0]
+end
+# Draw Graph
+require "gnuplot"
+ 
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new( gp ) do |plot|
+		plot.title  'sec2(x)'
+		plot.ylabel 'Y'
+		plot.xlabel 'X'
+ 		x = vv.map {|v| v[0]}
+		y = vv.map {|v| v[1]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.notitle
+		end
+		y = x.map {|v| Bs.sekibun(v)}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Integral"
+		end
+	end
+end

data/ext/pspline/example/excspline.rb

@@ -1,36 +1,57 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-n = 7
-j = 5
-
-#
-#puts "# Interpolation of the Bessel function"
-#
-#	Data points
-Ad = [ [0.0, 1.0    ],[0.8, 0.84629],[1.6, 0.45540],[2.0, 0.22389],
-	   [2.4, 0.00251],[3.2,-0.32019],[4.0,-0.39715], ]
-#	Additional data points
-C = [ [0.4, 0.96040],[1.2, 0.67113],[2.8,-0.18504],[3.6,-0.39177] ]
-Jbn = ARGV[0].to_i
-Dp = 10
-
-Bs = Bspline.new(Ad + C, n, j)
-
-vv = []
-for i in 0..6
-	(p, q) = Ad[i]
-	printf "%.2f, % f\n", p, q
-	vv.push p
-end
-printf "# value of interpolation points, Dp = %d", Dp
-if Jbn == 0
-	print "\n"
-else
-	printf ", Jbn = %d\n", Jbn
-end
-s = Bs.plot(vv, Dp, Jbn) do |u, v|
-	printf "% .2f % f\n", u, v[0]
-end
-#puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+n = 7
+j = 5
+
+#
+#puts "# Interpolation of the Bessel function"
+#
+#	Data points
+Ad = [ [0.0, 1.0    ],[0.8, 0.84629],[1.6, 0.45540],[2.0, 0.22389],
+	   [2.4, 0.00251],[3.2,-0.32019],[4.0,-0.39715], ]
+#	Additional data points
+C = [ [0.4, 0.96040],[1.2, 0.67113],[2.8,-0.18504],[3.6,-0.39177] ]
+Jbn = ARGV[0].to_i
+Dp = 10
+
+Bs = Bspline.new(Ad + C, n, j)
+
+vv = []
+for i in 0..6
+	(p, q) = Ad[i]
+	printf "%.2f, % f\n", p, q
+	vv.push p
+end
+printf "# value of interpolation points, Dp = %d", Dp
+if Jbn == 0
+	print "\n"
+else
+	printf ", Jbn = %d\n", Jbn
+end
+s = Bs.plot(vv, Dp, Jbn) do |u, v|
+	printf "% .2f % f\n", u, v[0]
+end
+# Draw Graph
+require "gnuplot"
+
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new( gp ) do |plot|
+		plot.title  'Bessel'
+		plot.ylabel 'Y'
+		plot.xlabel 'X'
+		x = vv.map {|v| v[0]}
+		y = vv.map {|v| v[1]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.notitle
+		end
+		y = x.map {|v| Bs.sekibun(v)}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Integral"
+		end
+	end
+end

data/ext/pspline/example/exdspline.rb

@@ -1,33 +1,55 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# Interpolation of the sin(x) function"
-
-XY = [ [0.0],[1.0],[2.0],[3.0],[4.0],[5.0],[6.283185] ]
-
-XY.each do |a|
-	a.push Math.sin(a[0])
-end
-
-C = [ [0.0, 1.0],[0.0, 0.0],[6.283185, 1.0],[6.283185, 0.0] ]
-D = [1, 2, 1, 2]
-Jbn = ARGV[0].to_i
-Dp = 10
-
-Bs = Bspline.new(XY + C, 7, 5, D);
-
-vv = []
-XY.each do |p|
-	printf "% f, % f\n", p[0], p[1]
-	vv.push p[0]
-end
-printf "# value of interpolation points, Dp = %d", Dp
-if Jbn == 0
-	print "\n"
-else
-	printf ", Jbn = %d\n", Jbn
-end
-s = Bs.plot(vv, 10, Jbn) do |a, b|
-	printf "% f % f\n", a, b[0]
-end
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# Interpolation of the sin(x) function"
+
+XY = [ [0.0],[1.0],[2.0],[3.0],[4.0],[5.0],[6.283185] ]
+
+XY.each do |a|
+	a.push Math.sin(a[0])
+end
+
+C = [ [0.0, 1.0],[0.0, 0.0],[6.283185, 1.0],[6.283185, 0.0] ]
+D = [1, 2, 1, 2]
+Jbn = ARGV[0].to_i
+Dp = 10
+
+Bs = Bspline.new(XY + C, 7, 5, D);
+
+vv = []
+XY.each do |p|
+	printf "% f, % f\n", p[0], p[1]
+	vv.push p[0]
+end
+printf "# value of interpolation points, Dp = %d", Dp
+if Jbn == 0
+	print "\n"
+else
+	printf ", Jbn = %d\n", Jbn
+end
+s = Bs.plot(vv, 10, Jbn) do |a, b|
+	printf "% f % f\n", a, b[0]
+end
+# Draw Graph
+require "gnuplot"
+
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new( gp ) do |plot|
+		plot.title  'sin(x)'
+		plot.ylabel 'Y'
+		plot.xlabel 'X'
+		x = vv.map {|v| v[0]}
+		y = vv.map {|v| v[1]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.notitle
+		end
+		y = x.map {|v| Bs.sekibun(v) - 1}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Integral"
+		end
+	end
+end

data/ext/pspline/example/exfspline.rb

@@ -1,44 +1,44 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# complex FFT data points"
-
-XY = [ [0, 1],[0, 1],[0, 1],[1, 0,],[1, 0,],[1, 0],[1, 0],[0, 1],[0, 1],[0, 1] ]
-
-fs = Cfft.new(XY)
-fs.real.each {|x| printf("% .2f ", x) }
-puts
-fs.imag.each {|y| printf("% .2f ", y) }
-puts
-
-printf "Real :"
--5.upto(5) {|t|
-	w = fs[t]
-	printf "% .2f ", w[0]
-}
-puts
-printf "Imag :"
--5.upto(5) {|t|
-	w = fs[t]
-	printf "% .2f ", w[1]
-}
-puts
-
-bs = fs.spline(1)
-vv = fs.axis
-
-puts "# Interpolation points"
-
-s = bs.plot(vv, 4) do |a, b|
-	printf "% .2f %10.7f %10.7f\n", a, b[0], b[1]
-end
-#STDERR.puts s
-
-st = fs.inverse
-printf "["
-st.each {|x| printf "% .1f ", x[0]}
-printf "]\n["
-st.each {|x| printf "% .1f ", x[1]}
-puts "]"
-
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# complex FFT data points"
+
+XY = [ [0, 1],[0, 1],[0, 1],[1, 0,],[1, 0,],[1, 0],[1, 0],[0, 1],[0, 1],[0, 1] ]
+
+fs = Cfft.new(XY)
+fs.real.each {|x| printf("% .2f ", x) }
+puts
+fs.image.each {|y| printf("% .2f ", y) }
+puts
+
+printf " Real:"
+-5.upto(5) {|t|
+	w = fs[t]
+	printf "% .2f ", w[0]
+}
+puts
+printf "Image:"
+-5.upto(5) {|t|
+	w = fs[t]
+	printf "% .2f ", w[1]
+}
+puts
+
+bs = fs.spline order:1
+vv = fs.axis
+
+puts "# Interpolation points"
+
+s = bs.plot(vv, 4) do |a, b|
+	printf "% .2f %10.7f %10.7f\n", a, b[0], b[1]
+end
+#STDERR.puts s
+
+st = fs.inverse
+printf "["
+st.each {|x| printf "% .1f ", x[0]}
+printf "]\n["
+st.each {|x| printf "% .1f ", x[1]}
+puts "]"
+