lunation 0.1.5 → 0.1.6

data/config/periodic_terms_moon_longitude_distance.yml

@@ -1,360 +1,60 @@
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 6288774
-  cosine_coefficient: -20905355
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 1274027
-  cosine_coefficient: -3699111
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 658314
-  cosine_coefficient: -2955968
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 213618
-  cosine_coefficient: -569925
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: -185116
-  cosine_coefficient: 48888
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 2
-  sine_coefficient: -114332
-  cosine_coefficient: -3149
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 58793
-  cosine_coefficient: 246158
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 57066
-  cosine_coefficient: -152138
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 53322
-  cosine_coefficient: -170733
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 45758
-  cosine_coefficient: -204586
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -40923
-  cosine_coefficient: -129620
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: -34720
-  cosine_coefficient: 108743
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -30383
-  cosine_coefficient: 104755
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: -2
-  sine_coefficient: 15327
-  cosine_coefficient: 10321
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 2
-  sine_coefficient: -12528
-  cosine_coefficient: null
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: -2
-  sine_coefficient: 10980
-  cosine_coefficient: 79661
-- moon_mean_elongation: 4
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 10675
-  cosine_coefficient: -34782
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 3
-  moon_argument_of_latitude: 0
-  sine_coefficient: 10034
-  cosine_coefficient: -23210
-- moon_mean_elongation: 4
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 8548
-  cosine_coefficient: -21636
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -7888
-  cosine_coefficient: 24208
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: -6766
-  cosine_coefficient: 30824
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -5163
-  cosine_coefficient: -8379
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 4987
-  cosine_coefficient: -16675
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 4036
-  cosine_coefficient: -12831
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 3994
-  cosine_coefficient: -10445
-- moon_mean_elongation: 4
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 3861
-  cosine_coefficient: -11650
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -3
-  moon_argument_of_latitude: 0
-  sine_coefficient: 3665
-  cosine_coefficient: 14403
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: -2689
-  cosine_coefficient: -7003
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 2
-  sine_coefficient: -2602
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 2390
-  cosine_coefficient: 10056
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -2348
-  cosine_coefficient: 6322
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -2
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 2236
-  cosine_coefficient: -9884
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: 0
-  sine_coefficient: -2120
-  cosine_coefficient: 5751
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 2
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: -2069
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -2
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 2048
-  cosine_coefficient: -4950
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: -2
-  sine_coefficient: -1773
-  cosine_coefficient: 4130
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 2
-  sine_coefficient: -1595
-  cosine_coefficient: null
-- moon_mean_elongation: 4
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 1215
-  cosine_coefficient: -3958
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: 2
-  sine_coefficient: -1110
-  cosine_coefficient: null
-- moon_mean_elongation: 3
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -892
-  cosine_coefficient: 3258
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -810
-  cosine_coefficient: 2616
-- moon_mean_elongation: 4
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 759
-  cosine_coefficient: -1897
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 2
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -713
-  cosine_coefficient: -2117
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 2
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -700
-  cosine_coefficient: 2354
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 691
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: -2
-  sine_coefficient: 596
-  cosine_coefficient: null
-- moon_mean_elongation: 4
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 549
-  cosine_coefficient: -1423
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 4
-  moon_argument_of_latitude: 0
-  sine_coefficient: 537
-  cosine_coefficient: -1117
-- moon_mean_elongation: 4
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: 0
-  sine_coefficient: 520
-  cosine_coefficient: -1571
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: -487
-  cosine_coefficient: -1739
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 0
-  moon_argument_of_latitude: -2
-  sine_coefficient: -399
-  cosine_coefficient: null
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: -2
-  sine_coefficient: -381
-  cosine_coefficient: -4421
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 351
-  cosine_coefficient: null
-- moon_mean_elongation: 3
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -2
-  moon_argument_of_latitude: 0
-  sine_coefficient: -340
-  cosine_coefficient: null
-- moon_mean_elongation: 4
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -3
-  moon_argument_of_latitude: 0
-  sine_coefficient: 330
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: -1
-  moon_mean_anomaly: 2
-  moon_argument_of_latitude: 0
-  sine_coefficient: 327
-  cosine_coefficient: null
-- moon_mean_elongation: 0
-  sun_mean_anomaly: 2
-  moon_mean_anomaly: 1
-  moon_argument_of_latitude: 0
-  sine_coefficient: -323
-  cosine_coefficient: 1165
-- moon_mean_elongation: 1
-  sun_mean_anomaly: 1
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: 0
-  sine_coefficient: 299
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: 3
-  moon_argument_of_latitude: 0
-  sine_coefficient: 294
-  cosine_coefficient: null
-- moon_mean_elongation: 2
-  sun_mean_anomaly: 0
-  moon_mean_anomaly: -1
-  moon_argument_of_latitude: -2
-  sine_coefficient: null
-  cosine_coefficient: 8752
+- [0, 0, 1, 0, 6288774, -20905355]
+- [2, 0, -1, 0, 1274027, -3699111]
+- [2, 0, 0, 0, 658314, -2955968]
+- [0, 0, 2, 0, 213618, -569925]
+- [0, 1, 0, 0, -185116, 48888]
+- [0, 0, 0, 2, -114332, -3149]
+- [2, 0, -2, 0, 58793, 246158]
+- [2, -1, -1, 0, 57066, -152138]
+- [2, 0, 1, 0, 53322, -170733]
+- [2, -1, 0, 0, 45758, -204586]
+- [0, 1, -1, 0, -40923, -129620]
+- [1, 0, 0, 0, -34720, 108743]
+- [0, 1, 1, 0, -30383, 104755]
+- [2, 0, 0, -2, 15327, 10321]
+- [0, 0, 1, 2, -12528, null]
+- [0, 0, 1, -2, 10980, 79661]
+- [4, 0, -1, 0, 10675, -34782]
+- [0, 0, 3, 0, 10034, -23210]
+- [4, 0, -2, 0, 8548, -21636]
+- [2, 1, -1, 0, -7888, 24208]
+- [2, 1, 0, 0, -6766, 30824]
+- [1, 0, -1, 0, -5163, -8379]
+- [1, 1, 0, 0, 4987, -16675]
+- [2, -1, 1, 0, 4036, -12831]
+- [2, 0, 2, 0, 3994, -10445]
+- [4, 0, 0, 0, 3861, -11650]
+- [2, 0, -3, 0, 3665, 14403]
+- [0, 1, -2, 0, -2689, -7003]
+- [2, 0, -1, 2, -2602, null]
+- [2, -1, -2, 0, 2390, 10056]
+- [1, 0, 1, 0, -2348, 6322]
+- [2, -2, 0, 0, 2236, -9884]
+- [0, 1, 2, 0, -2120, 5751]
+- [0, 2, 0, 0, -2069, null]
+- [2, -2, -1, 0, 2048, -4950]
+- [2, 0, 1, -2, -1773, 4130]
+- [2, 0, 0, 2, -1595, null]
+- [4, -1, -1, 0, 1215, -3958]
+- [0, 0, 2, 2, -1110, null]
+- [3, 0, -1, 0, -892, 3258]
+- [2, 1, 1, 0, -810, 2616]
+- [4, -1, -2, 0, 759, -1897]
+- [0, 2, -1, 0, -713, -2117]
+- [2, 2, -1, 0, -700, 2354]
+- [2, 1, -2, 0, 691, null]
+- [2, -1, 0, -2, 596, null]
+- [4, 0, 1, 0, 549, -1423]
+- [0, 0, 4, 0, 537, -1117]
+- [4, -1, 0, 0, 520, -1571]
+- [1, 0, -2, 0, -487, -1739]
+- [2, 1, 0, -2, -399, null]
+- [0, 0, 2, -2, -381, -4421]
+- [1, 1, 1, 0, 351, null]
+- [3, 0, -2, 0, -340, null]
+- [4, 0, -3, 0, 330, null]
+- [2, -1, 2, 0, 327, null]
+- [0, 2, 1, 0, -323, 1165]
+- [1, 1, -1, 0, 299, null]
+- [2, 0, 3, 0, 294, null]
+- [2, 0, -1, -2, null, 8752]

