discretefunctionselector.hh

#ifndef __DUNE_ACFEM_DISCRETE_FUNCTION_SELECTOR_HH__
#define __DUNE_ACFEM_DISCRETE_FUNCTION_SELECTOR_HH__
 
#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
#include <dune/fem/space/lagrange.hh>
#include <dune/fem/space/discontinuousgalerkin.hh>
#include <dune/fem/function/adaptivefunction.hh>
#include <dune/fem/function/blockvectorfunction.hh>
 
#include "solverselector.hh"
 
namespace Dune {
 
  namespace ACFem {
 
    template<class DiscreteFunctionSpaceType, SolverFamily::Family family>
    struct DiscreteFunctionSelector;
 
    template<class DiscreteFunctionSpaceType>
    struct DiscreteFunctionSelector<DiscreteFunctionSpaceType, SolverFamily::ISTL>
    {
      typedef Dune::Fem::ISTLBlockVectorDiscreteFunction<DiscreteFunctionSpaceType> DiscreteFunctionType;
    };
 
    template<class DiscreteFunctionSpaceType>
    struct DiscreteFunctionSelector<DiscreteFunctionSpaceType, SolverFamily::FEM>
    {
      typedef Dune::Fem::AdaptiveDiscreteFunction<DiscreteFunctionSpaceType> DiscreteFunctionType;
    };
 
    template<class DiscreteFunctionSpaceType>
    struct DiscreteFunctionSelector<DiscreteFunctionSpaceType, SolverFamily::PETSC>
    {
      typedef Dune::Fem::PetscDiscreteFunction<DiscreteFunctionSpaceType> DiscreteFunctionType;
    };
 
    template<class Grid, class FunctionSpace, unsigned PolOrder, SolverFamily::Family family = SolverFamily::family>
    struct LagrangeLeafFunctionSelector
    {
      typedef Grid GridType; 
      typedef FunctionSpace FunctionSpaceType; 
 
      enum { polynomialOrder = PolOrder };
 
      typedef Dune::Fem::AdaptiveLeafGridPart<GridType, Dune::InteriorBorder_Partition> GridPartType;
      typedef Dune::Fem::LagrangeDiscreteFunctionSpace<FunctionSpaceType, GridPartType, PolOrder> DiscreteFunctionSpaceType;
 
      typedef typename DiscreteFunctionSelector<DiscreteFunctionSpaceType, family>::DiscreteFunctionType DiscreteFunctionType;
    };
 
    template<class Grid, class FunctionSpace, unsigned PolOrder, SolverFamily::Family family = SolverFamily::family>
    struct DiscontinuousGalerkinLeafFunctionSelector
    {
      typedef Grid GridType;
      typedef FunctionSpace FunctionSpaceType;
 
      enum { polynomialOrder = PolOrder };
 
      typedef Dune::Fem::AdaptiveLeafGridPart<GridType, Dune::InteriorBorder_Partition> GridPartType;
      typedef Dune::Fem::DiscontinuousGalerkinSpace<FunctionSpaceType, GridPartType, PolOrder> DiscreteFunctionSpaceType;
 
      typedef typename DiscreteFunctionSelector<DiscreteFunctionSpaceType, family>::DiscreteFunctionType DiscreteFunctionType;
    };
 
    template<class Grid, class FunctionSpace, unsigned PolOrder, SolverFamily::Family family = SolverFamily::family>
    struct DiscontinuousLagrangeLeafFunctionSelector
    {
      typedef Grid GridType;
      typedef FunctionSpace FunctionSpaceType;
 
      enum { polynomialOrder = PolOrder };
 
      typedef Dune::Fem::AdaptiveLeafGridPart<GridType, Dune::InteriorBorder_Partition> GridPartType;
      typedef Dune::Fem::LagrangeDiscontinuousGalerkinSpace<FunctionSpaceType, GridPartType, PolOrder> DiscreteFunctionSpaceType;
 
      typedef typename DiscreteFunctionSelector<DiscreteFunctionSpaceType, family>::DiscreteFunctionType DiscreteFunctionType;
    };
 
 
  } // ACFem::
 
} // Dune::
 
#endif // __DUNE_ACFEM_DISCRETE_FUNCTION_SELECTOR_HH__