diff --git a/dune/gfe/riemanniantrsolver.hh b/dune/gfe/riemanniantrsolver.hh index 8a5c6f48ae4d4084dfa5c76c3fb670f496b6a34f..553badb7829cedc308fb6f188d52c5ce7d10a0a5 100644 --- a/dune/gfe/riemanniantrsolver.hh +++ b/dune/gfe/riemanniantrsolver.hh @@ -13,6 +13,8 @@ #include <dune/solvers/solvers/iterativesolver.hh> #include <dune/solvers/solvers/loopsolver.hh> +#include <dune/gfe/periodic1dpq1nodalbasis.hh> + #include "geodesicfeassembler.hh" #include <dune/grid/utility/globalindexset.hh> #include <dune/gfe/parallel/globalp1mapper.hh> @@ -32,6 +34,15 @@ struct MapperFactory<GridView, Dune::Functions::PQkNodalBasis<GridView,1> > typedef Dune::MultipleCodimMultipleGeomTypeMapper<GridView, Dune::MCMGVertexLayout> LocalMapper; }; +// This case is not going to actually work, but I need the specialization to make +// the sequential code compile. +template <typename GridView> +struct MapperFactory<GridView, Dune::Functions::Periodic1DPQ1NodalBasis<GridView> > +{ + typedef Dune::GlobalP1Mapper<GridView> GlobalMapper; + typedef Dune::MultipleCodimMultipleGeomTypeMapper<GridView, Dune::MCMGVertexLayout> LocalMapper; +}; + template <typename GridView> struct MapperFactory<GridView, Dune::Functions::PQkNodalBasis<GridView,2> > {