lunation 0.1.5 → 0.1.6

data/lib/lunation/calculation/nutation_and_obliquity.rb

@@ -5,11 +5,11 @@ module Lunation
       NUTATION_IN_OBLIQUITY_PERIODIC_TERMS_PATH = File.expand_path("../../../config/periodic_terms_earth_nutation.yml", __dir__)
       NUTATION_IN_OBLIQUITY_PERIODIC_TERMS = YAML.load_file(NUTATION_IN_OBLIQUITY_PERIODIC_TERMS_PATH).freeze
       ECLIPTIC_MEAN_OBLIQUITY_CONSTANTS = [21.448, -4680.93, -1.55, 1_999.25, -51.38, -249.67, -39.05, 7.12, 27.87, 5.79, 2.45].freeze
-      MOON_MEAN_ELONGATION_FROM_SUN_CONSTANTS = [297.85036, 445_267.111480, -0.0019142, 1 / 189_474.0].freeze
-      SUN_MEAN_ANOMALY_CONSTANTS = [357.52772, 35_999.050340, -0.0001603, -1 / 300_000.0].freeze
-      MOON_MEAN_ANOMALY_CONSTANTS = [134.96298, 477_198.867398, 0.0086972, 1 / 56_250.0].freeze
-      MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.27191, 483_202.017538, -0.0036825, 1 / 327_270.0].freeze
-      LONGITUDE_OF_THE_ASCENDING_NODE_CONSTANTS = [125.04452, -1934.136261, 0.0020708, 1 / 450_000.0].freeze
+      MOON_MEAN_ELONGATION_FROM_SUN_CONSTANTS = [297.85036, 445_267.111480, -0.0019142, 1.fdiv(189_474)].freeze
+      SUN_MEAN_ANOMALY_CONSTANTS = [357.52772, 35_999.050340, -0.0001603, -1.fdiv(300_000)].freeze
+      MOON_MEAN_ANOMALY_CONSTANTS = [134.96298, 477_198.867398, 0.0086972, 1.fdiv(56_250)].freeze
+      MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.27191, 483_202.017538, -0.0036825, 1.fdiv(327_270)].freeze
+      LONGITUDE_OF_THE_ASCENDING_NODE_CONSTANTS = [125.04452, -1934.136261, 0.0020708, 1.fdiv(450_000)].freeze
       # rubocop:enable Layout/LineLength
 
       # (D) Mean elongation of the moon from the sun (A.A. p. 144)
@@ -64,13 +64,13 @@ module Lunation
         @nutation_in_longitude ||= begin
           result = NUTATION_IN_OBLIQUITY_PERIODIC_TERMS.inject(0.0) do |acc, elem|
             argument = Angle.from_decimal_degrees(
-              elem["moon_mean_elongation"] * moon_mean_elongation_from_sun.decimal_degrees +
-              elem["sun_mean_anomaly"] * sun_mean_anomaly.decimal_degrees +
-              elem["moon_mean_anomaly"] * moon_mean_anomaly.decimal_degrees +
-              elem["moon_argument_of_latitude"] * moon_argument_of_latitude.decimal_degrees +
-              elem["longitude_of_ascending_node"] * longitude_of_ascending_node.decimal_degrees
+              elem[0] * moon_mean_elongation_from_sun.decimal_degrees +
+              elem[1] * sun_mean_anomaly.decimal_degrees +
+              elem[2] * moon_mean_anomaly.decimal_degrees +
+              elem[3] * moon_argument_of_latitude.decimal_degrees +
+              elem[4] * longitude_of_ascending_node.decimal_degrees
             )
-            coefficient = elem["sine_coefficient_first_term"] + elem["sine_coefficient_second_term"] * time
+            coefficient = elem[5] + elem[6] * time
             acc + coefficient * argument.sin
           end / 10_000.0
           Angle.from_decimal_arcseconds(result)
@@ -83,13 +83,13 @@ module Lunation
         @nutation_in_obliquity ||= begin
           result = NUTATION_IN_OBLIQUITY_PERIODIC_TERMS.inject(0.0) do |acc, elem|
             argument = Angle.from_decimal_degrees(
-              elem["moon_mean_elongation"] * moon_mean_elongation_from_sun.decimal_degrees +
-              elem["sun_mean_anomaly"] * sun_mean_anomaly.decimal_degrees +
-              elem["moon_mean_anomaly"] * moon_mean_anomaly.decimal_degrees +
-              elem["moon_argument_of_latitude"] * moon_argument_of_latitude.decimal_degrees +
-              elem["longitude_of_ascending_node"] * longitude_of_ascending_node.decimal_degrees
+              elem[0] * moon_mean_elongation_from_sun.decimal_degrees +
+              elem[1] * sun_mean_anomaly.decimal_degrees +
+              elem[2] * moon_mean_anomaly.decimal_degrees +
+              elem[3] * moon_argument_of_latitude.decimal_degrees +
+              elem[4] * longitude_of_ascending_node.decimal_degrees
             )
