Tweak accuracy and timestep to improve testing

Author: joaogvcarneiro

src/simulation/dynamics/linearTranslationalBodies/linearTranslationBodiesNDOF/_UnitTest/test_linearTranslationNDOFStateEffector.py

@@ -308,7 +308,7 @@ def translatingBodyNoInput(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy, err_msg="Orbital energy is not constant.")
     np.testing.assert_allclose(finalRotEnergy, initialRotEnergy, rtol=accuracy,
@@ -556,7 +556,7 @@ def translatingBodyLockAxis(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy,
                                err_msg="Orbital energy is not constant.")
@@ -797,7 +797,7 @@ def translatingBodyCommandedForce(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy,
                                err_msg="Orbital energy is not constant.")

src/simulation/dynamics/linearTranslationalBodies/linearTranslationBodiesOneDOF/_UnitTest/test_linearTranslationOneDOFStateEffector.py

@@ -99,7 +99,7 @@ def translatingBodyNoInput(show_plots):
     unitTestSim = SimulationBaseClass.SimBaseClass()
 
     # Create test thread
-    testProcessRate = macros.sec2nano(0.001)  # update process rate update time
+    testProcessRate = macros.sec2nano(0.0001)  # update process rate update time
     testProc = unitTestSim.CreateNewProcess(unitProcessName)
     testProc.addTask(unitTestSim.CreateNewTask(unitTaskName, testProcessRate))
 
@@ -255,7 +255,7 @@ def translatingBodyNoInput(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy, err_msg="Orbital energy is not constant.")
     np.testing.assert_allclose(finalRotEnergy, initialRotEnergy, rtol=accuracy, err_msg="Rotational energy is not constant.")
@@ -278,7 +278,7 @@ def translatingBodyLockFlag(show_plots):
     unitTestSim = SimulationBaseClass.SimBaseClass()
 
     # Create test thread
-    testProcessRate = macros.sec2nano(0.001)  # update process rate update time
+    testProcessRate = macros.sec2nano(0.0001)  # update process rate update time
     testProc = unitTestSim.CreateNewProcess(unitProcessName)
     testProc.addTask(unitTestSim.CreateNewTask(unitTaskName, testProcessRate))
 
@@ -440,7 +440,7 @@ def translatingBodyLockFlag(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy, err_msg="Orbital energy is not constant.")
     np.testing.assert_allclose(finalRotEnergy, initialRotEnergy, rtol=accuracy, err_msg="Rotational energy is not constant.")
@@ -463,7 +463,7 @@ def translatingBodyCommandedForce(show_plots, cmdForce):
     unitTestSim = SimulationBaseClass.SimBaseClass()
 
     # Create test thread
-    testProcessRate = macros.sec2nano(0.001)  # update process rate update time
+    testProcessRate = macros.sec2nano(0.0001)  # update process rate update time
     testProc = unitTestSim.CreateNewProcess(unitProcessName)
     testProc.addTask(unitTestSim.CreateNewTask(unitTaskName, testProcessRate))
 
@@ -615,7 +615,7 @@ def translatingBodyCommandedForce(show_plots, cmdForce):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy, err_msg="Orbital energy is not constant.")
     for i in range(3):
@@ -637,7 +637,7 @@ def translatingBodyRhoReference(show_plots, rhoRef):
     unitTestSim = SimulationBaseClass.SimBaseClass()
 
     # Create test thread
-    testProcessRate = macros.sec2nano(0.001)  # update process rate update time
+    testProcessRate = macros.sec2nano(0.0001)  # update process rate update time
     testProc = unitTestSim.CreateNewProcess(unitProcessName)
     testProc.addTask(unitTestSim.CreateNewTask(unitTaskName, testProcessRate))
 
@@ -726,7 +726,7 @@ def translatingBodyRhoReference(show_plots, rhoRef):
     unitTestSim.AddModelToTask(unitTaskName, rhoData)
 
     # Setup and run the simulation
-    stopTime = 25000 * testProcessRate
+    stopTime = 30000 * testProcessRate  # needs a little more time to settle to the correct angle
     unitTestSim.ConfigureStopTime(stopTime)
     unitTestSim.ExecuteSimulation()
 
@@ -791,7 +791,7 @@ def translatingBodyRhoReference(show_plots, rhoRef):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalAngle, rhoRef, atol=0.01, err_msg="Angle doesn't settle to reference angle.")
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy, err_msg="Orbital energy is not constant.")

src/simulation/dynamics/spinningBodies/spinningBodiesNDOF/_UnitTest/test_spinningBodyNDOFStateEffector.py

@@ -319,7 +319,7 @@ def spinningBodyNoInput(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy)
     np.testing.assert_allclose(finalRotEnergy, initialRotEnergy, rtol=accuracy)
@@ -578,7 +578,7 @@ def spinningBodyLockAxis(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy)
     np.testing.assert_allclose(finalRotEnergy, initialRotEnergy, rtol=accuracy)
@@ -836,7 +836,7 @@ def spinningBodyCommandedTorque(show_plots):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
 
     np.testing.assert_allclose(finalOrbEnergy, initialOrbEnergy, rtol=accuracy)
     for i in range(3):

src/simulation/dynamics/spinningBodies/spinningBodiesOneDOF/_UnitTest/test_spinningBodyOneDOFStateEffector.py

@@ -166,7 +166,7 @@ def spinningBody(show_plots, cmdTorque, lock, thetaRef):
     # Add energy and momentum variables to log
     scObjectLog = scObject.logger(["totOrbAngMomPntN_N", "totRotAngMomPntC_N", "totOrbEnergy", "totRotEnergy"])
     unitTestSim.AddModelToTask(unitTaskName, scObjectLog)
-    
+
     # Initialize the simulation
     unitTestSim.InitializeSimulation()
 
@@ -248,7 +248,7 @@ def spinningBody(show_plots, cmdTorque, lock, thetaRef):
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 1e-13
     finalOrbAngMom = numpy.delete(finalOrbAngMom, 0, axis=1)  # remove time column
     finalRotAngMom = numpy.delete(finalRotAngMom, 0, axis=1)  # remove time column
     finalRotEnergy = numpy.delete(finalRotEnergy, 0, axis=1)  # remove time column

src/simulation/dynamics/spinningBodies/spinningBodiesTwoDOF/_UnitTest/test_spinningBodyTwoDOFStateEffector.py

@@ -192,7 +192,7 @@ def spinningBody(show_plots, cmdTorque1, lock1, theta1Ref, cmdTorque2, lock2, th
     # Add energy and momentum variables to log
     scObjectLog = scObject.logger(["totOrbAngMomPntN_N", "totRotAngMomPntC_N", "totOrbEnergy", "totRotEnergy"])
     unitTestSim.AddModelToTask(unitTaskName, scObjectLog)
-    
+
     # Initialize the simulation
     unitTestSim.InitializeSimulation()
 
@@ -290,7 +290,7 @@ def spinningBody(show_plots, cmdTorque1, lock1, theta1Ref, cmdTorque2, lock2, th
     plt.close("all")
 
     # Testing setup
-    accuracy = 1e-12
+    accuracy = 3e-13
     finalOrbAngMom = numpy.delete(finalOrbAngMom, 0, axis=1)  # remove time column
     finalRotAngMom = numpy.delete(finalRotAngMom, 0, axis=1)  # remove time column
     finalRotEnergy = numpy.delete(finalRotEnergy, 0, axis=1)  # remove time column