astronoby 0.10.0 → 0.11.0

data/lib/astronoby/horizontal_transformation.rb

@@ -0,0 +1,112 @@
+# frozen_string_literal: true
+
+module Astronoby
+  # Transforms equatorial coordinates at J2000.0 into horizontal coordinates
+  # for one observer at one instant.
+  #
+  # Everything that depends only on the observer and the instant is computed
+  # once at construction: the rotation from the ICRS to the true equator and
+  # equinox of date (nutation * precession, with the ICRS frame bias), and the
+  # local apparent sidereal time. Converting a point then costs a matrix
+  # product and a few trigonometric functions, instead of rebuilding the
+  # reference frame chain for each point. It is meant for converting many
+  # points for the same observer and instant, such as the stars of a catalogue
+  # for a sky chart.
+  #
+  # It plays the role of the "star-independent astrometry parameters" of SOFA
+  # and ERFA, restricted to what Astronoby applies without an ephemeris. The
+  # result is the same as
+  # `DeepSkyObject#at(instant).observed_by(observer).horizontal`, but the
+  # points are treated as fixed directions: no aberration, no proper motion and
+  # no parallax. The topocentric correction is not applied either, as it is
+  # nil for objects at stellar distances. Atmospheric refraction is optional.
+  #
+  # @example Convert the stars of a catalogue for one observer and instant
+  #   transformation = Astronoby::HorizontalTransformation.new(
+  #     observer: observer,
+  #     instant: instant
+  #   )
+  #   horizontals = catalogue.map { |star| transformation.apply(star) }
+  class HorizontalTransformation
+    # @return [Astronoby::Observer] the observer
+    attr_reader :observer
+
+    # @return [Astronoby::Instant] the instant
+    attr_reader :instant
+
+    # @param observer [Astronoby::Observer] the observer
+    # @param instant [Astronoby::Instant] the instant of the observation
+    def initialize(observer:, instant:)
+      @observer = observer
+      @instant = instant
+      @rotation_rows = (
+        Nutation.matrix_for(instant) * Precession.matrix_for(instant)
+      ).to_a.freeze
+      @local_apparent_sidereal_time = LocalApparentSiderealTime.from_utc(
+        instant.to_time,
+        longitude: observer.longitude
+      ).time
+      freeze
+    end
+
+    # @param equatorial_coordinates [Astronoby::Coordinates::Equatorial]
+    #   equatorial coordinates at epoch J2000.0
+    # @param refraction [Boolean] whether to apply atmospheric refraction
+    #   correction (default: false)
+    # @return [Astronoby::Coordinates::Horizontal] horizontal coordinates
+    # @raise [Astronoby::IncompatibleArgumentsError] if the coordinates are not
+    #   referred to epoch J2000.0
+    def apply(equatorial_coordinates, refraction: false)
+      ensure_j2000!(equatorial_coordinates)
+
+      horizontal = apparent_equatorial(equatorial_coordinates).to_horizontal(
+        time: @instant.to_time,
+        observer: @observer
+      )
+
+      if refraction
+        Refraction.correct_horizontal_coordinates(coordinates: horizontal)
+      else
+        horizontal
+      end
+    end
+
+    private
+
+    # Equatorial coordinates of date, with the hour angle already known so that
+    # the sidereal time is not computed again.
+    def apparent_equatorial(coordinates)
+      declination = coordinates.declination
+      right_ascension = coordinates.right_ascension
+      cos_declination = declination.cos
+      x = cos_declination * right_ascension.cos
+      y = cos_declination * right_ascension.sin
+      z = declination.sin
+
+      row_x, row_y, row_z = @rotation_rows
+      apparent_x = row_x[0] * x + row_x[1] * y + row_x[2] * z
+      apparent_y = row_y[0] * x + row_y[1] * y + row_y[2] * z
+      apparent_z = row_z[0] * x + row_z[1] * y + row_z[2] * z
+
+      apparent_right_ascension =
+        Angle.from_radians(Math.atan2(apparent_y, apparent_x))
+
+      Coordinates::Equatorial.new(
+        right_ascension: apparent_right_ascension,
+        declination: Angle.asin(apparent_z.clamp(-1.0, 1.0)),
+        hour_angle:
+          Angle.from_hours(@local_apparent_sidereal_time) -
+            apparent_right_ascension,
+        epoch: @instant.tt
+      )
+    end
+
+    def ensure_j2000!(coordinates)
+      return if coordinates.epoch == JulianDate::J2000
+
+      raise IncompatibleArgumentsError,
+        "Expected equatorial coordinates at epoch J2000.0, got epoch " \
+        "#{coordinates.epoch}"
+    end
+  end
+end

data/lib/astronoby/instant.rb

@@ -1,6 +1,7 @@
 # frozen_string_literal: true
 
 require "date"
+require "horologium"
 
 module Astronoby
   # Represents a specific instant in time using Terrestrial Time (TT) as its
@@ -8,9 +9,15 @@ module Astronoby
   # time scales commonly used in astronomy:
   # - Terrestrial Time (TT)
   # - International Atomic Time (TAI)
+  # - Barycentric Dynamical Time (TDB)
   # - Universal Time Coordinated (UTC)
   # - Greenwich Mean Sidereal Time (GMST)
   #
+  # TAI and TDB are delegated to the horologium gem, which owns the scales
+  # that convert by definition or by model. UT1 is not: it rests on
+  # observations of the rotation of the Earth, and comes from
+  # {Astronoby::DeltaT}.
+  #
   # @example Create an instant from the current time
   #   instant = Astronoby::Instant.from_time(Time.now)
   #   instant.tai  # Get International Atomic Time
@@ -28,6 +35,16 @@ module Astronoby
     # midnight. This constant adds 0.5 days (12 hours) to make the conversion.
     DATETIME_JD_EPOCH_ADJUSTMENT = 0.5
 
+    # Astronoby pins the precision it asks horologium for rather than reading
+    # horologium's global default, so a host application configuring
+    # horologium for itself cannot change what Astronoby computes.
+    #
+    # Exact rather than standard: standard would round the Julian Date into a
+    # two-part float on the way in, and TAI is a fixed offset from TT, so
+    # there is nothing to round it for. It is also the cheaper of the two
+    # here, having no conversion to make in either direction.
+    PRECISION = :exact
+
     class << self
       # Creates a new Instant from a Terrestrial Time value
       #
@@ -44,7 +61,7 @@ module Astronoby
       # @param time [Time] a Time object to convert
       # @return [Astronoby::Instant] a new Instant object
       def from_time(time)
-        delta_t = Util::Time.terrestrial_universal_time_delta(time)
+        delta_t = DeltaT.at(time)
         terrestrial_time = time.utc.to_datetime.ajd +
           Rational(delta_t, Constants::SECONDS_PER_DAY)
         from_terrestrial_time(terrestrial_time)
@@ -55,7 +72,7 @@ module Astronoby
       # @param julian_date [Numeric] the Julian Date in UTC
       # @return [Astronoby::Instant] a new Instant object
       def from_utc_julian_date(julian_date)
-        delta_t = Util::Time.terrestrial_universal_time_delta(julian_date)
+        delta_t = DeltaT.at(julian_date)
         terrestrial_time = julian_date +
           Rational(delta_t, Constants::SECONDS_PER_DAY)
         from_terrestrial_time(terrestrial_time)
@@ -120,8 +137,7 @@ module Astronoby
     #
     # @return [Numeric] Delta T in seconds
     def delta_t
-      @memo[:delta_t] ||=
-        Util::Time.terrestrial_universal_time_delta(@terrestrial_time)
+      @memo[:delta_t] ||= DeltaT.at(@terrestrial_time)
     end
 
     # Get the Greenwich Mean Sidereal Time
@@ -156,21 +172,26 @@ module Astronoby
 
     # Get the International Atomic Time (TAI)
     #
-    # @return [Numeric] TAI as Julian Date
+    # Read as a Rational rather than a Float. A Julian Date is around 2.46
+    # million, so a single Float rounds to about 47 microseconds of the day,
+    # which would coarsen an instant that TT carries exactly.
+    #
+    # @return [Rational] TAI as Julian Date
     def tai
-      @memo[:tai] ||= @terrestrial_time -
-        Rational(Constants::TAI_TT_OFFSET, Constants::SECONDS_PER_DAY)
+      @memo[:tai] ||= scale_free_instant
+        .as(:julian_date, scale: :tai, as: :rational)
     end
 
     # Get the Barycentric Dynamical Time (TDB)
-    # Note: Currently approximated as equal to TT
     #
-    # @return [Numeric] TDB as Julian Date
+    # TDB - TT comes from the Fairhead and Bretagnon series ERFA uses. It
+    # peaks around 1.7 milliseconds, so a Float Julian Date would quantise
+    # most of the correction away; see {#tai} on the shape.
+    #
+    # @return [Rational] TDB as Julian Date
     def tdb
-      # This is technically false, there is a slight difference between TT and
-      # TDB. However, this difference is so small that currenly Astronoby
-      # doesn't support it and consider they are the same value.
-      @memo[:tdb] ||= @terrestrial_time
+      @memo[:tdb] ||= scale_free_instant
+        .as(:julian_date, scale: :tdb, as: :rational)
     end
 
     # Get the offset between TT and UTC for this instant
@@ -198,5 +219,21 @@ module Astronoby
       @terrestrial_time <=> other.terrestrial_time
     end
     alias_method :eql?, :==
+
+    private
+
+    # The instant as horologium holds it: a point on the timeline with no
+    # scale of its own. Building it costs more than reading a Julian Date back
+    # out of it, so it is built on the first scale conversion that needs it
+    # and never on the paths that only want TT.
+    #
+    # @return [Horologium::Instant]
+    def scale_free_instant
+      @memo[:scale_free_instant] ||= Horologium::Instant.from_julian_date(
+        @terrestrial_time,
+        scale: :tt,
+        precision: PRECISION
+      )
+    end
   end
 end

data/lib/astronoby/precession.rb

@@ -163,7 +163,7 @@ module Astronoby
 
     def t
       @t ||= Rational(
-        @instant.tdb - JulianDate::DEFAULT_EPOCH,
+        @instant.tt - JulianDate::DEFAULT_EPOCH,
         Constants::DAYS_PER_JULIAN_CENTURY
       )
     end

data/lib/astronoby/reference_frame.rb

@@ -39,13 +39,18 @@ module Astronoby
       @target_body = target_body
     end
 
+    # @return [Numeric] Julian Date epoch
+    def epoch
+      JulianDate::DEFAULT_EPOCH
+    end
+
     # @return [Astronoby::Coordinates::Equatorial] equatorial coordinates
-    #   derived from the position vector
+    #   derived from the position vector, referred to the frame's epoch
     def equatorial
       @equatorial ||= begin
-        return Coordinates::Equatorial.zero if distance.zero?
+        return Coordinates::Equatorial.zero(epoch: epoch) if distance.zero?
 
-        Coordinates::Equatorial.from_position_vector(@position)
+        Coordinates::Equatorial.from_position_vector(@position, epoch: epoch)
       end
     end
 

data/lib/astronoby/reference_frames/apparent.rb

@@ -6,29 +6,35 @@ module Astronoby
   class Apparent < ReferenceFrame
     # Builds an apparent frame from an astrometric frame by applying
     # aberration in GCRS, then rotating by nutation * precession (with ICRS
-    # frame bias).
+    # frame bias). Without Earth's geometric frame there is no observer
+    # velocity to derive the aberration from, so only the rotation is applied.
     #
     # @param instant [Astronoby::Instant] the time instant
     # @param target_astrometric [Astronoby::Astrometric] target's astrometric
     #   frame
-    # @param earth_geometric [Astronoby::Geometric] Earth's geometric frame
+    # @param earth_geometric [Astronoby::Geometric, nil] Earth's geometric
+    #   frame
     # @param target_body [Astronoby::Body, nil] the target body
     # @return [Astronoby::Apparent] a new apparent frame
     def self.build_from_astrometric(
       instant:,
       target_astrometric:,
-      earth_geometric:,
-      target_body:
+      target_body:,
+      earth_geometric: nil
     )
       position = target_astrometric.position
       velocity = target_astrometric.velocity
       precession_matrix = Precession.matrix_for(instant)
       nutation_matrix = Nutation.matrix_for(instant)
 
-      corrected_position = Aberration.new(
-        astrometric_position: position,
-        observer_velocity: earth_geometric.velocity
-      ).corrected_position
+      corrected_position = if earth_geometric
+        Aberration.new(
+          astrometric_position: position,
+          observer_velocity: earth_geometric.velocity
+        ).corrected_position
+      else
+        position
+      end
 
       # In theory, here we should also apply light deflection. However, so far
       # the deflection algorithm hasn't shown any significant changes to the
@@ -53,6 +59,12 @@ module Astronoby
       )
     end
 
+    # @return [Numeric] Julian Date epoch of the frame's coordinates, which
+    #   are referred to the true equator and equinox of date
+    def epoch
+      instant.tt
+    end
+
     # @return [Astronoby::Coordinates::Ecliptic] ecliptic coordinates at the
     #   current instant (true ecliptic and equinox of date)
     def ecliptic

data/lib/astronoby/reference_frames/mean_of_date.rb

@@ -37,6 +37,12 @@ module Astronoby
       )
     end
 
+    # @return [Numeric] Julian Date epoch of the frame's coordinates, which
+    #   are referred to the mean equator and equinox of date
+    def epoch
+      instant.tt
+    end
+
     # @return [Astronoby::Coordinates::Ecliptic] ecliptic coordinates at the
     #   current instant (mean equinox of date)
     def ecliptic

data/lib/astronoby/reference_frames/teme.rb

@@ -43,6 +43,12 @@ module Astronoby
       )
     end
 
+    # @return [Numeric] Julian Date epoch of the frame's coordinates, which
+    #   are referred to the true equator and mean equinox of date
+    def epoch
+      instant.tt
+    end
+
     # Converts TEME position and velocity to ECEF using the canonical
     # Vallado formulation with R₃(GMST).
     #

data/lib/astronoby/reference_frames/topocentric.rb

@@ -64,6 +64,12 @@ module Astronoby
       @observer = observer
     end
 
+    # @return [Numeric] Julian Date epoch of the frame's coordinates, which
+    #   are referred to the true equator and equinox of date
+    def epoch
+      instant.tt
+    end
+
     # @return [Astronoby::Coordinates::Ecliptic] ecliptic coordinates at the
     #   current instant (true ecliptic and equinox of date)
     def ecliptic

data/lib/astronoby/root_finder.rb

@@ -3,7 +3,7 @@
 module Astronoby
   class RootFinder
     MIN_SAMPLES_PER_PERIOD = 20
-    BISECTION_TOLERANCE_DAYS = 1e-5
+    BISECTION_TOLERANCE_DAYS = 1e-8
 
     # @param value_at [#call] callable mapping a Julian Date (Terrestrial Time)
     #   to a Float
@@ -22,21 +22,28 @@ module Astronoby
     # @return [Array<Float>] root times in Julian Date (Terrestrial Time),
     #   sorted by time
     def roots(start_jd, end_jd, accept: nil)
-      brackets = find_brackets(start_jd, end_jd)
-      located = brackets.map { |a, b| bisect(a, b) }
+      samples = collect_samples(start_jd, end_jd)
+      located = (
+        samples_on_root(samples) +
+          find_brackets(samples).map { |a, b, f_a| bisect(a, b, f_a) }
+      ).uniq.sort
       located.select { |jd| accept.nil? || accept.call(jd) }
     end
 
     private
 
-    def find_brackets(start_jd, end_jd)
-      samples = collect_samples(start_jd, end_jd)
+    def samples_on_root(samples)
+      samples.filter_map { |sample| sample[:jd] if sample[:value].zero? }
+    end
 
+    def find_brackets(samples)
       (0...samples.length - 1).filter_map do |i|
         current = samples[i]
         following = samples[i + 1]
 
-        [current[:jd], following[:jd]] if sign_change?(current, following)
+        if sign_change?(current, following)
+          [current[:jd], following[:jd], current[:value]]
+        end
       end
     end
 
@@ -61,11 +68,11 @@ module Astronoby
       [base_samples, MIN_SAMPLES_PER_PERIOD].max
     end
 
-    def bisect(a, b)
-      f_a = @value_at.call(a)
-
+    def bisect(a, b, f_a)
       while (b - a).abs > BISECTION_TOLERANCE_DAYS
         midpoint = (a + b) / 2
+        break if midpoint == a || midpoint == b
+
         f_midpoint = @value_at.call(midpoint)
         return midpoint if f_midpoint.zero?
 

data/lib/astronoby/time/greenwich_apparent_sidereal_time.rb

@@ -10,16 +10,18 @@ module Astronoby
     # @param utc [Time] the UTC time
     # @return [Astronoby::GreenwichApparentSiderealTime]
     def self.from_utc(utc)
-      gmst = GreenwichMeanSiderealTime.from_utc(utc)
-      instant = Instant.from_time(utc)
-      nutation = Nutation.new(instant: instant)
-      mean_obliquity = MeanObliquity.at(instant)
+      Astronoby.cache.fetch(CacheKey.generate_for_utc(:gast, utc)) do
+        gmst = GreenwichMeanSiderealTime.from_utc(utc)
+        instant = Instant.from_time(utc)
+        nutation = Nutation.new(instant: instant)
+        mean_obliquity = MeanObliquity.at(instant)
 
-      equation_of_equinoxes = nutation.nutation_in_longitude.hours *
-        mean_obliquity.cos
-      gast_time = normalize_time(gmst.time + equation_of_equinoxes)
+        equation_of_equinoxes = nutation.nutation_in_longitude.hours *
+          mean_obliquity.cos
+        gast_time = normalize_time(gmst.time + equation_of_equinoxes)
 
-      new(date: gmst.date, time: gast_time)
+        new(date: gmst.date, time: gast_time)
+      end
     end
 
     # @param date [Date] the calendar date

data/lib/astronoby/time/greenwich_mean_sidereal_time.rb

@@ -27,16 +27,18 @@ module Astronoby
     # @param utc [Time] the UTC time
     # @return [Astronoby::GreenwichMeanSiderealTime]
     def self.from_utc(utc)
-      date = utc.to_date
-      gmst_hours = begin
-        gmst_radians = IERS::GMST.at(utc)
-        gmst_radians * 12.0 / Math::PI
-      rescue IERS::OutOfRangeError
-        from_utc_polynomial(utc)
-      end
-      gmst_hours = normalize_time(gmst_hours)
+      Astronoby.cache.fetch(CacheKey.generate_for_utc(:gmst, utc)) do
+        date = utc.to_date
+        gmst_hours = begin
+          gmst_radians = IERS::GMST.at(utc)
+          gmst_radians * 12.0 / Math::PI
+        rescue IERS::OutOfRangeError
+          from_utc_polynomial(utc)
+        end
+        gmst_hours = normalize_time(gmst_hours)
 
-      new(date: date, time: gmst_hours)
+        new(date: date, time: gmst_hours)
+      end
     end
 
     # Polynomial fallback for dates outside the IERS EOP data range.

data/lib/astronoby/util/time.rb

@@ -1,7 +1,5 @@
 # frozen_string_literal: true
 
-require "iers"
-
 module Astronoby
   module Util
     module Time
@@ -40,26 +38,6 @@ module Astronoby
 
         ::Time.utc(date.year, date.month, date.day, hour, minute, second).round
       end
-
-      # @param instant [Numeric, Time, Date, DateTime]
-      # @return [Numeric] Delta T (TT - UT1) in seconds for the given instant
-      def self.terrestrial_universal_time_delta(instant)
-        case instant
-        when Numeric
-          IERS::DeltaT.at(jd: instant).delta_t
-        when ::Time, ::Date, ::DateTime
-          IERS::DeltaT.at(instant).delta_t
-        else
-          raise IncompatibleArgumentsError,
-            "Expected a Numeric, Time, Date or DateTime object, got #{instant.class}"
-        end
-      rescue IERS::OutOfRangeError => e
-        if e.available_range
-          IERS::DeltaT.at(mjd: e.available_range.end).delta_t
-        else
-          0
-        end
-      end
     end
   end
 end

data/lib/astronoby/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Astronoby
-  VERSION = "0.10.0"
+  VERSION = "0.11.0"
 end

data/lib/astronoby.rb

@@ -18,6 +18,7 @@ require "astronoby/body"
 require "astronoby/position"
 require "astronoby/ephem"
 require "astronoby/julian_date"
+require "astronoby/delta_t"
 require "astronoby/instant"
 require "astronoby/vector"
 require "astronoby/velocity"
@@ -47,6 +48,8 @@ require "astronoby/equinox_solstice"
 require "astronoby/errors"
 require "astronoby/extremum_finder"
 require "astronoby/root_finder"
+require "astronoby/horizon"
+require "astronoby/horizontal_transformation"
 require "astronoby/events/conjunction"
 require "astronoby/events/opposition"
 require "astronoby/events/greatest_elongation"
@@ -55,15 +58,20 @@ require "astronoby/events/greatest_elongation_calculator"
 require "astronoby/events/extremum_event"
 require "astronoby/events/extremum_calculator"
 require "astronoby/events/moon_phases"