-            coefficient = elem["cosine_coefficient_first_term"] + elem["cosine_coefficient_second_term"] * time
+            coefficient = elem[7] + elem[8] * time
             acc + coefficient * argument.cos
           end / 10_000.0
           Angle.from_decimal_arcseconds(result)

data/lib/lunation/calculation/sun_position.rb

@@ -58,7 +58,10 @@ module Lunation
       # UNIT: Astronomical Units (AU)
       def distance_between_earth_and_sun_in_astronomical_units
         @distance_between_earth_and_sun_in_astronomical_units ||= begin
-          result = 1.000001018 * (1 - earth_orbit_eccentricity**2) / (1 + earth_orbit_eccentricity * sun_anomaly.cos)
+          result = 1.000001018 *
+            (1 - earth_orbit_eccentricity**2) / (
+              (1 + earth_orbit_eccentricity * sun_anomaly.cos)
+            )
           result.round(7)
         end
       end
@@ -70,8 +73,8 @@ module Lunation
           (distance_between_earth_and_sun_in_astronomical_units * 149_597_870).floor
       end
 
-      # (Ω) Longitude of the ascending node of the Moon's mean orbit on the ecliptic (low precision)
-      #   A.A. p. 164
+      # (Ω) Longitude of the ascending node of the Moon's mean orbit on the ecliptic
+      #   (low precision) A.A. p. 164
       # UNIT: Angle
       def longitude_of_ascending_node_low_precision
         @longitude_of_ascending_node_low_precision ||=

data/lib/lunation/calculation/timekeeping.rb

@@ -7,14 +7,14 @@ module Lunation
       # rubocop:disable Layout/LineLength
       DELTA_T_500BC_500AD = [10_583.6, -1_014.41, 33.78311, -5.952053, -0.1798452, 0.022174192, 0.0090316521].freeze
       DELTA_T_500AD_1600AD = [1_574.2, -556.01, 71.23472, 0.319781, -0.8503463, -0.005050998, 0.0083572073].freeze
-      DELTA_T_1600AD_1700AD = [120, -0.9808, -0.01532, 1 / 7_129.0].freeze
-      DELTA_T_1700AD_1800AD = [8.83, 0.1603, -0.0059285, 0.00013336, -1 / 1_174_000.0].freeze
+      DELTA_T_1600AD_1700AD = [120, -0.9808, -0.01532, 1.fdiv(7_129)].freeze
+      DELTA_T_1700AD_1800AD = [8.83, 0.1603, -0.0059285, 0.00013336, -1.fdiv(1_174_000)].freeze
       DELTA_T_1800AD_1860AD = [13.72, -0.332447, 0.0068612, 0.0041116, -0.00037436, 0.0000121272, -0.0000001699, 0.000000000875].freeze
-      DELTA_T_1860AD_1900AD = [7.62, 0.5737, -0.251754, 0.01680668, -0.0004473624, 1 / 233_174.0].freeze
+      DELTA_T_1860AD_1900AD = [7.62, 0.5737, -0.251754, 0.01680668, -0.0004473624, 1.fdiv(233_174)].freeze
       DELTA_T_1900AD_1920AD = [-2.79, 1.494119, -0.0598939, 0.0061966, -0.000197].freeze
       DELTA_T_1920AD_1941AD = [21.20, 0.84493, -0.076100, 0.0020936].freeze
-      DELTA_T_1941AD_1961AD = [29.07, 0.407, -1 / 233.0, 1 / 2_547.0].freeze
-      DELTA_T_1961AD_1986AD = [45.45, 1.067, -1 / 260.0, -1 / 718.0].freeze
+      DELTA_T_1941AD_1961AD = [29.07, 0.407, -1.fdiv(233), 1.fdiv(2_547)].freeze
+      DELTA_T_1961AD_1986AD = [45.45, 1.067, -1.fdiv(260), -1.fdiv(718)].freeze
       DELTA_T_1986AD_2005AD = [63.86, 0.3345, -0.060374, 0.0017275, 0.000651814, 0.00002373599].freeze
       DELTA_T_2005AD_2050AD = [62.92, 0.32217, 0.005589].freeze
       # rubocop:enable Layout/LineLength
@@ -23,7 +23,7 @@ module Lunation
       # UNIT: centuries from the Epoch J2000.0
       def time
         @time ||=
-          ((julian_ephemeris_day - EPOCH_J2000_JDE) / JULIAN_CENTURY_IN_DAYS).round(12)
+          (julian_ephemeris_day - EPOCH_J2000_JDE).fdiv(JULIAN_CENTURY_IN_DAYS).round(12)
       end
 
       # ΔT Difference between Dynamical Time (TD) and Universal Time (UT, more commonly
