David Luff: Here is an update to the engine model. It now takes

[flightgear.git] / tests / test-up.cxx

diff --git a/tests/test-up.cxx b/tests/test-up.cxx
index 3623847841797e9fdfb4f3af1aa1e97989c5efab..96fb2c3411a237d8f5ed784d0bfa42bfe5bc5f8c 100644 (file)
--- a/tests/test-up.cxx
+++ b/tests/test-up.cxx
@@ -1,6 +1,8 @@
 // do some test relating to the concept of "up"
 
-#include <iostream>
+#include <simgear/compiler.h>
+
+#include STL_IOSTREAM
 
 #include <plib/sg.h>
 
@@ -28,18 +30,25 @@ int main() {
        Point3D pc = sgGeodToCart( pgd );
        cout << "  cartesian = " << pc << endl;
 
-       sgMat4 UP;
-       sgVec3 geod_up;
-       sgMakeRotMat4( UP, lon, 0.0, -lat );
-       sgSetVec3( geod_up, UP[0][0], UP[0][1], UP[0][2] );
+       sgdMat4 GEOD_UP;
+       sgdVec3 geod_up;
+       sgdMakeRotMat4( GEOD_UP, lon, 0.0, -lat );
+       sgdSetVec3( geod_up, GEOD_UP[0][0], GEOD_UP[0][1], GEOD_UP[0][2] );
        cout << "  geod up = " << geod_up[0] << ", " << geod_up[1] << ", "
             << geod_up[2] << endl;
 
+       sgdMat4 GEOC_UP;
+       sgdVec3 geoc_up;
+       sgdMakeRotMat4( GEOC_UP, lon, 0.0, -lat_geoc * RAD_TO_DEG );
+       sgdSetVec3( geoc_up, GEOC_UP[0][0], GEOC_UP[0][1], GEOC_UP[0][2] );
+       cout << "  geoc up = " << geoc_up[0] << ", " << geoc_up[1] << ", "
+            << geoc_up[2] << endl;
+
        double slope = geod_up[2] / geod_up[0];
        double intercept = pc.z() - slope * pc.x();
        cout << "  Z intercept (based on geodetic up) = " << intercept << endl;
 
-       slope = pc.z() / pc.x();
+       slope = geoc_up[2] / geoc_up[0];
        intercept = pc.z() - slope * pc.x();
        cout << "  Z intercept (based on geocentric up) = " << intercept << endl;