Fix bug in homotopy method

Do not set lambda=1 in each step when using global homotopy
OpenModelica · Jul 26, 2017 · 21ce932 · 21ce932
1 parent 441b181
commit 21ce932
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/SimulationRuntime/c/simulation/solver/nonlinearSystem.c b/SimulationRuntime/c/simulation/solver/nonlinearSystem.c
@@ -848,7 +848,9 @@ int solve_nonlinear_system(DATA *data, threadData_t *threadData, int sysNumber)
 #endif
   }
 
-  data->simulationInfo->lambda = 1.0;
+  if(data->callback->useHomotopy == 0)
+    data->simulationInfo->lambda = 1.0;
+
   /* SOLVE! */
   nonlinsys->solved = solveNLS(data, threadData, sysNumber);