esa / pykep

PyKEP is a scientific library providing basic tools for research in interplanetary trajectory design.
http://esa.github.io/pykep/
GNU General Public License v3.0
334 stars 96 forks source link

keplerian planets are singular in retrogate 0 inclination cases #162

Open darioizzo opened 1 year ago

darioizzo commented 1 year ago

position=[10000000, 0,0] 
velocity=[0, -8000, 0]
epoch=pk.epoch(0)
prova=pk.planet.keplerian(pk.epoch(0),position,velocity,pk.MU_EARTH,1.0,1.0,1.0,"prova")
prova.eph(pk.epoch(0))```

this returns nans! Note inclination here is 3.14, not zero as checked!