lunation 0.1.4 → 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

@@ -2,19 +2,32 @@ module Lunation
   class Calculation
     module EarthPositionVSOP87
       # rubocop:disable Layout/LineLength
-      PERIODIC_TERMS_B0 = YAML.load_file("config/periodic_terms_earth_position_b0.yml").freeze
-      PERIODIC_TERMS_B1 = YAML.load_file("config/periodic_terms_earth_position_b1.yml").freeze
-      PERIODIC_TERMS_L0 = YAML.load_file("config/periodic_terms_earth_position_l0.yml").freeze
-      PERIODIC_TERMS_L1 = YAML.load_file("config/periodic_terms_earth_position_l1.yml").freeze
-      PERIODIC_TERMS_L2 = YAML.load_file("config/periodic_terms_earth_position_l2.yml").freeze
-      PERIODIC_TERMS_L3 = YAML.load_file("config/periodic_terms_earth_position_l3.yml").freeze
-      PERIODIC_TERMS_L4 = YAML.load_file("config/periodic_terms_earth_position_l4.yml").freeze
-      PERIODIC_TERMS_L5 = YAML.load_file("config/periodic_terms_earth_position_l5.yml").freeze
-      PERIODIC_TERMS_R0 = YAML.load_file("config/periodic_terms_earth_position_r0.yml").freeze
-      PERIODIC_TERMS_R1 = YAML.load_file("config/periodic_terms_earth_position_r1.yml").freeze
-      PERIODIC_TERMS_R2 = YAML.load_file("config/periodic_terms_earth_position_r2.yml").freeze
-      PERIODIC_TERMS_R3 = YAML.load_file("config/periodic_terms_earth_position_r3.yml").freeze
-      PERIODIC_TERMS_R4 = YAML.load_file("config/periodic_terms_earth_position_r4.yml").freeze
+      PERIODIC_TERMS_B0_PATH = File.expand_path("../../../config/periodic_terms_earth_position_b0.yml", __dir__).freeze
+      PERIODIC_TERMS_B1_PATH = File.expand_path("../../../config/periodic_terms_earth_position_b1.yml", __dir__).freeze
+      PERIODIC_TERMS_L0_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l0.yml", __dir__).freeze
+      PERIODIC_TERMS_L1_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l1.yml", __dir__).freeze
+      PERIODIC_TERMS_L2_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l2.yml", __dir__).freeze
+      PERIODIC_TERMS_L3_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l3.yml", __dir__).freeze
+      PERIODIC_TERMS_L4_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l4.yml", __dir__).freeze
+      PERIODIC_TERMS_L5_PATH = File.expand_path("../../../config/periodic_terms_earth_position_l5.yml", __dir__).freeze
+      PERIODIC_TERMS_R0_PATH = File.expand_path("../../../config/periodic_terms_earth_position_r0.yml", __dir__).freeze
+      PERIODIC_TERMS_R1_PATH = File.expand_path("../../../config/periodic_terms_earth_position_r1.yml", __dir__).freeze
+      PERIODIC_TERMS_R2_PATH = File.expand_path("../../../config/periodic_terms_earth_position_r2.yml", __dir__).freeze
+      PERIODIC_TERMS_R3_PATH = File.expand_path("../../../config/periodic_terms_earth_position_r3.yml", __dir__).freeze
+      PERIODIC_TERMS_R4_PATH = File.expand_path("../../../config/periodic_terms_earth_position_r4.yml", __dir__).freeze
+      PERIODIC_TERMS_B0 = YAML.load_file(PERIODIC_TERMS_B0_PATH).freeze
+      PERIODIC_TERMS_B1 = YAML.load_file(PERIODIC_TERMS_B1_PATH).freeze
+      PERIODIC_TERMS_L0 = YAML.load_file(PERIODIC_TERMS_L0_PATH).freeze
+      PERIODIC_TERMS_L1 = YAML.load_file(PERIODIC_TERMS_L1_PATH).freeze
+      PERIODIC_TERMS_L2 = YAML.load_file(PERIODIC_TERMS_L2_PATH).freeze
+      PERIODIC_TERMS_L3 = YAML.load_file(PERIODIC_TERMS_L3_PATH).freeze
+      PERIODIC_TERMS_L4 = YAML.load_file(PERIODIC_TERMS_L4_PATH).freeze
+      PERIODIC_TERMS_L5 = YAML.load_file(PERIODIC_TERMS_L5_PATH).freeze
+      PERIODIC_TERMS_R0 = YAML.load_file(PERIODIC_TERMS_R0_PATH).freeze
+      PERIODIC_TERMS_R1 = YAML.load_file(PERIODIC_TERMS_R1_PATH).freeze
+      PERIODIC_TERMS_R2 = YAML.load_file(PERIODIC_TERMS_R2_PATH).freeze
+      PERIODIC_TERMS_R3 = YAML.load_file(PERIODIC_TERMS_R3_PATH).freeze
+      PERIODIC_TERMS_R4 = YAML.load_file(PERIODIC_TERMS_R4_PATH).freeze
       # rubocop:enable Layout/LineLength
 
       # (L) Ecliptical longitude of the earth (A.A. p. 219, 32.2)
@@ -45,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))