# Path: skyfield/constants.py
# AU_M = 149597870700             # per IAU 2012 Resolution B2
# AU_KM = 149597870.700
# ASEC360 = 1296000.0
# DAY_S = 86400.0
# ASEC2RAD = 4.848136811095359935899141e-6
# DEG2RAD = 0.017453292519943296
# RAD2DEG = 57.295779513082321
# C = 299792458.0                            # m/s
# ANGVEL = 7.2921150e-5           # radians/s
# ERAD = 6378136.6                # meters
# IERS_2010_INVERSE_EARTH_FLATTENING = 298.25642
# GS = 1.32712440017987e+20
# T0 = 2451545.0
# B1950 = 2433282.4235
# C_AUDAY = C * DAY_S / AU_M
# 
# Path: skyfield/functions.py
# def load_bundled_npy(filename):
#     """Load a binary NumPy array file that is bundled with Skyfield."""
#     data = get_data('skyfield', 'data/{0}'.format(filename))
#     return load(BytesIO(data))


from numpy import array, cos, dot, fmod, sin, outer, zeros
from .constants import ASEC2RAD, ASEC360, DEG2RAD, tau, T0
from .functions import load_bundled_npy
    a = next(fa) * t
    for fa_i in fa:
        a += fa_i
        a *= t
    a += fa0
    fmod(a, ASEC360, out=a)
    a *= ASEC2RAD
    if getattr(t, 'shape', ()):
        return a
    return a[:,0]

# Sine and cosine coefficients for t^1.

se1_0 = -0.87e-6
se1_1 = +0.00e-6

# Deprecated functions that third-party code might still call; several
# of our tests also still call them, to help keep them working.

def compute_nutation(t):
    """Deprecated: this is now a method on the Time object."""
    return t.N

def earth_tilt(t):
    """Deprecated: these are now computed separately on the Time object."""
    d_psi, d_eps = t._nutation_angles_radians
    mean_ob = t._mean_obliquity_radians
    true_ob = mean_ob + d_eps
    c_terms = equation_of_the_equinoxes_complimentary_terms(t.tt)
    eq_eq = d_psi * cos(mean_ob) + c_terms

# Path: skyfield/constants.py
# AU_M = 149597870700             # per IAU 2012 Resolution B2
# AU_KM = 149597870.700
# ASEC360 = 1296000.0
# DAY_S = 86400.0
# ASEC2RAD = 4.848136811095359935899141e-6
# DEG2RAD = 0.017453292519943296
# RAD2DEG = 57.295779513082321
# C = 299792458.0                            # m/s
# ANGVEL = 7.2921150e-5           # radians/s
# ERAD = 6378136.6                # meters
# IERS_2010_INVERSE_EARTH_FLATTENING = 298.25642
# GS = 1.32712440017987e+20
# T0 = 2451545.0
# B1950 = 2433282.4235
# C_AUDAY = C * DAY_S / AU_M
# 
# Path: skyfield/functions.py
# def load_bundled_npy(filename):
#     """Load a binary NumPy array file that is bundled with Skyfield."""
#     data = get_data('skyfield', 'data/{0}'.format(filename))
#     return load(BytesIO(data))


# Path: skyfield/nutationlib.py
from numpy import array, cos, dot, fmod, sin, outer, zeros
from .constants import ASEC2RAD, ASEC360, DEG2RAD, tau, T0
from .functions import load_bundled_npy

    a = next(fa) * t
    for fa_i in fa:
        a += fa_i
        a *= t
    a += fa0
    fmod(a, ASEC360, out=a)
    a *= ASEC2RAD
    if getattr(t, 'shape', ()):
        return a
    return a[:,0]

# Sine and cosine coefficients for t^1.

se1_0 = -0.87e-6
se1_1 = +0.00e-6

# Deprecated functions that third-party code might still call; several
# of our tests also still call them, to help keep them working.

def compute_nutation(t):
    """Deprecated: this is now a method on the Time object."""
    return t.N

def earth_tilt(t):
    """Deprecated: these are now computed separately on the Time object."""
    d_psi, d_eps = t._nutation_angles_radians
    mean_ob = t._mean_obliquity_radians
    true_ob = mean_ob + d_eps
    c_terms = equation_of_the_equinoxes_complimentary_terms(t.tt)
    eq_eq = d_psi * cos(mean_ob) + c_terms

    return (mean_ob / DEG2RAD, true_ob / DEG2RAD, eq_eq / ASEC2RAD / 15.0,