mapping_wrapper.h

Go to the documentation of this file.

#ifndef MAPPING_WRAPPER_H_
#define MAPPING_WRAPPER_H_

#include <dofs/dof_handler.h>
#include <multigrid/mg_dof_handler.h>
#include <hp/dof_handler.h>
#include <hp/mapping_collection.h>
#include <fe/mapping_q.h>

namespace DOpEWrapper
{ 
  template<int dim, template<int, int> class DH = dealii::DoFHandler>
    class Mapping
    {
      private:
        Mapping()
        {
        }

      ~Mapping()
      {
      }

    };

  /************************************************************************************/

  template<int dim>
    class Mapping<dim, dealii::DoFHandler > : public dealii::MappingQ<dim>
    {
      public:
        Mapping(const unsigned int p, const bool use_mapping_q_on_all_elements =
            false) :
            dealii::MappingQ<dim>(p, use_mapping_q_on_all_elements)
        {
        }

        Mapping(const dealii::MappingQ<dim> &mapping) :
            dealii::MappingQ<dim>(mapping)
        {
        }

      ~Mapping()
        {
        }

        const typename dealii::MappingQ<dim> &
        operator[](const unsigned int index) const
        {
          assert(index == 0);
          return *this;
        }

    };

  /************************************************************************************/

  template<int dim>
    class Mapping<dim, dealii::MGDoFHandler > : public dealii::MappingQ<dim>
    {
      public:
        Mapping(const unsigned int p, const bool use_mapping_q_on_all_elements =
            false)
            : dealii::MappingQ<dim>(p, use_mapping_q_on_all_elements)
        {
        }

        Mapping(const dealii::MappingQ<dim> &mapping)
            : dealii::MappingQ<dim>(mapping)
        {
        }

      ~Mapping()
        {
        }

        const typename dealii::MappingQ<dim> &
        operator[](const unsigned int index) const
        {
          assert(index == 0);
          return *this;
        }

    };

  /************************************************************************************/

  template<int dim>
    class Mapping<dim, dealii::hp::DoFHandler> : public dealii::hp::MappingCollection<
        dim>
    {
      public:
        Mapping() :
            dealii::hp::MappingCollection<dim>()
        {
        }

      ~Mapping()
        {
        }

        Mapping(const dealii::Mapping<dim>& mapping)
            : dealii::hp::MappingCollection<dim>(mapping)
        {
        }
        Mapping(const dealii::hp::MappingCollection<dim> & mapping_collection) :
            dealii::hp::MappingCollection<dim>(mapping_collection)
        {
        }
    };

  /************************************************************************************/

  template<int dim, template<int, int> class DH>
    struct StaticMappingQ1
    {
    };

  template<int dim>
    struct StaticMappingQ1<dim, dealii::DoFHandler>
    {
      public:
        static Mapping<dim, dealii::DoFHandler> mapping_q1;
    };

 template<int dim>
    struct StaticMappingQ1<dim, dealii::MGDoFHandler >
    {
      public:
        static Mapping<dim, dealii::MGDoFHandler > mapping_q1;
    };

  template<int dim>
    struct StaticMappingQ1<dim, dealii::hp::DoFHandler>
    {
      public:
        static Mapping<dim, dealii::hp::DoFHandler> mapping_q1;
    };
}

#endif /* MAPPING_WRAPPER_H_ */