In fgCalculatePlanet() pass a pointer to a structure to be modified, rather

[simgear.git] / Astro / planets.cxx

diff --git a/Astro/planets.cxx b/Astro/planets.cxx
index 7b4351639b439c656b85fab195a946838ea086ad..1fa7ce1b0cf8a9b88f44f8df0a9ac23485a36c27 100644 (file)
--- a/Astro/planets.cxx
+++ b/Astro/planets.cxx
@@ -45,11 +45,10 @@
 
 GLint planets = 0;
 
-struct CelestialCoord fgCalculatePlanet(struct OrbElements planet,
-                                        struct OrbElements theSun,
-                                        struct fgTIME t, int idx)
+void fgCalculatePlanet(struct OrbElements planet, struct OrbElements theSun,
+                      fgTIME t, int idx, struct CelestialCoord *result)
 {
-    struct CelestialCoord result;
+    // struct CelestialCoord result;
 
     fgSUNPOS SolarPosition;
 
@@ -86,8 +85,8 @@ struct CelestialCoord fgCalculatePlanet(struct OrbElements planet,
     ye = yg * cos(ecl) - zg * sin(ecl);
     ze = yg * sin(ecl) + zg * cos(ecl);
 
-    result.RightAscension = atan2(ye,xe);
-    result.Declination = atan2(ze, sqrt(xe*xe + ye*ye));
+    result->RightAscension = atan2(ye,xe);
+    result->Declination = atan2(ze, sqrt(xe*xe + ye*ye));
         
     /* Let's calculate the brightness of the planet */
     R = sqrt ( xg*xg + yg*yg + zg*zg);
@@ -97,53 +96,54 @@ struct CelestialCoord fgCalculatePlanet(struct OrbElements planet,
     switch(idx)
       {
       case 2: /* mercury */
-       result.magnitude = -0.36 + 5*log10( r*R ) + 0.027 * FV + 2.2E-13 * pow(FV, 6);
+       result->magnitude = -0.36 + 5*log10( r*R ) + 0.027 * FV + 2.2E-13 * pow(FV, 6);
        break;
       case 3: /*venus */
-       result.magnitude = -4.34 + 5*log10( r*R ) + 0.013 * FV + 4.2E-07 * pow(FV,3);
+       result->magnitude = -4.34 + 5*log10( r*R ) + 0.013 * FV + 4.2E-07 * pow(FV,3);
        break;
       case 4: /* mars */
-       result.magnitude = -1.51 + 5*log10( r*R ) + 0.016 * FV;
+       result->magnitude = -1.51 + 5*log10( r*R ) + 0.016 * FV;
        break;
       case 5: /* Jupiter */
-       result.magnitude = -9.25 + 5*log10( r*R ) + 0.014 * FV;
+       result->magnitude = -9.25 + 5*log10( r*R ) + 0.014 * FV;
        break;
       case 6: /* Saturn */
        
        ir = 0.4897394;
        Nr = 2.9585076 + 6.6672E-7*actTime;
        
-       B = asin (sin (result.Declination) * cos (ir) -
-                 cos (result.Declination) * sin (ir) *
-                 sin (result.RightAscension - Nr));
+       B = asin (sin (result->Declination) * cos (ir) -
+                 cos (result->Declination) * sin (ir) *
+                 sin (result->RightAscension - Nr));
        ring_magn = -2.6 * sin (fabs(B)) + 1.2 * pow(sin(B),2);
-       result.magnitude = -9.0 + 5*log10( r*R ) + 0.044 * FV + ring_magn;
+       result->magnitude = -9.0 + 5*log10( r*R ) + 0.044 * FV + ring_magn;
        break;
       case 7: /* Uranus */
-       result.magnitude = -7.15 + 5*log10( r*R) + 0.001 * FV;
+       result->magnitude = -7.15 + 5*log10( r*R) + 0.001 * FV;
        break;
       case 8:  /* Neptune */
-       result.magnitude = -6.90 + 5*log10 (r*R) + 0.001 *FV;
+       result->magnitude = -6.90 + 5*log10 (r*R) + 0.001 *FV;
        break;
       default:
        fgPrintf( FG_ASTRO, FG_ALERT, "index %d out of range !!!!\n", idx);
       }
     fgPrintf( FG_ASTRO, FG_DEBUG,
              "    Planet found at %f (ra), %f (dec)\n",
-             result.RightAscension, result.Declination);
+             result->RightAscension, result->Declination);
     fgPrintf( FG_ASTRO, FG_DEBUG,
              "      Geocentric dist     %f\n"
              "      Heliocentric dist   %f\n"
              "      Distance to the sun %f\n"
              "      Phase angle         %f\n"
-             "      Brightness          %f\n", R, r, s, FV, result.magnitude);
-    return result;
+             "      Brightness          %f\n", R, r, s, FV, result->magnitude);
+
+    // return result;
 }
 
 
 void fgPlanetsInit( void )
 {
-  struct fgLIGHT *l;
+  fgLIGHT *l;
   int i;
   struct CelestialCoord pltPos;
   double magnitude;
@@ -164,8 +164,8 @@ void fgPlanetsInit( void )
   /* Add the planets to all four display lists */
   for ( i = 2; i < 9; i++ ) {
     fgSolarSystemUpdate(&(pltOrbElements[i]), cur_time_params);
-    pltPos = fgCalculatePlanet(pltOrbElements[i], 
-                              pltOrbElements[0], cur_time_params, i);
+    fgCalculatePlanet(pltOrbElements[i], pltOrbElements[0], 
+                     cur_time_params, i, &pltPos);
     
     magnitude = (0.0 - pltPos.magnitude) / 5.0 + 1.0;
     
@@ -210,10 +210,23 @@ void fgPlanetsRender( void ) {
 
 
 /* $Log$
-/* Revision 1.2  1998/04/24 00:45:02  curt
-/* Wrapped "#include <config.h>" in "#ifdef HAVE_CONFIG_H"
-/* Fixed a bug when generating sky colors.
+/* Revision 1.6  1998/06/27 16:51:54  curt
+/* In fgCalculatePlanet() pass a pointer to a structure to be modified, rather
+/* than returning the entire structure.
 /*
+ * Revision 1.5  1998/04/28 01:19:01  curt
+ * Type-ified fgTIME and fgVIEW
+ *
+ * Revision 1.4  1998/04/26 05:10:01  curt
+ * "struct fgLIGHT" -> "fgLIGHT" because fgLIGHT is typedef'd.
+ *
+ * Revision 1.3  1998/04/25 22:06:25  curt
+ * Edited cvs log messages in source files ... bad bad bad!
+ *
+ * Revision 1.2  1998/04/24 00:45:02  curt
+ * Wrapped "#include <config.h>" in "#ifdef HAVE_CONFIG_H"
+ * Fixed a bug when generating sky colors.
+ *
  * Revision 1.1  1998/04/22 13:21:31  curt
  * C++ - ifing the code a bit.
  *
@@ -242,7 +255,7 @@ void fgPlanetsRender( void ) {
  * To version 0.29
  *
  * Revision 1.3  1998/01/27 00:47:47  curt
- * Incorporated Paul Bleisch's <bleisch@chromatic.com> new debug message
+ * Incorporated Paul Bleisch's <pbleisch@acm.org> new debug message
  * system and commandline/config file processing code.
  *
  * Revision 1.2  1998/01/19 19:26:59  curt