Examples of GillIntegrator

org.apache.commons.math.ode.GillIntegrator
This class implements the Gill fourth order Runge-Kutta integrator for Ordinary Differential Equations .
This method is an explicit Runge-Kutta method, its Butcher-array is the following one :
```
 0  |    0        0       0      0 1/2 |   1/2       0       0      0 1/2 | (q-1)/2  (2-q)/2    0      0 1  |    0       -q/2  (2+q)/2   0 |------------------------------- |   1/6    (2-q)/6 (2+q)/6  1/6 
```
where q = sqrt(2)
@see EulerIntegrator @see ClassicalRungeKuttaIntegrator @see MidpointIntegrator @see ThreeEighthesIntegrator @version $Revision: 620312 $ $Date: 2008-02-10 12:28:59 -0700 (Sun, 10 Feb 2008) $ @since 1.2
org.apache.commons.math.ode.nonstiff.GillIntegrator
This class implements the Gill fourth order Runge-Kutta integrator for Ordinary Differential Equations .
This method is an explicit Runge-Kutta method, its Butcher-array is the following one :
```
 0  |    0        0       0      0 1/2 |   1/2       0       0      0 1/2 | (q-1)/2  (2-q)/2    0      0 1  |    0       -q/2  (2+q)/2   0 |------------------------------- |   1/6    (2-q)/6 (2+q)/6  1/6 
```
where q = sqrt(2)
@see EulerIntegrator @see ClassicalRungeKuttaIntegrator @see MidpointIntegrator @see ThreeEighthesIntegrator @version $Revision: 990655 $ $Date: 2010-08-29 23:49:40 +0200 (dim. 29 août 2010) $ @since 1.2

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

    throws DerivativeException, IntegratorException {


    final TestProblem3 pb  = new TestProblem3(0.9);
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.0003;


    FirstOrderIntegrator integ = new GillIntegrator(step);
    integ.addStepHandler(new KeplerStepHandler(pb));
    integ.integrate(pb,
                    pb.getInitialTime(), pb.getInitialState(),
                    pb.getFinalTime(), new double[pb.getDimension()]);
  }

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

  }


  public void testUnstableDerivative()
  throws DerivativeException, IntegratorException {
    final StepProblem stepProblem = new StepProblem(0.0, 1.0, 2.0);
    FirstOrderIntegrator integ = new GillIntegrator(0.3);
    integ.addEventHandler(stepProblem, 1.0, 1.0e-12, 1000);
    double[] y = { Double.NaN };
    integ.integrate(stepProblem, 0.0, new double[] { 0.0 }, 10.0, y);
    assertEquals(8.0, y[0], 1.0e-12);
  }

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

  }


  public void testStepSize()
    throws DerivativeException, IntegratorException {
      final double step = 1.23456;
      FirstOrderIntegrator integ = new GillIntegrator(step);
      integ.addStepHandler(new StepHandler() {
          public void handleStep(StepInterpolator interpolator, boolean isLast) {
              if (! isLast) {
                  assertEquals(step,
                               interpolator.getCurrentTime() - interpolator.getPreviousTime(),
                               1.0e-12);
              }
          }
          public boolean requiresDenseOutput() {
              return false;
          }
          public void reset() {
          }
      });
      integ.integrate(new FirstOrderDifferentialEquations() {
          private static final long serialVersionUID = 0L;
          public void computeDerivatives(double t, double[] y, double[] dot) {
              dot[0] = 1.0;
          }
          public int getDimension() {

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

  @Test
  public void testDerivativesConsistency()
  throws DerivativeException, IntegratorException {
    TestProblem3 pb = new TestProblem3();
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.001;
    GillIntegrator integ = new GillIntegrator(step);
    StepInterpolatorTestUtils.checkDerivativesConsistency(integ, pb, 1.0e-10);
  }

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

    throws DerivativeException, IntegratorException,
           IOException, ClassNotFoundException {


    TestProblem3 pb = new TestProblem3(0.9);
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.0003;
    GillIntegrator integ = new GillIntegrator(step);
    integ.addStepHandler(new ContinuousOutputModel());
    integ.integrate(pb,
                    pb.getInitialTime(), pb.getInitialState(),
                    pb.getFinalTime(), new double[pb.getDimension()]);


    ByteArrayOutputStream bos = new ByteArrayOutputStream();
    ObjectOutputStream    oos = new ObjectOutputStream(bos);
    for (StepHandler handler : integ.getStepHandlers()) {
        oos.writeObject(handler);
    }


    assertTrue(bos.size () > 753000);
    assertTrue(bos.size () < 754000);

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

  @Test
  public void testDerivativesConsistency()
  throws DerivativeException, IntegratorException {
    TestProblem3 pb = new TestProblem3();
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.001;
    GillIntegrator integ = new GillIntegrator(step);
    StepInterpolatorTestUtils.checkDerivativesConsistency(integ, pb, 1.0e-10);
  }

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

    throws DerivativeException, IntegratorException,
           IOException, ClassNotFoundException {


    TestProblem3 pb = new TestProblem3(0.9);
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.0003;
    GillIntegrator integ = new GillIntegrator(step);
    integ.addStepHandler(new ContinuousOutputModel());
    integ.integrate(pb,
                    pb.getInitialTime(), pb.getInitialState(),
                    pb.getFinalTime(), new double[pb.getDimension()]);


    ByteArrayOutputStream bos = new ByteArrayOutputStream();
    ObjectOutputStream    oos = new ObjectOutputStream(bos);
    for (StepHandler handler : integ.getStepHandlers()) {
        oos.writeObject(handler);
    }


    assertTrue(bos.size () > 700000);
    assertTrue(bos.size () < 701000);

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

  }


  public void testDimensionCheck() {
    try  {
      TestProblem1 pb = new TestProblem1();
      new GillIntegrator(0.01).integrate(pb,
                                         0.0, new double[pb.getDimension()+10],
                                         1.0, new double[pb.getDimension()+10]);
        fail("an exception should have been thrown");
    } catch(DerivativeException de) {
      fail("wrong exception caught");

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator


        TestProblemAbstract pb = problems[k].copy();
        double step = (pb.getFinalTime() - pb.getInitialTime())
          * Math.pow(2.0, -i);


        FirstOrderIntegrator integ = new GillIntegrator(step);
        TestProblemHandler handler = new TestProblemHandler(pb, integ);
        integ.addStepHandler(handler);
        EventHandler[] functions = pb.getEventsHandlers();
        for (int l = 0; l < functions.length; ++l) {
          integ.addEventHandler(functions[l],
                                     Double.POSITIVE_INFINITY, 1.0e-6 * step, 1000);
        }
        double stopTime = integ.integrate(pb, pb.getInitialTime(), pb.getInitialState(),
                                          pb.getFinalTime(), new double[pb.getDimension()]);
        if (functions.length == 0) {
            assertEquals(pb.getFinalTime(), stopTime, 1.0e-10);
        }

View Full Code Here

Examples of org.apache.commons.math.ode.nonstiff.GillIntegrator

    throws DerivativeException, IntegratorException {


    TestProblem1 pb = new TestProblem1();
    double step = (pb.getFinalTime() - pb.getInitialTime()) * 0.001;


    FirstOrderIntegrator integ = new GillIntegrator(step);
    TestProblemHandler handler = new TestProblemHandler(pb, integ);
    integ.addStepHandler(handler);
    integ.integrate(pb, pb.getInitialTime(), pb.getInitialState(),
                    pb.getFinalTime(), new double[pb.getDimension()]);


    assertTrue(handler.getLastError() < 2.0e-13);
    assertTrue(handler.getMaximalValueError() < 4.0e-12);
    assertEquals(0, handler.getMaximalTimeError(), 1.0e-12);
    assertEquals("Gill", integ.getName());


  }

View Full Code Here

0 1 2 3

TOP

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.