pspline 5.0.5 → 5.1.0

data/ext/pspline/example/exbspline.rb

@@ -1,57 +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
-# 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
+#! /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 % .15f\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,57 +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
-# 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
+#! /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, % .5f\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 % .15f\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,55 +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
-# 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
+#! /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 % .15f\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.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 "]"
-
+#! /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 "]"
+

data/ext/pspline/example/exfspline1.rb

@@ -1,40 +1,40 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-printf "# real FFT data points\n"
-
-X = [ 0, 0, 0, 1, 1, 1, 1, 0, 0, 0 ]
-
-fs = Rfft.new(X);
-fs.each {|x| printf("% .2f ", x) }
-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:3
-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}
-puts "]"
-
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+printf "# real FFT data points\n"
+
+X = [ 0, 0, 0, 1, 1, 1, 1, 0, 0, 0 ]
+
+fs = Rfft.new(X);
+fs.each {|x| printf("% .2f ", x) }
+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:3
+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}
+puts "]"
+