pspline 5.0.5 → 5.1.0

data/ext/pspline/example/exqspline4.rb

@@ -1,35 +1,35 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# Parametric period interpolation (Closed curve)"
-
-XY = [ [ 0, 1.0, 0.0],[ 1, 0.923880, 0.382683],[ 2, 0.707107, 0.707107],[ 3, 0.382683, 0.923880],
-	   [ 4, 0.0, 1.0],[ 5,-0.382683, 0.923880],[ 6,-0.707107, 0.707107],[ 7,-0.923880, 0.382683],
-	   [ 8,-1.0, 0.0],[ 9,-0.923880,-0.382683],[10,-0.707107,-0.707107],[11,-0.382683,-0.923880],
-	   [12, 0.0,-1.0],[13, 0.382683,-0.923880],[14, 0.707107,-0.707107],[15, 0.923880,-0.382683],
-	   [16, 1.0, 0.0] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 8
-
-Qs = Bspline.new(XY, 16, 9, 1)
-
-vv = []
-XY.each do |p|
-	printf "% f, % f\n", p[1], p[2]
-	vv.push p[0]
-end
-
-Sint = Qs.line_integral(lambda{|x,y|Math.sqrt(x*x+y*y)*0.5}, vv)
-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|
-	printf "%6.3f % f % f % f\n", a, b[0], b[1], Sint[a]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# Parametric period interpolation (Closed curve)"
+
+XY = [ [ 0, 1.0, 0.0],[ 1, 0.923880, 0.382683],[ 2, 0.707107, 0.707107],[ 3, 0.382683, 0.923880],
+	   [ 4, 0.0, 1.0],[ 5,-0.382683, 0.923880],[ 6,-0.707107, 0.707107],[ 7,-0.923880, 0.382683],
+	   [ 8,-1.0, 0.0],[ 9,-0.923880,-0.382683],[10,-0.707107,-0.707107],[11,-0.382683,-0.923880],
+	   [12, 0.0,-1.0],[13, 0.382683,-0.923880],[14, 0.707107,-0.707107],[15, 0.923880,-0.382683],
+	   [16, 1.0, 0.0] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 8
+
+Qs = Bspline.new(XY, 16, 9, 1)
+
+vv = []
+XY.each do |p|
+	printf "% f, % f\n", p[1], p[2]
+	vv.push p[0]
+end
+
+Sint = Qs.line_integral(lambda{|x,y|Math.sqrt(x*x+y*y)*0.5}, vv)
+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|
+	printf "%6.3f % f % f % f\n", a, b[0], b[1], Sint[a]
+end
+#STDERR.puts s

data/ext/pspline/example/exrspline.rb

@@ -1,34 +1,34 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-n   = 7
-j   = 2
-s   = 2
-
-puts "# Riesenfeld interpolation"
-
-X = [0, 1, 2, 3, 4, 5, 6]
-Y = [ [ 0.0, 0.0],[ 3.0, 0.0],[ 3.0, 1.0],[ 0.5, 2.0],
-	  [ 0.5, 6.0],[ 5.0, 6.0],[ 5.0, 9.0] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Rs = Pspline.new(Y, [X], [n], [j], [s])
-	
-Y.each do |p|
-	printf "% f, % f\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 = Rs.plot(vv, Dp, [Jbn]) do |a, b|
-	printf "% .2f, <% .7f, % .7f>\n", a[0], b[0], b[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+n   = 7
+j   = 2
+s   = 2
+
+puts "# Riesenfeld interpolation"
+
+X = [0, 1, 2, 3, 4, 5, 6]
+Y = [ [ 0.0, 0.0],[ 3.0, 0.0],[ 3.0, 1.0],[ 0.5, 2.0],
+	  [ 0.5, 6.0],[ 5.0, 6.0],[ 5.0, 9.0] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Rs = Pspline.new(Y, [X], [n], [j], [s])
+	
+Y.each do |p|
+	printf "% f, % f\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 = Rs.plot(vv, Dp, [Jbn]) do |a, b|
+	printf "% .2f, <% .7f, % .7f>\n", a[0], b[0], b[1]
+end
+#STDERR.puts s

data/ext/pspline/example/exrspline1.rb

@@ -1,34 +1,34 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-n   = 7
-j   = 2
-s   = 2
-
-puts "# Riesenfeld interpolation"
-
-XY = [ [0, 0.0, 0.0],[1, 3.0, 0.0],[2, 3.0, 1.0],[3, 0.5, 2.0],
-	   [4, 0.5, 6.0],[5, 5.0, 6.0],[6, 5.0, 9.0] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Rs = Bspline.new(XY, n, j, s)
-
-vv = []	
-XY.each do |p|
-	printf "% f, % f\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 = Rs.plot(vv, Dp, Jbn) do |a, b|
-	printf "% .2f % .7f % .7f\n", a, b[0], b[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+n   = 7
+j   = 2
+s   = 2
+
+puts "# Riesenfeld interpolation"
+
+XY = [ [0, 0.0, 0.0],[1, 3.0, 0.0],[2, 3.0, 1.0],[3, 0.5, 2.0],
+	   [4, 0.5, 6.0],[5, 5.0, 6.0],[6, 5.0, 9.0] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Rs = Bspline.new(XY, n, j, s)
+
+vv = []	
+XY.each do |p|
+	printf "% f, % f\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 = Rs.plot(vv, Dp, Jbn) do |a, b|
+	printf "% .2f % .7f % .7f\n", a, b[0], b[1]
+end
+#STDERR.puts s

data/ext/pspline/example/exrspline2.rb

@@ -1,44 +1,44 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-nn = [4, 6]
-jj = [3, 2]
-ss = [3, 2]
-
-puts "# Reisenfeld interpolation"
-
-S = [0, 1, 2, 3, 4]
-T = [0, 1, 2, 3, 4, 5]
-
-XYZ = [ [[ 0.0, 5.0, 9.0],[ 0.0, 5.0, 6.0],[ 0.0, 0.5, 6.0],[ 0.0, 0.5, 2.0],[ 0.0, 3.0, 1.0],[ 0.0, 3.0, 0.0]],
-		[[ 5.0,-2.5, 6.5],[ 5.0,-2.5, 3.5],[ 0.5,-0.25,5.75],[0.5,-0.25,1.75],[3.0,-1.5,-0.5],[ 3.0,-1.5,-1.5]],
-		[[ 0.0,-5.0, 9.0],[ 0.0,-5.0, 6.0],[ 0.0,-0.5, 6.0],[ 0.0,-0.5, 2.0],[ 0.0,-3.0, 1.0],[ 0.0,-3.0, 0.0]],
-		[[-5.0, 2.5,11.5],[-5.0, 2.5, 8.5],[-0.5,0.25,6.25],[-0.5,0.25,2.25],[-3.0, 1.5, 2.5],[-3.0, 1.5, 1.5]] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 8
-
-Rs = Pspline.new(XYZ, [S, T], nn, jj, ss)
-
-for i in 0..3
-	for j in 0..5
-		printf("<%d,%d>", S[i], T[j])
-		u = XYZ[i][j]
-		printf(" % .7f % .7f %10.7f\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 = Rs.plot(vv, Dp, [Jbn,0]) do |a, b|
-	if (a[1] == 0.0)
-		printf("%.3f %.3f % .7f % .7f %10.7f\n", a[0], a[1], b[0], b[1], b[2])
-	end
-end
-#STDERR.puts N
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+nn = [4, 6]
+jj = [3, 2]
+ss = [3, 2]
+
+puts "# Reisenfeld interpolation"
+
+S = [0, 1, 2, 3, 4]
+T = [0, 1, 2, 3, 4, 5]
+
+XYZ = [ [[ 0.0, 5.0, 9.0],[ 0.0, 5.0, 6.0],[ 0.0, 0.5, 6.0],[ 0.0, 0.5, 2.0],[ 0.0, 3.0, 1.0],[ 0.0, 3.0, 0.0]],
+		[[ 5.0,-2.5, 6.5],[ 5.0,-2.5, 3.5],[ 0.5,-0.25,5.75],[0.5,-0.25,1.75],[3.0,-1.5,-0.5],[ 3.0,-1.5,-1.5]],
+		[[ 0.0,-5.0, 9.0],[ 0.0,-5.0, 6.0],[ 0.0,-0.5, 6.0],[ 0.0,-0.5, 2.0],[ 0.0,-3.0, 1.0],[ 0.0,-3.0, 0.0]],
+		[[-5.0, 2.5,11.5],[-5.0, 2.5, 8.5],[-0.5,0.25,6.25],[-0.5,0.25,2.25],[-3.0, 1.5, 2.5],[-3.0, 1.5, 1.5]] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 8
+
+Rs = Pspline.new(XYZ, [S, T], nn, jj, ss)
+
+for i in 0..3
+	for j in 0..5
+		printf("<%d,%d>", S[i], T[j])
+		u = XYZ[i][j]
+		printf(" % .7f % .7f %10.7f\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 = Rs.plot(vv, Dp, [Jbn,0]) do |a, b|
+	if (a[1] == 0.0)
+		printf("%.3f %.3f % .7f % .7f %10.7f\n", a[0], a[1], b[0], b[1], b[2])
+	end
+end
+#STDERR.puts N

data/ext/pspline/example/exsspline.rb

@@ -1,35 +1,35 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-n = 8
-j = 3
-s = 3
-
-puts "# Riesenfeld interpolation （Closed curve）"
-
-X = [0, 1, 2, 3, 4, 5, 6, 7, 8]
-
-Y = [ [ 4.0, 0.0],[ 3.0, 3.0],[ 0.0, 4.0],[-3.0, 3.0],
-	  [-4.0, 0.0],[-3.0,-3.0],[ 0.0,-4.0],[ 3.0,-3.0] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Sp = Pspline.new(Y, [X], [n], [j], [s])
-
-Y.each do |p|
-	printf "% f, % f\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 = Sp.plot(vv, Dp, [Jbn]) do |a, b|
-	printf "% .2f, <% .7f, % .7f>\n", a[0], b[0], b[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+n = 8
+j = 3
+s = 3
+
+puts "# Riesenfeld interpolation （Closed curve）"
+
+X = [0, 1, 2, 3, 4, 5, 6, 7, 8]
+
+Y = [ [ 4.0, 0.0],[ 3.0, 3.0],[ 0.0, 4.0],[-3.0, 3.0],
+	  [-4.0, 0.0],[-3.0,-3.0],[ 0.0,-4.0],[ 3.0,-3.0] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Sp = Pspline.new(Y, [X], [n], [j], [s])
+
+Y.each do |p|
+	printf "% f, % f\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 = Sp.plot(vv, Dp, [Jbn]) do |a, b|
+	printf "% .2f, <% .7f, % .7f>\n", a[0], b[0], b[1]
+end
+#STDERR.puts s

data/ext/pspline/example/exsspline1.rb

@@ -1,35 +1,35 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-puts "# Riesenfeld interpolation （Closed curve interpolation）"
-
-n = 8
-j = 3
-s = 3
-
-XY = [ [0, 4.0, 0.0],[1, 3.0, 3.0],[2, 0.0, 4.0],[3,-3.0, 3.0],
-	   [4,-4.0, 0.0],[5,-3.0,-3.0],[6, 0.0,-4.0],[7, 3.0,-3.0],
-	   [8, 4.0, 0.0] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Sp = Bspline.new(XY, n, j, s)
-
-vv = []
-XY.each do |p|
-	printf "% f, % f\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 = Sp.plot(vv, Dp, Jbn) do |a, b|
-	printf "% .2f % .7f % .7f\n", a, b[0], b[1]
-end
-#STDERR.puts s
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+puts "# Riesenfeld interpolation （Closed curve interpolation）"
+
+n = 8
+j = 3
+s = 3
+
+XY = [ [0, 4.0, 0.0],[1, 3.0, 3.0],[2, 0.0, 4.0],[3,-3.0, 3.0],
+	   [4,-4.0, 0.0],[5,-3.0,-3.0],[6, 0.0,-4.0],[7, 3.0,-3.0],
+	   [8, 4.0, 0.0] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Sp = Bspline.new(XY, n, j, s)
+
+vv = []
+XY.each do |p|
+	printf "% f, % f\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 = Sp.plot(vv, Dp, Jbn) do |a, b|
+	printf "% .2f % .7f % .7f\n", a, b[0], b[1]
+end
+#STDERR.puts s

data/ext/pspline/example/extspline.rb

@@ -1,54 +1,54 @@
-#! /usr/local/bin/ruby
-require 'pspline'
-include PSPLINE
-
-#
-#puts "# Interpolation of the Jacobi function"
-#
-XY = [ [ 0.0,    0.0    ],[ 0.5474, 0.5    ],[ 1.2837, 0.86603],[ 2.7681, 1.0    ],
-	   [ 4.2525, 0.86603],[ 4.9888, 0.5    ],[ 5.5362, 0.0    ],[ 6.0836,-0.5    ],
-	   [ 6.8199,-0.86603],[ 8.3043,-1.0    ],[ 9.7887,-0.86603],[10.5250,-0.5    ],
-	   [11.0724, 0.0    ] ]
-
-Jbn = ARGV[0].to_i
-
-Dp = 10
-
-Bs = Bspline.new(XY, 12, 5, 1)
-
-vv = []
-XY.each do |p|
-	printf "%7.4f, % 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 |u, v|
-	printf "%8.5f % f\n", u, v[0]
-end
-# Draw Graph
-require "gnuplot"
- 
-Gnuplot.open do |gp|
-	Gnuplot::Plot.new( gp ) do |plot|
-		plot.title  'Jacobi'
-		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.value(v, 1)}
-		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
-			ds.with = "lines"
-			ds.title = "Differential"
-		end
-	end
-end
+#! /usr/local/bin/ruby
+require 'pspline'
+include PSPLINE
+
+#
+#puts "# Interpolation of the Jacobi function"
+#
+XY = [ [" 0.0",   " 0.0"    ],[" 0.5474"," 0.5"],[" 1.2837"," 0.86603"],[" 2.7681"," 1.0"],
+	   [" 4.2525"," 0.86603"],[" 4.9888"," 0.5"],[" 5.5362"," 0.0"    ],[" 6.0836","-0.5"],
+	   [" 6.8199","-0.86603"],[" 8.3043","-1.0"],[" 9.7887","-0.86603"],["10.5250","-0.5"],
+	   ["11.0724"," 0.0"    ] ]
+
+Jbn = ARGV[0].to_i
+
+Dp = 10
+
+Bs = Bspline.new(XY, 12, 5, 1)
+
+vv = []
+XY.each do |p|
+	printf "%7.4f, % .5f\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 |u, v|
+	printf "%8.5f % .15f\n", u, v[0]
+end
+# Draw Graph
+require "gnuplot"
+ 
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new( gp ) do |plot|
+		plot.title  'Jacobi'
+		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.value(v, 1)}
+		plot.data << Gnuplot::DataSet.new( [x, y] ) do |ds|
+			ds.with = "lines"
+			ds.title = "Differential"
+		end
+	end
+end