equationoftime 4.1.7 → 4.1.8

data/lib/eot/constants.rb

@@ -1,13 +1,17 @@
-# class Eot file = constants.rb
+##
+# class Eot file = constants.rb:
 # has the constants used and more
+
 class Eot
-  # Array result for time_julian_century default = [0.0, 0.0, 0.0, 0.0, 0.0]
-  # A2000       = [0.0, 0.0, 0.0, 0.0, 0.0]
 
+  ##
   # Arc seconds in a degree = 3_600.0
+
   ARCSEC      = 3_600.0
 
+  ##
   # Arc seconds in a degree = 3_600.0
+
   ASD         = 3_600.0
 
   # Arc seconds in an hour = 240.0
@@ -25,31 +29,46 @@ class Eot
   # 2Pi = 6.283185307179586476925287
   # D2PI        = 6.283185307179586476925287
 
+  ##
   # from desktop calculator DAS2R = 4.8481368110953599358991410235795e-6
+
   DAS2R       = 4.8481368110953599358991410235795e-6
 
   # Astronomical unit (m) = 149597870e3
   # DAU         = 149597870e3
 
+  ##  
   # Hours in a day = 24.0
+
   DAY_HOURS   = 24.0
 
+  ##  
   # Minutes in a day = 1_440.0
+
   DAY_MINUTES = 1_440.0
 
+  ##
   # Seconds in a day = 86_400.0
+
   DAY_SECONDS = 86_400.0
 
+  ##
   # Seconds in a day = 86_400.0
+
   DAYSEC      = 86_400.0
 
+  ##
   # Micro Seconds in a day = 86_400_000_000.0
+
   DAY_USECS   = 86_400_000_000.0
 
+  ##
   # Speed of light (AU per day) = DAYSEC / AULT
   # DC = DAYSEC / AULT
 
+  ##
   # from desktop calculator D2R = 0.017453292519943295769236907684886
+
   D2R          = 0.017453292519943295769236907684886
 
   # dint(A) - truncate to nearest whole number towards zero (double)
@@ -61,11 +80,15 @@ class Eot
   # dsign(A,B) - magnitude of A with sign of B (double)
   # dsign(A,B) = ((B)<0.0?-fabs(A):fabs(A))
 
+  ##
   # Reference epoch (J2000.0), Julian Date
   # Default Julian Number = 2451545.0
+
   DJ00        = 2_451_545.0
 
+  ##
   # Days per Julian century = 36525.0
+
   DJC         = 36_525.0
 
   # Days per Julian millennium  = 365250.0
@@ -94,7 +117,9 @@ class Eot
   # Seconds of time to radians = 7.272205216643039903848712e-5
   # DS2R        = 7.272205216643039903848712e-5
 
+  ##
   # Default DateTime = DateTime.new( 2000, 01, 01, 12, 00, 00, "+00:00" )
+
   DT2000      = DateTime.new(2000, 01, 01, 12, 00, 00, '+00:00')
 
   # arc seconds degrees to radians = PI / 180.0 / ARCSEC
@@ -124,20 +149,32 @@ class Eot
   # 1858, 11, 17, 0.0 midnight start of calendar reform = 2400000.5
   # MJD0        = 2400000.5
 
+  ##
   # from desktop calculator PI = 3.1415926535897932384626433832795
+
   PI          = 3.1415926535897932384626433832795
 
+  ##
   # 2Pi from Math module = Math::PI * 2.0
+
   P2          = PI * 2.0
 
+  ##
   # from desktop calculator R2D = 57.295779513082320876798154814105
+
   R2D         = 57.295779513082320876798154814105
 
+  ##
   # from desktop calculator RTD = 0.015915494309189533576888376337251
+
   RTD         = 0.015915494309189533576888376337251
 
+  ##
   # from desktop calculator Sidereal minutes =
   # 4.0 / 1.0027379093507953456536618754278
+  # probably 3.989 would be close enough because
+  # it's not that constant anyway 
+
   SM          = 4.0 / 1.0027379093507953456536618754278
 
   # Schwarzschild radius of the Sun (au) =

data/lib/eot/times.rb

@@ -1,9 +1,17 @@
-# class Eot file = times.rb
+##
+# class Eot file = times.rb:
 # methods calculating times
+
 class Eot
+
+  ##
   # From times.rb:
