Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
freebsd
GitHub Repository: freebsd/freebsd-ports-gnome
 Path: blob/main/astro/libnova/files/patch-src_precession.c
16147 views
1
--- src/precession.c.orig	2011-07-28 11:02:50 UTC

2
+++ src/precession.c

3
@@ -82,21 +82,21 @@ void ln_get_equ_prec (struct ln_equ_posn

4
 	theta = ln_deg_to_rad (theta); 

5
 

6
 	/* calc A,B,C equ 20.4 */

7
-	A = cosl (mean_dec) * sinl (mean_ra + zeta);

8
-	B = cosl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) - sinl (theta) * sinl (mean_dec);

9
-	C = sinl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) + cosl (theta) * sinl (mean_dec);

10
+	A = cos (mean_dec) * sin (mean_ra + zeta);

11
+	B = cos (theta) * cos (mean_dec) * cos (mean_ra + zeta) - sin (theta) * sin (mean_dec);

12
+	C = sin (theta) * cos (mean_dec) * cos (mean_ra + zeta) + cos (theta) * sin (mean_dec);

13
 	

14
-	ra = atan2l (A,B) + eta;

15
+	ra = atan2 (A,B) + eta;

16
 	

17
 	/* check for object near celestial pole */

18
 	if (mean_dec > (0.4 * M_PI) || mean_dec < (-0.4 * M_PI)) {

19
 		/* close to pole */

20
-		dec = acosl (sqrt(A * A + B * B));

21
+		dec = acos (sqrt(A * A + B * B));

22
 		if (mean_dec < 0.)

23
 		  dec *= -1; /* 0 <= acos() <= PI */

24
 	} else {

25
 		/* not close to pole */

26
-		dec = asinl (C);

27
+		dec = asin (C);

28
 	}

29
 

30
 	/* change to degrees */

31
@@ -140,21 +140,21 @@ void ln_get_equ_prec2 (struct ln_equ_pos

32
 	theta = ln_deg_to_rad (theta); 

33
 

34
 	/* calc A,B,C equ 20.4 */

35
-	A = cosl (mean_dec) * sinl (mean_ra + zeta);

36
-	B = cosl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) - sinl (theta) * sinl (mean_dec);

37
-	C = sinl (theta) * cosl (mean_dec) * cosl (mean_ra + zeta) + cosl (theta) * sinl (mean_dec);

38
+	A = cos (mean_dec) * sin (mean_ra + zeta);

39
+	B = cos (theta) * cos (mean_dec) * cos (mean_ra + zeta) - sin (theta) * sin (mean_dec);

40
+	C = sin (theta) * cos (mean_dec) * cos (mean_ra + zeta) + cos (theta) * sin (mean_dec);

41
 	

42
-	ra = atan2l (A,B) + eta;

43
+	ra = atan2 (A,B) + eta;

44
 	

45
 	/* check for object near celestial pole */

46
 	if (mean_dec > (0.4 * M_PI) || mean_dec < (-0.4 * M_PI)) {

47
 		/* close to pole */

48
-		dec = acosl (sqrt(A * A + B * B));

49
+		dec = acos (sqrt(A * A + B * B));

50
 		if (mean_dec < 0.)

51
 		  dec *= -1; /* 0 <= acos() <= PI */

52
 	} else {

53
 		/* not close to pole */

54
-		dec = asinl (C);

55
+		dec = asin (C);

56
 	}

57
 

58
 	/* change to degrees */

59


60