fedealwrapper.h

Go to the documentation of this file.

#ifndef  _MY_FEDealWrapper_
#define _MY_FEDealWrapper_
#include "finiteelementinterface.h"

class FEDealWrapper : public FiniteElementInterface
{
 private:
   ArrayOfVecDouble m_support_points;
 protected:
  dealii::FiniteElement<3> *_dealFE;
  std::vector<VecIndex > _face_to_cell_map;

 public:
  FEDealWrapper()
    :_dealFE(NULL)
  {
  }

 public:
  void wrap(dealii::FiniteElement<3> &fe)
  {
    _dealFE = &fe;
    //Build the maps between the local dofs of the finite element
    //with the local dofs of the faces
    _face_to_cell_map.resize(6);

    for (unsigned i=0;i<6;i++) //for all faces
    {
      //initialize the map of that face
      _face_to_cell_map[i].resize(fe.n_dofs_per_face());

      //for each dof j on the face
      for (unsigned j=0;j<fe.n_dofs_per_face();j++)
        _face_to_cell_map[i][j]=fe.face_to_cell_index(j,i);
    }

    m_support_points=_dealFE->get_unit_support_points();

  }
  virtual unsigned n_fields(){return _dealFE->n_components();}

  ~FEDealWrapper(){}


  virtual Index n_dofs_per_face()
  {
    return _dealFE->n_dofs_per_face();
  }
  
  virtual const VecIndex& face_to_cell_local_map(Index faceId)
  {
    return _face_to_cell_map[faceId];
  }



  virtual Index n_local_dofs()
  {
    return _dealFE->n_dofs_per_cell();
  }


  
  virtual Index local_dof_field(unsigned local_dof)
  {
    return _dealFE->system_to_component_index(local_dof).first;
  }

  virtual double shape_value(unsigned i,const VecDouble &p)
  {
    Point3D pt(p(0),p(1),p(2));
    return _dealFE->shape_value(i,pt);
  }

  virtual void shape_grad(unsigned i,const VecDouble  &pt,VecDouble &vGrad)
  {
    Point3D p(pt(0),pt(1),pt(2));
    _dealFE->shape_grad(i,p).unroll(vGrad);
  }

  virtual const ArrayOfVecDouble& get_support_points()
  {
    return m_support_points;
  }


    
  
};





#endif

---