RubyGems - lunation - Versions diffs - 0.1.5 → 0.1.6 - Mend

lunation 0.1.5 → 0.1.6

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (16) hide show

checksums.yaml +4 -4
data/.rubocop.yml +18 -0
data/CHANGELOG.md +8 -0
data/config/periodic_terms_earth_nutation.yml +63 -567
data/config/periodic_terms_moon_latitude.yml +60 -300
data/config/periodic_terms_moon_longitude_distance.yml +60 -360
data/lib/lunation/calculation/angle.rb +9 -6
data/lib/lunation/calculation/earth_position_vsop87.rb +1 -1
data/lib/lunation/calculation/moon_illuminated_fraction.rb +3 -2
data/lib/lunation/calculation/moon_position.rb +47 -45
data/lib/lunation/calculation/nutation_and_obliquity.rb +17 -17
data/lib/lunation/calculation/sun_position.rb +6 -3
data/lib/lunation/calculation/timekeeping.rb +14 -13
data/lib/lunation/calculation.rb +5 -1
data/lib/lunation/version.rb +1 -1
metadata +3 -3

data/lib/lunation/calculation/nutation_and_obliquity.rb CHANGED Viewed

@@ -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 CHANGED Viewed

@@ -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 CHANGED Viewed

@@ -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 CHANGED Viewed

@@ -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 CHANGED Viewed

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

metadata CHANGED Viewed

@@ -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: