Skip to content
Snippets Groups Projects
Commit f43674bc authored by graeser's avatar graeser
Browse files

Add PQ1 postfix to global assembler

parent f199c76b
Branches
No related tags found
No related merge requests found
...@@ -98,7 +98,7 @@ int main(int argc, char** argv) ...@@ -98,7 +98,7 @@ int main(int argc, char** argv)
auto rhsFunction = [](auto x) { auto rhsFunction = [](auto x) {
return x.two_norm() < .5; return x.two_norm() < .5;
}; };
assemblePoissonProblem(gridView, A, rhs, rhsFunction); assemblePoissonProblemPQ1(gridView, A, rhs, rhsFunction);
x.resize(rhs.size(), 0); x.resize(rhs.size(), 0);
......
...@@ -130,7 +130,7 @@ void assembleLocalRHS( ...@@ -130,7 +130,7 @@ void assembleLocalRHS(
template<class GridView, class RHSFunction> template<class GridView, class RHSFunction>
void assemblePoissonProblem( void assemblePoissonProblemPQ1(
const GridView& gridView, const GridView& gridView,
Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1>>& matrix, Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1>>& matrix,
Dune::BlockVector<Dune::FieldVector<double,1>>& rhs, Dune::BlockVector<Dune::FieldVector<double,1>>& rhs,
......
...@@ -146,7 +146,7 @@ int main(int argc, char** argv) ...@@ -146,7 +146,7 @@ int main(int argc, char** argv)
auto rhsFunction = [](auto x) { auto rhsFunction = [](auto x) {
return (x.two_norm() < .5)*100.0; return (x.two_norm() < .5)*100.0;
}; };
assemblePoissonProblem(gridView, A, rhs, rhsFunction); assemblePoissonProblemPQ1(gridView, A, rhs, rhsFunction);
x.resize(rhs.size(), 0); x.resize(rhs.size(), 0);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment