コード例 #1
0
ファイル: hindu.py プロジェクト: rob-blackbourn/PyFinance
def hindu_solar_longitude_at_or_after(lam, tee):
    """Return the moment of the first time at or after moment, tee
    when Hindu solar longitude will be lam degrees."""
    tau = tee + (HinduSolarDate.SIDEREAL_YEAR * (1 / 360) * mod(lam - HinduDate.solar_longitude(tee), 360))
    a = max(tee, tau - 5)
    b = tau +5
    return invert_angular(HinduDate.solar_longitude, lam, a, b)
コード例 #2
0
ファイル: hindu.py プロジェクト: rob-blackbourn/PyFinance
 def lunar_day_at_or_after(cls, k, tee):
     """Return the time lunar_day (tithi) number, k, begins at or after
     moment, tee.  k can be fractional (for karanas)."""
     phase = (k - 1) * 12
     tau   = tee + ((1/360) * mod(phase - cls.lunar_phase(tee), 360) * cls.SYNODIC_MONTH)
     a = max(tee, tau - 2)
     b = tau + 2
     return invert_angular(cls.lunar_phase, phase, a, b)
コード例 #3
0
ファイル: solar.py プロジェクト: rob-blackbourn/PyFinance
 def solar_longitude_after(cls, lam, tee):
     """Return the moment UT of the first time at or after moment, tee,
     when the solar longitude will be lam degrees."""
     rate = cls.MEAN_TROPICAL_YEAR / 360
     tau = tee + rate * mod(lam - cls.solar_longitude(tee), 360)
     a = max(tee, tau - 5)
     b = tau + 5
     return invert_angular(cls.solar_longitude, lam, a, b)
コード例 #4
0
ファイル: lunar.py プロジェクト: rob-blackbourn/PyFinance
 def lunar_phase_at_or_after(cls, phi, tee):
     """Return the moment UT of the next time at or after moment, tee,
     when the lunar_phase is phi degrees."""
     tau = (tee +
            (cls.MEAN_SYNODIC_MONTH    *
             (1/360) *
             mod(phi - cls.lunar_phase(tee), 360)))
     a = max(tee, tau - 2)
     b = tau + 2
     return invert_angular(cls.lunar_phase, phi, a, b)
コード例 #5
0
ファイル: lunar.py プロジェクト: rob-blackbourn/PyFinance
 def lunar_phase_at_or_before(cls, phi, tee):
     """Return the moment UT of the last time at or before moment, tee,
     when the lunar_phase was phi degrees."""
     tau = (tee -
            (cls.MEAN_SYNODIC_MONTH  *
             (1/360) *
             mod(cls.lunar_phase(tee) - phi, 360)))
     a = tau - 2
     b = min(tee, tau +2)
     return invert_angular(cls.lunar_phase, phi, a, b)