diff options
author | delphij <delphij@FreeBSD.org> | 2006-08-18 18:15:16 +0000 |
---|---|---|
committer | delphij <delphij@FreeBSD.org> | 2006-08-18 18:15:16 +0000 |
commit | c637bf0e7ee1e34bd8ed0c9146f44bc0998f9140 (patch) | |
tree | 0ba4e7bbb0ecdc894a23f4f0483d93fd56a13384 /astro/libnova/files/patch-src_precession.c | |
parent | f24e7b3daa43166950eb13892167f56548817b7f (diff) | |
download | FreeBSD-ports-c637bf0e7ee1e34bd8ed0c9146f44bc0998f9140.zip FreeBSD-ports-c637bf0e7ee1e34bd8ed0c9146f44bc0998f9140.tar.gz |
Update to 0.12.0.
Submitted by: Frank W. Josellis <frank dynamical-systems org>
PR: ports/102119
Diffstat (limited to 'astro/libnova/files/patch-src_precession.c')
-rw-r--r-- | astro/libnova/files/patch-src_precession.c | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/astro/libnova/files/patch-src_precession.c b/astro/libnova/files/patch-src_precession.c new file mode 100644 index 0000000..0f7784f --- /dev/null +++ b/astro/libnova/files/patch-src_precession.c @@ -0,0 +1,58 @@ +--- src/precession.c.orig Tue Jan 10 23:35:08 2006 ++++ src/precession.c Tue Aug 15 18:03:50 2006 +@@ -56,21 +56,21 @@ + theta = ln_deg_to_rad (theta); + + /* calc A,B,C equ 20.4 */ +- A = cosl (mean_dec) * sinl (mean_ra + zeta); +- B = cosl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) - sinl (theta) * sinl (mean_dec); +- C = sinl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) + cosl (theta) * sinl (mean_dec); ++ A = cos (mean_dec) * sin (mean_ra + zeta); ++ B = cos (theta) * cos (mean_dec) * cos (mean_ra + zeta) - sin (theta) * sin (mean_dec); ++ C = sin (theta) * cos (mean_dec) * cos (mean_ra + zeta) + cos (theta) * sin (mean_dec); + +- ra = atan2l (A,B) + eta; ++ ra = atan2 (A,B) + eta; + + /* check for object near celestial pole */ + if (mean_dec > (0.4 * M_PI) || mean_dec < (-0.4 * M_PI)) { + /* close to pole */ +- dec = acosl (sqrt(A * A + B * B)); ++ dec = acos (sqrt(A * A + B * B)); + if (mean_dec < 0.) + dec *= -1; /* 0 <= acos() <= PI */ + } else { + /* not close to pole */ +- dec = asinl (C); ++ dec = asin (C); + } + + /* change to degrees */ +@@ -114,21 +114,21 @@ + theta = ln_deg_to_rad (theta); + + /* calc A,B,C equ 20.4 */ +- A = cosl (mean_dec) * sinl (mean_ra + zeta); +- B = cosl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) - sinl (theta) * sinl (mean_dec); +- C = sinl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) + cosl (theta) * sinl (mean_dec); ++ A = cos (mean_dec) * sin (mean_ra + zeta); ++ B = cos (theta) * cos (mean_dec) * cos (mean_ra + zeta) - sin (theta) * sin (mean_dec); ++ C = sin (theta) * cos (mean_dec) * cos (mean_ra + zeta) + cos (theta) * sin (mean_dec); + +- ra = atan2l (A,B) + eta; ++ ra = atan2 (A,B) + eta; + + /* check for object near celestial pole */ + if (mean_dec > (0.4 * M_PI) || mean_dec < (-0.4 * M_PI)) { + /* close to pole */ +- dec = acosl (sqrt(A * A + B * B)); ++ dec = acos (sqrt(A * A + B * B)); + if (mean_dec < 0.) + dec *= -1; /* 0 <= acos() <= PI */ + } else { + /* not close to pole */ +- dec = asinl (C); ++ dec = asin (C); + } + + /* change to degrees */ |