data/lib/lunation/calculation/angle.rb

@@ -40,7 +40,7 @@ module Lunation
       end
 
       def decimal_hours
-        @decimal_hours ||= decimal_degrees / 15.0
+        @decimal_hours ||= decimal_degrees.fdiv(15)
       end
 
       def decimal_minutes
@@ -95,7 +95,9 @@ module Lunation
 
         def from_decimal_degrees(decimal_degrees, normalize: true)
           new(
-            decimal_degrees: (normalize ? decimal_degrees % 360 : decimal_degrees).round(9),
+            decimal_degrees: (
+                normalize ? decimal_degrees % 360 : decimal_degrees
+              ).round(9),
             radians: decimal_degrees_to_radians(decimal_degrees, normalize: normalize)
           )
         end
@@ -149,7 +151,8 @@ module Lunation
         end
 
         private def decimal_degrees_to_radians(decimal_degrees, normalize: true)
-          result = (normalize ? decimal_degrees % 360 : decimal_degrees) * Math::PI / 180.0
+          result = (normalize ? decimal_degrees % 360 : decimal_degrees) *
+            Math::PI / 180.0
           result.round(9)
         end
 
@@ -159,7 +162,7 @@ module Lunation
           decimal_arcseconds,
           normalize: true
         )
-          result = degrees + arcminutes / 60.0 + decimal_arcseconds / 3600.0
+          result = degrees + arcminutes.fdiv(60) + decimal_arcseconds.fdiv(3600)
           (normalize ? result % 360 : result).round(9)
         end
 
