Return to trapezoidal.c CVS log

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

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

1.1       root        1: #include "map.h"
                      2: 
                      3: static struct coord stdpar0, stdpar1;
                      4: static double k;
                      5: static double yeq;
                      6: 
                      7: static int
                      8: Xtrapezoidal(struct place *place, double *x, double *y)
                      9: {
                     10:        *y = yeq + place->nlat.l;
                     11:        *x = *y*k*place->wlon.l;
                     12:        return 1;
                     13: }
                     14: 
                     15: proj
                     16: trapezoidal(double par0, double par1)
                     17: {
                     18:        if(fabs(fabs(par0)-fabs(par1))<.1)
                     19:                return rectangular(par0);
                     20:        deg2rad(par0,&stdpar0);
                     21:        deg2rad(par1,&stdpar1);
                     22:        if(fabs(par1-par0) < .1)
                     23:                k = stdpar1.s;
                     24:        else
                     25:                k = (stdpar1.c-stdpar0.c)/(stdpar0.l-stdpar1.l);
                     26:        yeq = -stdpar1.l - stdpar1.c/k;
                     27:        return Xtrapezoidal;
                     28: }

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