Changeset 1681 in Sophya

Timestamp:

Oct 11, 2001, 2:40:41 PM (24 years ago)

Author:

cmv

Message:

correction d un bug d adressage t_abs[-1] cmv 11/10/01

File:

• trunk/SophyaExt/XephemAstroLib/vsop87.c (modified) (2 diffs)

trunk/SophyaExt/XephemAstroLib/vsop87.c

@@ -152 +152 @@
 
             /* precision threshold */
-            p = q/(t_abs[alpha] + alpha * t_abs[alpha-1] * 1e-4 + 1e-35);
+            if(alpha>0) p = t_abs[alpha-1]; else p=0.;
+            p = q/(t_abs[alpha] + alpha * p * 1e-4 + 1e-35);
 #if VSOP_SPHERICAL
             if (cooidx == 2)    /* scale by semimajor axis for radius */
@@ -211 +212 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2001-10-11 12:40:41 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};