Return to guyou.c CVS log

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

Annotation of researchv10no/cmd/map/libmap/guyou.c, revision 1.1.1.1

1.1       root        1: #include "map.h"
                      2: 
                      3: static struct place gywhem, gyehem;
                      4: static struct coord gytwist;
                      5: static double gyconst, gykc, gyside;
                      6: 
                      7: 
                      8: static void
                      9: dosquare(double z1, double z2, double *x, double *y)
                     10: {
                     11:        double w1,w2;
                     12:        w1 = z1 -1;
                     13:        if(fabs(w1*w1+z2*z2)>.000001) {
                     14:                cdiv(z1+1,z2,w1,z2,&w1,&w2);
                     15:                w1 *= gyconst;
                     16:                w2 *= gyconst;
                     17:                if(w1<0)
                     18:                        w1 = 0;
                     19:                elco2(w1,w2,gykc,1.,1.,x,y);
                     20:        } else {
                     21:                *x = gyside;
                     22:                *y = 0;
                     23:        }
                     24: }
                     25: 
                     26: int
                     27: Xguyou(struct place *place, double *x, double *y)
                     28: {
                     29:        int ew;         /*which hemisphere*/
                     30:        double z1,z2;
                     31:        struct place pl;
                     32:        ew = place->wlon.l<0;
                     33:        copyplace(place,&pl);
                     34:        norm(&pl,ew?&gyehem:&gywhem,&gytwist);
                     35:        Xstereographic(&pl,&z1,&z2);
                     36:        dosquare(z1/2,z2/2,x,y);
                     37:        if(!ew)
                     38:                *x -= gyside;
                     39:        return(1);
                     40: }
                     41: 
                     42: proj
                     43: guyou(void)
                     44: {
                     45:        double junk;
                     46:        gykc = 1/(3+2*sqrt(2.));
                     47:        gyconst = -(1+sqrt(2.));
                     48:        elco2(-gyconst,0.,gykc,1.,1.,&gyside,&junk);
                     49:        gyside *= 2;
                     50:        latlon(0.,90.,&gywhem);
                     51:        latlon(0.,-90.,&gyehem);
                     52:        deg2rad(0.,&gytwist);
                     53:        return(Xguyou);
                     54: }
                     55: 
                     56: int
                     57: guycut(struct place *g, struct place *og, double *cutlon)
                     58: {
                     59:        int c;
                     60:        c = picut(g,og,cutlon);
                     61:        if(c!=1)
                     62:                return(c);
                     63:        *cutlon = 0.;
                     64:        if(g->nlat.c<.7071||og->nlat.c<.7071)
                     65:                return(ckcut(g,og,0.));
                     66:        return(1);
                     67: }
                     68: 
                     69: static int
                     70: Xsquare(struct place *place, double *x, double *y)
                     71: {
                     72:        double z1,z2;
                     73:        double r, theta;
                     74:        struct place p;
                     75:        copyplace(place,&p);
                     76:        if(place->nlat.l<0) {
                     77:                p.nlat.l = -p.nlat.l;
                     78:                p.nlat.s = -p.nlat.s;
                     79:        }
                     80:        if(p.nlat.l<FUZZ && fabs(p.wlon.l)>PI-FUZZ){
                     81:                *y = -gyside/2;
                     82:                *x = p.wlon.l>0?0:gyside;
                     83:                return(1);
                     84:        }
                     85:        Xstereographic(&p,&z1,&z2);
                     86:        r = sqrt(sqrt(hypot(z1,z2)/2));
                     87:        theta = atan2(z1,-z2)/4;
                     88:        dosquare(r*sin(theta),-r*cos(theta),x,y);
                     89:        if(place->nlat.l<0)
                     90:                *y = -gyside - *y;
                     91:        return(1);
                     92: }
                     93: 
                     94: proj
                     95: square(void)
                     96: {
                     97:        guyou();
                     98:        return(Xsquare);
                     99: }

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