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/ |