pspline 5.0.4 → 5.0.5

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.real.each {|x| printf("% .2f ", x) }
-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(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 "]"
+

data/ext/pspline/example/exfspline2.rb

@@ -0,0 +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

data/ext/pspline/example/exfspline3.rb

@@ -0,0 +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

data/ext/pspline/example/exmspline.rb

@@ -0,0 +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

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

data/ext/pspline/example/expspline1.rb

@@ -1,29 +1,29 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# Parametric spline interpolation"
-
-XY = [ [0, 1.0, 0.0], [1, 0.0, 1.0], [2,-1.0, 0.0], [3, 0.0,-1.0],
-	   [4, 2.0, 0.0], [5, 0.0, 2.0], [6,-2.0, 0.0], [7, 0.0,-2.0] ]
-Jbn = ARGV[0].to_i
-Dp = 20
-
-Qs = Bspline.new(XY, 8, 2)
-
-vv = []
-XY.each do |p|
-	printf "<% .1f % .1f>\n", p[1], p[2]
-	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 = Qs.plot(vv, Dp, Jbn) do |a, b|
-	c = Qs.sekibun(a);
-	printf "% .2f % .7f % .7f % .7f % .7f\n", a, b[0], b[1], c[0], c[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# Parametric spline interpolation"
+
+XY = [ [0, 1.0, 0.0], [1, 0.0, 1.0], [2,-1.0, 0.0], [3, 0.0,-1.0],
+	   [4, 2.0, 0.0], [5, 0.0, 2.0], [6,-2.0, 0.0], [7, 0.0,-2.0] ]
+Jbn = ARGV[0].to_i
+Dp = 20
+
+Qs = Bspline.new(XY, 8, 2)
+
+vv = []
+XY.each do |p|
+	printf "<% .1f % .1f>\n", p[1], p[2]
+	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 = Qs.plot(vv, Dp, Jbn) do |a, b|
+	c = Qs.sekibun(a);
+	printf "% .2f % .7f % .7f % .7f % .7f\n", a, b[0], b[1], c[0], c[1]
+end
+#STDERR.puts s

data/ext/pspline/example/expspline2.rb

@@ -1,47 +1,47 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-# Parametric interpolation
-
-nn = [5, 3]
-jj = [3, 2]
-ss = [0, 0]
-
-puts "# 2 variable 3 dimension parametric spline interpolation"
-
-S = [0, 1, 2, 3, 4]
-T = [0, 1, 2]
-
-XYZ = [ [[ 2.0,-1.0,-2.0], [ 2.0,-2.0,-1.0], [ 2.0,-2.0,-3.0]],
-		[[ 1.0, 0.0,-1.0], [ 1.0,-1.0, 0.0], [ 1.0,-1.0,-2.0]],
-		[[ 0.0, 2.0, 0.0], [ 0.0, 0.0, 2.0], [ 0.0, 0.0,-2.0]],
-		[[-1.0, 2.0, 1.0], [-1.0, 1.0, 2.0], [-1.0, 1.0, 0.0]],
-		[[-2.0, 3.0, 2.0], [-2.0, 2.0, 3.0], [-2.0, 2.0, 1.0]] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Ps = Pspline.new(XYZ, [S, T], nn, jj, ss)
-
-for i in 0..4
-	for j in 0..2
-		printf "%d,%d", S[i], T[j]
-		u = XYZ[i][j]
-		printf " <% .1f % .1f % .1f>\n", u[0], u[1], u[2]
-	end
-end
-printf "# value of interpolation points, Dp = %d", Dp
-if Jbn == 0
-	print "\n"
-else
-	printf ", Jbn = %d\n", Jbn
-end
-vv = [S, T]
-N = Ps.plot(vv, Dp, [Jbn,0]) do |a, b|
-	if (a[0] == 0.1)
-		printf "% .2f % .2f % f % f % f\n", a[0], a[1], b[0], b[1], b[2]
-	end
-end
-#STDERR.puts N
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+# Parametric interpolation
+
+nn = [5, 3]
+jj = [3, 2]
+ss = [0, 0]
+
+puts "# 2 variable 3 dimension parametric spline interpolation"
+
+S = [0, 1, 2, 3, 4]
+T = [0, 1, 2]
+
+XYZ = [ [[ 2.0,-1.0,-2.0], [ 2.0,-2.0,-1.0], [ 2.0,-2.0,-3.0]],
+		[[ 1.0, 0.0,-1.0], [ 1.0,-1.0, 0.0], [ 1.0,-1.0,-2.0]],
+		[[ 0.0, 2.0, 0.0], [ 0.0, 0.0, 2.0], [ 0.0, 0.0,-2.0]],
+		[[-1.0, 2.0, 1.0], [-1.0, 1.0, 2.0], [-1.0, 1.0, 0.0]],
+		[[-2.0, 3.0, 2.0], [-2.0, 2.0, 3.0], [-2.0, 2.0, 1.0]] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Ps = Pspline.new(XYZ, [S, T], nn, jj, ss)
+
+for i in 0..4
+	for j in 0..2
+		printf "%d,%d", S[i], T[j]
+		u = XYZ[i][j]
+		printf " <% .1f % .1f % .1f>\n", u[0], u[1], u[2]
+	end
+end
+printf "# value of interpolation points, Dp = %d", Dp
+if Jbn == 0
+	print "\n"
+else
+	printf ", Jbn = %d\n", Jbn
+end
+vv = [S, T]
+N = Ps.plot(vv, Dp, [Jbn,0]) do |a, b|
+	if (a[0] == 0.1)
+		printf "% .2f % .2f % f % f % f\n", a[0], a[1], b[0], b[1], b[2]
+	end
+end
+#STDERR.puts N