# Path: pyLIMA/microlorbitalmotion.py
# def orbital_motion_shifts(orbital_motion_model, time, pyLIMA_parameters):
# def orbital_motion_2D_trajectory_shift(to_om, time, dalpha_dt):
# def orbital_motion_2D_separation_shift(to_om, time, ds_dt):
# def orbital_motion_circular(to_om, v_para, v_perp, v_radial, separation_0, time):
# def orbital_motion_keplerian(to_om, v_para, v_perp, v_radial, separation_0, separation_z, mass, rE, time):
# def orbital_motion_keplerian_direct(to_om, a_true,period,eccentricity,inclination,omega_node,omega_periastron,t_periastron,
#                                   rE, time):
# def eccentric_anomaly_function(time, ellipticity, t_periastron, speed):
# def mean_to_eccentric_anomaly(params, mean_anomaly, ellipticity):
# def mean_to_eccentric_anomaly_prime(params, mean_anomaly, ellipticity):
# def mean_to_eccentric_anomaly_prime2(params, mean_anomaly, ellipticity):
#     N = 2 * np.pi / period * 1 / 365.25
#     N = 2*np.pi/period/365.25


import numpy as np
import unittest.mock as mock
import collections
from pyLIMA import microlorbitalmotion
# -*- coding: utf-8 -*-
"""
Created on Mon Dec  7 13:46:13 2015

@author: ebachelet
"""



def test_orbital_motion_shifts():
    model = ['2D', 2458000]

    pyLIMA_parameters = collections.namedtuple('params', ['dsdt', 'dalphadt'])
    pyLIMA_parameters.dsdt = 0.0001
    pyLIMA_parameters.dalphadt = 0.0001

    time = np.arange(2457950., 2458050., 10)


# Path: pyLIMA/microlorbitalmotion.py
# def orbital_motion_shifts(orbital_motion_model, time, pyLIMA_parameters):
# def orbital_motion_2D_trajectory_shift(to_om, time, dalpha_dt):
# def orbital_motion_2D_separation_shift(to_om, time, ds_dt):
# def orbital_motion_circular(to_om, v_para, v_perp, v_radial, separation_0, time):
# def orbital_motion_keplerian(to_om, v_para, v_perp, v_radial, separation_0, separation_z, mass, rE, time):
# def orbital_motion_keplerian_direct(to_om, a_true,period,eccentricity,inclination,omega_node,omega_periastron,t_periastron,
#                                   rE, time):
# def eccentric_anomaly_function(time, ellipticity, t_periastron, speed):
# def mean_to_eccentric_anomaly(params, mean_anomaly, ellipticity):
# def mean_to_eccentric_anomaly_prime(params, mean_anomaly, ellipticity):
# def mean_to_eccentric_anomaly_prime2(params, mean_anomaly, ellipticity):
#     N = 2 * np.pi / period * 1 / 365.25
#     N = 2*np.pi/period/365.25


# Path: pyLIMA/test/test_microlorbitalmotion.py
import numpy as np
import unittest.mock as mock
import collections
from pyLIMA import microlorbitalmotion

# -*- coding: utf-8 -*-
"""
Created on Mon Dec  7 13:46:13 2015

@author: ebachelet
"""



def test_orbital_motion_shifts():
    model = ['2D', 2458000]

    pyLIMA_parameters = collections.namedtuple('params', ['dsdt', 'dalphadt'])
    pyLIMA_parameters.dsdt = 0.0001
    pyLIMA_parameters.dalphadt = 0.0001

    time = np.arange(2457950., 2458050., 10)


    ds, dalpha = microlorbitalmotion.orbital_motion_shifts(model, time, pyLIMA_parameters)