Fix scaling of Jacobian of new LinearSolver

The previous version was comparing with 1.0 instead of assigning it.

SimulationRuntime/cpp/Solver/LinearSolver/LinearSolver.cpp

@@ -206,18 +206,14 @@ void LinearSolver::solve()
 
     memcpy(_A, Atemp, _dimSys*_dimSys*sizeof(double));
 
+    // scale Jacobian
+    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]);
       }
     }
 
-    for (int i=0; i<_dimSys; i++) {
-      if (_fNominal[i] == 0.0) {
-        _fNominal[i]==1.0;// if the row contains only zeros, there is no need to scale that row.
-      }
-    }
-
     LOGGER_WRITE_VECTOR("fNominal", _fNominal, _dimSys, LC_LS, LL_DEBUG);
 
     for (int j = 0, idx = 0; j < _dimSys; j++)
@@ -227,8 +223,6 @@ void LinearSolver::solve()
     for (int i = 0; i < _dimSys; i++)
       _b[i] /= _fNominal[i];
 
-
-
     if (_generateoutput) {
       std::cout << std::endl;
       std::cout << "We solve a linear system with coefficient matrix" << std::endl;