+require "astronoby/events/altitude_crossing"
+require "astronoby/events/altitude_extreme"
+require "astronoby/events/altitude_crossing_calculator"
 require "astronoby/events/rise_transit_set_event"
 require "astronoby/events/rise_transit_set_events"
 require "astronoby/events/rise_transit_set_calculator"
 require "astronoby/events/twilight_event"
 require "astronoby/events/twilight_events"
 require "astronoby/events/twilight_calculator"
-require "astronoby/events/eclipse_phase"
+require "astronoby/events/lunar_eclipse_phase"
+require "astronoby/events/lunar_eclipse_geometry"
 require "astronoby/events/lunar_eclipse"
 require "astronoby/events/lunar_eclipse_calculator"
+require "astronoby/events/lunar_eclipse_visibility"
 require "astronoby/geocentric_parallax"
 require "astronoby/libration"
 require "astronoby/mean_obliquity"

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.10.0
+  version: 0.11.0
 platform: ruby
 authors:
 - Rémy Hannequin
@@ -24,7 +24,7 @@ dependencies:
       - !ruby/object:Gem::Version
         version: '0.5'
 - !ruby/object:Gem::Dependency
-  name: iers
+  name: horologium
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
@@ -37,6 +37,26 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: '0.1'
+- !ruby/object:Gem::Dependency
+  name: iers
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.3.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.3.1
 - !ruby/object:Gem::Dependency
   name: matrix
   requirement: !ruby/object:Gem::Requirement
@@ -127,14 +147,20 @@ dependencies:
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '0.22'
+        version: '1.0'
+    - - "<"
+      - !ruby/object:Gem::Version
+        version: '1.2'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '0.22'
+        version: '1.0'
+    - - "<"
+      - !ruby/object:Gem::Version
+        version: '1.2'
 - !ruby/object:Gem::Dependency
   name: standard
   requirement: !ruby/object:Gem::Requirement
@@ -167,6 +193,7 @@ files:
 - Rakefile
 - UPGRADING.md
 - docs/README.md
+- docs/altitude_crossings.md
 - docs/angles.md
 - docs/configuration.md
 - docs/coordinates.md
@@ -225,21 +252,27 @@ files:
 - lib/astronoby/data/constellations/sorted_declinations.dat
 - lib/astronoby/data/constellations/sorted_right_ascensions.dat
 - lib/astronoby/deflection.rb
+- lib/astronoby/delta_t.rb
 - lib/astronoby/distance.rb
 - lib/astronoby/duration.rb
 - lib/astronoby/earth_rotation.rb
 - lib/astronoby/ephem.rb
 - lib/astronoby/equinox_solstice.rb
 - lib/astronoby/errors.rb
+- lib/astronoby/events/altitude_crossing.rb
+- lib/astronoby/events/altitude_crossing_calculator.rb
+- lib/astronoby/events/altitude_extreme.rb
 - lib/astronoby/events/conjunction.rb
 - lib/astronoby/events/conjunction_opposition_calculator.rb
-- lib/astronoby/events/eclipse_phase.rb
 - lib/astronoby/events/extremum_calculator.rb
 - lib/astronoby/events/extremum_event.rb
 - lib/astronoby/events/greatest_elongation.rb
 - lib/astronoby/events/greatest_elongation_calculator.rb
 - lib/astronoby/events/lunar_eclipse.rb
 - lib/astronoby/events/lunar_eclipse_calculator.rb
+- lib/astronoby/events/lunar_eclipse_geometry.rb
+- lib/astronoby/events/lunar_eclipse_phase.rb
+- lib/astronoby/events/lunar_eclipse_visibility.rb
 - lib/astronoby/events/moon_phases.rb
 - lib/astronoby/events/opposition.rb
 - lib/astronoby/events/rise_transit_set_calculator.rb
@@ -250,6 +283,8 @@ files:
 - lib/astronoby/events/twilight_events.rb
 - lib/astronoby/extremum_finder.rb
 - lib/astronoby/geocentric_parallax.rb
+- lib/astronoby/horizon.rb
+- lib/astronoby/horizontal_transformation.rb
 - lib/astronoby/instant.rb
 - lib/astronoby/julian_date.rb
 - lib/astronoby/libration.rb