Package com.nr.ran

Examples of com.nr.ran.Ran.doub()

40
41
42
43
44
45
46
47
48
49
50
      }
    }
    Bilin_interp z = new Bilin_interp(x1,x2,y);
    Ran myran = new Ran(17);
    for (i=0;i<N;i++) {
      xx1=myran.doub();
      xx2=myran.doub();
      yy[i]=z.interp(xx1,xx2);  // interpolated values
      zz[i]=xx1+xx2;        // Actual values
    }
    sbeps=1.e-15;
View Full Code Here
41
42
43
44
45
46
47
48
49
50
51
    }
    Bilin_interp z = new Bilin_interp(x1,x2,y);
    Ran myran = new Ran(17);
    for (i=0;i<N;i++) {
      xx1=myran.doub();
      xx2=myran.doub();
      yy[i]=z.interp(xx1,xx2);  // interpolated values
      zz[i]=xx1+xx2;        // Actual values
    }
    sbeps=1.e-15;
    System.out.printf("     Bilin_interp: Max. actual error:    %f\n", maxel(vecsub(zz,yy)));
View Full Code Here
65
66
67
68
69
70
71
72
73
74
75
    // Round trip test
    Ran myran = new Ran(17);
    sbeps=1.e-8;
    for (i=0;i<M;i++) {
      c=10.0*myran.doub();
      d=10.0*myran.doub();
      u=myran.doub();
//      System.out.printf(i << " %f\n", u;
      r=Beta.betai(c,d,u);
//      System.out.println(" %f\n", r;
View Full Code Here
66
67
68
69
70
71
72
73
74
75
76
    // Round trip test
    Ran myran = new Ran(17);
    sbeps=1.e-8;
    for (i=0;i<M;i++) {
      c=10.0*myran.doub();
      d=10.0*myran.doub();
      u=myran.doub();
//      System.out.printf(i << " %f\n", u;
      r=Beta.betai(c,d,u);
//      System.out.println(" %f\n", r;
      uu=Beta.invbetai(r,c,d);
View Full Code Here
67
68
69
70
71
72
73
74
75
76
77
    Ran myran = new Ran(17);
    sbeps=1.e-8;
    for (i=0;i<M;i++) {
      c=10.0*myran.doub();
      d=10.0*myran.doub();
      u=myran.doub();
//      System.out.printf(i << " %f\n", u;
      r=Beta.betai(c,d,u);
//      System.out.println(" %f\n", r;
      uu=Beta.invbetai(r,c,d);
//      System.out.println(" %f\n", uu);
View Full Code Here
50
51
52
53
54
55
56
57
58
59
60
    Ran myran = new Ran(17);
   
    // Test rj(x,x,x,x) = 1/x^(3/2)
    for (i=0;i<N;i++) {
      x=myran.doub();
      f1[i]=rj(x,x,x,x);
      f2[i]=1.0/pow(x,3.0/2.0);
    }

    System.out.printf("rj: Maximum discrepancy = %f\n", maxel(vecsub(f1,f2)));
View Full Code Here
66
67
68
69
70
71
72
73
74
75
76
     
    }

    // Symmetry test
    for (i=0;i<N;i++) {
      x=10.0*myran.doub();
      y=10.0*myran.doub();
      z=10.0*myran.doub();
      p=10.0*myran.doub();

      f1[i]=rj(x,y,z,p);
View Full Code Here
67
68
69
70
71
72
73
74
75
76
77
    }

    // Symmetry test
    for (i=0;i<N;i++) {
      x=10.0*myran.doub();
      y=10.0*myran.doub();
      z=10.0*myran.doub();
      p=10.0*myran.doub();

      f1[i]=rj(x,y,z,p);
      f2[i]=rj(y,x,z,p);
View Full Code Here
68
69
70
71
72
73
74
75
76
77
78
    // Symmetry test
    for (i=0;i<N;i++) {
      x=10.0*myran.doub();
      y=10.0*myran.doub();
      z=10.0*myran.doub();
      p=10.0*myran.doub();

      f1[i]=rj(x,y,z,p);
      f2[i]=rj(y,x,z,p);
      f3[i]=rj(x,z,y,p);
View Full Code Here
69
70
71
72
73
74
75
76
77
78
79
    // Symmetry test
    for (i=0;i<N;i++) {
      x=10.0*myran.doub();
      y=10.0*myran.doub();
      z=10.0*myran.doub();
      p=10.0*myran.doub();

      f1[i]=rj(x,y,z,p);
      f2[i]=rj(y,x,z,p);
      f3[i]=rj(x,z,y,p);
    }
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.