suntrack 0.0.3 → 0.0.4

data/test/test_point3d.rb

@@ -0,0 +1,97 @@
+require 'test/unit'
+require 'suntrack'
+
+# Test cases for Suntrack::Point3D.
+
+class Point3DTest < Test::Unit::TestCase
+
+  # Confirm that we get expected results for basic transformations
+  def test_ecliptic_Cartesian_1
+
+    # Ecliptic point: radius, ecliptic latitude, ecliptic longitude
+    # By setting the ecliptic latitude to zero, the point will be in the
+    # xy plane
+    pt = Suntrack::Point3D.new(1,0,30)
+    pt.polar_to_cartesian!
+    assert_equal pt.x.round(6), 0.866025
+    assert_equal pt.y.round(6), 0.5
+    assert_equal pt.z.round(6), 0.0
+  end
+
+  # Show the effect of precessing an ecliptic point between two epochs.
+  def test_precessed_Cartesian_ecliptic
+
+    # Spherical/Ecliptic/1950
+    pt = Suntrack::Point3D.new(1,0,30)
+    pt.polar_to_cartesian!
+    # Cartesian/Ecliptic/1950
+    pt.precess_ecliptic_cartesian!(Suntrack::RAstro::B1950_EPOCH,Suntrack::RAstro::J2000_EPOCH)
+    # Cartesian/Ecliptic/2000
+    assert_equal pt.x.round(6), 0.859866
+    assert_equal pt.y.round(6), 0.510519
+    assert_equal pt.z.round(6), 0.000067
+  end
+
+  # Show the effect of precessing an ecliptic point between two epochs, and
+  # show how much the original polar coordinates are affected.
+  def test_precessed_Cartesian_ecliptic_to_polar
+
+    # Spherical/Ecliptic/1950
+    pt = Suntrack::Point3D.new(1,0,30)
+    pt.polar_to_cartesian!
+    # Cartesian/Ecliptic/1950
+    pt.precess_ecliptic_cartesian!(Suntrack::RAstro::B1950_EPOCH,Suntrack::RAstro::J2000_EPOCH)
+    # Cartesian/Ecliptic/2000
+    d = Suntrack::Point3D.new(pt.x,pt.y,pt.z)
+    d.cartesian_to_polar!
+    # Polar/Ecliptic/2000
+    assert_equal d.x.round(6), 1.0
+    assert_equal d.y.round(6), 0.003811
+    assert_equal d.z.round(6), 30.698411
+  end
+
+  # The two approaches should give the same answer.
+  def test_consistency_of_conversions
+
+    # Spherical/Ecliptic/1950
+    pt = Suntrack::Point3D.new(1,0,30)
+    pt.polar_to_cartesian!
+    # Cartesian/Ecliptic/1950
+    pt.precess_ecliptic_cartesian!(Suntrack::RAstro::B1950_EPOCH,Suntrack::RAstro::J2000_EPOCH)
+    # Cartesian/Ecliptic/2000
+    d = Suntrack::Point3D.new(pt.x,pt.y,pt.z)
+    d.cartesian_to_polar!
+    # Polar/Ecliptic/2000
+
+    pt.ecliptic_to_equatorial!(Suntrack::RAstro::J2000_EPOCH)
+    # Cartesian/Equatorial/2000
+    pt.equatorial_to_ecliptic!(Suntrack::RAstro::J2000_EPOCH)
+    # Cartesian/Ecliptic/2000
+    pt.cartesian_to_polar!
+    # Polar/Ecliptic/2000
+    assert_equal pt.x, d.x
+    assert_equal pt.y.round(6), d.y.round(6)
+    assert_equal pt.z, d.z
+  end
+
+  # Montebruck/Pfleger test case on page 31
+  def test_MP_case
+
+    # Equatorial point r = 1, delta = 0, alpha = 0, in B1950
+    # delta = declination, alpha = right ascension
+    pt = Suntrack::Point3D.new(1,0,0)
+    pt.polar_to_cartesian!
+    pt.precess_equatorial_cartesian!(Suntrack::RAstro::B1950_EPOCH,Suntrack::RAstro::J2000_EPOCH)
+    pt.equatorial_to_ecliptic!(Suntrack::RAstro::J2000_EPOCH)
+    pt.cartesian_to_polar!
+    assert_equal pt.x.round(6), 1.0
+    dms = Suntrack::RAstro.get_dms(pt.y)
+    assert_equal dms.x, 0
+    assert_equal dms.y, 0
+    assert_equal dms.z.round(6), 2.335464
+    dms = Suntrack::RAstro.get_dms(pt.z)
+    assert_equal dms.x, 0
+    assert_equal dms.y, 41
+    assert_equal dms.z.round(6), 54.280965
+  end
+end

data/test/test_sirius.rb

@@ -0,0 +1,24 @@
+require 'test/unit'
+require 'suntrack'
+
+# Test cases for Suntrack::Point3D and Suntrack::RAstro, to track the
+# altitude/azimuth of Sirius.
+
+class SiriusTest < Test::Unit::TestCase
+  def test_somerset_nj_1
+    x = DateTime.new(2012,1,16,0,0,0)
+    # y component is elevation
+    # z component is azimuth
+    r = Suntrack.sirius_location(x,40.5,74.5)
+    assert_equal r.y.round(3),9.598
+    assert_equal r.z.round(3),238.104
+  end
+  def test_somerset_nj_2
+    x = DateTime.new(2012,1,16,4,0,0)
+    # y component is elevation
+    # z component is azimuth
+    r = Suntrack.sirius_location(x,40.5,74.5)
+    assert_equal r.y.round(3),32.765
+    assert_equal r.z.round(3),0.885
+  end
+end

