Examples of interp()

Examples of com.nr.interp.RBF_interp.interp()

92
93
94
95
96
97
98
99
100
101
102
    RBF_gauss gauss = new RBF_gauss(scale);
    RBF_interp myRBFgf = new RBF_interp (pts,y,gauss,false);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFgf.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.RBF_interp.interp()

109
110
111
112
113
114
115
116
117
118
119
    RBF_inversemultiquadric inversemultiquadric = new RBF_inversemultiquadric(scale);
    RBF_interp myRBFimqf =new RBF_interp(pts,y,inversemultiquadric,false);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFimqf.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.RBF_interp.interp()

126
127
128
129
130
131
132
133
134
135
136
    System.out.println("Testing RBF_interp with multiquadric function");
    RBF_interp myRBFmqt = new RBF_interp(pts,y,multiquadric,true);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFmqt.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.RBF_interp.interp()

142
143
144
145
146
147
148
149
150
151
152
    scale=2.0;
    RBF_interp myRBFtpt =new RBF_interp(pts,y,thinplate,true);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFtpt.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.RBF_interp.interp()

158
159
160
161
162
163
164
165
166
167
168
    scale=5.0;
    RBF_interp myRBFgt = new RBF_interp(pts,y,gauss,true);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFgt.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.RBF_interp.interp()

174
175
176
177
178
179
180
181
182
183
184
    scale=2.0;
    RBF_interp myRBFimqt = new RBF_interp(pts,y,inversemultiquadric,true);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=myRBFimqt.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.Rational_interp.interp()

38
39
40
41
42
43
44
45
46
47
48
    }
    Ran myran = new Ran(17);
    Rational_interp z = new Rational_interp(x,y,3);
    for (i=0;i<N;i++) {
      xx[i]=myran.doub();
      yy[i]=z.interp(xx[i]);      // interpolated values
      dyy[i]=z.dy;          // Estimated error
      zz[i]=(1.0-xx[i]*(0.5-0.1*xx[i]))/(1.0+(xx[i]-2.0)*(xx[i]-2.0)); // Actual values
    }
    System.out.printf("     Rat_interp: Max. estimated error: %f\n", maxel(dyy));
    System.out.printf("     Rat_interp: Max. actual error:    %f\n", maxel(vecsub(zz,yy)));
View Full Code Here

Examples of com.nr.interp.Shep_interp.interp()

51
52
53
54
55
56
57
58
59
60
61
    }
    Shep_interp shep =new Shep_interp(pts,y,p);
    for (i=0;i<M;i++) {
      ppt[0]=pt[i][0];
      ppt[1]=pt[i][1];
      estim[i]=shep.interp(ppt);
    }
    System.out.printf("     Discrepancy: %f\n", maxel(vecsub(actual,estim)));
    localflag = maxel(vecsub(actual,estim)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
View Full Code Here

Examples of com.nr.interp.Spline2D_interp.interp()

35
36
37
38
39
40
41
42
43
44
45
    System.out.println("Testing Spline2D_interp");
    for(i=0;i<4;i++)
      for(j=0;j<4;j++)
        y[i][j]=xx1[i]+xx2[j];
    Spline2D_interp z = new Spline2D_interp(x1,x2,y);
    localflag = abs(z.interp(2.5,2.5)-5.0) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
      fail("*** Spline2D_interp: .");
    }

View Full Code Here

Examples of com.nr.interp.Spline_interp.interp()

44
45
46
47
48
49
50
51
52
53
54
    ypn=-1.0;
    Spline_interp z = new Spline_interp(x,y,yp1,ypn);
    Ran myran = new Ran(17);
    for(i=0;i<N;i++) {
      xx[i]=pi*myran.doub();
      yy[i]=z.interp(xx[i]);
      zz[i]=sin(xx[i]);
    }
    System.out.printf("     Spline_interp: Max. actual error:    %f\n", maxel(vecsub(zz,yy)));
    sbeps=1.e-5*maxel(y);
    localflag = maxel(vecsub(zz,yy)) > sbeps;
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. 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.