Return to newspline CVS log

Up to [Research Unix] / researchv9 / jtools / src / Jpic

Annotation of researchv9/jtools/src/Jpic/newspline, revision 1.1.1.1

1.1       root        1: From theo Wed Aug  6 16:12 EDT 1980
                      2: Brian:
                      3: This is the efficient procedure for splines.
                      4: It requires points from 0 to N-1. To keep the
                      5: same convention as before we must have
                      6: P(0) = P(1)  and P(N-2) = P(N-1).
                      7: spline(){
                      8:        int i,j,ILEN=10,I,Xs,Ys;
                      9:        double X1, X2, X3, Y1, Y2, Y3;
                     10:        double xi,u,LEN=10.0;
                     11:        I=0;
                     12:        for(i=0;i< N-2; i++){
                     13:                xi=i;
                     14:                X1 = (X[i+2]+X[i])/2 - X[i+1];
                     15:                X2 = X[i+1]-X[i];
                     16:                X3 = (X[i+1]+X[i])/2;
                     17:                Y1 = (Y[i+2]+Y[i])/2 - Y[i+1];
                     18:                Y2 = Y[i+1]-Y[i];
                     19:                Y3 = (Y[i+1]+Y[i])/2;
                     20:                for(j=0;j<ILEN;j++){
                     21:                        u = I; u = u/LEN - xi;
                     22:                        Xs = (X1*u + X2)*u + X3;
                     23:                        Ys = (Y1*u + Y2)*u + Y3;
                     24:                        if(!I) move(Xs,Ys); else cont(Xs,Ys);
                     25:                        I++;
                     26:                        }
                     27:                }
                     28:        }
                     29: 
                     30: 

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