extract_curves 0.0.1-i586-linux

data/CVS/Entries

@@ -0,0 +1,4 @@
+/extract_curves.gemspec/1.1.1.1/Tue Nov 15 17:32:08 2005//
+D/bin////
+D/ruby_ext////
+D/ruby_libs////

data/CVS/Repository

@@ -0,0 +1 @@
+extract-curves

data/CVS/Root

@@ -0,0 +1 @@
+:ext:pavpen@rubyforge.org:/var/cvs/extract-curves

data/bin/CVS/Entries

@@ -0,0 +1,5 @@
+/ec_rect2polar.rb/1.1.1.1/Tue Nov 15 17:32:11 2005//
+/ec_rev_lines.rb/1.1.1.1/Tue Nov 15 17:32:11 2005//
+/ec_sph_area.rb/1.1.1.1/Tue Nov 15 17:32:11 2005//
+/extract_curves.rb/1.1.1.1/Tue Nov 15 17:32:11 2005//
+D/tests////

data/bin/CVS/Repository

@@ -0,0 +1 @@
+extract-curves/bin

data/bin/CVS/Root

@@ -0,0 +1 @@
+:ext:pavpen@rubyforge.org:/var/cvs/extract-curves

data/bin/ec_rect2polar.rb

@@ -0,0 +1,22 @@
+#!/usr/bin/ruby
+
+require 'complex'
+require 'pav/floatsio'
+
+x0, y0, kr, ktheta, format_str = *ARGV
+unless x0 && y0
+	puts "Usage: #{File.basename($0)} x_0 y_0 [k_r [k_theta [format_str]]]"
+	puts "  (x_0, y_0) = Projection of the zero pole."
+	puts "  k_r, k_theta = Coefficients for the linear (r, theta) pixel " +
+		"coordinates."
+	exit 1
+end
+x0 = x0.to_f; y0 = y0.to_f
+kr ? kr = kr.to_f : kr = 1.0
+ktheta ? ktheta = ktheta.to_f : ktheta = 1.0
+format_str = "%.29g\t%.29g" unless format_str
+
+FloatsInp.new($stdin).each_rec { |rec| x, y = *rec
+	r, theta = *Complex(x-x0, y-y0).polar
+	puts format_str % [theta * ktheta, r * kr]
+}

data/bin/ec_rev_lines.rb

@@ -0,0 +1,5 @@
+#!/usr/bin/ruby -w
+
+lns = $stdin.readlines
+lns.reverse!
+print lns

data/bin/ec_sph_area.rb

@@ -0,0 +1,30 @@
+#!/usr/bin/ruby
+
+require 'pav/floatsio'
+
+format_str = *ARGV
+if ARGV.include?("--usage") || ARGV.include?("--help")
+	puts "Usage: #{File.basename($0)} [format_str]"
+	puts "  Calculates the area of a spherical polygon given by "
+	puts "  (azimuth; angular distance from the pole) points in steradians."
+	exit 1
+end
+format_str = "%.30g" unless format_str
+
+area = 0.0
+EPS = 1e-300
+inp = FloatsInp.new($stdin)
+az0, d0 = inp.get_rec
+inp.each_rec { |rec| az, d = *rec
+	#puts "d0=#{d0.inspect}, az0=#{az0.inspect}, d=#{d.inspect}, az=#{
+	#	az.inspect}"
+	if (d - d0).abs < EPS
+		term = Math.cos( (d+d0)/2.0 )
+	else
+		term = (Math.sin(d) - Math.sin(d0)) / (d - d0)
+	end
+	area += (az - az0) * (1.0 - term)
+	d0 = d; az0 = az
+}
+
+puts format_str % area