Add test coverage.

Author: Chris Birkhold

tests/unit_tests/FGAtmosphereTest.h

@@ -150,8 +150,24 @@ class FGAtmosphereTest : public CxxTest::TestSuite
     const double nu0 = mu0/rho0;
 
     for(double h=-1000.0; h<10000; h+= 1000) {
-      double T = T0 + 0.1*h;
-      double P = P0 + 1.0*h;
+      const double T = T0 + 0.1*h;
+      const double P = P0 + 1.0*h;
+      const double rho = P/(R*T);
+      const double a = sqrt(gama*R*T);
+
+      FGAtmosphere::AltitudeDependentValues adv;
+      atm.GetAltitudeDependentValues(adv, h);
+      TS_ASSERT_DELTA(adv.Temperature, T, epsilon);
+      TS_ASSERT_DELTA(adv.Pressure, P, epsilon);
+      TS_ASSERT_DELTA(adv.Density, rho, epsilon);
+      TS_ASSERT_DELTA(adv.SoundSpeed, a, epsilon);
+
+      FGAtmosphere::AltitudeDependentValues adv_0;
+      atm.GetAltitudeDependentValues(adv, 0.0);
+      TS_ASSERT_EQUALS(adv_0.Temperature, T0);
+      TS_ASSERT_EQUALS(adv_0.Pressure, P0);
+      TS_ASSERT_DELTA(adv_0.Density, rho0, epsilon);
+      TS_ASSERT_DELTA(adv_0.SoundSpeed, a0, epsilon);
 
       TS_ASSERT_DELTA(atm.GetTemperature(h), T, epsilon);
       TS_ASSERT_EQUALS(atm.GetTemperature(0.0), T0);
@@ -160,11 +176,9 @@ class FGAtmosphereTest : public CxxTest::TestSuite
       TS_ASSERT_DELTA(atm.GetPressure(h), P, epsilon);
       TS_ASSERT_EQUALS(atm.GetPressure(0.0), P0);
 
-      double rho = P/(R*T);
       TS_ASSERT_DELTA(atm.GetDensity(h), rho, epsilon);
       TS_ASSERT_DELTA(atm.GetDensity(0.0), rho0, epsilon);
 
-      double a = sqrt(gama*R*T);
       TS_ASSERT_DELTA(atm.GetSoundSpeed(h), a, epsilon);
       TS_ASSERT_DELTA(atm.GetSoundSpeed(0.0), a0, epsilon);
 

tests/unit_tests/FGColumnVector3Test.h

@@ -577,6 +577,24 @@ class FGColumnVector3Test : public CxxTest::TestSuite
     TS_ASSERT_EQUALS(Z(3), 1.0);
   }
 
+  void testSqrMagnitude(void) {
+    JSBSim::FGColumnVector3 v0;
+    JSBSim::FGColumnVector3 v(3.0, 4.0, 0.0);
+
+    TS_ASSERT_EQUALS(v0.SqrMagnitude(), 0.0);
+    TS_ASSERT_EQUALS(v.SqrMagnitude(), 25.0);
+    TS_ASSERT_EQUALS(DotProduct(v,v), v.SqrMagnitude());
+
+    // Verify that the operands are not modified
+    TS_ASSERT_EQUALS(v0(1), 0.0);
+    TS_ASSERT_EQUALS(v0(2), 0.0);
+    TS_ASSERT_EQUALS(v0(3), 0.0);
+
+    TS_ASSERT_EQUALS(v(1), 3.0);
+    TS_ASSERT_EQUALS(v(2), 4.0);
+    TS_ASSERT_EQUALS(v(3), 0.0);
+  }
+
   void testMagnitude(void) {
     JSBSim::FGColumnVector3 v0;
     JSBSim::FGColumnVector3 v(3.0, 4.0, 0.0);

tests/unit_tests/FGGroundCallbackTest.h

@@ -55,6 +55,8 @@ class FGGroundCallbackTest : public CxxTest::TestSuite
         double lat_rad = lat*M_PI/180.;
         loc = FGLocation(lon_rad, lat_rad, RadiusReference);
         double agl = cb->GetAGLevel(loc, contact, normal, v, w);
+        double agl_2 = cb->GetAGLevel(loc);
+        TS_ASSERT_EQUALS(agl, agl_2);
         TS_ASSERT_DELTA(0.0, agl, 1e-8);
         TS_ASSERT_VECTOR_EQUALS(v, zero);
         TS_ASSERT_VECTOR_EQUALS(w, zero);
@@ -88,6 +90,8 @@ class FGGroundCallbackTest : public CxxTest::TestSuite
         double lat_rad = lat*M_PI/180.;
         loc = FGLocation(lon_rad, lat_rad, RadiusReference+h);
         double agl = cb->GetAGLevel(loc, contact, normal, v, w);
+        double agl_2 = cb->GetAGLevel(loc);
+        TS_ASSERT_EQUALS(agl, agl_2);
         TS_ASSERT_DELTA(h/agl, 1.0, epsilon*100.);
         TS_ASSERT_VECTOR_EQUALS(v, zero);
         TS_ASSERT_VECTOR_EQUALS(w, zero);
@@ -129,6 +133,8 @@ class FGGroundCallbackTest : public CxxTest::TestSuite
         double lat_rad = lat*M_PI/180.;
         loc = FGLocation(lon_rad, lat_rad, RadiusReference+h);
         double agl = cb->GetAGLevel(loc, contact, normal, v, w);
+        double agl_2 = cb->GetAGLevel(loc);
+        TS_ASSERT_EQUALS(agl, agl_2);
         TS_ASSERT_DELTA((h-elevation)/agl, 1.0, epsilon*100.);
         TS_ASSERT_VECTOR_EQUALS(v, zero);
         TS_ASSERT_VECTOR_EQUALS(w, zero);
@@ -195,6 +201,8 @@ class FGGroundCallbackTest : public CxxTest::TestSuite
         double lat_rad = lat*M_PI/180.;
         loc.SetPositionGeodetic(lon_rad, lat_rad, h);
         double agl = cb->GetAGLevel(loc, contact, normal, v, w);
+        double agl_2 = cb->GetAGLevel(loc);
+        TS_ASSERT_EQUALS(agl, agl_2);
         TS_ASSERT_DELTA(h, agl, 1e-8);
         TS_ASSERT_VECTOR_EQUALS(v, zero);
         TS_ASSERT_VECTOR_EQUALS(w, zero);

tests/unit_tests/FGLocationTest.h

@@ -115,6 +115,9 @@ class  FGLocationTest : public CxxTest::TestSuite
 
     lv3.SetEllipse(1., 1.);
     CheckLocation(lv3, v);
+
+    JSBSim::FGLocation lv4(v, 1., 1.);
+    CheckLocation(lv3, v);
   }
 
   void testCopyConstructor() {