| Server IP : 82.208.35.60 / Your IP : 216.73.216.238 Web Server : Apache/2.4.55 (FreeBSD) OpenSSL/1.1.1q-freebsd PHP/7.3.31 System : FreeBSD server7.d2m.cz 12.4-RELEASE-p9 FreeBSD 12.4-RELEASE-p9 GENERIC amd64 User : studiokobylisy_cz ( 1008) PHP Version : 7.3.31 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : OFF | Perl : OFF | Python : OFF | Sudo : OFF | Pkexec : OFF Directory : /usr/ports/astro/libnova/files/ |
Upload File : |
--- src/precession.c.orig 2011-07-28 11:02:50 UTC
+++ src/precession.c
@@ -82,21 +82,21 @@ void ln_get_equ_prec (struct ln_equ_posn
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 */
@@ -140,21 +140,21 @@ void ln_get_equ_prec2 (struct ln_equ_pos
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 */