@@ -169,7 +172,7 @@ module Lunation
           decimal_seconds,
           normalize: true
         )
-          result = hours * 15.0 + minutes / 4.0 + decimal_seconds / 240.0
+          result = hours * 15.0 + minutes.fdiv(4.0) + decimal_seconds.fdiv(240.0)
           (normalize ? result % 360 : result).round(9)
         end
 
@@ -177,7 +180,7 @@ module Lunation
           decimal_arcseconds,
           normalize: true
         )
-          result = decimal_arcseconds / 3600.0
+          result = decimal_arcseconds.fdiv(3600.0)
           (normalize ? result % 360 : result).round(9)
         end
       end

data/lib/lunation/calculation/earth_position_vsop87.rb

@@ -58,7 +58,7 @@ module Lunation
           second_series = reduce_periodic_terms(PERIODIC_TERMS_B1)
 
           Angle.from_radians(
-            (first_series + second_series * time_millennia) / 100_000_000.0,
+            (first_series + second_series * time_millennia).fdiv(100_000_000),
             normalize: false
           )
         end

data/lib/lunation/calculation/moon_illuminated_fraction.rb

@@ -5,14 +5,15 @@ module Lunation
       # UNIT: fraction (decimal)
       def moon_illuminated_fraction
         @moon_illuminated_fraction ||=
-          ((1 + moon_phase_angle.cos) / 2.0).round(4)
+          (1 + moon_phase_angle.cos).fdiv(2).round(4)
       end
 
       # (i) phase angle of the moon (48.3, A.A. p. 346)
       # UNIT: Angle
       def moon_phase_angle
         @moon_phase_angle ||= begin
-          numerator = distance_between_earth_and_sun_in_kilometers * moon_elongation_from_sun.sin
+          numerator = distance_between_earth_and_sun_in_kilometers *
+            moon_elongation_from_sun.sin
           denominator = distance_between_earth_and_moon -
             distance_between_earth_and_sun_in_kilometers * moon_elongation_from_sun.cos
           Angle.from_radians(Math.atan2(numerator, denominator))

data/lib/lunation/calculation/moon_position.rb

@@ -4,14 +4,14 @@ module Lunation
   class Calculation
     module MoonPosition
       # rubocop:disable Layout/LineLength
