astronoby 0.7.0 → 0.8.0

data/lib/astronoby/data/constellations/sorted_declinations.dat

@@ -0,0 +1,202 @@
+-85.0
+-82.5
+-76.0
+-75.0
+-70.0
+-67.5
+-65.0
+-64.0
+-63.583333333333336
+-61.0
+-60.0
+-59.0
+-58.5
+-58.0
+-57.5
+-57.0
+-56.5
+-55.0
+-54.5
+-54.0
+-53.5
+-53.166666666666664
+-53.0
+-52.5
+-51.5
+-51.0
+-50.75
+-50.0
+-49.0
+-48.166666666666664
+-48.0
+-46.5
+-46.0
+-45.5
+-44.0
+-43.0
+-42.0
+-40.0
+-39.75
+-39.583333333333336
+-37.0
+-36.75
+-36.0
+-35.0
+-33.0
+-31.166666666666668
+-30.0
+-29.5
+-29.166666666666668
+-28.0
+-27.25
+-26.5
+-25.5
+-24.583333333333332
+-24.5
+-24.383333333333333
+-24.0
+-22.0
+-20.0
+-19.25
+-19.0
+-18.25
+-17.0
+-16.0
+-15.0
+-14.5
+-12.033333333333333
+-11.666666666666666
+-11.0
+-10.0
+-9.0
+-8.0
+-7.0
+-6.0
+-4.0
+-3.25
+-1.75
+0.0
+1.5
+1.75
+2.0
+2.75
+3.0
+4.0
+4.5
+5.5
+6.0
+6.25
+7.0
+7.5
+8.0
+8.5
+9.916666666666666
+10.0
+11.0
+11.833333333333334
+12.0
+12.5
+12.833333333333334
+13.5
+14.0
+14.333333333333334
+15.0
+15.5
+15.75
+16.0
+16.166666666666668
+17.5
+18.0
+18.5
+19.0
+19.166666666666668
+19.5
+20.0
+20.5
+21.0
+21.083333333333332
+21.25
+21.5
+22.0
+22.833333333333332
+23.5
+23.75
+25.0
+25.5
+26.0
+27.0
+27.25
+27.5
+28.0
+28.5
+29.0
+30.0
+30.666666666666668
+30.75
+31.333333333333332
+32.0
+32.083333333333336
+33.0
+33.5
+34.0
+34.5
+35.0
+35.5
+36.0
+36.5
+36.75
+39.75
+40.0
+42.0
+43.5
+43.75
+44.0
+44.5
+45.0
+46.0
+47.0
+47.5
+48.0
+48.5
+50.0
+50.5
+51.5
+52.5
+52.75
+53.0
+54.0
+54.833333333333336
+55.0
+55.5
+56.0
+56.25
+57.0
+57.5
+58.0
+58.5
+59.083333333333336
+59.5
+60.0
+60.916666666666664
+61.5
+62.0
+63.0
+64.0
+66.0
+66.5
+67.0
+68.0
+70.0
+73.5
+75.0
+77.0
+80.0
+82.0
+85.0
+86.0
+86.16666666666667
+86.5
+88.0
+89.0
+89.5
+90.0

data/lib/astronoby/data/constellations/sorted_right_ascensions.dat

