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

More porting from the dune-fufem bases to the dune-functions bases

parent 1c959fb9
No related branches found
No related tags found
No related merge requests found
...@@ -6,7 +6,7 @@ ...@@ -6,7 +6,7 @@
#include <dune/istl/io.hh> #include <dune/istl/io.hh>
#include <dune/fufem/functionspacebases/p1nodalbasis.hh> #include <dune/functions/functionspacebases/pqknodalbasis.hh>
#include <dune/fufem/functionspacebases/dunefunctionsbasis.hh> #include <dune/fufem/functionspacebases/dunefunctionsbasis.hh>
#include <dune/fufem/assemblers/operatorassembler.hh> #include <dune/fufem/assemblers/operatorassembler.hh>
#include <dune/fufem/assemblers/localassemblers/laplaceassembler.hh> #include <dune/fufem/assemblers/localassemblers/laplaceassembler.hh>
...@@ -197,11 +197,11 @@ setup(const GridType& grid, ...@@ -197,11 +197,11 @@ setup(const GridType& grid,
{ {
if (numLevels>1) { if (numLevels>1) {
typedef typename TruncatedCompressedMGTransfer<CorrectionType0>::TransferOperatorType TransferOperatorType; typedef typename TruncatedCompressedMGTransfer<CorrectionType0>::TransferOperatorType TransferOperatorType;
P1NodalBasis<typename GridType::LeafGridView,double> p1Basis(grid_->leafGridView()); DuneFunctionsBasis<Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView,1> > p1Basis(grid_->leafGridView());
TransferOperatorType pkToP1TransferMatrix; TransferOperatorType pkToP1TransferMatrix;
assembleBasisInterpolationMatrix<TransferOperatorType, assembleBasisInterpolationMatrix<TransferOperatorType,
P1NodalBasis<typename GridType::LeafGridView,double>, DuneFunctionsBasis<Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView,1> >,
FufemBasis0>(pkToP1TransferMatrix,p1Basis,assembler->basis0_); FufemBasis0>(pkToP1TransferMatrix,p1Basis,assembler->basis0_);
mmgStep0->mgTransfer_.back() = new TruncatedCompressedMGTransfer<CorrectionType0>; mmgStep0->mgTransfer_.back() = new TruncatedCompressedMGTransfer<CorrectionType0>;
...@@ -237,11 +237,11 @@ setup(const GridType& grid, ...@@ -237,11 +237,11 @@ setup(const GridType& grid,
{ {
if (numLevels>1) { if (numLevels>1) {
typedef typename TruncatedCompressedMGTransfer<CorrectionType1>::TransferOperatorType TransferOperatorType; typedef typename TruncatedCompressedMGTransfer<CorrectionType1>::TransferOperatorType TransferOperatorType;
P1NodalBasis<typename GridType::LeafGridView,double> p1Basis(grid_->leafGridView()); DuneFunctionsBasis<Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView,1> > p1Basis(grid_->leafGridView());
TransferOperatorType pkToP1TransferMatrix; TransferOperatorType pkToP1TransferMatrix;
assembleBasisInterpolationMatrix<TransferOperatorType, assembleBasisInterpolationMatrix<TransferOperatorType,
P1NodalBasis<typename GridType::LeafGridView,double>, DuneFunctionsBasis<Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView,1> >,
FufemBasis1>(pkToP1TransferMatrix,p1Basis,assembler->basis1_); FufemBasis1>(pkToP1TransferMatrix,p1Basis,assembler->basis1_);
mmgStep0->mgTransfer_.back() = new TruncatedCompressedMGTransfer<CorrectionType1>; mmgStep0->mgTransfer_.back() = new TruncatedCompressedMGTransfer<CorrectionType1>;
......
...@@ -16,10 +16,6 @@ ...@@ -16,10 +16,6 @@
#include <dune/solvers/solvers/loopsolver.hh> #include <dune/solvers/solvers/loopsolver.hh>
#include <dune/solvers/iterationsteps/mmgstep.hh> #include <dune/solvers/iterationsteps/mmgstep.hh>
#include <dune/fufem/functionspacebases/p1nodalbasis.hh>
#include <dune/fufem/functionspacebases/p2nodalbasis.hh>
#include <dune/fufem/functionspacebases/p3nodalbasis.hh>
#include <dune/gfe/mixedgfeassembler.hh> #include <dune/gfe/mixedgfeassembler.hh>
/** \brief Riemannian trust-region solver for geodesic finite-element problems */ /** \brief Riemannian trust-region solver for geodesic finite-element problems */
......
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