Commit a7648318 authored by Lukas Riedel's avatar Lukas Riedel

switched to P0 FEM for CG space. no apparent effect

parent fabdf001
......@@ -40,13 +40,14 @@ template<typename GridView, typename RF>
struct GridFunctionSpaceHelper<GridView,RF,0,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,0>;
using FEM = typename Dune::PDELab::P0LocalFiniteElementMap<DF,RF,dim>;
/// Constraints type of the GFS
using CON = Dune::PDELab::P0ParallelConstraints;
/// GFS type
......@@ -59,7 +60,9 @@ public:
static Type create (const GridView& gv)
{
ES es(gv);
auto fem = std::make_shared<FEM>(es);
Dune::GeometryType geo;
geo.makeSimplex(dim);
auto fem = std::make_shared<FEM>(geo);
auto con = std::make_shared<CON>();
return Type(es,fem,con);
}
......@@ -77,7 +80,7 @@ public:
/// Entity set of the GFS
using ES = Dune::PDELab::OverlappingEntitySet<GridView>;
/// FiniteElementMap type of GFS
using FEM = typename Dune::PDELab::QkDGLocalFiniteElementMap<DF,RF,0,dim>;
using FEM = typename Dune::PDELab::P0LocalFiniteElementMap<DF,RF,dim>;
/// Constraints type of the GFS
using CON = Dune::PDELab::P0ParallelConstraints;
/// GFS type
......@@ -90,7 +93,9 @@ public:
static Type create (const GridView& gv)
{
ES es(gv);
auto fem = std::make_shared<FEM>();
Dune::GeometryType geo;
geo.makeCube(dim);
auto fem = std::make_shared<FEM>(geo);
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