public class oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator 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_01;
private static final double E1_06;
private static final double E1_07;
private static final double E1_08;
private static final double E1_09;
private static final double E1_10;
private static final double E1_11;
private static final double E1_12;
private static final double E2_01;
private static final double E2_06;
private static final double E2_07;
private static final double E2_08;
private static final double E2_09;
private static final double E2_10;
private static final double E2_11;
private static final double E2_12;
public void <init>(double, double, double, double)
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator v;
double[] v, v;
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator v;
double[][] v;
double v, v, v, v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_C>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[][] STATIC_A>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_B>;
v = new oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator;
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator: 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 8 (5, 3)", 1, v, v, v, v, v, v, v, v);
return;
}
public void <init>(double, double, double[], double[])
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator v;
double[] v, v, v, v;
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator v;
double[][] v;
double v, v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double[];
v := @parameter: double[];
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_C>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[][] STATIC_A>;
v = <oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_B>;
v = new oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator;
specialinvoke v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853StepInterpolator: 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 8 (5, 3)", 1, v, v, v, v, v, v, v, v);
return;
}
public int getOrder()
{
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator v;
v := @this: oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator;
return 8;
}
protected double estimateError(double[][], double[], double[], double)
{
byte v;
double[][] v;
double[] v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator 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, 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.DormandPrince853Integrator;
v := @parameter: double[][];
v := @parameter: double[];
v := @parameter: double[];
v := @parameter: double;
v = 0.0;
v = 0.0;
v = 0;
label:
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: int mainSetDimension>;
if v >= v goto label;
v = v[0];
v = v[v];
v = 0.01312004499419488 * v;
v = v[5];
v = v[v];
v = -1.2251564463762044 * v;
v = v + v;
v = v[6];
v = v[v];
v = -0.4957589496572502 * v;
v = v + v;
v = v[7];
v = v[v];
v = 1.6643771824549864 * v;
v = v + v;
v = v[8];
v = v[v];
v = -0.35032884874997366 * v;
v = v + v;
v = v[9];
v = v[v];
v = 0.3341791187130175 * v;
v = v + v;
v = v[10];
v = v[v];
v = 0.08192320648511571 * v;
v = v + v;
v = v[11];
v = v[v];
v = -0.022355307863886294 * v;
v = v + v;
v = v[0];
v = v[v];
v = -0.18980075407240762 * v;
v = v[5];
v = v[v];
v = 4.450312892752409 * v;
v = v + v;
v = v[6];
v = v[v];
v = 1.8915178993145003 * v;
v = v + v;
v = v[7];
v = v[v];
v = -5.801203960010585 * v;
v = v + v;
v = v[8];
v = v[v];
v = -0.42268232132379197 * v;
v = v + v;
v = v[9];
v = v[v];
v = -0.1521609496625161 * v;
v = v + v;
v = v[10];
v = v[v];
v = 0.20136540080403034 * v;
v = v + v;
v = v[11];
v = v[v];
v = 0.022651792198360825 * 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.DormandPrince853Integrator: double[] vecAbsoluteTolerance>;
if v != null goto label;
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double scalAbsoluteTolerance>;
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double scalRelativeTolerance>;
v = v * v;
v = v + v;
goto label;
label:
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] vecAbsoluteTolerance>;
v = v[v];
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: 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 * v;
v = v + v;
v = v + 1;
goto label;
label:
v = 0.01 * v;
v = v + v;
v = v cmpg 0.0;
if v > 0 goto label;
v = 1.0;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double abs(double)>(v);
v = v * v;
v = v.<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: int mainSetDimension>;
v = v * v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(v);
v = v / v;
return v;
}
static void <clinit>()
{
double[][] v;
double[] v, v, v, v, v, v, v, v, v, v, v, 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, 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, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v = newarray (double)[12];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 2.0 * v;
v = 12.0 - v;
v = v / 135.0;
v[0] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 45.0;
v[1] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 30.0;
v[2] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 + v;
v = v / 30.0;
v[3] = v;
v[4] = 0.3333333333333333;
v[5] = 0.25;
v[6] = 0.3076923076923077;
v[7] = 0.6512820512820513;
v[8] = 0.6;
v[9] = 0.8571428571428571;
v[10] = 1.0;
v[11] = 1.0;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_C> = v;
v = newarray (double[])[12];
v = newarray (double)[1];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 2.0 * v;
v = 12.0 - v;
v = v / 135.0;
v[0] = v;
v[0] = v;
v = newarray (double)[2];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 180.0;
v[0] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 60.0;
v[1] = v;
v[1] = v;
v = newarray (double)[3];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 120.0;
v[0] = v;
v[1] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 6.0 - v;
v = v / 40.0;
v[2] = v;
v[2] = v;
v = newarray (double)[4];
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 107.0 * v;
v = 462.0 + v;
v = v / 3000.0;
v[0] = v;
v[1] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 197.0 * v;
v = -402.0 - v;
v = v / 1000.0;
v[2] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 73.0 * v;
v = 168.0 + v;
v = v / 375.0;
v[3] = v;
v[3] = v;
v = newarray (double)[5];
v[0] = 0.037037037037037035;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 16.0 + v;
v = v / 108.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 16.0 - v;
v = v / 108.0;
v[4] = v;
v[4] = v;
v = newarray (double)[6];
v[0] = 0.037109375;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 23.0 * v;
v = 118.0 + v;
v = v / 1024.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 23.0 * v;
v = 118.0 - v;
v = v / 1024.0;
v[4] = v;
v[5] = -0.017578125;
v[5] = v;
v = newarray (double)[7];
v[0] = 0.03709200011850479;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 4784.0 * v;
v = 51544.0 + v;
v = v / 371293.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 4784.0 * v;
v = 51544.0 - v;
v = v / 371293.0;
v[4] = v;
v[5] = -0.015319437748624402;
v[6] = 0.008273789163814023;
v[6] = v;
v = newarray (double)[8];
v[0] = 0.6241109587160757;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 3.18801444819E11 * v;
v = -1.324889724104E12 - v;
v = v / 6.265569375E11;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 3.18801444819E11 * v;
v = -1.324889724104E12 + v;
v = v / 6.265569375E11;
v[4] = v;
v[5] = 27.59209969944671;
v[6] = 20.154067550477894;
v[7] = -43.48988418106996;
v[7] = v;
v = newarray (double)[9];
v[0] = 0.47766253643826434;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 1137963.0 * v;
v = -4521408.0 - v;
v = v / 2937500.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 1137963.0 * v;
v = -4521408.0 + v;
v = v / 2937500.0;
v[4] = v;
v[5] = 21.230051448181193;
v[6] = 15.279233632882423;
v[7] = -33.28821096898486;
v[8] = -0.020331201708508627;
v[8] = v;
v = newarray (double)[10];
v[0] = -0.9371424300859873;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 94326.0 * v;
v = 354216.0 + v;
v = v / 112847.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 94326.0 * v;
v = 354216.0 - v;
v = v / 112847.0;
v[4] = v;
v[5] = -8.149787010746927;
v[6] = -18.52006565999696;
v[7] = 22.739487099350505;
v[8] = 2.4936055526796523;
v[9] = -3.0467644718982196;
v[9] = v;
v = newarray (double)[11];
v[0] = 2.273310147516538;
v[1] = 0.0;
v[2] = 0.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 960905.0 * v;
v = -3457480.0 - v;
v = v / 551636.0;
v[3] = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double sqrt(double)>(6.0);
v = 960905.0 * v;
v = -3457480.0 + v;
v = v / 551636.0;
v[4] = v;
v[5] = -17.9589318631188;
v[6] = 27.94888452941996;
v[7] = -2.8589982771350235;
v[8] = -8.87285693353063;
v[9] = 12.360567175794303;
v[10] = 0.6433927460157636;
v[10] = v;
v = newarray (double)[12];
v[0] = 0.054293734116568765;
v[1] = 0.0;
v[2] = 0.0;
v[3] = 0.0;
v[4] = 0.0;
v[5] = 4.450312892752409;
v[6] = 1.8915178993145003;
v[7] = -5.801203960010585;
v[8] = 0.3111643669578199;
v[9] = -0.1521609496625161;
v[10] = 0.20136540080403034;
v[11] = 0.04471061572777259;
v[11] = v;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[][] STATIC_A> = v;
v = newarray (double)[13];
v[0] = 0.054293734116568765;
v[1] = 0.0;
v[2] = 0.0;
v[3] = 0.0;
v[4] = 0.0;
v[5] = 4.450312892752409;
v[6] = 1.8915178993145003;
v[7] = -5.801203960010585;
v[8] = 0.3111643669578199;
v[9] = -0.1521609496625161;
v[10] = 0.20136540080403034;
v[11] = 0.04471061572777259;
v[12] = 0.0;
<oadd.org.apache.commons.math.ode.nonstiff.DormandPrince853Integrator: double[] STATIC_B> = v;
return;
}
}