guard atol again nominal <= 0

perhaps we should guard nominal value git-svn-id: https://openmodelica.org/svn/OpenModelica/trunk@21334 f25d12d1-65f4-0310-ae8a-bbce733d8d8e
OpenModelica · Jun 30, 2014 · be7e0df · be7e0df
1 parent 65e0b71
commit be7e0df
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/SimulationRuntime/c/simulation/solver/dassl.c b/SimulationRuntime/c/simulation/solver/dassl.c
@@ -254,7 +254,7 @@ int dassl_initial(DATA* data, SOLVER_INFO* solverInfo, DASSL_DATA *dasslData)
   for(i=0;i<data->modelData.nStates;++i)
   {
     dasslData->rtol[i] = data->simulationInfo.tolerance;
-    dasslData->atol[i] = data->simulationInfo.tolerance * data->modelData.realVarsData[i].attribute.nominal;
+    dasslData->atol[i] = data->simulationInfo.tolerance * fmax(fabs(data->modelData.realVarsData[i].attribute.nominal),1e-32);
   }
 
   /* setup internal ring buffer for dassl */