Ticket #2539: stepSize.diff

SimulationRuntime/c/simulation/solver/perform_simulation.c

@@ -96 +96 @@
   SIMULATION_INFO *simInfo = &(data->simulationInfo);
 
   solverInfo->currentTime = simInfo->startTime;
+  
+  unsigned int __currStepNo = 0;
 
   if(measure_time_flag)
   {
@@ -151 +153 @@
       rotateRingBuffer(data->simulationData, 1, (void**) data->localData);
 
       /***** Calculation next step size *****/
-      /* Calculate new step size after an event */
       if(solverInfo->didEventStep == 1)
       {
-        if((solverInfo->currentTime - solverInfo->laststep) + DBL_EPSILON > simInfo->stepSize)
-          solverInfo->currentStepSize = simInfo->stepSize;
-        else
-          solverInfo->currentStepSize = simInfo->stepSize - (solverInfo->currentTime - solverInfo->laststep);
-        
         infoStreamPrint(LOG_SOLVER, 0, "offset value for the next step: %.10f", (solverInfo->currentTime - solverInfo->laststep));
       }
       else
       {
-        solverInfo->currentStepSize = simInfo->stepSize;
+        __currStepNo++;
       }
-
-      /* adjust final step? */
-      if(solverInfo->currentTime + solverInfo->currentStepSize > simInfo->stopTime) {
-        solverInfo->currentStepSize = simInfo->stopTime - solverInfo->currentTime;
-      }
+      solverInfo->currentStepSize = (double)(__currStepNo*(simInfo->stopTime-simInfo->startTime))/(simInfo->numSteps) + simInfo->startTime - solverInfo->currentTime;
       /***** End calculation next step size *****/
 
       /* check for next time event */

testsuite/openmodelica/cruntime/simoptions/testOutputIntervalDASSL.mos

@@ -56 +56 @@
 //     simulationOptions = "startTime = 0.0, stopTime = 2.0, numberOfIntervals = 6, tolerance = 0.000001, method = 'dassl', fileNamePrefix = 'testModel', options = '', outputFormat = 'mat', variableFilter = '.*', measureTime = false, cflags = '', simflags = '-noEventEmit'",
 //     messages = ""
 // end SimulationResult;
-// 8
+// 7
 // record SimulationResult
 //     resultFile = "testModel_res.mat",
 //     simulationOptions = "startTime = 0.0, stopTime = 3.0, numberOfIntervals = 500, tolerance = 0.000001, method = 'dassl', fileNamePrefix = 'testModel', options = '', outputFormat = 'mat', variableFilter = '.*', measureTime = false, cflags = '', simflags = '-noEventEmit'",

testsuite/openmodelica/cruntime/simoptions/testOutputIntervalEuler.mos

@@ -56 +56 @@
 //     simulationOptions = "startTime = 0.0, stopTime = 2.0, numberOfIntervals = 6, tolerance = 0.000001, method = 'euler', fileNamePrefix = 'testModel', options = '', outputFormat = 'mat', variableFilter = '.*', measureTime = false, cflags = '', simflags = '-noEventEmit'",
 //     messages = ""
 // end SimulationResult;
-// 8
+// 7
 // record SimulationResult
 //     resultFile = "testModel_res.mat",
 //     simulationOptions = "startTime = 0.0, stopTime = 3.0, numberOfIntervals = 500, tolerance = 0.000001, method = 'euler', fileNamePrefix = 'testModel', options = '', outputFormat = 'mat', variableFilter = '.*', measureTime = false, cflags = '', simflags = '-noEventEmit'",

testsuite/openmodelica/cruntime/simoptions/testOutputIntervalRK.mos

@@ -62 +62 @@
 //     simulationOptions = "startTime = 0.0, stopTime = 2.0, numberOfIntervals = 6, tolerance = 0.000001, method = 'rungekutta', fileNamePrefix = 'testModel', options = '', outputFormat = 'mat', variableFilter = '.*', measureTime = false, cflags = '', simflags = '-noEventEmit'",
 //     messages = ""
 // end SimulationResult;
-// 8
+// 7
 // endResult