-      LONGITUDE_AND_DISTANCE_OF_MOON_PERIODIC_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_longitude_distance.yml", __dir__).freeze
-      LATITUDE_OF_MOON_PERIODIC_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_latitude.yml", __dir__).freeze
-      LONGITUDE_AND_DISTANCE_OF_MOON_PERIODIC_TERMS = YAML.load_file(LONGITUDE_AND_DISTANCE_OF_MOON_PERIODIC_TERMS_PATH).freeze
-      LATITUDE_OF_MOON_PERIODIC_TERMS = YAML.load_file(LATITUDE_OF_MOON_PERIODIC_TERMS_PATH).freeze
-      MOON_MEAN_LONGITUDE_CONSTANTS = [218.3164477, 481_267.88123421, -0.0015786, 1 / 538_841.0, -1 / 65_194_000.0].freeze
-      MOON_MEAN_ELONGATION_CONSTANTS = [297.8501921, 445_267.1114034, -0.0018819, 1 / 545_868.0, -1 / 113_065_000.0].freeze
-      MOON_MEAN_ANOMALY_CONSTANTS = [134.9633964, 477_198.8675055, 0.0087414, 1 / 69_699.0, -1 / 14_712_000.0].freeze
-      MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.2720950, 483_202.0175233, -0.0036539, -1 / 3_526_000.0, 1 / 863_310_000.0].freeze
+      MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_longitude_distance.yml", __dir__).freeze
+      MOON_LATITUDE_PERIOD_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_latitude.yml", __dir__).freeze
+      MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS = YAML.load_file(MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS_PATH).freeze
+      MOON_LATITUDE_PERIOD_TERMS = YAML.load_file(MOON_LATITUDE_PERIOD_TERMS_PATH).freeze
+      MOON_MEAN_LONGITUDE_CONSTANTS = [218.3164477, 481_267.88123421, -0.0015786, 1.fdiv(538_841), -1.fdiv(65_194_000)].freeze
+      MOON_MEAN_ELONGATION_CONSTANTS = [297.8501921, 445_267.1114034, -0.0018819, 1.fdiv(545_868), -1.fdiv(113_065_000)].freeze
+      MOON_MEAN_ANOMALY_CONSTANTS = [134.9633964, 477_198.8675055, 0.0087414, 1.fdiv(69_699), -1.fdiv(14_712_000)].freeze
+      MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.2720950, 483_202.0175233, -0.0036539, -1.fdiv(3_526_000), 1.fdiv(863_310_000)].freeze
       CORRECTION_ECCENTRICITY_OF_EARTH_CONSTANTS = [1, -0.002516, -0.0000074].freeze
       FIXED_DISTANCE_BETWEEN_EARTH_AND_MOON = 385_000.56
       RADIUS_OF_EARTH = 6378.14
@@ -77,22 +77,22 @@ module Lunation
       # UNIT: decimal degrees
       def moon_heliocentric_longitude
         @moon_heliocentric_longitude ||= begin
-          result = LONGITUDE_AND_DISTANCE_OF_MOON_PERIODIC_TERMS.inject(0.0) do |acc, elem|
+          result = MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS.inject(0.0) do |acc, elem|
             sine_argument = Angle.from_decimal_degrees(
-              elem["moon_mean_elongation"] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
-              elem["sun_mean_anomaly"] * sun_mean_anomaly2.decimal_degrees +
-              elem["moon_mean_anomaly"] * moon_mean_anomaly_high_precision.decimal_degrees +
-              elem["moon_argument_of_latitude"] * moon_argument_of_latitude_high_precision.decimal_degrees
+              elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
+              elem[1] * sun_mean_anomaly2.decimal_degrees +
+              elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
+              elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
             )
 
-            if elem["sine_coefficient"].nil?
+            if elem[4].nil?
               next acc
-            elsif [1, -1].include?(elem["sun_mean_anomaly"])
-              acc + elem["sine_coefficient"] * correction_eccentricity_of_earth * sine_argument.sin
-            elsif [-2, 2].include?(elem["sun_mean_anomaly"])
-              acc + elem["sine_coefficient"] * correction_eccentricity_of_earth**2 * sine_argument.sin
+            elsif [1, -1].include?(elem[1])
+              acc + elem[4] * correction_eccentricity_of_earth * sine_argument.sin
+            elsif [-2, 2].include?(elem[1])
+              acc + elem[4] * correction_eccentricity_of_earth**2 * sine_argument.sin
             else
-              acc + elem["sine_coefficient"] * sine_argument.sin
+              acc + elem[4] * sine_argument.sin
             end
           end + 3958 * correction_venus.sin +
             1962 * (moon_mean_longitude - moon_argument_of_latitude_high_precision).sin +
@@ -105,22 +105,22 @@ module Lunation
       # UNIT: decimal degrees
       def moon_heliocentric_latitude
         @moon_heliocentric_latitude ||= begin
-          result = LATITUDE_OF_MOON_PERIODIC_TERMS.inject(0.0) do |acc, elem|
+          result = MOON_LATITUDE_PERIOD_TERMS.inject(0.0) do |acc, elem|
             sine_argument = Angle.from_decimal_degrees(
-              elem["moon_mean_elongation"] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
-              elem["sun_mean_anomaly"] * sun_mean_anomaly2.decimal_degrees +
-              elem["moon_mean_anomaly"] * moon_mean_anomaly_high_precision.decimal_degrees +
-              elem["moon_argument_of_latitude"] * moon_argument_of_latitude_high_precision.decimal_degrees
+              elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
+              elem[1] * sun_mean_anomaly2.decimal_degrees +
+              elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
+              elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
             )
 
-            if elem["sine_coefficient"].nil?
+            if elem[4].nil?
               next acc
