Skip to content
Snippets Groups Projects
Commit beb8c92b authored by Dr. Felix Tobias Schindler's avatar Dr. Felix Tobias Schindler
Browse files

[examples.elliptic.continuousgalerkin] updated constraints

parent d84b573a
No related branches found
No related tags found
No related merge requests found
......@@ -231,7 +231,9 @@ int main(int argc, char** argv)
*boundaryInfo, space.mapper().maxNumDofs(), space.mapper().maxNumDofs());
rhsVector->backend() =
forceVector->backend() + neumannVector->backend() - systemMatrix->backend() * dirichletVector->backend();
systemAssembler.applyConstraints(dirichletConstraints, *systemMatrix, *rhsVector);
systemAssembler.addLocalConstraints(dirichletConstraints, *systemMatrix);
systemAssembler.addLocalConstraints(dirichletConstraints, *rhsVector);
systemAssembler.applyConstraints();
info << "done (took " << timer.elapsed() << " sec)" << std::endl;
info << "solving linear system (of size " << systemMatrix->rows() << "x" << systemMatrix->cols() << ")"
......
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