Return to gilbert.c CVS log

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

Annotation of researchv8dc/cmd/map/libmap/gilbert.c, revision 1.1.1.1

1.1       root        1: #include "map.h"
                      2: 
                      3: Xgilbert(place,x,y)
                      4: struct place *place;
                      5: float *x,*y;
                      6: {
                      7:        float z1,z2;
                      8:        double w1,w2,t1,t2;
                      9:        struct place p;
                     10:        copyplace(place,&p);
                     11:        if(place->nlat.l<0) {
                     12:                p.nlat.l = -p.nlat.l;
                     13:                p.nlat.s = -p.nlat.s;
                     14:        }
                     15:        Xstereographic(&p,&z1,&z2);
                     16:        csqrt(-z2/2,z1/2,&w1,&w2);
                     17:        cdiv(w1-1,w2,w1+1,w2,&t1,&t2);
                     18:        *y = -t1;
                     19:        *x = t2;
                     20:        if(place->nlat.l<0)
                     21:                *y = -*y;
                     22:        return(1);
                     23: }
                     24: 
                     25: int (*gilbert())()
                     26: {
                     27:        return(Xgilbert);
                     28: }

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