Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/reduce.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/reduce.c (modified) (4 diffs)

trunk/SophyaExt/XephemAstroLib/reduce.c

@@ -1 +1 @@
 #include <math.h>
 
-#include "P_.h"
+#include <stdio.h>
+
 #include "astro.h"
 
-/* convert those orbital elements that change from epoch mjd0 to epoch mjd.
+/* convert those orbital elements that change from epoch mj0 to epoch mj.
  */
 void
-reduce_elements (mjd0, mjd, inc0, ap0, om0, inc, ap, om)
-double mjd0;    /* initial epoch */
-double mjd;     /* desired epoch */
-double inc0;    /* initial inclination, rads */
-double ap0;     /* initial argument of perihelion, as an mjd */
-double om0;     /* initial long of ascending node, rads */
-double *inc;    /* desired inclination, rads */
-double *ap;     /* desired epoch of perihelion, as an mjd */
-double *om;     /* desired long of ascending node, rads */
+reduce_elements (
+double mj0,     /* initial epoch */
+double mj,      /* desired epoch */
+double inc0,    /* initial inclination, rads */
+double ap0,     /* initial argument of perihelion, as an mj */
+double om0,     /* initial long of ascending node, rads */
+double *inc,    /* resultant inclination, rads */
+double *ap,     /* resultant arg of perihelion, as an mj */
+double *om)     /* resultant long of ascending node, rads */
 {
         double t0, t1;
@@ -26 +27 @@
         double seta, ceta;
 
-        if (fabs(mjd - mjd0) < 1e-5) {
+        if (fabs(mj - mj0) < 1e-5) {
             /* sin(eta) blows for inc < 10 degrees -- anyway, no need */
             *inc = inc0;
@@ -34 +35 @@
         }
 
-        t0 = mjd0/365250.0;
-        t1 = mjd/365250.0;
+        t0 = mj0/365250.0;
+        t1 = mj/365250.0;
 
         tt = t1-t0;
@@ -75 +76 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2004-06-15 16:52:40 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};