diff --git a/src/orbits.cpp b/src/orbits.cpp index d393e52..489125f 100644 --- a/src/orbits.cpp +++ b/src/orbits.cpp @@ -40,6 +40,8 @@ ellipseInitAB(double a, double b) return ep; } +// FIXME: we should avoid calling ellipseInitAB, and recalculate the properties +// to avoid floating point errors in the known quantity 'e' EllipseParameters ellipseInitAE(double a, double e) {