Return to guyou.c CVS log

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

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

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

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