vsop87.c   vsop87.c 
/************************************************************************ /************************************************************************
The PLANETARY SOLUTION VSOP87 by Bretagnon P. and Francou G. can be found a t The PLANETARY SOLUTION VSOP87 by Bretagnon P. and Francou G. can be found a t
ftp://ftp.imcce.fr/pub/ephem/planets/vsop87 ftp://ftp.imcce.fr/pub/ephem/planets/vsop87
I (Johannes Gajdosik) have just taken the data obtained from above I (Johannes Gajdosik) have just taken the data obtained from above
(VSOP87.mer,...,VSOP87.nep) and rearranged it into this piece of software. (VSOP87.mer,...,VSOP87.nep) and rearranged it into this piece of software.
I can neigther allow nor forbid the usage of VSOP87. I can neither allow nor forbid the usage of VSOP87.
The copyright notice below covers not the work of Bretagnon P. and Francou G. The copyright notice below covers not the work of Bretagnon P. and Francou G.
but just my work, that is the compilation of the VSOP87 data but just my work, that is the compilation of the VSOP87 data
into the software supplied in this file. into the software supplied in this file.
Copyright (c) 2006 Johannes Gajdosik Copyright (c) 2006 Johannes Gajdosik
Permission is hereby granted, free of charge, to any person obtaining a Permission is hereby granted, free of charge, to any person obtaining a
copy of this software and associated documentation files (the "Software"), copy of this software and associated documentation files (the "Software"),
to deal in the Software without restriction, including without limitation to deal in the Software without restriction, including without limitation
the rights to use, copy, modify, merge, publish, distribute, sublicense, the rights to use, copy, modify, merge, publish, distribute, sublicense,
skipping to change at line 137317 skipping to change at line 137317
/* /*
for (i=0;i<8;i++) { for (i=0;i<8;i++) {
printf("%f %d %15.10f %15.10f %15.10f %15.10f %15.10f %15.10f\n", printf("%f %d %15.10f %15.10f %15.10f %15.10f %15.10f %15.10f\n",
2451545.0+t*365250.0,i, 2451545.0+t*365250.0,i,
elem[i*6+0],elem[i*6+1],elem[i*6+2], elem[i*6+0],elem[i*6+1],elem[i*6+2],
elem[i*6+3],elem[i*6+3],elem[i*6+4]); elem[i*6+3],elem[i*6+3],elem[i*6+4]);
} }
*/ */
} }
/* dirty caching in static variables */ /* dirty caching in static variables
If you ever want to allow parallel execution,
make a struct from these and malloc such structs and add pointer argumen
ts to the calls as needed.
*/
#define VSOP87_DIM (8*6) #define VSOP87_DIM (8*6)
static double t_0 = -1e100; static double t_0 = -1e100;
static double t_1 = -1e100; static double t_1 = -1e100;
static double t_2 = -1e100; static double t_2 = -1e100;
static double vsop87_elem_0[VSOP87_DIM]; static double vsop87_elem_0[VSOP87_DIM];
static double vsop87_elem_1[VSOP87_DIM]; static double vsop87_elem_1[VSOP87_DIM];
static double vsop87_elem_2[VSOP87_DIM]; static double vsop87_elem_2[VSOP87_DIM];
/* 10 days: */ /* 10 days: */
#define DELTA_T (10.0/365250.0) #define DELTA_T (10.0/365250.0)
 End of changes. 2 change blocks. 
2 lines changed or deleted 6 lines changed or added

This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/