Commit b0386663 authored by Santiago Ospina's avatar Santiago Ospina

Use P0 finite element for finite volume method (fixed blocking is removed...

Use P0 finite element for finite volume method (fixed blocking is removed since it does nothing for P0 and one GFS)
Signed-off-by: Santiago Ospina De Los Ríos's avatarSantiago Ospina <santiago.ospina@iup.uni-heidelberg.de>
parent adb73315
......@@ -59,7 +59,7 @@ struct TransportSimulationTraits : public BaseTraits
using GFSaturation = GFSaturationType;
/// GFS Helper
using GFSHelper = Dune::Dorie::GridFunctionSpaceHelper<GV,RF,order,
using GFSHelper = Dune::Dorie::LinearSolverGridFunctionSpaceHelper<GV,RF,
BaseTraits::GridGeometryType>;
/// Problem GFS
using GFS = typename GFSHelper::Type;
......
......@@ -74,7 +74,7 @@ public:
using CON = Dune::PDELab::P0ParallelConstraints;
/// GFS type
using Type = typename Dune::PDELab::GridFunctionSpace<ES,FEM,CON,
Dune::PDELab::ISTL::VectorBackend<Dune::PDELab::ISTL::Blocking::fixed, 1>>;
Dune::PDELab::ISTL::VectorBackend<>>;
/// GFS Constraints Container type
using CC = typename Type::template ConstraintsContainer<RF>::Type;
......@@ -105,7 +105,7 @@ public:
using CON = Dune::PDELab::P0ParallelConstraints;
/// GFS type
using Type = typename Dune::PDELab::GridFunctionSpace<ES,FEM,CON,
Dune::PDELab::ISTL::VectorBackend<Dune::PDELab::ISTL::Blocking::fixed, 1>>;
Dune::PDELab::ISTL::VectorBackend<>>;
/// GFS Constraints Container type
using CC = typename Type::template ConstraintsContainer<RF>::Type;
......
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