astronoby 0.3.0 → 0.5.0

data/lib/astronoby/angle.rb

@@ -1,23 +1,10 @@
 # frozen_string_literal: true
 
-require "bigdecimal/math"
-
 module Astronoby
   class Angle
     include Comparable
 
-    PRECISION = 14
-    PI = BigMath.PI(PRECISION)
-    PI_IN_DEGREES = BigDecimal("180")
-
-    FULL_CIRCLE_IN_RADIANS = (2 * PI)
-
-    RADIAN_PER_HOUR = PI / BigDecimal("12")
-    MINUTES_PER_DEGREE = BigDecimal("60")
-    MINUTES_PER_HOUR = BigDecimal("60")
-    SECONDS_PER_MINUTE = BigDecimal("60")
-    SECONDS_PER_HOUR = MINUTES_PER_HOUR * SECONDS_PER_MINUTE
-
+    MIN_PRECISION = 10
     FORMATS = %i[dms hms].freeze
 
     class << self
@@ -26,28 +13,32 @@ module Astronoby
       end
 
       def from_radians(radians)
-        normalized_radians = radians.remainder(FULL_CIRCLE_IN_RADIANS)
+        normalized_radians = radians.remainder(Constants::RADIANS_PER_CIRCLE)
         new(normalized_radians)
       end
 
       def from_degrees(degrees)
-        radians = degrees / PI_IN_DEGREES * PI
+        radians = degrees / Constants::PI_IN_DEGREES * Math::PI
         from_radians(radians)
       end
 
       def from_hours(hours)
-        radians = hours * RADIAN_PER_HOUR
+        radians = hours * Constants::RADIAN_PER_HOUR
         from_radians(radians)
       end
 
       def from_hms(hour, minute, second)
-        hours = hour + minute / MINUTES_PER_HOUR + second / SECONDS_PER_HOUR
+        hours = hour +
+          minute / Constants::MINUTES_PER_HOUR +
+          second / Constants::SECONDS_PER_HOUR
         from_hours(hours)
       end
 
       def from_dms(degree, minute, second)
         sign = degree.negative? ? -1 : 1
-        degrees = degree.abs + minute / MINUTES_PER_HOUR + second / SECONDS_PER_HOUR
+        degrees = degree.abs +
+          minute / Constants::MINUTES_PER_HOUR +
+          second / Constants::SECONDS_PER_HOUR
         from_degrees(sign * degrees)
       end
 
@@ -70,20 +61,16 @@ module Astronoby
     attr_reader :radians
 
     def initialize(radians)
-      @radians = if radians.is_a?(Integer) || radians.is_a?(BigDecimal)
-        BigDecimal(radians)
-      else
-        BigDecimal(radians, PRECISION)
-      end
+      @radians = radians
       freeze
     end
 
     def degrees
-      @radians * PI_IN_DEGREES / PI
+      @radians * Constants::PI_IN_DEGREES / Math::PI
     end
 
     def hours
-      @radians / RADIAN_PER_HOUR
+      @radians / Constants::RADIAN_PER_HOUR
     end
 
     def +(other)
@@ -94,6 +81,10 @@ module Astronoby
       self.class.from_radians(@radians - other.radians)
     end
 
+    def -@
+      self.class.from_radians(-@radians)
+    end
+
     def sin
       Math.sin(radians)
     end
@@ -144,31 +135,32 @@ module Astronoby
       sign = deg.negative? ? "-" : "+"
       absolute_degrees = deg.abs
       degrees = absolute_degrees.floor
-      decimal_minutes = MINUTES_PER_DEGREE * (absolute_degrees - degrees)
+      decimal_minutes = Constants::MINUTES_PER_DEGREE *
+        (absolute_degrees - degrees)
       absolute_decimal_minutes = (
-        MINUTES_PER_DEGREE * (absolute_degrees - degrees)
+        Constants::MINUTES_PER_DEGREE * (absolute_degrees - degrees)
       ).abs
       minutes = decimal_minutes.floor