data/test/test_sun_points.rb

@@ -0,0 +1,64 @@
+require 'test/unit'
+require 'suntrack'
+require 'date'
+
+# Test cases for Suntrack::Point3D and Suntrack::RAstro, tracking the position
+# of the Sun, in declination and right ascension, as a function of time.
+# These results do not depend on latitude and longitude. For a given moment
+# in time, declination/right ascension of the Sun is a specific point on the
+# sky; latitude and longitude are then used along with an appropriate
+# coordinate system to compute altitude and azimuth.
+
+class SuntrackTest < Test::Unit::TestCase
+
+  def test_point_one
+
+    # August 7th 1997
+    # http://www.stargazing.net/kepler/sun.html
+    pointOne = DateTime.new(1997,8,7,11,0,0)
+    mjd = Suntrack::RAstro.to_mjd(pointOne)
+    assert_equal mjd.round(3), 50667.458
+    jd = Suntrack::RAstro.jd(mjd)
+    assert_equal jd.round(3), 2450667.958
+    tm = Suntrack::RAstro.get_time(jd)
+    assert_equal tm.round(4), -0.0240
+    dec = Suntrack::RAstro.sun_position(tm)
+    # Declination
+    assert_equal dec.y.round(5), 16.34011
+    # Right ascension
+    assert_equal dec.z.round(5), 9.16339
+  end
+
+  def test_point_two
+
+    pointTwo = DateTime.new(2004,1,3,11,0,0)
+    mjd = Suntrack::RAstro.to_mjd(pointTwo)
+    assert_equal mjd.round(3), 53007.458
+    jd = Suntrack::RAstro.jd(mjd)
+    assert_equal jd.round(3), 2453007.958
+    tm = Suntrack::RAstro.get_time(jd)
+    assert_equal tm.round(4), 0.0401
+    dec = Suntrack::RAstro.sun_position(tm)
+    # Declination
+    assert_equal dec.y.round(5), -22.86037
+    # Right ascension
+    assert_equal dec.z.round(5), 18.89909
+  end
+
+  def test_point_three
+
+    pointTwo = DateTime.new(2012,9,7,19,0,0)
+    mjd = Suntrack::RAstro.to_mjd(pointTwo)
+    assert_equal mjd.round(3), 56177.792
+    jd = Suntrack::RAstro.jd(mjd)
+    assert_equal jd.round(3), 2456178.292
+    tm = Suntrack::RAstro.get_time(jd)
+    assert_equal tm.round(4), 0.1269
+    dec = Suntrack::RAstro.sun_position(tm)
+    # Declination
+    assert_equal dec.y.round(5), 5.69534
+    # Right ascension
+    assert_equal dec.z.round(5), 11.11341
+  end
+
+end

data/test/test_suntrack.rb

@@ -0,0 +1,35 @@
+require 'test/unit'
+require 'suntrack'
+
+# Test cases for Suntrack::Point3D and RAstro, to track the position of the
+# Sun as a function of time, in altitude and azimuth.
+
+class SuntrackTest < Test::Unit::TestCase
+
+  def test_somerset_nj
+    x = DateTime.new(2012,9,7,19,0,0)
+    # y component is elevation
+    # z component is azimuth
+    r = Suntrack.sun_location(x,40.5,74.5)
+    assert_equal r.y.round(3),45.453
+    assert_equal r.z.round(3),47.026
+  end
+
+  def test_miami_fl
+    x = DateTime.new(2012,9,7,19,30,0)
+    # y component is elevation
+    # z component is azimuth
+    r = Suntrack.sun_location(x,25.76,80.21)
+    assert_equal r.y.round(3),52.751
+    assert_equal r.z.round(3),63.075
+  end
+
+  def test_san_francisco
+    x = DateTime.new(2012,9,7,19,45,0)
+    # y component is elevation
+    # z component is azimuth
+    r = Suntrack.sun_location(x,37.766,122.42)
+    assert_equal r.y.round(3),57.513
+    assert_equal r.z.round(3),10.446
+  end
+end

data/test/test_user_define.rb

@@ -0,0 +1,18 @@
+require 'test/unit'
+require 'suntrack'
+
+# Test cases for Suntrack::Point3D and Suntrack::RAstro, to track the
+# altitude/azimuth of a user-provided star
+
+class UserProvidedTest < Test::Unit::TestCase
+  def test_vega_somerset_nj_1
+    x = DateTime.new(2012,1,16,0,0,0)
+    # y component is elevation
+    # z component is azimuth
+    Suntrack.vega_declination=38.784
+    Suntrack.vega_ra=18.616
+    r = Suntrack.vega_location(x,40.5,74.5)
+    assert_equal r.y.round(3),5.773
+    assert_equal r.z.round(3),222.132
+  end
+end

metadata

@@ -2,7 +2,7 @@
 name: suntrack
 version: !ruby/object:Gem::Version 
   prerelease: 
-  version: 0.0.3
+  version: 0.0.4
 platform: ruby
 authors: 
 - Joel M. Gottlieb
@@ -33,6 +33,11 @@ files:
 - doc/Suntrack.html
 - doc/class_list.html
 - doc/_index.html
+- test/test_sun_points.rb
+- test/test_suntrack.rb
+- test/test_sirius.rb
+- test/test_user_define.rb
+- test/test_point3d.rb
 homepage: 
 licenses: []
 
@@ -59,7 +64,7 @@ rubyforge_project:
 rubygems_version: 1.8.21
 signing_key: 
 specification_version: 3
-summary: Suntrack
+summary: Suntrack offers methods for obtaining star positions given time, latitude and longitude
 test_files: []
 
 has_rdoc: