Use L1.2 model for Jupiter Moons

Instead of L1. This seems to improve the precision, fixing mutual eclipse of the moons bug (Stellarium#222). I haven't tested much, and the code might be slower than the previous L1 version that was using interpolation.
StarryLeo · Jul 26, 2018 · 748b13f · 748b13f
1 parent ba436fd
commit 748b13f
Show file tree

Hide file tree

Showing 4 changed files with 784 additions and 4 deletions.
diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
@@ -138,6 +138,8 @@ SET(stellarium_lib_SRCS
      core/planetsephems/gust86.h
      core/planetsephems/l1.c
      core/planetsephems/l1.h
+     core/planetsephems/l12.c
+     core/planetsephems/l12.h
      core/planetsephems/marssat.c
      core/planetsephems/marssat.h
      core/planetsephems/pluto.c

diff --git a/src/core/planetsephems/EphemWrapper.cpp b/src/core/planetsephems/EphemWrapper.cpp
@@ -25,6 +25,7 @@ Foundation, Inc., 51 Franklin Street, Suite 500, Boston, MA  02110-1335, USA.
 #include "elp82b.h"
 #include "marssat.h"
 #include "l1.h"
+#include "l12.h"
 #include "tass17.h"
 #include "gust86.h"
 #include "de431.hpp"
@@ -355,25 +356,25 @@ void get_deimos_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void*
 void get_io_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void* unused)
 {
 	Q_UNUSED(unused);
-	GetL1Coor(jd, L1_IO, xyz, xyzdot);
+	GetL12Coor(jd, L1_IO, xyz, xyzdot);
 }
 
 void get_europa_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void* unused)
 {
 	Q_UNUSED(unused);
-	GetL1Coor(jd, L1_EUROPA, xyz, xyzdot);
+	GetL12Coor(jd, L1_EUROPA, xyz, xyzdot);
 }
 
 void get_ganymede_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void* unused)
 {
 	Q_UNUSED(unused);
-	GetL1Coor(jd, L1_GANYMEDE, xyz, xyzdot);
+	GetL12Coor(jd, L1_GANYMEDE, xyz, xyzdot);
 }
 
 void get_callisto_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void* unused)
 {
 	Q_UNUSED(unused);
-	GetL1Coor(jd, L1_CALLISTO, xyz, xyzdot);
+	GetL12Coor(jd, L1_CALLISTO, xyz, xyzdot);
 }
 
 void get_mimas_parent_coordsv(double jd, double xyz[3], double xyzdot[3], void* unused)