major optimization updates

[debian/openrocket] / src / net / sf / openrocket / simulation / RK4SimulationStepper.java

diff --git a/src/net/sf/openrocket/simulation/RK4SimulationStepper.java b/src/net/sf/openrocket/simulation/RK4SimulationStepper.java
index 8ef74c8581c40791ffab558f99231ad74f8a1ed7..8d5e27dbf94264479d1d27bea78074e420a57e0a 100644 (file)
--- a/src/net/sf/openrocket/simulation/RK4SimulationStepper.java
+++ b/src/net/sf/openrocket/simulation/RK4SimulationStepper.java
@@ -21,6 +21,9 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
        
        private static final LogHelper log = Application.getLogger();
        
+       /** Random value with which to XOR the random seed value */
+       private static final int SEED_RANDOMIZATION = 0x23E3A01F;
+       
 
        /**
         * A recommended reasonably accurate time step.
@@ -50,8 +53,7 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
        private static final double MIN_TIME_STEP = 0.001;
        
 
-
-       private final Random random = new Random();
+       private Random random;
        
        
 
@@ -59,8 +61,6 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
        @Override
        public RK4SimulationStatus initialize(SimulationStatus original) {
                
-               log.info("Performing RK4SimulationStepper initialization");
-               
                RK4SimulationStatus status = new RK4SimulationStatus();
                
                status.copyFrom(original);
@@ -72,6 +72,8 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                                Math.cos(sim.getLaunchRodAngle())
                                ));
                
+               this.random = new Random(original.getSimulationConditions().getRandomSeed() ^ SEED_RANDOMIZATION);
+               
                return status;
        }
        
@@ -155,11 +157,11 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                
                double minTimeStep = status.getSimulationConditions().getTimeStep() / 20;
                if (store.timestep < minTimeStep) {
-                       log.debug("Too small time step " + store.timestep + " (limiting factor " + limitingValue + "), using " +
+                       log.verbose("Too small time step " + store.timestep + " (limiting factor " + limitingValue + "), using " +
                                        minTimeStep + " instead.");
                        store.timestep = minTimeStep;
                } else {
-                       log.debug("Selected time step " + store.timestep + " (limiting factor " + limitingValue + ")");
+                       log.verbose("Selected time step " + store.timestep + " (limiting factor " + limitingValue + ")");
                }
                checkNaN(store.timestep);
                
@@ -169,13 +171,13 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                 * diminished by it affecting only 1/6th of the total, so it's an acceptable error.
                 */
                double thrustEstimate = store.thrustForce;
-               store.thrustForce = calculateThrust(status, store.timestep, store.longitudalAcceleration,
+               store.thrustForce = calculateThrust(status, store.timestep, store.longitudinalAcceleration,
                                store.atmosphericConditions, true);
                double thrustDiff = Math.abs(store.thrustForce - thrustEstimate);
                // Log if difference over 1%, recompute if over 10%
                if (thrustDiff > 0.01 * thrustEstimate) {
                        if (thrustDiff > 0.1 * thrustEstimate + 0.001) {
-                               log.info("Thrust estimate differs from correct value by " +
+                               log.debug("Thrust estimate differs from correct value by " +
                                                (Math.rint(1000 * (thrustDiff + 0.000001) / thrustEstimate) / 10.0) + "%," +
                                                " estimate=" + thrustEstimate +
                                                " correct=" + store.thrustForce +
@@ -183,7 +185,7 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                                                ", recomputing k1 parameters");
                                k1 = computeParameters(status, store);
                        } else {
-                               log.debug("Thrust estimate differs from correct value by " +
+                               log.verbose("Thrust estimate differs from correct value by " +
                                                (Math.rint(1000 * (thrustDiff + 0.000001) / thrustEstimate) / 10.0) + "%," +
                                                " estimate=" + thrustEstimate +
                                                " correct=" + store.thrustForce +
@@ -357,8 +359,8 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                        double momZ = store.forces.getCroll() * dynP * refArea * refLength;
                        
                        // Compute acceleration in rocket coordinates
-                       store.angularAcceleration = new Coordinate(momX / store.massData.getLongitudalInertia(),
-                                               momY / store.massData.getLongitudalInertia(), momZ / store.massData.getRotationalInertia());
+                       store.angularAcceleration = new Coordinate(momX / store.massData.getLongitudinalInertia(),
+                                               momY / store.massData.getLongitudinalInertia(), momZ / store.massData.getRotationalInertia());
                        
                        store.rollAcceleration = store.angularAcceleration.z;
                        // TODO: LOW: This should be hypot, but does it matter?
@@ -576,6 +578,8 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                }
                if (store.massData != null) {
                        data.setValue(FlightDataType.TYPE_MASS, store.massData.getCG().weight);
+                       data.setValue(FlightDataType.TYPE_LONGITUDINAL_INERTIA, store.massData.getLongitudinalInertia());
+                       data.setValue(FlightDataType.TYPE_ROTATIONAL_INERTIA, store.massData.getRotationalInertia());
                }
                
                data.setValue(FlightDataType.TYPE_THRUST_FORCE, store.thrustForce);
@@ -667,7 +671,7 @@ public class RK4SimulationStepper extends AbstractSimulationStepper {
                
                public FlightConditions flightConditions;
                
-               public double longitudalAcceleration = Double.NaN;
+               public double longitudinalAcceleration = Double.NaN;
                
                public MassData massData;