Return to elliptic.c CVS log

Up to [Research Unix] / researchv8dc / cmd / map / libmap

Annotation of researchv8dc/cmd/map/libmap/elliptic.c, revision 1.1

1.1     ! root        1: #include "map.h"
        !             2: 
        !             3: struct coord center;
        !             4: 
        !             5: Xelliptic(p,x,y)
        !             6: struct place *p;
        !             7: float *x,*y;
        !             8: {
        !             9:        float r1,r2;
        !            10:        r1 = acos(p->nlat.c*(p->wlon.c*center.c
        !            11:                - p->wlon.s*center.s));
        !            12:        r2 = acos(p->nlat.c*(p->wlon.c*center.c
        !            13:                + p->wlon.s*center.s));
        !            14:        *x = -(r1*r1 - r2*r2)/(4*center.l);
        !            15:        *y = (r1*r1+r2*r2)/2 - (center.l*center.l+*x**x);
        !            16:        *y = sqrt(*y);
        !            17:        if(p->nlat.l<0)
        !            18:                *y = -*y;
        !            19:        return(1);
        !            20: }
        !            21: 
        !            22: int (*elliptic(l))()
        !            23: float l;
        !            24: {
        !            25:        int Xazequidistant();
        !            26:        l = fabs(l);
        !            27:        if(l>89)
        !            28:                return(0);
        !            29:        if(l<1)
        !            30:                return(Xazequidistant);
        !            31:        deg2rad(l,&center);
        !            32:        return(Xelliptic);
        !            33: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.