marsdate 1.1.7 → 2.0.1

data/test/timescale_test.rb

@@ -0,0 +1,124 @@
+#!/usr/bin/env ruby
+# Tests for MarsDateTime::TimeScale (NASA/Allison / Mars24).
+# Run: ruby test/timescale_test.rb
+
+require 'minitest/autorun'
+require_relative '../lib/marsdate/timescale'
+require 'date'
+
+class TimeScaleTest < Minitest::Test
+  TS = MarsDateTime::TimeScale
+
+  def test_tai_utc_table
+    assert_in_delta 10.0, TS.tai_utc(Date.new(1972, 1, 1)), 1e-9
+    assert_in_delta 32.0, TS.tai_utc(Date.new(1999, 6, 1)), 1e-9
+    assert_in_delta 32.0, TS.tai_utc(Date.new(2000, 1, 6)), 1e-9
+    assert_in_delta 37.0, TS.tai_utc(Date.new(2017, 1, 1)), 1e-9
+    assert_in_delta 37.0, TS.tai_utc(Date.new(2026, 8, 13)), 1e-9
+    assert_in_delta 64.184, TS.tt_utc(Date.new(2000, 1, 6)), 1e-9
+    assert_in_delta 69.184, TS.tt_utc(Date.new(2017, 1, 1)), 1e-9
+  end
+
+  def test_mars24_2000_jan_6
+    # Mars24 worked example: 2000-01-06 00:00:00 UTC → MST 23:59:39
+    dt = DateTime.new(2000, 1, 6, 0, 0, 0)
+    assert_in_delta 2_451_549.5, TS.jd_ut(dt), 1e-9
+    assert_in_delta 2_451_549.500742407, TS.jd_tt(dt), 2e-6
+    h, m, s = TS.mtc_hms(dt)
+    hours = hms_to_hours(h, m, s)
+    assert_in_delta hms_to_hours(23, 59, 39), hours, 1.5 / 3600.0
+    assert_equal 23, h
+    assert_equal 59, m
+  end
+
+  def test_offset_independent
+    utc = DateTime.new(2000, 1, 1, 12, 0, 0)
+    est = utc.new_offset(Rational(-5, 24))
+    assert_in_delta TS.msd(utc), TS.msd(est), 1e-12
+    assert_in_delta TS.mtc_hours(utc), TS.mtc_hours(est), 1e-12
+  end
+
+  def test_wikipedia_2026_aug_12
+    # Wikipedia Timekeeping_on_Mars snapshot: 12 Aug 2026 00:25:08 UTC
+    # JDTT = 2461264.51825, MSD = 54251.08514, MTC = 02:02:36
+    dt = DateTime.new(2026, 8, 12, 0, 25, 8)
+    assert_in_delta 2_461_264.51825, TS.jd_tt(dt), 5e-5
+    assert_in_delta 54_251.08514, TS.msd(dt), 5e-5
+    h, m, s = TS.mtc_hms(dt)
+    assert_in_delta hms_to_hours(2, 2, 36), hms_to_hours(h, m, s), 2.0 / 3600.0
+  end
+
+  def test_mxt_vs_mtc
+    dt = DateTime.new(2000, 1, 6, 0, 0, 0)
+    frac = TS.sol_fraction(dt)
+    assert_operator frac, :>=, 0
+    assert_operator frac, :<, 1
+    si = TS.mxt_si_seconds(dt)
+    assert_in_delta frac * TS::SOL_SI_SECONDS, si, 1e-6
+    mtc_h = TS.mtc_hours(dt)
+    mxt_h = si / 3600.0
+    assert_in_delta mtc_h * TS::SOL_DAYS, mxt_h, 1e-6
+    assert_includes 88_775.0..88_776.0, TS::SOL_SI_SECONDS
+  end
+
+  def test_mtc_rounding_stays_below_24_hours
+    h, m, s = TS.hours_to_hms(24.0.prev_float, limit: 24)
+
+    assert_equal [23, 59], [h, m]
+    assert_operator s, :<, 60
+  end
+
+  def test_ls_j2000_region
+    dt = DateTime.new(2000, 1, 6, 0, 0, 0)
+    ls = TS.ls_deg(TS.jd_tt(dt))
+    assert_in_delta 277.19, ls, 0.05
+  end
+
+  def test_earth_from_msd_round_trip
+    dt = DateTime.new(2000, 1, 6, 0, 0, 0)
+    back = TS.earth_from_msd(TS.msd(dt))
+    assert_in_delta 0, (back - dt).to_f, 1e-9
+
+    dt = DateTime.new(2026, 8, 13, 0, 6, 56)
+    back = TS.earth_from_msd(TS.msd(dt))
+    assert_in_delta 0, (back - dt).to_f, 1e-8
+  end
+
+  def test_year1_caption_treats_tt_as_ut
+    r = TS.epoch_report
+    assert_match(/no year-1/, r[:caveat])
+    assert_equal(-665773, r[:epoch_msd])
+    cap = r[:epoch_earth_if_tt_were_ut]
+    assert_equal 1, cap.year
+    assert_equal 1, cap.month
+    assert_equal 23, cap.day
+    assert_equal 13, cap.hour
+    assert_equal 40, cap.min
+    eq = r[:ls0_earth_if_tt_were_ut]
+    assert_equal 1, eq.year
+    assert_equal 1, eq.month
+    assert_equal 24, eq.day
+    assert_equal 5, eq.hour
+    assert_equal 53, eq.min
+  end
+
+  def test_year1_equinox_sol
+    hint = DateTime.new(1, 1, 22, 12, 0, 0)
+    eq = TS.vernal_equinox_jd_tt(hint.ajd.to_f)
+    ls = TS.ls_deg(eq)
+    dist = [ls, (ls - 360).abs].min
+    assert_operator dist, :<, 0.01, "year-1 Ls=0 search: Ls=#{ls} at JD_TT=#{eq}"
+    m = TS.msd(jd_tt: eq)
+    assert_operator m, :<, 0, "year-1 MSD should be before 1873 (got #{m})"
+    frac = m % 1.0
+    frac += 1 if frac < 0
+    assert_operator frac, :>=, 0
+    assert_operator frac, :<, 1
+  end
+
+  private
+
+  def hms_to_hours(h, m, s)
+    h + m / 60.0 + s / 3600.0
+  end
+end

