search

exoplanet-dev / exoplanet

Fast & scalable MCMC for all your exoplanet needs!
https://docs.exoplanet.codes
MIT License
206 stars 52 forks source link

Issue with contact points for some orbits #32

Closed dfm closed 5 years ago

dfm commented 5 years ago

Sometimes the start point is larger than the end point. For now, I'm just catching that and bailing but we should work out where this is coming from.

dfm commented 5 years ago

Here's an example where this failure happens:

orbit = xo.orbits.KeplerianOrbit(period=3.456, ecc=0.6, omega=-1.5)
t = np.linspace(-0.1, 0.1, 1000)
u = [0.3, 0.2]
light_curve = xo.StarryLightCurve(u).get_light_curve(
    orbit=orbit, r=0.1, t=t, texp=0.02).eval()