Package com.nr.ode

Examples of com.nr.ode.Odeint

40
41
42
43
44
45
46
47
48
49
50
51
    y[0]=1.0; y[1]=1.0; y[2]=0.0;
    Output out=new Output(20);
    rhs_StepperRoss d = new rhs_StepperRoss();
    StepperRoss s = new StepperRoss();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
//      System.out.printf(setprecision(12) << yout[i]);
    }
View Full Code Here
39
40
41
42
43
44
45
46
47
48
49
50
    y[0]=1.0; y[1]=1.0; y[2]=0.0;
    Output out=new Output(20);
    rhs_StepperSie d = new rhs_StepperSie();
    StepperSie s= new StepperSie();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
//      System.out.printf(setprecision(12) << yout[i]);
    }
View Full Code Here
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
   
    Output out=new Output(20);
    x1=0.0;x2=1.0;
    rhs_StepperStoerm d = new rhs_StepperStoerm();
    StepperStoerm s1= new StepperStoerm();
    Odeint ode1 = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s1);
    ode1.integrate();
    yexp[0]=bess.airy_ai(1.0);
    yexp[1]=bess.airy_bi(1.0);

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
      System.out.printf("%f  %f\n", yout[i],yexp[i]);
    }

    sbeps = 1.e-8;
    System.out.println(maxel(vecsub(yout,yexp)));
    localflag = maxel(vecsub(yout,yexp)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
      fail("*** StepperStoerm: Inaccurate integration #1");
     
    }

    y[0]=1.0/pow(3.0,2.0/3.0)/exp(gammln(2.0/3.0));     // Ai(x)
    y[nvar]=-1.0/pow(3.0,1.0/3.0)/exp(gammln(1.0/3.0));
    dydx[0]=0.0;
    y[1]=1.0/pow(3.0,1.0/6.0)/exp(gammln(2.0/3.0));     // Bi(x)
    y[nvar+1]=pow(3.0,1.0/6.0)/exp(gammln(1.0/3.0));
    dydx[1]=0.0;

    x1=0.0;x2=-1.0;
    StepperStoerm s2 = new StepperStoerm();
    Odeint ode2 = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s2);
    ode2.integrate();
    yexp[0]=bess.airy_ai(-1.0);
    yexp[1]=bess.airy_bi(-1.0);

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
View Full Code Here
44
45
46
47
48
49
50
51
52
53
54
55
      yexp[i]=bess.jn(i,x2);
    }
    Output out = new Output(20);
    rhs_StepperDopr5 d = new rhs_StepperDopr5();
    StepperDopr5 s = new StepperDopr5();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
      System.out.printf("%f  %f\n", yout[i],yexp[i]);
    }
View Full Code Here
67
68
69
70
71
72
73
74
75
76
    yy[2]=y[1].re();
    yy[3]=y[1].im();
    Hypderiv d = new Hypderiv(a,b,c,z0,dz);
    Output out = new Output();
    StepperBS s = new StepperBS();
    Odeint ode =new Odeint(yy,0.0,1.0,atol,rtol,0.1,0.0,out,d,s);
    ode.integrate();
    y[0]=new Complex(yy[0],yy[1]);
    return y[0];
  }
View Full Code Here
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
      v2[i] = v[n2+i];
    h1=(x2-x1)/100.0;
    y=buildVector(load1(x1,v));
    Output out = new Output();
    StepperDopr853 s1 = new StepperDopr853();
    Odeint integ1 = new Odeint(y,x1,xf,atol,rtol,h1,hmin,out,d,s1);
    integ1.integrate();
    f1=buildVector(score(xf,y));
    y=buildVector(load2(x2,v2));
   
    StepperDopr853 s2 = new StepperDopr853();
    Odeint integ2 = new Odeint(y,x2,xf,atol,rtol,h1,hmin,out,d,s2);
    integ2.integrate();
    f2=buildVector(score(xf,y));
    for (int i=0;i<nvar;i++) f1[i] -= f2[i];
    return f1;
  }
View Full Code Here
39
40
41
42
43
44
45
46
47
  public double[] funk (final double[] v) {
    h1=(x2-x1)/100.0;
    y=load(x1,v);
    Output out = new Output();
    StepperDopr853 s = new StepperDopr853();
    Odeint integ = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d, s);
    integ.integrate();
    return score(x2,y);
  }
View Full Code Here
44
45
46
47
48
49
50
51
52
53
54
55
      yexp[i]=bess.jn(i,x2);
    }
    Output out=new Output(20);
    rhs_StepperDopr853 d = new rhs_StepperDopr853();
    StepperDopr853 s = new StepperDopr853();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
      System.out.printf("%f  %f\n", yout[i],yexp[i]);
    }
View Full Code Here
TOP

Related Classes of com.nr.ode.Odeint

  • com.nr.bvp.Shoot
  • com.nr.bvp.Shootf
  • com.nr.sf.Hypergeo
  • com.nr.test.test_chapter17.Test_StepperDopr5
  • com.nr.test.test_chapter17.Test_StepperDopr853
  • com.nr.test.test_chapter17.Test_StepperRoss
  • com.nr.test.test_chapter17.Test_StepperSie
  • com.nr.test.test_chapter17.Test_StepperStoerm
    • 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.