PageRenderTime 25ms CodeModel.GetById 15ms RepoModel.GetById 0ms app.codeStats 0ms

/Proj4/aasincos.c

http://github.com/route-me/route-me
C | 36 lines | 33 code | 2 blank | 1 comment | 5 complexity | 824de4698703f9f1afb798a92399238e MD5 | raw file
    

  1. /* arc sin, cosine, tan2 and sqrt that will NOT fail */
    2. 

  2. #ifndef lint
    3. 

  3. static const char SCCSID[]="@(#)aasincos.c	4.6	93/12/12	GIE	REL";
    4. 

  4. #endif
    5. 

  5. #include "projects.h"
    6. 

  6. #define ONE_TOL	 1.00000000000001
    7. 

  7. #define TOL	0.000000001
    8. 

  8. #define ATOL 1e-50
    9. 

  9. 	double
    10. 

  10. aasin(double v) {
    11. 

  11. 	double av;
    12. 

  12. 

  13. 	if ((av = fabs(v)) >= 1.) {
    14. 

  14. 		if (av > ONE_TOL)
    15. 

  15. 			pj_errno = -19;
    16. 

  16. 		return (v < 0. ? -HALFPI : HALFPI);
    17. 

  17. 	}
    18. 

  18. 	return asin(v);
    19. 

  19. }
    20. 

  20. 	double
    21. 

  21. aacos(double v) {
    22. 

  22. 	double av;
    23. 

  23. 

  24. 	if ((av = fabs(v)) >= 1.) {
    25. 

  25. 		if (av > ONE_TOL)
    26. 

  26. 			pj_errno = -19;
    27. 

  27. 		return (v < 0. ? PI : 0.);
    28. 

  28. 	}
    29. 

  29. 	return acos(v);
    30. 

  30. }
    31. 

  31. 	double
    32. 

  32. asqrt(double v) { return ((v <= 0) ? 0. : sqrt(v)); }
    33. 

  33. 	double
    34. 

  34. aatan2(double n, double d) {
    35. 

  35. 	return ((fabs(n) < ATOL && fabs(d) < ATOL) ? 0. : atan2(n,d));
    36. 

  36. }
    37. 

  37.