Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
test for output for derivative function of ThreeBodyExample
  • Loading branch information
jic13003 committed Dec 5, 2015
1 parent 8d7f84a commit c920ba0
Show file tree
Hide file tree
Showing 2 changed files with 32 additions and 19 deletions.
37 changes: 19 additions & 18 deletions src/jat/examples/ThreeBodyExample/ThreeBodyAPL.java
Expand Up @@ -23,24 +23,25 @@ public class ThreeBodyAPL extends ThreeBody implements FirstOrderDifferentialEqu
return 6;
}
public double[] randv() {
double[] randv = new double[17];
randv[1] = 1.0;
randv[2] = 1.0;
randv[3] = 1.0;
randv[4] = 1.0;
randv[5] = 1.0;
randv[6] = 1.0;
randv[7] = 1.0;
randv[8] = 1.0;
randv[9] = 1.0;
randv[10] = 1.0;
randv[11] = 1.0;
randv[12] = 1.0;
randv[13] = 1.0;
randv[14] = 1.0;
randv[15] = 1.0;
randv[16] = 1.0;
randv[17] = 1.0;
double[] randv = new double[18];
randv[0] = 2.0;
randv[1] = 6.0;
randv[2] = 7.0;
randv[3] = 4.0;
randv[4] = 8.0;
randv[5] = 7.0;
randv[6] = 5.0;
randv[7] = 8.0;
randv[8] = 3.0;
randv[9] = 6.0;
randv[10] = 8.0;
randv[11] = 9.0;
randv[12] = 0.0;
randv[13] = 7.0;
randv[14] = 4.0;
randv[15] = 5.0;
randv[16] = 7.0;
randv[17] = 8.0;
return randv;
}
}
14 changes: 13 additions & 1 deletion src/jat/examples/ThreeBodyExample/ThreeBodyExample.java
Expand Up @@ -46,7 +46,19 @@ public class ThreeBodyExample {

// set the initial time to zero
double t0 = 0.0;
//test

// propagate the orbit
//FirstOrderIntegrator dp853 = new DormandPrince853Integrator(1.0e-8, 100.0, 1.0e-10, 1.0e-10);
//dp853.addStepHandler(sat.stepHandler);
// double[] y = new double[] { 7000.0, 0, 0, .0, 8, 0 }; // initial
// state
//test loop

for (int i = 0; i< derivs.length; i++)
{
System.out.println(derivs[i]);
}

}

}

0 comments on commit c920ba0

Please sign in to comment.