metadata

@@ -1,15 +1,42 @@
 --- !ruby/object:Gem::Specification
 name: marsdate
 version: !ruby/object:Gem::Version
-  version: 1.1.7
+  version: 2.0.1
 platform: ruby
 authors:
 - Hal Fulton
-autorequire: 
 bindir: bin
 cert_chain: []
-date: 2019-10-31 00:00:00.000000000 Z
-dependencies: []
+date: 1980-01-02 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '6.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '6.0'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '13.4'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '13.4'
 description: |2
                     This is a library for handling dates and times on Mars
                     for the Martian Common Era calendar (created by Hal Fulton).
@@ -20,15 +47,27 @@ executables:
 extensions: []
 extra_rdoc_files: []
 files:
+- CHANGELOG.md
+- LICENSE.txt
+- README.md
 - bin/marsdate
 - lib/marsdate.rb
-- test/spec/marsdate_spec.rb
+- lib/marsdate/calendar.rb
+- lib/marsdate/timescale.rb
+- lib/marsdate/version.rb
+- marsdate.gemspec
+- test/calendar_test.rb
+- test/cli_test.rb
 - test/test.rb
+- test/timescale_test.rb
 homepage: https://github.com/Hal9000/marsdate
 licenses:
-- Ruby
-metadata: {}
-post_install_message: 
+- MIT
+metadata:
+  bug_tracker_uri: https://github.com/Hal9000/marsdate/issues
+  changelog_uri: https://github.com/Hal9000/marsdate/blob/master/CHANGELOG.md
+  source_code_uri: https://github.com/Hal9000/marsdate
+post_install_message: "\n Success! Run executable 'marsdate' for help.\n "
 rdoc_options: []
 require_paths:
 - lib
@@ -36,15 +75,14 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: '0'
+      version: '2.7'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.0.4
-signing_key: 
+rubygems_version: 3.6.9
 specification_version: 4
 summary: Date/time library for Mars
 test_files: []

data/test/spec/marsdate_spec.rb

