180 lines
7.0 KiB
C
180 lines
7.0 KiB
C
|
|
/************************************************************
|
|
$Header: /home/dieter/sweph/RCS/swephlib.h,v 1.74 2008/06/16 10:07:20 dieter Exp $
|
|
|
|
Authors: Dieter Koch and Alois Treindl, Astrodienst Zurich
|
|
|
|
************************************************************/
|
|
/* Copyright (C) 1997 - 2008 Astrodienst AG, Switzerland. All rights reserved.
|
|
|
|
License conditions
|
|
------------------
|
|
|
|
This file is part of Swiss Ephemeris.
|
|
|
|
Swiss Ephemeris is distributed with NO WARRANTY OF ANY KIND. No author
|
|
or distributor accepts any responsibility for the consequences of using it,
|
|
or for whether it serves any particular purpose or works at all, unless he
|
|
or she says so in writing.
|
|
|
|
Swiss Ephemeris is made available by its authors under a dual licensing
|
|
system. The software developer, who uses any part of Swiss Ephemeris
|
|
in his or her software, must choose between one of the two license models,
|
|
which are
|
|
a) GNU public license version 2 or later
|
|
b) Swiss Ephemeris Professional License
|
|
|
|
The choice must be made before the software developer distributes software
|
|
containing parts of Swiss Ephemeris to others, and before any public
|
|
service using the developed software is activated.
|
|
|
|
If the developer choses the GNU GPL software license, he or she must fulfill
|
|
the conditions of that license, which includes the obligation to place his
|
|
or her whole software project under the GNU GPL or a compatible license.
|
|
See http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
|
|
|
|
If the developer choses the Swiss Ephemeris Professional license,
|
|
he must follow the instructions as found in http://www.astro.com/swisseph/
|
|
and purchase the Swiss Ephemeris Professional Edition from Astrodienst
|
|
and sign the corresponding license contract.
|
|
|
|
The License grants you the right to use, copy, modify and redistribute
|
|
Swiss Ephemeris, but only under certain conditions described in the License.
|
|
Among other things, the License requires that the copyright notices and
|
|
this notice be preserved on all copies.
|
|
|
|
Authors of the Swiss Ephemeris: Dieter Koch and Alois Treindl
|
|
|
|
The authors of Swiss Ephemeris have no control or influence over any of
|
|
the derived works, i.e. over software or services created by other
|
|
programmers which use Swiss Ephemeris functions.
|
|
|
|
The names of the authors or of the copyright holder (Astrodienst) must not
|
|
be used for promoting any software, product or service which uses or contains
|
|
the Swiss Ephemeris. This copyright notice is the ONLY place where the
|
|
names of the authors can legally appear, except in cases where they have
|
|
given special permission in writing.
|
|
|
|
The trademarks 'Swiss Ephemeris' and 'Swiss Ephemeris inside' may be used
|
|
for promoting such software, products or services.
|
|
*/
|
|
|
|
|
|
/* Set TRUE, to include Herring's (1987) corrections to IAU 1980
|
|
* nutation series. AA (1996) neglects them. */
|
|
#define NUT_CORR_1987 FALSE
|
|
|
|
/* Precession coefficients for remote past and future.
|
|
* One of the following four defines must be true.
|
|
*/
|
|
#define PREC_VONDRAK_2011 TRUE
|
|
#define PREC_WILLIAMS_1994 FALSE
|
|
#define PREC_SIMON_1994 FALSE
|
|
#define PREC_LASKAR_1986 FALSE
|
|
#define PREC_BRETAGNON_2003 FALSE
|
|
/* IAU precession 1976 or 2003 for recent centuries.
|
|
* only one of the following two defines may be TRUE */
|
|
#define PREC_IAU_1976 FALSE
|
|
#define PREC_IAU_2003 FALSE /* precession model P03 */
|
|
#define PREC_IAU_1976_CTIES 2.0 /* J2000 +/- two centuries */
|
|
/* we use P03 for whole ephemeris */
|
|
#define PREC_IAU_2003_CTIES 75.0 /* J2000 +/- 75 centuries */
|
|
|
|
/* choose between the following nutation models */
|
|
#define NUT_IAU_1980 FALSE
|
|
#define NUT_IAU_2000A FALSE /* very time consuming ! */
|
|
#define NUT_IAU_2000B TRUE /* fast, but precision of milli-arcsec */
|
|
|
|
/* coordinate transformation */
|
|
extern void swi_coortrf(double *xpo, double *xpn, double eps);
|
|
|
|
/* coordinate transformation */
|
|
extern void swi_coortrf2(double *xpo, double *xpn, double sineps, double coseps);
|
|
|
|
/* cartesian to polar coordinates */
|
|
extern void swi_cartpol(double *x, double *l);
|
|
|
|
/* cartesian to polar coordinates with velocity */
|
|
extern void swi_cartpol_sp(double *x, double *l);
|
|
extern void swi_polcart_sp(double *l, double *x);
|
|
|
|
/* polar to cartesian coordinates */
|
|
extern void swi_polcart(double *l, double *x);
|
|
|
|
/* GCRS to J2000 */
|
|
extern void swi_bias(double *x, int32 iflag, AS_BOOL backward);
|
|
/* GCRS to FK5 */
|
|
extern void swi_icrs2fk5(double *x, int32 iflag, AS_BOOL backward);
|
|
|
|
/* precession */
|
|
extern int swi_precess(double *R, double J, int direction );
|
|
extern void swi_precess_speed(double *xx, double t, int direction);
|
|
|
|
/* from sweph.c, light deflection, aberration, etc. */
|
|
extern void swi_deflect_light(double *xx, double dt, int32 iflag);
|
|
extern void swi_aberr_light(double *xx, double *xe, int32 iflag);
|
|
extern int swi_plan_for_osc_elem(int32 iflag, double tjd, double *xx);
|
|
extern int swi_trop_ra2sid_lon(double *xin, double *xout, double *xoutr, int32 iflag, char *serr);
|
|
extern int swi_trop_ra2sid_lon_sosy(double *xin, double *xout, double *xoutr, int32 iflag, char *serr);
|
|
extern int swi_get_observer(double tjd, int32 iflag,
|
|
AS_BOOL do_save, double *xobs, char *serr);
|
|
extern void swi_force_app_pos_etc();
|
|
|
|
/* obliquity of ecliptic */
|
|
extern void swi_check_ecliptic(double tjd);
|
|
extern double swi_epsiln(double J);
|
|
extern void swi_ldp_peps(double J, double *dpre, double *deps);
|
|
|
|
/* nutation */
|
|
extern void swi_check_nutation(double tjd, int32 iflag);
|
|
extern int swi_nutation(double J, double *nutlo);
|
|
extern void swi_nutate(double *xx, int32 iflag, AS_BOOL backward);
|
|
|
|
extern void swi_mean_lunar_elements(double tjd,
|
|
double *node, double *dnode,
|
|
double *peri, double *dperi);
|
|
/* */
|
|
extern double swi_mod2PI(double x);
|
|
|
|
/* evaluation of chebyshew series and derivative */
|
|
extern double swi_echeb(double x, double *coef, int ncf);
|
|
extern double swi_edcheb(double x, double *coef, int ncf);
|
|
|
|
/* cross product of vectors */
|
|
extern void swi_cross_prod(double *a, double *b, double *x);
|
|
/* dot product of vecotrs */
|
|
extern double swi_dot_prod_unit(double *x, double *y);
|
|
|
|
extern double swi_angnorm(double x);
|
|
|
|
/* generation of SWISSEPH file names */
|
|
extern void swi_gen_filename(double tjd, int ipli, char *fname);
|
|
|
|
/* cyclic redundancy checksum (CRC), 32 bit */
|
|
extern uint32 swi_crc32(unsigned char *buf, int len);
|
|
|
|
extern int swi_cutstr(char *s, char *cutlist, char *cpos[], int nmax);
|
|
extern char *swi_right_trim(char *s);
|
|
|
|
extern double swi_kepler(double E, double M, double ecce);
|
|
|
|
extern char *swi_get_fict_name(int32 ipl, char *s);
|
|
|
|
extern void swi_FK4_FK5(double *xp, double tjd);
|
|
|
|
extern char *swi_strcpy(char *to, char *from);
|
|
extern char *swi_strncpy(char *to, char *from, size_t n);
|
|
|
|
#ifdef TRACE
|
|
# define TRACE_COUNT_MAX 10000
|
|
extern FILE *swi_fp_trace_c;
|
|
extern FILE *swi_fp_trace_out;
|
|
extern void swi_open_trace(char *serr);
|
|
extern int32 swi_trace_count;
|
|
static char *fname_trace_c = "swetrace.c";
|
|
static char *fname_trace_out = "swetrace.txt";
|
|
#ifdef FORCE_IFLAG
|
|
static char *fname_force_flg = "force.flg";
|
|
#endif
|
|
#endif /* TRACE */
|