Statistics
| Revision:

svn-gvsig-desktop / tags / v1_1_2_Build_1041 / libraries / libjni-proj4 / src / PJ_august.c @ 33839

History | View | Annotate | Download (639 Bytes)

1 
#ifndef lint
2 
static const char SCCSID[]="@(#)PJ_august.c        4.1        94/02/15        GIE        REL";
3 
#endif
4 
#define PJ_LIB__
5 
#include        <projects.h>
6 
PROJ_HEAD(august, "August Epicycloidal") "\n\tMisc Sph, no inv.";
7 
#define M 1.333333333333333
8 
FORWARD(s_forward); /* spheroid */
9 
        double t, c1, c, x1, x12, y1, y12;
10 

    		
  

  

    
      11
    
    
      
        t = tan(.5 * lp.phi);


    		
  

  

    
      12
    
    
      
        c1 = sqrt(1. - t * t);


    		
  

  

    
      13
    
    
      
        c = 1. + c1 * cos(lp.lam *= .5);


    		
  

  

    
      14
    
    
      
        x1 = sin(lp.lam) *  c1 / c;


    		
  

  

    
      15
    
    
      
        y1 =  t / c;


    		
  

  

    
      16
    
    
      
        xy.x = M * x1 * (3. + (x12 = x1 * x1) - 3. * (y12 = y1 *  y1));


    		
  

  

    
      17
    
    
      
        xy.y = M * y1 * (3. + 3. * x12 - y12);


    		
  

  

    
      18
    
    
      
        return (xy);


    		
  

  

    
      19
    
    
      
}


    		
  

  

    
      20
    
    
      
FREEUP; if (P) pj_dalloc(P); }


    		
  

  

    
      21
    
    
      
ENTRY0(august) P->inv = 0; P->fwd = s_forward; P->es = 0.; ENDENTRY(P)