@@ -1,382 +0,0 @@
-$: << "lib"
-require 'marsdate'
-
-require 'rubygems'
-require 'minitest'
-require 'minitest/spec'
-require 'minitest/autorun'
-
-require 'pp'
-
-def my_assert(expr)
-  if !expr && block_given?
-    puts "\n--- DEBUG: "
-    yield 
-  end
-  assert(expr)
-end
-
-def my_assert_equal(e1, e2)
-  if e1 != e2 && block_given?
-    puts "\n--- DEBUG: "
-    yield 
-  end
-  assert_equal(e1, e2)
-end
-
-## experimenting...
-
-describe "Mars date 1/1/1" do
-  md = MarsDateTime.new(1,1,1)
-  it "should compute backwards to same result" do
-    assert_equal([1,1,1,0,0,0],md.ymshms)
-  end
-  it "should be a Sunday" do
-    assert_equal(0,md.dow)
-  end
-end
-
-describe "Mars date 2/3/4" do
-  md = MarsDateTime.new(2,3,4)
-  it "should compute backwards to same result" do
-    assert_equal([2,3,4,0,0,0],md.ymshms)
-  end
-end
-
-describe "Mars date 5/4/24" do
-  md = MarsDateTime.new(5,4,24)
-  it "should compute backwards to same result" do
-    assert_equal([5,4,24,0,0,0],md.ymshms)
-  end
-end
-
-describe "Mars date 1054/20/20" do
-  md = MarsDateTime.new(1054,20,20)
-  it "should compute backwards to same result" do
-    assert_equal([1054,20,20,0,0,0],md.ymshms)
-  end
-end
-
-describe "Earth date 1/1/22" do
-  ed = DateTime.new(1,1,22)
-  md = MarsDateTime.new(ed)
-  it "should be approximately M-date 1/1/1" do
-    assert_equal([1,1,1],md.ymshms[0..2])
-  end
-  it "should yield a Martian Sunday" do
-    assert_equal(md.day_of_week,"Sunday")
-  end
-end
-
-describe "Earth date May 23, 1961" do
-  e = DateTime.new(1961,5,23)
-  m = MarsDateTime.new(e)
-  it "should yield a Martian Thursday" do
-    assert_equal(m.day_of_week,"Thursday")
-  end
-end
-
-describe "Earth date May 24, 1961" do
-  e = DateTime.new(1961,5,24)
-  m = MarsDateTime.new(e)
-  it "should yield a Martian Friday" do
-    assert_equal(m.day_of_week,"Friday")
-  end
-end
-
-describe "A Martian date/time (y,m,s,h,m,s)" do
-  m = MarsDateTime.new(1043,2,15,12,34,45)
-  it "should set all its accessors properly" do
-    assert_equal(1043,m.myear)
-    assert_equal(2,m.month)
-    assert_equal(15,m.sol)
-    assert_equal(12,m.hr)
-    assert_equal(34,m.min)
-    assert_equal(45,m.sec)
-    assert_equal(4,m.dow)
-    assert_equal("Thursday",m.day_of_week)
-    assert_equal(43,m.year_sol)
-    assert_equal(696715,m.epoch_sol)
-  end
-end
-
-describe "The 2007 Martian equinox" do
-  e = DateTime.new(2007,12,9,17,20,0).new_offset(-5)
-  m = MarsDateTime.new(e)
-  m11 = MarsDateTime.new(1068,1,1)
-  it "should fall on or near 1/1 (MCE)" do
-    diff = (m11-m).abs
-    flag = diff < 1.0
-    assert(flag, "Difference: #{diff}")
-  end
-end
-
-describe "Constructing without parameters" do
-  m = MarsDateTime.new
-  e = DateTime.now
-  m2e = m.earth_date
-  it "should default to today" do
-    diff = (e-m2e).abs     # m - e doesn't work!
-    assert(diff < 1.0)
-  end
-end
-
-describe "Constructing with a DateTime" do
-  it "should convert from Gregorian to MCE" do
-    e = DateTime.new(1902,8,12,7,0,0)
-    m = MarsDateTime.new(e)
-    diff = (m.earth_date - e).abs   # days as a float
-    x = assert(diff < 1.0)
-  end
-end
-
-describe "A Martian date/time (933,1,2)" do
-  m = MarsDateTime.new(933,1,2)
-  it "should set the three obvious accessors properly" do
-    assert_equal(933,m.myear)
-    assert_equal(1,m.month)
-    assert_equal(2,m.sol)
-  end
-end
-
-describe "A Martian date/time (933,1,1)" do
-  m = MarsDateTime.new(933,1,1)
-  it "should set the three obvious accessors properly" do
-    assert_equal(933,m.myear)
-    assert_equal(1,m.month)
-    assert_equal(1,m.sol)
-  end
-end
-
-describe "The Martian date 24/25" do
-  describe "in a leap year" do
-    year = 1067
-    it "should exist" do
-      assert( MarsDateTime.leap?(year) )
-      MarsDateTime.new(year,24,25)   # may raise exception - no #wont_raise
-    end
-  end
-  describe "in a non-leap year" do
-    year = 1066
-    it "should not exist" do
-      assert( ! MarsDateTime.leap?(year) )
-      lambda { MarsDateTime.new(year,24,25) }.must_raise(RuntimeError) 
-    end
-  end
-end
-
-describe "The earth_date convertor" do
-  m = MarsDateTime.new(1,1,1)
-  it "should convert to a Gregorian date" do
-    e = m.earth_date
-    assert_equal(1,e.year)
-    assert_equal(1,e.month)
-    assert((e.day-22).abs <= 1, "e - 22 = #{e.day - 22}")
-  end
-end
-
-describe "An Earthly date" do
-  e1 = DateTime.new(1961,5,31)
-  it "should (approximately) make the conversion roundtrip" do
-    m1 = MarsDateTime.new(e1)
-    e2 = m1.earth_date
-    diff = e2 - e1
-    assert(diff.abs <= 0.01, "Difference = #{diff}")
-  end
-end
-
-describe "A Martian date" do
-  m1 = MarsDateTime.new(1067,1,1)
-  it "should (approximately) make the conversion roundtrip" do
-    e1 = m1.earth_date
-    m2 = MarsDateTime.new(e1)
-    diff = m2 - m1
-    assert(diff.abs < 0.01)
-  end
-end
-
-describe "An arbitrary Martian date" do
-  m1 = MarsDateTime.new(1069,15,24)
-  m2 = MarsDateTime.new(933,6,4)
-  m3 = MarsDateTime.new(1055,14,1)
-  it "should format well with strftime" do
-    assert_equal("Mon",m1.strftime("%a"))
-    assert_equal("Monday",m1.strftime("%A"))
-    assert_equal("Aug",m1.strftime("%b"))
-    assert_equal("M-August",m1.strftime("%B"))
-    assert_equal("24",m1.strftime("%d"))
-    assert_equal("24",m1.strftime("%e"))
-    assert_equal("1069-15-24",m1.strftime("%F"))
-    assert_equal("416",m1.strftime("%j"))
-    assert_equal("15",m1.strftime("%m"))
-    assert_equal("0",m1.strftime("%s"))
-    assert_equal("2",m1.strftime("%u"))
-    assert_equal("60",m1.strftime("%U"))
-    assert_equal("1",m1.strftime("%w"))
-    assert_equal("1069/15/24",m1.strftime("%x"))
-    assert_equal("1069",m1.strftime("%Y"))
-    assert_equal("\n",m1.strftime("%n"))
-    assert_equal("\t",m1.strftime("%t"))
-    assert_equal("%",m1.strftime("%%"))
-
-    assert_equal("Wed",m2.strftime("%a"))
-    assert_equal("Wednesday",m2.strftime("%A"))
-    assert_equal("Leo",m2.strftime("%b"))
-    assert_equal("Leo",m2.strftime("%B"))
-    assert_equal("04",m2.strftime("%d"))
-    assert_equal(" 4",m2.strftime("%e"))
-    assert_equal("933-06-04",m2.strftime("%F"))
-    assert_equal("144",m2.strftime("%j"))
-    assert_equal("06",m2.strftime("%m"))
-    assert_equal("0",m2.strftime("%s"))
-    assert_equal("4",m2.strftime("%u"))
-    assert_equal("21",m2.strftime("%U"))
-    assert_equal("3",m2.strftime("%w"))
-    assert_equal("933/06/04",m2.strftime("%x"))
-    assert_equal("933",m2.strftime("%Y"))
-
-    assert_equal("Sag",m3.strftime("%b"))
-    assert_equal("Sagittarius",m3.strftime("%B"))
-  end
-end
-
-describe "An Earthly date with h:m:s" do
-  e1 = DateTime.new(1976,7,4,16,30,0)  # July 4, 1976  16:30
-  m1 = MarsDateTime.new(e1)
-  it "should make the conversion roundtrip" do
-    30.times do
-      m1 = MarsDateTime.new(e1)
-      e2 = m1.earth_date
-      diff = e2 - e1
-      assert(diff.abs <= 0.01) { STDERR.puts "    diff = #{diff}" }
-      e1 += 1
-    end
-  end
-end
-
- describe "Martian dates" do
-   m1 = MarsDateTime.new(1068,14,22)
-   m2 = MarsDateTime.new(1068,14,21)
-   e1 = DateTime.new(2002,4,19)
-   e2 = DateTime.new(2012,9,29)
-   it "should compare with each other" do
-     assert(m1 > m2)
-     assert(m2 < m1)
-   end
-   it "should compare with Earth dates" do
-     assert(m1 > e1)
-     assert(m2 < e2)
-   end
-   it "should honor equality (within roundoff)" do
-     assert(m1 == m1)
-     assert(m1 != m2)
-   end
- end
-
-describe "The 'now' and 'today' class methods" do
-  t1 = MarsDateTime.now
-  t0 = MarsDateTime.today
-  it "should be less than a day apart" do
-    diff = t1 - t0
-    assert(diff < 1.0)
-  end
-end
-
-describe "In M-year 1043, each month" do
-  md = (1..24).to_a.map {|mm| MarsDateTime.new(1043,mm,1) }
-  it "should start on a Thursday" do
-    md.each {|m| assert(m.day_of_week == "Thursday") }
-  end
-end
-
-describe "Converting 1961/5/31 to Martian" do
-  e = DateTime.new(1961,5,31) 
-  m = MarsDateTime.new(e)
-  it "should yield Friday M-April 9, 1043 MCE" do
-    assert(m.year == 1043, "bad year: #{m.year}")
-    assert(m.month == 7, "bad month: #{m.month}")
-    assert(m.sol == 9, "bad sol: #{m.sol}")
-  end
-  it "should give a Martian Saturday" do
-    assert(m.day_of_week == "Friday", "got: #{m.day_of_week}")
-  end
-end
-
-describe "The 1609 Martian equinox" do
-  e = DateTime.new(1609,3,12).new_offset(-5)
-  m = MarsDateTime.new(e)
-  m11 = MarsDateTime.new(856,1,1)
-  it "should fall on or near 1/1 (MCE)" do
-    diff = (m11-m).abs
-    assert(diff < 1.0, "Difference = #{diff}")
-  end
-end
-
-describe "The 1902 Martian equinox" do
-  e = DateTime.new(1902,8,12,7,0,0).new_offset(-5)
-  m = MarsDateTime.new(e)
-  m11 = MarsDateTime.new(1012,1,1)
-  it "should fall on or near 1/1 (MCE)" do
-    diff = (m11-m).abs
-    assert(diff < 1.0, "Difference: #{diff}")
-  end
-end
-
-describe "A Martian date" do
-  puts "***** Entering problem test..."
-  m = MarsDateTime.new(1068,1,1)
-  p m
-  describe "plus 7 days" do
-    m2 = m + 7
-    p m2
-    it "should be the same day of the week" do
-      assert_equal(m.day_of_week, m2.day_of_week)
-    end
-  end
-  describe "plus 14 days" do
-    m3 = m + 14
-    p m3
-    it "should be the same day of the week" do
-      assert_equal(m.day_of_week, m3.day_of_week)
-    end
-  end
-  describe "plus 21 days" do
-    m4 = m + 21
-    p m4
-    it "should be the same day of the week" do
-      assert_equal(m.day_of_week, m4.day_of_week)
-    end
-  end
-  puts "***** Exiting problem test"
-end
-
-describe "Martian April 1, 1043" do
-  m = MarsDateTime.new(1043,7,1)
-  it "should be a Thursday" do
-    assert_equal(m.day_of_week,"Thursday")
-  end
-end
- 
-describe "A Martian date" do
-  md = MarsDateTime.new(1062,20,20)
-  m2 = MarsDateTime.new(1062,20,17) 
-  it "should allow subtraction of another date" do
-    m = md - m2
-    assert( (m - 3.0) < 0.01)   # sols
-  end
-  it "should allow subtraction of an Earth date" do
-    e = DateTime.new(1998, 3, 12, 16, 45, 0)
-    m = md - e
-    assert( (m - 3.0) < 0.01)   # sols
-  end
-  it "should allow subtraction of a Fixnum" do
-    m = md - 3                 # sols
-    assert( (m - m2) < 0.01)
-  end
-  it "should allow subtraction of a Float" do
-    m = md - 3.0               # sols
-    assert( (m - m2) < 0.01)
-  end
-end
-
-