-      seconds = SECONDS_PER_MINUTE * (
+      seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Dms.new(sign, degrees, minutes, seconds.to_f.floor(4))
+      Dms.new(sign, degrees, minutes, seconds.floor(4))
     end
 
     def to_hms(hrs)
       absolute_hours = hrs.abs
       hours = absolute_hours.floor
-      decimal_minutes = MINUTES_PER_HOUR * (absolute_hours - hours)
+      decimal_minutes = Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
       absolute_decimal_minutes = (
-        MINUTES_PER_HOUR * (absolute_hours - hours)
+        Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
       ).abs
       minutes = decimal_minutes.floor
-      seconds = SECONDS_PER_MINUTE * (
+      seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Hms.new(hours, minutes, seconds.to_f.floor(4))
+      Hms.new(hours, minutes, seconds.floor(4))
     end
   end
 end

data/lib/astronoby/astronomical_models/ephemeride_lunaire_parisienne.rb

@@ -0,0 +1,143 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class EphemerideLunaireParisienne
+    # Lunar theory: Éphéméride Lunaire Parisienne
+    # By Jean Chapront, Michelle Chapront-Touzé
+    # https://en.wikipedia.org/wiki/Ephemeride_Lunaire_Parisienne
+
+    DEGREES_UNIT = 10**-6
+
+    # @return [Array] Periodic terms for the Moon's longitude and distance
+    def self.periodic_terms_for_moon_longitude_and_distance
+      [
+        [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, 0],
+        [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, 0],
+        [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, 0],
+        [2, -2, -1, 0, 2048, -4950],
+        [2, 0, 1, -2, -1773, 4130],
+        [2, 0, 0, 2, -1595, 0],
+        [4, -1, -1, 0, 1215, -3958],
+        [0, 0, 2, 2, -1110, 0],
+        [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, 0],
+        [2, -1, 0, -2, 596, 0],
+        [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, 0],
+        [0, 0, 2, -2, -381, -4421],
+        [1, 1, 1, 0, 351, 0],
+        [3, 0, -2, 0, -340, 0],
+        [4, 0, -3, 0, 330, 0],
+        [2, -1, 2, 0, 327, 0],
+        [0, 2, 1, 0, -323, 1165],
+        [1, 1, -1, 0, 299, 0],
+        [2, 0, 3, 0, 294, 0],
+        [2, 0, -1, -2, 0, 8752]
+      ]
+    end
+
+    # @return [Array] Periodic terms for the Moon's latitude
+    def self.periodic_terms_for_moon_latitude
+      [
+        [0, 0, 0, 1, 5128122],
+        [0, 0, 1, 1, 280602],
+        [0, 0, 1, -1, 277693],
+        [2, 0, 0, -1, 173237],
+        [2, 0, -1, 1, 55413],
+        [2, 0, -1, -1, 46271],
+        [2, 0, 0, 1, 32573],
+        [0, 0, 2, 1, 17198],
+        [2, 0, 1, -1, 9266],
+        [0, 0, 2, -1, 8822],
+        [2, -1, 0, -1, 8216],
+        [2, 0, -2, -1, 4324],
+        [2, 0, 1, 1, 4200],
+        [2, 1, 0, -1, -3359],
+        [2, -1, -1, 1, 2463],
+        [2, -1, 0, 1, 2211],
+        [2, -1, -1, -1, 2065],
+        [0, 1, -1, -1, -1870],
+        [4, 0, -1, -1, 1828],
+        [0, 1, 0, 1, -1794],
+        [0, 0, 0, 3, -1749],
+        [0, 1, -1, 1, -1565],
+        [1, 0, 0, 1, -1491],
+        [0, 1, 1, 1, -1475],
+        [0, 1, 1, -1, -1410],
+        [0, 1, 0, -1, -1344],
+        [1, 0, 0, -1, -1335],
+        [0, 0, 3, 1, 1107],
+        [4, 0, 0, -1, 1021],
+        [4, 0, -1, 1, 833],
+        [0, 0, 1, -3, 777],
+        [4, 0, -2, 1, 671],
+        [2, 0, 0, -3, 607],
+        [2, 0, 2, -1, 596],
+        [2, -1, 1, -1, 491],
+        [2, 0, -2, 1, -451],
+        [0, 0, 3, -1, 439],
+        [2, 0, 2, 1, 422],
+        [2, 0, -3, -1, 421],
+        [2, 1, -1, 1, -366],
+        [2, 1, 0, 1, -351],
+        [4, 0, 0, 1, 331],
+        [2, -1, 1, 1, 315],
+        [2, -2, 0, -1, 302],
+        [0, 0, 1, 3, -283],
+        [2, 1, 1, -1, -229],
+        [1, 1, 0, -1, 223],
+        [1, 1, 0, 1, 223],
+        [0, 1, -2, -1, -220],
+        [2, 1, -1, -1, -220],
+        [1, 0, 1, 1, -185],
+        [2, -1, -2, -1, 181],
+        [0, 1, 2, 1, -177],
+        [4, 0, -2, -1, 176],
+        [4, -1, -1, -1, 166],
+        [1, 0, 1, -1, -164],
+        [4, 0, 1, -1, 132],
+        [1, 0, -1, -1, -119],
+        [4, -1, 0, -1, 115],
+        [2, -2, 0, 1, 107]
+      ]
+    end
+  end
+end

data/lib/astronoby/astronomical_models/moon_phases_periodic_terms.rb

@@ -0,0 +1,249 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class MoonPhasesPeriodicTerms
+    # Lunar Solution ELP 2000-82B (Chapront-Touze+, 1988)
+
+    # Source:
+    #  Title: Astronomical Algorithms
+    #  Author: Jean Meeus
+    #  Edition: 2nd edition
+    #  Chapter: 49 - Phases of the Moon
+
+    # @param julian_centuries [Float] Julian centuries
+    # @param time [Float] Time
+    # @param eccentricity_correction [Float] Eccentricity correction
+    # @param moon_mean_anomaly [Float] Moon mean anomaly
+    # @param sun_mean_anomaly [Float] Sun mean anomaly
+    # @param moon_argument_of_latitude [Float] Moon argument of latitude
+    # @param longitude_of_the_ascending_node [Float] Longitude of the ascending node
+    def initialize(
+      julian_centuries:,
+      time:,
+      eccentricity_correction:,
+      moon_mean_anomaly:,
+      sun_mean_anomaly:,
+      moon_argument_of_latitude:,
+      longitude_of_the_ascending_node:
+
+    )
+      @julian_centuries = julian_centuries
+      @time = time
+      @eccentricity_correction = eccentricity_correction
+      @moon_mean_anomaly = moon_mean_anomaly
+      @sun_mean_anomaly = sun_mean_anomaly
+      @moon_argument_of_latitude = moon_argument_of_latitude
+      @longitude_of_the_ascending_node = longitude_of_the_ascending_node
+    end
+
+    # @return [Float] New moon correction
+    def new_moon_correction
+      ecc = @eccentricity_correction
+      mma = @moon_mean_anomaly.radians
+      sma = @sun_mean_anomaly.radians
+      maol = @moon_argument_of_latitude.radians
+      lotan = @longitude_of_the_ascending_node.radians
+      [
+        [-0.40720, mma],
+        [0.17241 * ecc, sma],
+        [0.01608, 2 * mma],
+        [0.01039, 2 * maol],
+        [0.00739 * ecc, mma - sma],
+        [-0.00514 * ecc, mma + sma],
+        [0.00208 * ecc * ecc, 2 * sma],
+        [-0.00111, mma - 2 * maol],
+        [-0.00057, mma + 2 * maol],
+        [0.00056 * ecc, 2 * mma + sma],
+        [-0.00042, 3 * mma],
+        [0.00042 * ecc, sma + 2 * maol],
+        [0.00038 * ecc, sma - 2 * maol],
+        [-0.00024 * ecc, 2 * mma - sma],
+        [-0.00017, lotan],
+        [-0.00007, mma + 2 * sma],
+        [0.00004, 2 * mma - 2 * maol],
+        [0.00004, 3 * sma],
+        [0.00003, mma + sma - 2 * maol],
+        [0.00003, 2 * mma + 2 * maol],
+        [-0.00003, mma + sma + 2 * maol],
+        [0.00003, mma - sma + 2 * maol],
+        [-0.00002, mma - sma - 2 * maol],
+        [-0.00002, 3 * mma + sma],
+        [0.00002, 4 * mma]
+      ].map { _1.first * Math.sin(_1.last) }.sum
+    end
+
+    # @return [Float] First quarter correction
+    def first_quarter_correction
+      first_and_last_quarter_correction + first_and_last_quarter_final_correction
+    end
+
+    # @return [Float] Full moon correction
+    def full_moon_correction
+      ecc = @eccentricity_correction
+      mma = @moon_mean_anomaly.radians
+      sma = @sun_mean_anomaly.radians
+      maol = @moon_argument_of_latitude.radians
+      lotan = @longitude_of_the_ascending_node.radians
+      [
+        [-0.40614, mma],
+        [0.17302 * ecc, sma],
+        [0.01614, 2 * mma],
+        [0.01043, 2 * maol],
+        [0.00734 * ecc, mma - sma],
+        [-0.00515 * ecc, mma + sma],
+        [0.00209 * ecc * ecc, 2 * sma],
+        [-0.00111, mma - 2 * maol],
+        [-0.00057, mma + 2 * maol],
+        [0.00056 * ecc, 2 * mma + sma],
+        [-0.00042, 3 * mma],
+        [0.00042 * ecc, sma + 2 * maol],
+        [0.00038 * ecc, sma - 2 * maol],
+        [-0.00024 * ecc, 2 * mma - sma],
+        [-0.00017, lotan],
+        [-0.00007, mma + 2 * sma],
+        [0.00004, 2 * mma - 2 * maol],
+        [0.00004, 3 * sma],
+        [0.00003, mma + sma - 2 * maol],
+        [0.00003, 2 * mma + 2 * maol],
+        [-0.00003, mma + sma + 2 * maol],
+        [0.00003, mma - sma + 2 * maol],
+        [-0.00002, mma - sma - 2 * maol],
+        [-0.00002, 3 * mma + sma],
+        [0.00002, 4 * mma]
+      ].map { _1.first * Math.sin(_1.last) }.sum
+    end
+
+    # @return [Float] Last quarter correction
+    def last_quarter_correction
+      first_and_last_quarter_correction - first_and_last_quarter_final_correction
+    end
+
+    # @return [Float] Additional corrections
+    def additional_corrections
+      [
+        [0.000325, a1],
+        [0.000165, a2],
+        [0.000164, a3],
+        [0.000126, a4],
+        [0.000110, a5],
+        [0.000062, a6],
+        [0.000060, a7],
+        [0.000056, a8],
+        [0.000047, a9],
+        [0.000042, a10],
+        [0.000040, a11],
+        [0.000037, a12],
+        [0.000035, a13],
+        [0.000023, a14]
+      ].map { _1.first * _1.last.sin }.sum
+    end
+
+    private
+
+    def a1
+      Angle.from_degrees(
+        (
+          299.77 +
+            0.107408 * @time -
+            0.009173 * @julian_centuries**2
+        ) % 360
+      )
+    end
+
+    def a2
+      Angle.from_degrees (251.88 + 0.016321 * @time) % 360
+    end
+
+    def a3
+      Angle.from_degrees (251.83 + 26.651886 * @time) % 360
+    end
+
+    def a4
+      Angle.from_degrees (349.42 + 36.412478 * @time) % 360
+    end
+
+    def a5
+      Angle.from_degrees (84.66 + 18.206239 * @time) % 360
+    end
+
+    def a6
+      Angle.from_degrees (141.74 + 53.303771 * @time) % 360
+    end
+
+    def a7
+      Angle.from_degrees (207.14 + 2.453732 * @time) % 360
+    end
+
+    def a8
+      Angle.from_degrees (154.84 + 7.306860 * @time) % 360
+    end
+
+    def a9
+      Angle.from_degrees (34.52 + 27.261239 * @time) % 360
+    end
+
+    def a10
+      Angle.from_degrees (207.19 + 0.121824 * @time) % 360
+    end
+
+    def a11
+      Angle.from_degrees (291.34 + 1.844379 * @time) % 360
+    end
+
+    def a12
+      Angle.from_degrees (161.72 + 24.198154 * @time) % 360
+    end
+
+    def a13
+      Angle.from_degrees (239.56 + 25.513099 * @time) % 360
+    end
+
+    def a14
+      Angle.from_degrees (331.55 + 3.592518 * @time) % 360
+    end
+
+    def first_and_last_quarter_final_correction
+      0.00306 -
+        0.00038 * @eccentricity_correction * @sun_mean_anomaly.cos +
+        0.00026 * @moon_mean_anomaly.cos -
+        0.00002 * (@moon_mean_anomaly - @sun_mean_anomaly).cos +
+        0.00002 * (@moon_mean_anomaly + @sun_mean_anomaly).cos +
+        0.00002 * Math.cos(2 * @moon_argument_of_latitude.radians)
+    end
+
+    def first_and_last_quarter_correction
+      ecc = @eccentricity_correction
+      mma = @moon_mean_anomaly.radians
+      sma = @sun_mean_anomaly.radians
+      maol = @moon_argument_of_latitude.radians
+      lotan = @longitude_of_the_ascending_node.radians
+      [
+        [-0.62801, mma],
+        [0.17172 * ecc, sma],
+        [-0.01183, mma + sma],
+        [0.00862, 2 * mma],
+        [0.00804, 2 * maol],
+        [0.00454 * ecc, mma - sma],
+        [0.00204 * ecc**2, 2 * sma],
+        [-0.00180, mma - 2 * maol],
+        [-0.00070, mma + 2 * maol],
+        [-0.00040, 3 * mma],
+        [-0.00034 * ecc, 2 * mma - sma],
+        [0.00032 * ecc, sma + 2 * maol],
+        [0.00032, sma - 2 * maol],
+        [-0.00028 * ecc**2, mma + 2 * sma],
+        [0.00027 * ecc, 2 * mma + sma],
+        [-0.00017, lotan],
+        [-0.00005, mma - sma - 2 * maol],
+        [0.00004, 2 * mma + 2 * maol],
+        [-0.00004, mma + sma + 2 * maol],
+        [0.00004, mma - 2 * sma],
+        [0.00003, mma + sma - 2 * maol],
+        [0.00003, 3 * sma],
+        [0.00002, 2 * mma - 2 * maol],
+        [0.00002, mma - sma + 2 * maol],
+        [-0.00002, 3 * mma + sma]
+      ].map { _1.first * Math.sin(_1.last) }.sum
+    end
+  end
+end