Skip to content
Snippets Groups Projects
Commit d5cb47ac authored by Oliver Sander's avatar Oliver Sander
Browse files

Adapt to recent changes in dune-functions

parent 69a96e14
No related branches found
No related tags found
No related merge requests found
...@@ -77,7 +77,7 @@ getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const ...@@ -77,7 +77,7 @@ getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
nb.resize(n, n); nb.resize(n, n);
// A view on the FE basis on a single element // A view on the FE basis on a single element
typename Basis::LocalView localView(&basis_); auto localView = basis_.localView();
auto localIndexSet = basisIndexSet_.localIndexSet(); auto localIndexSet = basisIndexSet_.localIndexSet();
ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>(); ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>();
...@@ -129,7 +129,7 @@ assembleGradientAndHessian(const std::vector<TargetSpace>& sol, ...@@ -129,7 +129,7 @@ assembleGradientAndHessian(const std::vector<TargetSpace>& sol,
gradient = 0; gradient = 0;
// A view on the FE basis on a single element // A view on the FE basis on a single element
typename Basis::LocalView localView(&basis_); auto localView = basis_.localView();
auto localIndexSet = basisIndexSet_.localIndexSet(); auto localIndexSet = basisIndexSet_.localIndexSet();
ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>(); ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>();
...@@ -186,7 +186,7 @@ assembleGradient(const std::vector<TargetSpace>& sol, ...@@ -186,7 +186,7 @@ assembleGradient(const std::vector<TargetSpace>& sol,
grad = 0; grad = 0;
// A view on the FE basis on a single element // A view on the FE basis on a single element
typename Basis::LocalView localView(&basis_); auto localView = basis_.localView();
auto localIndexSet = basisIndexSet_.localIndexSet(); auto localIndexSet = basisIndexSet_.localIndexSet();
ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>(); ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>();
...@@ -231,7 +231,7 @@ computeEnergy(const std::vector<TargetSpace>& sol) const ...@@ -231,7 +231,7 @@ computeEnergy(const std::vector<TargetSpace>& sol) const
DUNE_THROW(Dune::Exception, "Coefficient vector doesn't match the function space basis!"); DUNE_THROW(Dune::Exception, "Coefficient vector doesn't match the function space basis!");
// A view on the FE basis on a single element // A view on the FE basis on a single element
typename Basis::LocalView localView(&basis_); auto localView = basis_.localView();
auto localIndexSet = basisIndexSet_.localIndexSet(); auto localIndexSet = basisIndexSet_.localIndexSet();
ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>(); ElementIterator it = basis_.gridView().template begin<0,Dune::Interior_Partition>();
......
...@@ -149,9 +149,9 @@ int main (int argc, char *argv[]) try ...@@ -149,9 +149,9 @@ int main (int argc, char *argv[]) try
GridView gridView = grid->leafGridView(); GridView gridView = grid->leafGridView();
#ifdef SECOND_ORDER #ifdef SECOND_ORDER
typedef Dune::Functions::PQKNodalBasis<typename GridType::LeafGridView, 2> FEBasis; typedef Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView, 2> FEBasis;
#else #else
typedef Dune::Functions::PQKNodalBasis<typename GridType::LeafGridView, 1> FEBasis; typedef Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView, 1> FEBasis;
#endif #endif
FEBasis feBasis(gridView); FEBasis feBasis(gridView);
......
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