+
   # Pass in an AJD number
   # Returns a DateTime object
+  # only DateTime#jd() to work with so
+  # add a half day to make it work.
+
   def ajd_to_datetime(ajd)
     DateTime.jd(ajd + 0.5)
   end

data/lib/eot/version.rb

@@ -2,5 +2,5 @@
 # contains gem version
 class Eot
   # current version
-  VERSION = '4.1.7'
+  VERSION = '4.1.8'
 end

data/test/eot/aliased_angles_spec.rb

@@ -0,0 +1,293 @@
+# aliased_angles_spec.rb
+gem 'minitest'
+require 'minitest/autorun'
+
+lib = File.expand_path('../../../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'eot'
+
+aliased_angles = Eot.new
+
+describe 'tests ajd of 2456885.0 ' do
+
+  before(:each) do
+    aliased_angles.ajd  = 2_456_885.0
+    ajd = aliased_angles.ajd
+    aliased_angles.ma_ta_set
+    # somtimes need date to check values somewhere else
+    aliased_angles.date = aliased_angles.ajd_to_datetime(ajd)
+  end
+
+  it 'expected   2_456_885.0 for aliased_angles.ajd 'do
+    assert_equal(2_456_885.0, aliased_angles.ajd)
+  end
+
+  it 'expected   "2014-08-15T12:00:00+00:00" from \
+      aliased_angles.date.to_s ' do
+    assert_equal('2014-08-15T12:00:00+00:00', aliased_angles.date.to_s)
+  end
+
+  it 'expected   3.8508003966038915 from aliased_angles.ma 'do
+    assert_equal(3.8508003966038915, aliased_angles.ma)
+  end
+
+  it 'expected   2.4887103398436143 from \
+      aliased_angles.apparent_longitude()? ' do
+    assert_equal(2.4887103398436143, aliased_angles.apparent_longitude)
+  end
+
+  it 'expected   -0.7943361570447028 from \
+      aliased_angles.cosine_apparent_longitude()? ' do
+    assert_equal(-0.7943361570447028, \
+                 aliased_angles.cosine_apparent_longitude)
+  end
+
+  it 'expected   -0.7943772759574919 from \
+      aliased_angles.cosine_true_longitude()? ' do
+    assert_equal(-0.7943772759574919, aliased_angles.cosine_true_longitude)
+  end
+
+  it 'expected   0.9175115346811911 from \
+      aliased_angles.cosine_true_obliquity()? ' do
+    assert_equal(0.9175115346811911, aliased_angles.cosine_true_obliquity)
+  end
+
+  it 'expected   0.24401410218543554 from \
+      aliased_angles.declination()? ' do
+    assert_equal(0.24401410218543554, aliased_angles.declination)
+  end
+
+  it 'expected   -0.04103082558803539 from \
+      aliased_angles.delta_t_ecliptic()? ' do
+    assert_equal(-0.04103082558803539, aliased_angles.delta_t_ecliptic)
+  end
+
+  it 'expected   0.021413249720702462 from \
+      aliased_angles.delta_t_elliptic()? ' do
+    assert_equal(0.021413249720702462, aliased_angles.delta_t_elliptic)
+  end
+
+  it 'expected   0.016702468499021204 from \
+      aliased_angles.eccentricity_earth_orbit()? ' do
+    assert_equal(0.016702468499021204, \
+                 aliased_angles.eccentricity_earth_orbit)
+  end
+
+  it 'expected   -0.021413249720702462 from \
+      aliased_angles.equation_of_center()? ' do
+    assert_equal(-0.021413249720702462, aliased_angles.equation_of_center)
+  end
+
+  it 'expected   2.5101912804141424 from \
+      aliased_angles.geometric_mean_longitude()? ' do
+    assert_equal(2.5101912804141424, \
+                 aliased_angles.geometric_mean_longitude)
+  end
+
+  it 'expected   1.5857841877939605 from \
+      aliased_angles.horizon_angle(1)? ' do
+    assert_equal(1.5857841877939605, aliased_angles.horizon_angle(1))
+  end
+
+  it 'expected   3.8508003966038915 from \
+      aliased_angles.mean_anomaly()? ' do
+    assert_equal(3.8508003966038915, aliased_angles.mean_anomaly)
+  end
+
+  it 'expected   2.510089864980358 from \
+      aliased_angles.mean_longitude_aries()? ' do
+    assert_equal(2.510089864980358, aliased_angles.mean_longitude_aries)
+  end
+
+  it 'expected   0.40905940254265843 from \
+      aliased_angles.mean_obliquity()? ' do
+    assert_equal(0.40905940254265843, aliased_angles.mean_obliquity)
+  end
+
+  it 'expected   0.40905940254265843 from \
+      aliased_angles.mean_obliquity_of_ecliptic()? ' do
+    assert_equal(0.40905940254265843, \
+                 aliased_angles.mean_obliquity_of_ecliptic)
+  end
+
+  it 'expected   0.40901870461547685 from \
+      aliased_angles.obliquity_correction()? ' do
+    assert_equal(0.40901870461547685, aliased_angles.obliquity_correction)
+  end
+
+  it 'expected   2.5297411654316497 from \
+      aliased_angles.right_ascension()? ' do
+    assert_equal(2.5297411654316497, aliased_angles.right_ascension)
+  end
+
+  it 'expected   0.6074784519729512 from \
+      aliased_angles.sine_apparent_longitude()? ' do
+    assert_equal(0.6074784519729512, aliased_angles.sine_apparent_longitude)
+  end
+
+  it 'expected   0.6074246812917259 from \
+      aliased_angles.sine_true_longitude()? ' do
+    assert_equal(0.6074246812917259, aliased_angles.sine_true_longitude)
+  end
+
+  it 'expected   3.8293871468831893 from \
+      aliased_angles.true_anomaly()? ' do
+    assert_equal(3.8293871468831893, aliased_angles.true_anomaly)
+  end
+
+  it 'expected   2.48877803069344 from \
+      aliased_angles.true_longitude()? ' do
+    assert_equal(2.48877803069344, aliased_angles.true_longitude)
+  end
+
+  it 'expected   2.5101242776531474 from \
+      aliased_angles.true_longitude_aries()? ' do
+    assert_equal(2.5101242776531474, aliased_angles.true_longitude_aries)
+  end
+
+  it 'expected   0.40901870461547685 from \
+      aliased_angles.true_obliquity()? ' do
+    assert_equal(0.40901870461547685, aliased_angles.true_obliquity)
+  end
+end
+
+describe 'tests ajd of 2455055.5 ' do
+
+  before(:each) do
+    aliased_angles.ajd                   = 2_455_055.0
+    ajd = aliased_angles.ajd
+    aliased_angles.ma_ta_set
+    # check date for this ajd when needed.
+    aliased_angles.date = aliased_angles.ajd_to_datetime(ajd)
+  end
+
+  it 'expected   2_455_055.0, from aliased_angles.' do
+    assert_equal(2_455_055.0, aliased_angles.ajd)
+  end
+
+  it 'expected   "2009-08-11T12:00:00+00:00" from \
+      aliased_angles.date.to_s ' do
+    assert_equal('2009-08-11T12:00:00+00:00', aliased_angles.date.to_s)
+  end
+
+  it 'expected   3.7871218188949207, from aliased_angles.ma ' do
+    assert_equal(3.7871218188949207, aliased_angles.ma)
+  end
+
+  it 'expected   3.7871218188949207 from \
+      aliased_angles.ma  from Eot_angles.mean_anomaly() ' do
+    assert_equal(3.7871218188949207, aliased_angles.mean_anomaly)
+  end
+
+  it 'expected   2.4252140645725033 from \
+      aliased_angles.apparent_longitude()? ' do
+    assert_equal(2.4252140645725033, aliased_angles.apparent_longitude)
+  end
+
+  it 'expected   -0.7541886969975007 from \
+      aliased_angles.cosine_apparent_longitude()? ' do
+    assert_equal(-0.7541886969975007, \
+                 aliased_angles.cosine_apparent_longitude)
+  end
+
+  it 'expected   -0.7542060769936684 from \
+      aliased_angles.cosine_true_longitude()? ' do
+    assert_equal(-0.7542060769936684, aliased_angles.cosine_true_longitude)
+  end
+
+  it 'expected   0.9174818088112336 from \
+      aliased_angles.cosine_true_obliquity()? ' do
+    assert_equal(0.9174818088112336, aliased_angles.cosine_true_obliquity)
+  end
+
+  it 'expected   0.2642691272294404 from \
+      aliased_angles.declination()? ' do
+    assert_equal(0.2642691272294404, aliased_angles.declination)
+  end
+
+  it 'expected   -0.04234904897476355 from \
+      aliased_angles.delta_t_ecliptic()? ' do
+    assert_equal(-0.04234904897476355, aliased_angles.delta_t_ecliptic)
+  end
+
+  it 'expected    0.019768413456709915 from \
+      aliased_angles.delta_t_elliptic()? ' do
+    assert_equal(0.019768413456709915, aliased_angles.delta_t_elliptic)
+  end
+
+  it 'expected   0.016704576164208475 from \
+      aliased_angles.eccentricity_earth_orbit()? ' do
+    assert_equal(0.016704576164208475, \
+                 aliased_angles.eccentricity_earth_orbit)
+  end
+
+  it 'expected   -0.019768413456709915 from \
+      aliased_angles.equation_of_center()? ' do
+    assert_equal(-0.019768413456709915, aliased_angles.equation_of_center)
+  end
+
+  it 'expected   2.445008945789877 from \
+      aliased_angles.geometric_mean_longitude()? ' do
+    assert_equal(2.445008945789877, aliased_angles.geometric_mean_longitude)
+  end
+
+  it 'expected   1.585863261753274 from \
+      aliased_angles.horizon_angle()? ' do
+    assert_equal(1.585863261753274, aliased_angles.horizon_angle(1))
+  end
+
+  it 'expected   3.7871218188949207 from \
+      aliased_angles.mean_anomaly()? ' do
+    assert_equal(3.7871218188949207, aliased_angles.mean_anomaly)
+  end
+
+  it 'expected   2.444907382260759 from \
+      aliased_angles.mean_longitude_aries()? ' do
+    assert_equal(2.444907382260759, aliased_angles.mean_longitude_aries)
+  end
+
+  it 'expected   0.4090707793981491 from \
+      aliased_angles.mean_obliquity()? ' do
+    assert_equal(0.4090707793981491, aliased_angles.mean_obliquity)
+  end
+
+  it 'expected   0.4090934409048494 from \
+      aliased_angles.obliquity_correction()? ' do
+    assert_equal(0.4090934409048494, aliased_angles.obliquity_correction)
+  end
+
+  it 'expected   2.467563113547267 from \
+      aliased_angles.right_ascension()? ' do
+    assert_equal(2.467563113547267, aliased_angles.right_ascension)
+  end
+
+  it 'expected   0.6566577566139093 from \
+      aliased_angles.sine_apparent_longitude()? ' do
+    assert_equal(0.6566577566139093, aliased_angles.sine_apparent_longitude)
+  end
+
+  it 'expected   0.6566377946979757 from \
+      aliased_angles.sine_true_longitude()? ' do
+    assert_equal(0.6566377946979757, aliased_angles.sine_true_longitude)
+  end
+
+  it 'expected   3.767353405438211 from aliased_angles.true_anomaly()? ' do
+    assert_equal(3.767353405438211, aliased_angles.true_anomaly)
+  end
+
+  it 'expected   2.4252405323331674 from \
+      aliased_angles.true_longitude()? ' do
+    assert_equal(2.4252405323331674, aliased_angles.true_longitude)
+  end
+
+  it 'expected   2.4449774607872907 from \
+      aliased_angles.true_longitude_aries()? ' do
+    assert_equal(2.4449774607872907, aliased_angles.true_longitude_aries)
+  end
+
+  it 'expected   0.4090934409048494 from \
+      aliased_angles.true_obliquity()? ' do
+    assert_equal(0.4090934409048494, aliased_angles.true_obliquity)
+  end
+
+end

data/test/eot/aliased_displays_spec.rb

@@ -0,0 +1,127 @@
+# aliased_displays_spec.rb
+gem 'minitest'
+require 'minitest/autorun'
+
+lib = File.expand_path('../../../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'eot'
+
+aliased_displays = Eot.new
+
+describe 'aliased_displays using ajd of 2456885.0' do
+
+  before(:each) do
+    aliased_displays.ajd  = 2_456_885.0
+    ajd = aliased_displays.ajd
+    aliased_displays.ma_ta_set
+    # check date for this ajd when needed.
+    aliased_displays.date = aliased_displays.ajd_to_datetime(ajd)
+
+  end
+
+  it 'expected   2456885.0 , from aliased_displays.' do
+    assert_equal 2_456_885.0, aliased_displays.ajd
+  end
+
+  it 'expected   "2014-08-15T12:00:00+00:00", \
+      from aliased_displays.date.to_s' do
+    assert_equal '2014-08-15T12:00:00+00:00', aliased_displays.date.to_s
+  end
+
+  it 'expected   3.8508003966038915, from aliased_displays.' do
+    assert_equal 3.8508003966038915, aliased_displays.ma
+  end
+
+  it 'expected   "+142:35:33.356" \
+      from aliased_displays.apparent_longitude_string()? ' do
+    assert_equal '+142:35:33.356', aliased_displays.apparent_longitude_string
+  end
+
+  it 'expected   "+013:58:51.522" \
+      from aliased_displays.declination_string()? ' do
+    assert_equal '+013:58:51.522', aliased_displays.declination_string
+  end
+
+  it 'expected   "-04m, 29.2s" from \
+      aliased_displays.display_equation_of_time()? ' do
+    assert_equal '-04m, 29.2s', aliased_displays.display_equation_of_time
+  end
+
+  it 'expected   "12:00:00.000" \
+      from aliased_displays.display_time_string()? ' do
+    assert_equal '12:00:00.000', aliased_displays.display_time_string
+    assert_equal '12:00:00.000', aliased_displays.display_time_string(nil)
+    assert_equal '12:00:00.000', aliased_displays.display_time_string(0)
+  end
+
+  it 'expected   "12:00:00.000" from \
+      aliased_displays.display_time_string(aliased_displays.date)? ' do
+    assert_equal '12:00:00.000', \
+                 aliased_displays.display_time_string(aliased_displays.date)
+  end
+
+  it 'expected   "2000-01-01" from aliased_displays.jd_to_date_string()? ' do
+    assert_equal '2000-01-01', aliased_displays.jd_to_date_string
+    assert_equal '2000-01-01', aliased_displays.jd_to_date_string(nil)
+    assert_equal '2000-01-01', aliased_displays.jd_to_date_string(0)
+  end
+
+  it 'expected   "2014-08-15" \
+      from aliased_displays.jd_to_date_string(aliased_displays.ajd)? ' do
+    assert_equal '2014-08-15', \
+                 aliased_displays.jd_to_date_string(aliased_displays.ajd)
+  end
+
+  it 'expected   "12:00:00" from \
+      aliased_displays.julian_period_day_fraction_to_time()? ' do
+    assert_equal '12:00:00', \
+                 aliased_displays.julian_period_day_fraction_to_time
+    assert_equal '12:00:00', \
+                 aliased_displays.julian_period_day_fraction_to_time(nil)
+    assert_equal '12:00:00', \
+                 aliased_displays.julian_period_day_fraction_to_time(0)
+  end
+
+  it 'expected   "+220:38:04.598" from \
+      aliased_displays.mean_anomaly_string()? ' do
+    assert_equal '+220:38:04.598', aliased_displays.mean_anomaly_string
+  end
+
+  it 'expected   "+144:56:36.571" from \
+      aliased_displays.right_ascension_string()? ' do
+    assert_equal '+144:56:36.571', aliased_displays.right_ascension_string
+  end
+
+  it 'expected   "+219:24:27.798" from \
+      aliased_displays.true_anomaly_string()? ' do
+    assert_equal '+219:24:27.798', aliased_displays.true_anomaly_string
+  end
+
+  it 'expected   "+142:35:47.318" from \
+      aliased_displays.true_longitude_string()? ' do
+    assert_equal '+142:35:47.318', aliased_displays.true_longitude_string
+  end
+
+  it 'expected   "+023:26:06.164" from \
+      aliased_displays.true_obliquity_string()? ' do
+    assert_equal '+023:26:06.164', aliased_displays.true_obliquity_string
+  end
+
+end
+
+describe 'Eot aliased displays explicit values' do
+
+  it 'expected   "16:40:40.800" from \
+      aliased_displays.display_time_string(16.6780)? ' do
+    assert_equal '16:40:40.800', \
+                 aliased_displays.display_time_string(16.6780)
+  end
+
+  it 'expected   "17:59:16.800" from \
+      aliased_displays.display_time_string(17988)? ' do
+    aliased_displays.date = Date.today.to_s
+    assert_equal '17:59:16.800', \
+                 aliased_displays.display_time_string(17.988)
+  end
+
+end