/Proj4/pj_auth.c

http://github.com/route-me/route-me · C · 34 lines · 32 code · 1 blank · 1 comment · 1 complexity · 7aa5fcdf573190d6e88be0f01dbce4ec MD5 · raw file 

    

  1. /* determine latitude from authalic latitude */
    2. 

  2. #ifndef lint
    3. 

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

  4. #endif
    5. 

  5. #include "projects.h"
    6. 

  6. # define P00 .33333333333333333333
    7. 

  7. # define P01 .17222222222222222222
    8. 

  8. # define P02 .10257936507936507936
    9. 

  9. # define P10 .06388888888888888888
    10. 

  10. # define P11 .06640211640211640211
    11. 

  11. # define P20 .01641501294219154443
    12. 

  12. #define APA_SIZE 3
    13. 

  13. 	double *
    14. 

  14. pj_authset(double es) {
    15. 

  15. 	double t, *APA;
    16. 

  16. 

  17. 	APA = (double *)pj_malloc(APA_SIZE * sizeof(double));
    18. 

  18. 	if (APA) {
    19. 

  19. 		APA[0] = es * P00;
    20. 

  20. 		t = es * es;
    21. 

  21. 		APA[0] += t * P01;
    22. 

  22. 		APA[1] = t * P10;
    23. 

  23. 		t *= es;
    24. 

  24. 		APA[0] += t * P02;
    25. 

  25. 		APA[1] += t * P11;
    26. 

  26. 		APA[2] = t * P20;
    27. 

  27. 	}
    28. 

  28. 	return APA;
    29. 

  29. }
    30. 

  30. 	double
    31. 

  31. pj_authlat(double beta, double *APA) {
    32. 

  32. 	double t = beta+beta;
    33. 

  33. 	return(beta + APA[0] * sin(t) + APA[1] * sin(t+t) + APA[2] * sin(t+t+t));
    34. 

  34. }
    35.