@@ -0,0 +1,237 @@
+0.06666666666666667
+0.14166666666666666
+0.16666666666666666
+0.3333333333333333
+0.7166666666666667
+0.75
+0.85
+0.8666666666666667
+1.1166666666666667
+1.3333333333333333
+1.3666666666666667
+1.4083333333333334
+1.5833333333333333
+1.6666666666666667
+1.7
+1.8333333333333333
+1.9083333333333334
+1.9166666666666667
+2.0
+2.0416666666666665
+2.1666666666666665
+2.3333333333333335
+2.4166666666666665
+2.433333333333333
+2.5166666666666666
+2.566666666666667
+2.65
+2.6666666666666665
+2.716666666666667
+3.0
+3.1
+3.1666666666666665
+3.2
+3.283333333333333
+3.3333333333333335
+3.3666666666666667
+3.4166666666666665
+3.5
+3.5083333333333333
+3.5833333333333335
+3.75
+3.8333333333333335
+3.8666666666666667
+4.0
+4.083333333333333
+4.266666666666667
+4.333333333333333
+4.5
+4.583333333333333
+4.616666666666666
+4.666666666666667
+4.691666666666666
+4.7
+4.75
+4.833333333333333
+4.916666666666667
+4.966666666666667
+5.0
+5.083333333333333
+5.333333333333333
+5.5
+5.6
+5.7
+5.766666666666667
+5.833333333333333
+5.883333333333334
+6.0
+6.1
+6.116666666666666
+6.166666666666667
+6.216666666666667
+6.241666666666666
+6.308333333333334
+6.5
+6.533333333333333
+6.583333333333333
+6.8
+6.833333333333333
+6.933333333333334
+7.0
+7.016666666666667
+7.2
+7.366666666666666
+7.5
+7.666666666666667
+7.75
+7.808333333333334
+7.883333333333334
+7.925
+7.966666666666667
+8.0
+8.083333333333334
+8.166666666666666
+8.366666666666667
+8.416666666666666
+8.45
+8.583333333333334
+8.833333333333334
+9.033333333333333
+9.083333333333334
+9.166666666666666
+9.25
+9.366666666666667
+9.583333333333334
+9.75
+9.883333333333333
+10.166666666666666
+10.25
+10.5
+10.583333333333334
+10.666666666666666
+10.75
+10.783333333333333
+10.833333333333334
+11.0
+11.25
+11.333333333333334
+11.5
+11.516666666666667
+11.833333333333334
+11.866666666666667
+12.0
+12.083333333333334
+12.25
+12.333333333333334
+12.583333333333334
+12.833333333333334
+13.0
+13.25
+13.5
+13.583333333333334
+13.666666666666666
+13.958333333333334
+14.0
+14.033333333333333
+14.166666666666666
+14.25
+14.416666666666666
+14.5
+14.533333333333333
+14.666666666666666
+14.75
+14.916666666666666
+15.05
+15.083333333333334
+15.166666666666666
+15.183333333333334
+15.25
+15.333333333333334
+15.433333333333334
+15.666666666666666
+15.75
+15.916666666666666
+16.0
+16.033333333333335
+16.083333333333332
+16.166666666666668
+16.266666666666666
+16.333333333333332
+16.375
+16.420833333333334
+16.533333333333335
+16.583333333333332
+16.75
+16.833333333333332
+17.0
+17.166666666666668
+17.25
+17.5
+17.583333333333332
+17.6
+17.666666666666668
+17.833333333333332
+17.966666666666665
+18.0
+18.175
+18.233333333333334
+18.25
+18.366666666666667
+18.425
+18.583333333333332
+18.662222222222223
+18.866666666666667
+19.0
+19.083333333333332
+19.166666666666668
+19.25
+19.258333333333333
+19.358333333333334
+19.4
+19.416666666666668
+19.666666666666668
+19.766666666666666
+19.833333333333332
+20.0
+20.141666666666666
+20.166666666666668
+20.25
+20.3
+20.333333333333332
+20.416666666666668
+20.533333333333335
+20.536666666666665
+20.566666666666666
+20.6
+20.666666666666668
+20.833333333333332
+20.875
+20.916666666666668
+21.0
+21.05
+21.116666666666667
+21.25
+21.333333333333332
+21.416666666666668
+21.466666666666665
+21.666666666666668
+21.733333333333334
+21.866666666666667
+21.875
+21.908333333333335
+21.966666666666665
+22.0
+22.133333333333333
+22.316666666666666
+22.75
+22.816666666666666
+22.866666666666667
+23.0
+23.166666666666668
+23.333333333333332
+23.5
+23.583333333333332
+23.75
+23.833333333333332
+23.916666666666668
+24.0

data/lib/astronoby/equinox_solstice.rb

@@ -110,7 +110,7 @@ module Astronoby
     private
 
     def compute
