Commit fd9be4b2 authored by Lukas Riedel's avatar Lukas Riedel

Switch to OPBLocalFiniteElementMap for simplex grids

parent 7b1dcc23
#ifndef DUNE_DORIE_UTIL_HH
#define DUNE_DORIE_UTIL_HH
#include <dune/pdelab/finiteelementmap/opbfem.hh>
namespace Dune{
namespace Dorie{
......@@ -114,13 +116,14 @@ template<typename GridView, typename RF, int order>
struct GridFunctionSpaceHelper<GridView,RF,order,Dune::GeometryType::BasicType::simplex>
{
private:
static constexpr int dim = GridView::dimension;
using DF = typename GridView::Grid::ctype;
public:
/// Entity set of the GFS
using ES = Dune::PDELab::OverlappingEntitySet<GridView>;
/// FiniteElementMap type of GFS
using FEM = typename Dune::PDELab::PkLocalFiniteElementMap<ES,DF,RF,order>;
using FEM = typename Dune::PDELab::OPBLocalFiniteElementMap<DF, RF, order, dim, Dune::GeometryType::BasicType::simplex>;
/// Constraints type of the GFS
using CON = Dune::PDELab::P0ParallelConstraints;
/// GFS type
......@@ -133,7 +136,7 @@ public:
static Type create (const GridView& gv)
{
ES es(gv);
auto fem = std::make_shared<FEM>(es);
auto fem = std::make_shared<FEM>();
auto con = std::make_shared<CON>();
return Type(es,fem,con);
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment