Changeset 1719 in Sophya for trunk/SophyaExt/XephemAstroLib/eq_ecl.c

Timestamp:

Oct 22, 2001, 2:08:28 PM (24 years ago)

Author:

cmv

Message:

Adapted to version 3.5 xephem cmv 22/10/2001

File:

• trunk/SophyaExt/XephemAstroLib/eq_ecl.c (modified) (3 diffs)

trunk/SophyaExt/XephemAstroLib/eq_ecl.c

@@ -51 +51 @@
         static double lastmjd = -10000; /* last mjd calculated */
         static double seps, ceps;       /* sin and cos of mean obliquity */
-        double sx, cx, sy, cy, ty;
+        double sx, cx, sy, cy, ty, sq;
 
         if (mjd != lastmjd) {
@@ -67 +67 @@
         cx = cos(x);
         sx = sin(x);
-        /* modifs CMV *q = asin((sy*ceps)-(cy*seps*sx*sw)); */
-        *q = (sy*ceps)-(cy*seps*sx*sw);
-        if(*q<-1.) *q = -PI/2.; else if(*q>1.) *q = PI/2.; else *q = asin(*q);
+        sq = (sy*ceps)-(cy*seps*sx*sw);
+        if (sq < -1) sq = -1;
+        if (sq >  1) sq =  1;
+        *q = asin(sq);
         *p = atan(((sx*ceps)+(ty*seps*sw))/cx);
         if (cx<0) *p += PI;             /* account for atan quad ambiguity */
@@ -76 +77 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};