-      t = (julian_day - Epoch::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
+      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
       w = Angle.from_degrees(35999.373 * t) - Angle.from_degrees(2.47)
       delta = 1 +
         0.0334 * w.cos +
@@ -123,7 +123,7 @@ module Astronoby
       delta_days = 0.00001 * s / delta
       epoch = julian_day + delta_days
 
-      Epoch.to_utc(epoch)
+      Instant.from_terrestrial_time(epoch).to_time
     end
 
     def julian_day

data/lib/astronoby/events/moon_phases.rb

@@ -15,17 +15,19 @@ module Astronoby
       # @param month [Integer] Requested month
       # @return [Array<Astronoby::MoonPhase>] List of Moon phases
       def self.phases_for(year:, month:)
-        [
-          MoonPhase.first_quarter(new(year, month, :first_quarter, -0.75).time),
-          MoonPhase.full_moon(new(year, month, :full_moon, -0.5).time),
-          MoonPhase.last_quarter(new(year, month, :last_quarter, -0.25).time),
-          MoonPhase.new_moon(new(year, month, :new_moon, 0).time),
-          MoonPhase.first_quarter(new(year, month, :first_quarter, 0.25).time),
-          MoonPhase.full_moon(new(year, month, :full_moon, 0.5).time),
-          MoonPhase.last_quarter(new(year, month, :last_quarter, 0.75).time),
-          MoonPhase.new_moon(new(year, month, :new_moon, 1).time),
-          MoonPhase.first_quarter(new(year, month, :first_quarter, 1.25).time)
-        ].select { _1.time.month == month }
+        Astronoby.cache.fetch([:moon_phases, year, month]) do
+          [
+            MoonPhase.first_quarter(new(year, month, :first_quarter, -0.75).time),
+            MoonPhase.full_moon(new(year, month, :full_moon, -0.5).time),
+            MoonPhase.last_quarter(new(year, month, :last_quarter, -0.25).time),
+            MoonPhase.new_moon(new(year, month, :new_moon, 0).time),
+            MoonPhase.first_quarter(new(year, month, :first_quarter, 0.25).time),
+            MoonPhase.full_moon(new(year, month, :full_moon, 0.5).time),
+            MoonPhase.last_quarter(new(year, month, :last_quarter, 0.75).time),
+            MoonPhase.new_moon(new(year, month, :new_moon, 1).time),
+            MoonPhase.first_quarter(new(year, month, :first_quarter, 1.25).time)
+          ].select { _1.time.month == month }
+        end
       end
 
       # @param year [Integer] Requested year
@@ -43,13 +45,12 @@ module Astronoby
       def time
         correction = moon_phases_periodic_terms
           .public_send(:"#{@phase}_correction")
-        terrestrial_time = Epoch.to_utc(
+        terrestrial_time = Instant.from_terrestrial_time(
           julian_ephemeris_day +
             correction +
             moon_phases_periodic_terms.additional_corrections
         )
-        delta = Util::Time.terrestrial_universal_time_delta(terrestrial_time)
-        (terrestrial_time - delta).round
+        terrestrial_time.to_time.round
       end
 
       private

data/lib/astronoby/events/rise_transit_set_calculator.rb

@@ -237,13 +237,14 @@ module Astronoby
         )
 
         # Calculate time differences
-        time_differences = new_instants.each_with_index.map do |instant, i|
-          Instant.from_terrestrial_time(instant.tt - old_instants[i].tt)
+        time_differences_in_days = new_instants.each_with_index.map do |instant, i|
+          instant.tt - old_instants[i].tt
         end
 
         # Calculate hour angle rate (radians per day)
         hour_angle_rates = hour_angle_changes.each_with_index.map do |angle, i|
-          angle.radians / time_differences[i].tt
+          denominator = time_differences_in_days[i]
+          angle.radians / denominator
         end
 
         # Store current values for next iteration
@@ -255,7 +256,8 @@ module Astronoby
           .each_with_index
           .map do |angle, i|
             ratio = angle.radians / hour_angle_rates[i]
-            [ratio.nan? ? 0 : ratio, MIN_TIME_ADJUSTMENT].max
+            time_adjustment = (ratio.nan? || ratio.infinite?) ? 0 : ratio
+            [time_adjustment, MIN_TIME_ADJUSTMENT].max
           end
 
         # Apply time adjustments
@@ -268,10 +270,31 @@ module Astronoby
     end
 
     def calculate_positions_at_instants(instants)
-      instants.map do |instant|
-        @body
-          .new(instant: instant, ephem: @ephem)
-          .observed_by(@observer)
+      if Astronoby.configuration.cache_enabled?
+        instants.map do |instant|
+          cache_key = CacheKey.generate(
+            :observed_by,
+            instant,
+            @body.to_s,
+            @observer.hash
+          )
+          Astronoby.cache.fetch(cache_key) do
+            @body
+              .new(instant: instant, ephem: @ephem)
+              .observed_by(@observer)
+          end
+        end
+      else
+        @positions_cache ||= {}
+        precision = Astronoby.configuration.cache_precision(:observed_by)
+        instants.map do |instant|
+          rounded_instant = Instant.from_terrestrial_time(
+            instant.tt.round(precision)
+          )
+          @positions_cache[rounded_instant] ||= @body
+            .new(instant: rounded_instant, ephem: @ephem)
+            .observed_by(@observer)
+        end
       end
     end
 
@@ -347,6 +370,7 @@ module Astronoby
       @sample_instants = nil
       @start_instant = nil
       @end_instant = nil
+      @positions_cache = nil
     end
   end
 end