public class oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator extends oadd.org.apache.commons.math.ode.nonstiff.EmbeddedRungeKuttaIntegrator
{
private static final java.lang.String METHOD_NAME;
private static final double[] STATIC_C;
private static final double[][] STATIC_A;
private static final double[] STATIC_B;
private static final double E1;
private static final double E3;
private static final double E4;
private static final double E5;
private static final double E6;
private static final double E7;
public void <init>(double, double, double, double)
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator v;
double[] v, v;
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator v;
double[][] v;
double v, v, v, v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_C>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[][] STATIC_A>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_B>;
v = new oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator;
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator: void <init>()>();
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.EmbeddedRungeKuttaIntegrator: void <init>(java.lang.String,boolean,double[],double[][],double[],oadd.org.apache.commons.math.ode.nonstiff.RungeKuttaStepInterpolator,double,double,double,double)>("Dormand-Prince 5(4)", 1, v, v, v, v, v, v, v, v);
return;
}
public void <init>(double, double, double[], double[])
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator v;
double[] v, v, v, v;
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator v;
double[][] v;
double v, v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double[];
v := @parameter: double[];
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_C>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[][] STATIC_A>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_B>;
v = new oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator;
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54StepInterpolator: void <init>()>();
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.EmbeddedRungeKuttaIntegrator: void <init>(java.lang.String,boolean,double[],double[][],double[],oadd.org.apache.commons.math.ode.nonstiff.RungeKuttaStepInterpolator,double,double,double[],double[])>("Dormand-Prince 5(4)", 1, v, v, v, v, v, v, v, v);
return;
}
public int getOrder()
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator;
return 5;
}
protected double estimateError(double[][], double[], double[], double)
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator v;
double[][] v;
double[] v, v, v, v, v, v, v, v, v, v, v;
int v, v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator;
v := @parameter: double[][];
v := @parameter: double[];
v := @parameter: double[];
v := @parameter: double;
v = 0.0;
v = 0;
label:
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: int mainSetDimension>;
if v >= v goto label;
v = v[0];
v = v[v];
v = 0.0012326388888888888 * v;
v = v[2];
v = v[v];
v = -0.0042527702905061394 * v;
v = v + v;
v = v[3];
v = v[v];
v = 0.03697916666666667 * v;
v = v + v;
v = v[4];
v = v[v];
v = -0.05086379716981132 * v;
v = v + v;
v = v[5];
v = v[v];
v = 0.0419047619047619 * v;
v = v + v;
v = v[6];
v = v[v];
v = -0.025 * v;
v = v + v;
v = v[v];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double abs(double)>(v);
v = v[v];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double abs(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double max(double,double)>(v, v);
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] vecAbsoluteTolerance>;
if v != null goto label;
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double scalAbsoluteTolerance>;
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double scalRelativeTolerance>;
v = v * v;
v = v + v;
goto label;
label:
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] vecAbsoluteTolerance>;
v = v[v];
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] vecRelativeTolerance>;
v = v[v];
v = v * v;
v = v + v;
label:
v = v * v;
v = v / v;
v = v * v;
v = v + v;
v = v + 1;
goto label;
label:
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: int mainSetDimension>;
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(v);
return v;
}
static void <clinit>()
{
double[] v, v, v, v, v, v, v, v;
double[][] v;
v = newarray (double)[6];
v[0] = 0.2;
v[1] = 0.3;
v[2] = 0.8;
v[3] = 0.8888888888888888;
v[4] = 1.0;
v[5] = 1.0;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_C> = v;
v = newarray (double[])[6];
v = newarray (double)[1];
v[0] = 0.2;
v[0] = v;
v = newarray (double)[2];
v[0] = 0.075;
v[1] = 0.225;
v[1] = v;
v = newarray (double)[3];
v[0] = 0.9777777777777777;
v[1] = -3.7333333333333334;
v[2] = 3.5555555555555554;
v[2] = v;
v = newarray (double)[4];
v[0] = 2.9525986892242035;
v[1] = -11.595793324188385;
v[2] = 9.822892851699436;
v[3] = -0.2908093278463649;
v[3] = v;
v = newarray (double)[5];
v[0] = 2.8462752525252526;
v[1] = -10.757575757575758;
v[2] = 8.906422717743473;
v[3] = 0.2784090909090909;
v[4] = -0.2735313036020583;
v[4] = v;
v = newarray (double)[6];
v[0] = 0.09114583333333333;
v[1] = 0.0;
v[2] = 0.44923629829290207;
v[3] = 0.6510416666666666;
v[4] = -0.322376179245283;
v[5] = 0.13095238095238096;
v[5] = v;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[][] STATIC_A> = v;
v = newarray (double)[7];
v[0] = 0.09114583333333333;
v[1] = 0.0;
v[2] = 0.44923629829290207;
v[3] = 0.6510416666666666;
v[4] = -0.322376179245283;
v[5] = 0.13095238095238096;
v[6] = 0.0;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince54Integrator: double[] STATIC_B> = v;
return;
}
}