Skip to content
Snippets Groups Projects
Commit 9cb09946 authored by Klaus Böhnlein's avatar Klaus Böhnlein
Browse files

bugfix: remove if(instrumented_) when setting up h1SemiNorm.

	The h1SemiNorm was only setup when the instrumented-flag was true.
	These are however independent functionalities.
parent f4040907
No related branches found
No related tags found
1 merge request!150Add option to choose a different Regularization-Norm for the Riemannian Newton with Hessian modification
......@@ -109,14 +109,12 @@ setup(const GridType& grid,
ScalarMatrixType,
LocalMapper,
LocalMapper> matrixComm(*globalMapper_, grid_->leafGridView(), localMapper, localMapper, 0);
if (instrumented_) {
auto A = std::make_shared<ScalarMatrixType>(matrixComm.reduceAdd(localA));
auto A = std::make_shared<ScalarMatrixType>(matrixComm.reduceAdd(localA));
#else
if (instrumented_) {
auto A = std::make_shared<ScalarMatrixType>(localA);
auto A = std::make_shared<ScalarMatrixType>(localA);
#endif
h1SemiNorm_ = std::make_shared<H1SemiNorm<CorrectionType> >(A);
}
h1SemiNorm_ = std::make_shared<H1SemiNorm<CorrectionType> >(A);
//////////////////////////////////////////////////////////////////
// Create the inner solver using a cholmod solver
//////////////////////////////////////////////////////////////////
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment