diff --git a/dune/solvers/solvers/loopsolver.cc b/dune/solvers/solvers/loopsolver.cc index 5dee67692f51919d127260289b42730182cdbacf..35f3b22e9f77c8acbd98ccbdbed9cac6c14127f1 100644 --- a/dune/solvers/solvers/loopsolver.cc +++ b/dune/solvers/solvers/loopsolver.cc @@ -87,7 +87,7 @@ void ::LoopSolver<VectorType, BitVectorType>::solve() oldSolution -= iterationStep_->getSol(); double normOfCorrection; - double normOfError; + double normOfError=std::numeric_limits<double>::quiet_NaN(); double convRate; normOfCorrection = this->errorNorm_->operator()(oldSolution);