Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/twobody.c

Timestamp:

Jun 15, 2004, 6:54:12 PM (21 years ago)

Author:

cmv

Message:

nouvelle version de xephem/libastro (3.6) cmv 15/6/04

File:

• trunk/SophyaExt/XephemAstroLib/twobody.c (modified) (6 diffs)

trunk/SophyaExt/XephemAstroLib/twobody.c

@@ -22 +22 @@
  
 /* Constants used when solving Kepler's equation */
+#undef  EPSILON
 #define EPSILON   3E-8
+#undef  INFINITY
 #define INFINITY  1E+10
  
 /* Math constants */
+#undef  PI
 #define PI      3.14159265358979323846
 #define RADEG   ( 180.0 / PI )
@@ -125 +128 @@
  
  
-void vrc( double *v, double *r, double tp, double e, double q )
+/* return 0 if ok, else -1 */
+int vrc( double *v, double *r, double tp, double e, double q )
 /*
  * Elliptic, hyperbolic and near-parabolic orbits:
@@ -145 +149 @@
         *v = 0.0;
         *r = q;
-        return;
+        return 0;
     }
  
@@ -176 +180 @@
                 printf( "\nNear-parabolic orbit: inaccurate result."
                         "\n  e = %f, lambda = %f, w = %f", e, lambda, w );
-                exit(1);
+                return -1;
             }
             else
@@ -202 +206 @@
         *v = 2.0 * atand(w);
         *r = q * (1+w2) / ( 1.0 + w2*lambda );
-        return;  /* Near-parabolic orbit */
+        return 0;  /* Near-parabolic orbit */
     }
  
@@ -232 +236 @@
         *r = q * (1.0+e) / ( 1.0 + e*cosd(*v) );
     }
+    return 0;
  
 } /* vrc */
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: twobody.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.1 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: twobody.c,v $ $Date: 2004-06-15 16:52:41 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};