Package com.nr.interp

Examples of com.nr.interp.Curve_interp


    }
    target[0]=2.5*vec[0];
    target[1]=2.5*vec[1];
    target[2]=2.5*vec[2];
//    System.out.printf(target[0] << " " << target[1] << " " << target[2] << endl;
    Curve_interp myCurve = new Curve_interp(ptsin);
    double[] f=myCurve.interp(2.5/(NPTS-1))
//    System.out.printf(f[0]/target[0] << " " << f[1]/target[1] << " " << f[2]/target[2] << endl;
    sbeps=1.e-3;
    System.out.printf("     Discrepancy (straight line): %f\n", maxel(vecsub(f,target)));
    localflag = maxel(vecsub(f,target)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
      fail("*** Curve_interp: Does not work for data on a 3D straight line.");
     
    }

    // Test closed circle
    double a[]={0.0,vec[2],-vec[1]}// vec x (1,0,0)
    double norm=sqrt(SQR(vec[2])+SQR(vec[1]))// normalize a[]
    for (i=0;i<NDIM;i++) a[i] /= norm;
    double b[]={vec[1]*a[2]-vec[2]*a[1],vec[2]*a[0]-vec[0]*a[2],vec[0]*a[1]-vec[1]*a[0]}; // perpendicular to a[] and vec[]
    double pi=acos(-1.0);
    double theta;
    for (i=0;i<NPTS;i++) {
      theta=i*2.0*pi/NPTS;
      ptsin[i][0]=a[0]*cos(theta)+b[0]*sin(theta);
      ptsin[i][1]=a[1]*cos(theta)+b[1]*sin(theta);
      ptsin[i][2]=a[2]*cos(theta)+b[2]*sin(theta);
    }
    theta=2.5*2.0*pi/NPTS;
    target[0]=a[0]*cos(theta)+b[0]*sin(theta);
    target[1]=a[1]*cos(theta)+b[1]*sin(theta);
    target[2]=a[2]*cos(theta)+b[2]*sin(theta);
//    System.out.printf(target[0] << " " << target[1] << " " << target[2] << endl;
    close=true;
    Curve_interp myCurve2 = new Curve_interp(ptsin,close);
    f=myCurve2.interp(2.5/NPTS);
//    System.out.printf(f[0]/target[0] << " " << f[1]/target[1] << " " << f[2]/target[2] << endl;
    sbeps=1.e-4;
    System.out.printf("     Discrepancy (circle): %f\n", maxel(vecsub(f,target)));
    localflag = maxel(vecsub(f,target)) > sbeps;
    globalflag = globalflag || localflag;
View Full Code Here

TOP

Related Classes of com.nr.interp.Curve_interp

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.