pspline 5.0.5 → 5.1.0

data/ext/pspline/example/exfspline2.rb

@@ -1,68 +1,68 @@
-#! /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 = XY.fft_complex_forward
-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:3
-vv = fs.axis
-
-s = bs.plot(vv, 10)
-
-puts "# Inverse data"
-
-st = fs.inverse
-printf "["
-st.each {|x| printf "% .1f ", x[0]}
-printf "]\n["
-st.each {|x| printf "% .1f ", x[1]}
-puts "]"
-
-# 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][0]}
-		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
-			ds.with = "lines"
-			ds.title = "Real"
-		end
-		z = vv.map {|v| v[1][1]}
-		plot.data << Gnuplot::DataSet.new( [x, z] ) do |ds|
-			ds.with = "lines"
-			ds.title = "Image"
-		end
-		p = vv.map {|v| a = v[1][0]; b = v[1][1]; Math.sqrt(a*a+b*b)}
-		plot.data << Gnuplot::DataSet.new( [x, p] ) do |ds|
-			ds.with = "lines"
-			ds.linewidth = 2
-			ds.title = "Power"
-		end
-	end
-end
+#! /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 = XY.fft_complex_forward
+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:3
+vv = fs.axis
+
+s = bs.plot(vv, 10)
+
+puts "# Inverse data"
+
+st = fs.inverse
+printf "["
+st.each {|x| printf "% .1f ", x[0]}
+printf "]\n["
+st.each {|x| printf "% .1f ", x[1]}
+puts "]"
+
+# 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][0]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Real"
+		end
+		z = vv.map {|v| v[1][1]}
+		plot.data << Gnuplot::DataSet.new( [x, z] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Image"
+		end
+		p = vv.map {|v| a = v[1][0]; b = v[1][1]; Math.sqrt(a*a+b*b)}
+		plot.data << Gnuplot::DataSet.new( [x, p] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.title = "Power"
+		end
+	end
+end

data/ext/pspline/example/exfspline3.rb

@@ -1,64 +1,64 @@
-#! /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 = X.fft_real_forward
-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, 10)
-
-st = fs.inverse
-printf "["
-st.each {|x| printf "% .1f ", x}
-puts "]"
-
-# 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][0]}
-		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
-			ds.with = "lines"
-			ds.title = "Real"
-		end
-		z = vv.map {|v| v[1][1]}
-		plot.data << Gnuplot::DataSet.new( [x, z] ) do |ds|
-			ds.with = "lines"
-			ds.title = "Image"
-		end
-		p = vv.map {|v| a = v[1][0]; b = v[1][1]; Math.sqrt(a*a+b*b)}
-		plot.data << Gnuplot::DataSet.new( [x, p] ) do |ds|
-			ds.with = "lines"
-			ds.linewidth = 2
-			ds.title = "Power"
-		end
-	end
-end
+#! /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 = X.fft_real_forward
+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, 10)
+
+st = fs.inverse
+printf "["
+st.each {|x| printf "% .1f ", x}
+puts "]"
+
+# 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][0]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Real"
+		end
+		z = vv.map {|v| v[1][1]}
+		plot.data << Gnuplot::DataSet.new( [x, z] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Image"
+		end
+		p = vv.map {|v| a = v[1][0]; b = v[1][1]; Math.sqrt(a*a+b*b)}
+		plot.data << Gnuplot::DataSet.new( [x, p] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.title = "Power"
+		end
+	end
+end

data/ext/pspline/example/exmspline.rb

@@ -1,68 +1,68 @@
-#! /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]]
-
-Bd = [[-4.0, -4.0],[-3.6, -3.6],[-3.2, -3.2],[-2.8, -2.8],
-	  [-2.4, -2.4],[-2.0, -2.0],[-1.6, -1.6],[-1.2, -1.2],
-	  [-0.8, -0.8],[-0.4, -0.4],[ 0.4,  0.4],[ 0.8,  0.8],
-	  [ 1.2,  1.2],[ 1.6,  1.6],[ 2.0,  2.0],[ 2.4,  2.4],
-	  [ 2.8,  2.8],[ 3.2,  3.2],[ 3.6,  3.6],[ 4.0,  4.0]]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-As = Bspline.new(Ad, 20, 5)
-Bs = Bspline.new(Bd, 20, 2)
-Ps = As * Bs
-
-puts "# Interpolation of multiple Bspline 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 = Ps.plot(vv, Dp, Jbn) do |a,b|
-	c = As.sekibun(a[0])
-	d = Ps.sekibun(a[0])
-	printf "% .2f % f % f % f\n", a[0], b[0], c, d
-end
-# STDERR.puts s
-require "gnuplot"
- 
-Gnuplot.open do |gp|
-	Gnuplot::Plot.new( gp ) do |plot|
-		plot.title  'Multiple Bspline'
-		plot.ylabel 'Y'
-		plot.xlabel 'X'
- 		x = vv.map {|u,v| u[0]}
-		y = vv.map {|u,v| v[0]}
-		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
-			ds.with = "lines"
-			ds.linewidth = 2
-			ds.notitle
-		end
-		y = vv.map {|u,v| Ps[[u[0], 0.5*u[0]]]}
-		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
-			ds.with = "lines"
-			ds.linewidth = 2
-			ds.notitle
-		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]]
+
+Bd = [[-4.0, -4.0],[-3.6, -3.6],[-3.2, -3.2],[-2.8, -2.8],
+	  [-2.4, -2.4],[-2.0, -2.0],[-1.6, -1.6],[-1.2, -1.2],
+	  [-0.8, -0.8],[-0.4, -0.4],[ 0.4,  0.4],[ 0.8,  0.8],
+	  [ 1.2,  1.2],[ 1.6,  1.6],[ 2.0,  2.0],[ 2.4,  2.4],
+	  [ 2.8,  2.8],[ 3.2,  3.2],[ 3.6,  3.6],[ 4.0,  4.0]]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+As = Bspline.new(Ad, 20, 5)
+Bs = Bspline.new(Bd, 20, 2)
+Ps = As * Bs
+
+puts "# Interpolation of multiple Bspline 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 = Ps.plot(vv, Dp, Jbn) do |a,b|
+	c = As.sekibun(a[0])
+	d = Ps.sekibun(a[0])
+	printf "% .2f % f % f % f\n", a[0], b[0], c, d
+end
+# STDERR.puts s
+require "gnuplot"
+ 
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new( gp ) do |plot|
+		plot.title  'Multiple Bspline'
+		plot.ylabel 'Y'
+		plot.xlabel 'X'
+ 		x = vv.map {|u,v| u[0]}
+		y = vv.map {|u,v| v[0]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.notitle
+		end
+		y = vv.map {|u,v| Ps[[u[0], 0.5*u[0]]]}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.linewidth = 2
+			ds.notitle
+		end
+	end
+end

data/ext/pspline/example/expspline.rb

@@ -1,29 +1,29 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# Parametric spline interpolation"
-
-X = [0,1,2,3,4,5,6,7]
-Y = [ [ 1.0, 0.0], [ 0.0, 1.0], [-1.0, 0.0], [ 0.0,-1.0],
-	  [ 2.0, 0.0], [ 0.0, 2.0], [-2.0, 0.0], [ 0.0,-2.0] ]
-Jbn = ARGV[0].to_i
-Dp = 20
-
-Ps = Pspline.new(Y, [X], [8], [2], [0])
-
-Y.each do |p|
-	printf "% .2f % .2f\n", p[0], p[1]
-end
-printf "# value of interpolation points, Dp = %d", Dp
-if Jbn == 0
-	print "\n"
-else
-	printf ", Jbn = %d\n", Jbn
-end
-vv = [X]
-s = Ps.plot(vv, Dp, [Jbn]) do |a, b|
-	c = Ps.sekibun(a[0])
-	printf "% .2f % .7f % .7f % .7f % .7f\n", a[0], b[0], b[1], c[0], c[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# Parametric spline interpolation"
+
+X = [0,1,2,3,4,5,6,7]
+Y = [ [ 1.0, 0.0], [ 0.0, 1.0], [-1.0, 0.0], [ 0.0,-1.0],
+	  [ 2.0, 0.0], [ 0.0, 2.0], [-2.0, 0.0], [ 0.0,-2.0] ]
+Jbn = ARGV[0].to_i
+Dp = 20
+
+Ps = Pspline.new(Y, [X], [8], [2], [0])
+
+Y.each do |p|
+	printf "% .2f % .2f\n", p[0], p[1]
+end
+printf "# value of interpolation points, Dp = %d", Dp
+if Jbn == 0
+	print "\n"
+else
+	printf ", Jbn = %d\n", Jbn
+end
+vv = [X]
+s = Ps.plot(vv, Dp, [Jbn]) do |a, b|
+	c = Ps.sekibun(a[0])
+	printf "% .2f % .7f % .7f % .7f % .7f\n", a[0], b[0], b[1], c[0], c[1]
+end
+#STDERR.puts s