Minor code updates for testing

tudat-team · Jan 12, 2024 · 1be70b5 · 1be70b5
1 parent 515423e
commit 1be70b5
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/astro/gravitation/polyhedronGravityField.cpp b/src/astro/gravitation/polyhedronGravityField.cpp
@@ -19,7 +19,7 @@ namespace gravitation
 
 void PolyhedronGravityCache::update (const Eigen::Vector3d& currentBodyFixedPosition)
 {
-    if ( currentBodyFixedPosition != currentBodyFixedPosition_ )
+    if ( currentBodyFixedPosition_.hasNaN( ) || currentBodyFixedPosition != currentBodyFixedPosition_ )
     {
         currentBodyFixedPosition_ = currentBodyFixedPosition;
 

diff --git a/tests/src/astro/propagators/unitTestDependentVariableOutput.cpp b/tests/src/astro/propagators/unitTestDependentVariableOutput.cpp
@@ -1670,7 +1670,7 @@ BOOST_AUTO_TEST_CASE( test_GravitationalPotentialAndLaplacianSaving )
                         moonGravitationalPotential,
                         moonGravityModel->getGravitationalPotential( moonBodyFixedCartesianPosition ),
                         1e-15 );*/
-                std::cout<<"Moon potential "<<moonGravitationalPotential<<" "<<moonGravityModel->getGravitationalPotential( moonBodyFixedCartesianPosition )<<std::endl;
+                std::cout<<"Moon potential "<<moonGravitationalPotential<<" "<<moonGravityModel->getGravitationalPotential( moonBodyFixedCartesianPosition ) - moonGravitationalPotential<<std::endl;
                 BOOST_CHECK( std::fabs( moonGravitationalPotential - moonGravityModel->getGravitationalPotential( moonBodyFixedCartesianPosition ) )
                             < std::fabs( std::min( moonGravitationalPotential, moonGravityModel->getGravitationalPotential( moonBodyFixedCartesianPosition ) ) * 1e-15 ) );