vsop87.h   vsop87.h 
/************************************************************************ /************************************************************************
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 43 skipping to change at line 43
This is my implementation of the VSOP87 planetary solution. This is my implementation of the VSOP87 planetary solution.
I tried to optimize for speed by rearranging the terms so that I tried to optimize for speed by rearranging the terms so that
for a given touple of (a[0],a[1],...,a[11]) the values for a given touple of (a[0],a[1],...,a[11]) the values
(cos,sin)(a[0]*lambda[0](T)+...a[11]*lambda[11](T)) (cos,sin)(a[0]*lambda[0](T)+...a[11]*lambda[11](T))
have only to be calculated once. have only to be calculated once.
Furthermore I used the addition formulas Furthermore I used the addition formulas
(cos,sin)(x+y) = ... (cos,sin)(x+y) = ...
so that for given T the functions cos and sin have only to be called 12 tim es. so that for given T the functions cos and sin have only to be called 12 tim es.
ATTENTION! Due to static caching this solution is not reentrant and cannot
be parallelized to run in several threads.
****************************************************************/ ****************************************************************/
#ifndef _VSOP87_H_ #ifndef _VSOP87_H_
#define _VSOP87_H_ #define _VSOP87_H_
#ifdef __cplusplus #ifdef __cplusplus
extern "C" { extern "C" {
#endif #endif
void GetVsop87Coor(double jd,int body,double *xyz); void GetVsop87Coor(double jd,int body,double *xyz);
 End of changes. 2 change blocks. 
1 lines changed or deleted 4 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/