@@ -36,12 +36,12 @@ module Lunation
         @delta_t ||= begin
           result = case datetime.year
           when -1_999...-500
-            elapsed_years = (@decimal_year - 1_820) / 100
+            elapsed_years = (@decimal_year - 1_820).fdiv(100)
             -20 + 32 * elapsed_years**2
           when -500...500
-            Horner.compute(@decimal_year / 100.0, DELTA_T_500BC_500AD)
+            Horner.compute(@decimal_year.fdiv(100), DELTA_T_500BC_500AD)
           when 500...1_600
-            Horner.compute((@decimal_year - 1_000.0) / 100.0, DELTA_T_500AD_1600AD)
+            Horner.compute((@decimal_year - 1_000.0).fdiv(100), DELTA_T_500AD_1600AD)
           when 1_600...1_700
             Horner.compute(@decimal_year - 1_600.0, DELTA_T_1600AD_1700AD)
           when 1_700...1_800
@@ -63,7 +63,8 @@ module Lunation
           when 2_005...2_050
             Horner.compute(@decimal_year - 2_000, DELTA_T_2005AD_2050AD)
           when 2_050...2_150
-            -20 + 32 * ((@decimal_year - 1_820) / 100)**2 - 0.5628 * (2_150 - @decimal_year)
+            -20 + 32 * ((@decimal_year - 1_820) / 100)**2 - 0.5628 *
+              (2_150 - @decimal_year)
           when 2_150..3_000
             -20 + 32 * @decimal_year**2
           end
@@ -74,7 +75,7 @@ module Lunation
       # (TD) Dynamical time (A.A. p. 77)
       # UNIT: ISO 8601 Date and time with offset
       def dynamical_time
-        @dynamical_time ||= datetime + delta_t.round.to_f / SECONDS_PER_DAY
+        @dynamical_time ||= datetime + delta_t.round.fdiv(SECONDS_PER_DAY)
       end
 
       # (JDE) Julian ephemeris day (A.A. p. 59)
@@ -86,13 +87,13 @@ module Lunation
       # (t) time, measured in Julian millennia from the epoch J2000.0 (32.1, A.A. p. 218)
       # UNIT: millennia from the Epoch J2000.0
       def time_millennia
-        @time_millennia ||= (time / 10.0).round(12)
+        @time_millennia ||= time.fdiv(10).round(12)
       end
 
       # (U) Time measured in units of 10_000 Julian years from J2000.0 (A.A. p. 147)
       # UNIT: 10_000 years (myriads) from the Epoch J2000.0
       def time_myriads
-        @time_myriads ||= time / 100.0
+        @time_myriads ||= time.fdiv(100)
       end
     end
   end

data/lib/lunation/calculation.rb

@@ -25,11 +25,14 @@ module Lunation
 
     def to_h
       {
+        # rubocop:disable Layout/LineLength
         corrected_obliquity_of_ecliptic: corrected_obliquity_of_ecliptic,
         correction_eccentricity_of_earth: correction_eccentricity_of_earth,
         correction_jupiter: correction_jupiter,
         correction_latitude: correction_latitude,
         correction_venus: correction_venus,
+        datetime: datetime,
+        decimal_year: decimal_year,
         delta_t: delta_t,
         distance_between_earth_and_moon: distance_between_earth_and_moon,
         distance_between_earth_and_sun_in_astronomical_units: distance_between_earth_and_sun_in_astronomical_units,
@@ -78,11 +81,12 @@ module Lunation
         time_millennia: time_millennia,
         time_myriads: time_myriads,
         time: time,
+        # rubocop:enable Layout/LineLength
       }
     end
 
     def to_s
-      <<-HEREDOC
+      <<~HEREDOC
         DATE AND TIME (UT):                           #{datetime}
         ----------------------------------------------------------------------------------
         DECIMAL YEAR:                                 #{decimal_year}

data/lib/lunation/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Lunation
-  VERSION = "0.1.5"
+  VERSION = "0.1.6"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: lunation
 version: !ruby/object:Gem::Version
-  version: 0.1.5
+  version: 0.1.6
 platform: ruby
 authors:
 - Ivo Kalverboer
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2024-12-30 00:00:00.000000000 Z
+date: 2024-12-31 00:00:00.000000000 Z
 dependencies: []
 description: Lunation offers a Ruby implementation of Meeus's Astronomical Algorithms.
 email:
@@ -59,7 +59,7 @@ licenses:
 metadata:
   homepage_uri: https://www.github.com/valerius/lunation
   source_code_uri: https://www.github.com/valerius/lunation
-  changelog_uri: https://www.github.com/valerius/lunation/CHANGELOG.md
+  changelog_uri: https://github.com/valerius/lunation/blob/main/CHANGELOG.md
 post_install_message: 
 rdoc_options: []
 require_paths: