Examples of jn()


Examples of com.nr.sf.Bessjy.jn()

          zz2[j]=bess.jnu(nnu,xx);
        } else if (i == 1) {
          zz1[j]=bjy.j1(xx);
          zz2[j]=bess.jnu(nnu,xx);
        } else {
          zz1[j]=bjy.jn(i,xx);
          zz2[j]=bess.jnu(nnu,xx);
        }
      }
      localflag = maxel(vecsub(zz1,zz2)) > sbeps;
      globalflag = globalflag || localflag;
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

    // Test Bessjy (jn)
    System.out.println("Testing Bessjy (jn)");

    sbeps=5.e-14;
    for (i=0;i<N;i++) zz5[i]=bess.jn(n5[i],x5[i]);
    for (i=0;i<N;i++) {
//      System.out.printf(fabs((zz5[i]-yy5[i])/yy5[i]));
      localflag = localflag || abs((zz5[i]-yy5[i])/yy5[i]) > sbeps;
    }
    globalflag = globalflag || localflag;
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

      y=x[0]*(2.0*myran.doub()-1.0);
      Weights.weights(y,x,c);
      // Calculate function
      f=f1=f2=0;
      for (j=0;j<N;j++) {
        f += c[j][0]*bess.jn(5,x[j]);
        f1 += c[j][1]*bess.jn(5,x[j]);
        f2 += c[j][2]*bess.jn(5,x[j]);
      }
//      System.out.printf(y);
//      System.out.printf(abs(bess.jn(5,y)-f) << " ";
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

      Weights.weights(y,x,c);
      // Calculate function
      f=f1=f2=0;
      for (j=0;j<N;j++) {
        f += c[j][0]*bess.jn(5,x[j]);
        f1 += c[j][1]*bess.jn(5,x[j]);
        f2 += c[j][2]*bess.jn(5,x[j]);
      }
//      System.out.printf(y);
//      System.out.printf(abs(bess.jn(5,y)-f) << " ";
//      System.out.printf(abs(0.5*(bess.jn(4,y)-bess.jn(6,y))-f1) << " ";
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

      // Calculate function
      f=f1=f2=0;
      for (j=0;j<N;j++) {
        f += c[j][0]*bess.jn(5,x[j]);
        f1 += c[j][1]*bess.jn(5,x[j]);
        f2 += c[j][2]*bess.jn(5,x[j]);
      }
//      System.out.printf(y);
//      System.out.printf(abs(bess.jn(5,y)-f) << " ";
//      System.out.printf(abs(0.5*(bess.jn(4,y)-bess.jn(6,y))-f1) << " ";
//      System.out.printf(abs(0.25*(bess.jn(3,y)-2.0*bess.jn(5,y)+bess.jn(7,y))-f2));
 
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

//      System.out.printf(abs(bess.jn(5,y)-f) << " ";
//      System.out.printf(abs(0.5*(bess.jn(4,y)-bess.jn(6,y))-f1) << " ";
//      System.out.printf(abs(0.25*(bess.jn(3,y)-2.0*bess.jn(5,y)+bess.jn(7,y))-f2));

      sbeps=1.e-3;
      localflag = abs(bess.jn(5,y)-f) > sbeps;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** weights: Evaluation of function was inaccurate");
       
      }
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

        fail("*** weights: Evaluation of function was inaccurate");
       
      }

      sbeps=1.e-3;
      localflag = abs(0.5*(bess.jn(4,y)-bess.jn(6,y))-f1) > sbeps;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** weights: Evaluation of first derivative was inaccurate");
       
      }
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

        fail("*** weights: Evaluation of function was inaccurate");
       
      }

      sbeps=1.e-3;
      localflag = abs(0.5*(bess.jn(4,y)-bess.jn(6,y))-f1) > sbeps;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** weights: Evaluation of first derivative was inaccurate");
       
      }
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

        fail("*** weights: Evaluation of first derivative was inaccurate");
       
      }

      sbeps=1.e-3;
      localflag = abs(0.25*(bess.jn(3,y)-2.0*bess.jn(5,y)+bess.jn(7,y))-f2) > sbeps;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** weights: Evaluation of second derivative was inaccurate");
       
      }
View Full Code Here

Examples of com.nr.sf.Bessjy.jn()

        fail("*** weights: Evaluation of first derivative was inaccurate");
       
      }

      sbeps=1.e-3;
      localflag = abs(0.25*(bess.jn(3,y)-2.0*bess.jn(5,y)+bess.jn(7,y))-f2) > sbeps;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** weights: Evaluation of second derivative was inaccurate");
       
      }
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.