Fix bug in homotopy method

ptaeuber · OpenModelica-Hudson · commit 21ce9323ad26 · 2017-07-26T16:49:01.000+02:00
Do not set lambda=1 in each step when using global homotopy
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);
 

Original file line number	Diff line number	Diff line change
`@@ -848,7 +848,9 @@ int solve_nonlinear_system(DATA data, threadData_t threadData, int sysNumber)`
`848`	`848`	`#endif`
`849`	`849`	`}`
`850`	`850`
`851`		`- data->simulationInfo->lambda = 1.0;`
	`851`	`+ if(data->callback->useHomotopy == 0)`
	`852`	`+ data->simulationInfo->lambda = 1.0;`
	`853`	`+`
`852`	`854`	`/* SOLVE! */`
`853`	`855`	`nonlinsys->solved = solveNLS(data, threadData, sysNumber);`
`854`	`856`