-            elsif [1, -1].include?(elem["sun_mean_anomaly"])
-              acc + elem["sine_coefficient"] * correction_eccentricity_of_earth * sine_argument.sin
-            elsif [-2, 2].include?(elem["sun_mean_anomaly"])
-              acc + elem["sine_coefficient"] * correction_eccentricity_of_earth**2 * sine_argument.sin
+            elsif [1, -1].include?(elem[1])
+              acc + elem[4] * correction_eccentricity_of_earth * sine_argument.sin
+            elsif [-2, 2].include?(elem[1])
+              acc + elem[4] * correction_eccentricity_of_earth**2 * sine_argument.sin
             else
-              acc + elem["sine_coefficient"] * sine_argument.sin
+              acc + elem[4] * sine_argument.sin
             end
           end - 2235 * moon_mean_longitude.sin +
             382 * correction_latitude.sin +
@@ -136,22 +136,22 @@ module Lunation
       # UNIT: 1000 km
       def moon_heliocentric_distance
         @moon_heliocentric_distance ||= begin
-          result = LONGITUDE_AND_DISTANCE_OF_MOON_PERIODIC_TERMS.inject(0.0) do |acc, elem|
+          result = MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS.inject(0.0) do |acc, elem|
             cosine_argument = Angle.from_decimal_degrees(
-              elem["moon_mean_elongation"] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
-              elem["sun_mean_anomaly"] * sun_mean_anomaly2.decimal_degrees +
-              elem["moon_mean_anomaly"] * moon_mean_anomaly_high_precision.decimal_degrees +
-              elem["moon_argument_of_latitude"] * moon_argument_of_latitude_high_precision.decimal_degrees
+              elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
+              elem[1] * sun_mean_anomaly2.decimal_degrees +
+              elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
+              elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
             )
 
-            if elem["cosine_coefficient"].nil?
+            if elem[5].nil?
               next acc
-            elsif [1, -1].include?(elem["sun_mean_anomaly"])
-              acc + elem["cosine_coefficient"] * correction_eccentricity_of_earth * cosine_argument.cos
-            elsif [-2, 2].include?(elem["sun_mean_anomaly"])
-              acc + elem["cosine_coefficient"] * correction_eccentricity_of_earth**2 * cosine_argument.cos
+            elsif [1, -1].include?(elem[1])
+              acc + elem[5] * correction_eccentricity_of_earth * cosine_argument.cos
+            elsif [-2, 2].include?(elem[1])
+              acc + elem[5] * correction_eccentricity_of_earth**2 * cosine_argument.cos
             else
-              acc + elem["cosine_coefficient"] * cosine_argument.cos
+              acc + elem[5] * cosine_argument.cos
             end
           end
           result.round
@@ -162,7 +162,8 @@ module Lunation
       # UNIT: Angle
       def moon_ecliptic_longitude
         @moon_ecliptic_longitude ||= Angle.from_decimal_degrees(
-          moon_mean_longitude.decimal_degrees + moon_heliocentric_longitude / 1_000_000.0
+          moon_mean_longitude.decimal_degrees +
+            moon_heliocentric_longitude.fdiv(1_000_000)
         )
       end
 
@@ -170,7 +171,7 @@ module Lunation
       # UNIT: Angle
       def moon_ecliptic_latitude
         @moon_ecliptic_latitude ||= Angle.from_decimal_degrees(
-          moon_heliocentric_latitude / 1_000_000.0,
+          moon_heliocentric_latitude.fdiv(1_000_000),
           normalize: false
         )
       end
@@ -180,7 +181,7 @@ module Lunation
       def distance_between_earth_and_moon
         @distance_between_earth_and_moon ||= begin
           result = FIXED_DISTANCE_BETWEEN_EARTH_AND_MOON +
-            (moon_heliocentric_distance / 1_000.0)
+            moon_heliocentric_distance.fdiv(1_000)
           result.round(1)
         end
       end
@@ -196,7 +197,8 @@ module Lunation
       # (apparent λ) Moon apparent longitude (A.A. p. 343)
       # UNIT: Angle
       def moon_apparent_ecliptic_longitude
-        @moon_apparent_ecliptic_longitude ||= moon_ecliptic_longitude + nutation_in_longitude
+        @moon_apparent_ecliptic_longitude ||=
+          moon_ecliptic_longitude + nutation_in_longitude
       end
 
       # (α) geocentric (apparent) right ascension of the moon (13.3 A.A. p. 93)