Fix initialization of nominal values in DgesvSolver, ticket:4459

OpenModelica · Jul 1, 2017 · 0d0414c · 0d0414c
1 parent 61870e3
commit 0d0414c
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/SimulationRuntime/cpp/Solver/Dgesv/DgesvSolver.cpp b/SimulationRuntime/cpp/Solver/Dgesv/DgesvSolver.cpp
@@ -138,6 +138,7 @@ void DgesvSolver::solve()
 
   memcpy(_A, Atemp, _dimSys*_dimSys*sizeof(double));
 
+  std::fill(_fNominal, _fNominal + _dimSys, 1e-6);
   for (int j = 0, idx = 0; j < _dimSys; j++)
     for (int i = 0; i < _dimSys; i++, idx++)
       _fNominal[i] = std::max(std::abs(Atemp[idx]), _fNominal[i]);