Skip to content
Snippets Groups Projects
Commit 75b35cea authored by Elias Pipping's avatar Elias Pipping
Browse files

MMGStep: Restore precision

parent b32be161
No related branches found
No related tags found
No related merge requests found
...@@ -316,8 +316,12 @@ void MonotoneMGStep<MatrixType, VectorType>::iterate() ...@@ -316,8 +316,12 @@ void MonotoneMGStep<MatrixType, VectorType>::iterate()
} }
// Debug: output energy // Debug: output energy
if (level==(int) this->numLevels()-1 && this->verbosity_==NumProc::FULL) if (level==(int) this->numLevels()-1 && this->verbosity_==NumProc::FULL) {
std::streamsize const oldPrecision = std::cout.precision()
std::cout << "Total energy: " std::cout << "Total energy: "
<< std::setprecision(10) << computeEnergy(*mat[level], *x[level], rhs[level]) << std::endl; << std::setprecision(10)
<< computeEnergy(*mat[level], *x[level], rhs[level]) << std::endl
<< std::setprecision(oldPrecision);
}
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment