BaseDerivativesDataOperators.hpp

Go to the documentation of this file.

/** \file BaseDerivativesDataOperators.hpp
  * \brief Base derivatives data operators
 
*/
 
 
 
#ifndef __BASE_DERIVATIVES_DATA_OPERATORS_HPP__
#define __BASE_DERIVATIVES_DATA_OPERATORS_HPP__
 
namespace MoFEM {
 
struct OpBaseDerivativesBase : public ForcesAndSourcesCore::UserDataOperator {
 
  OpBaseDerivativesBase(boost::shared_ptr<MatrixDouble> base_mass_ptr,
                        boost::shared_ptr<EntitiesFieldData> data_l2,
                        const FieldApproximationBase b, const FieldSpace s,
                        int verb = QUIET, Sev sev = Sev::verbose);
 
protected:
  FieldApproximationBase base;
  int verbosity;
  Sev severityLevel;
 
  boost::shared_ptr<MatrixDouble> baseMassPtr;
  boost::shared_ptr<EntitiesFieldData> dataL2;
 
  using GetOrderFun = boost::function<int()>;
  MoFEMErrorCode calculateBase(GetOrderFun get_order);
  MoFEMErrorCode calculateMass();
  
};
 
template <int BASE_DIM>
struct OpBaseDerivativesMass;
 
template <> struct OpBaseDerivativesMass<1> : public OpBaseDerivativesBase  {
 
  using OpBaseDerivativesBase::OpBaseDerivativesBase;
 
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
};
 
template <> struct OpBaseDerivativesMass<3> : public OpBaseDerivativesMass<1> {
  using OpBaseDerivativesMass<1>::OpBaseDerivativesMass;
};
 
template <int DIM> struct OpBaseDerivativesSetHOInvJacobian;
 
template <>
struct OpBaseDerivativesSetHOInvJacobian<2>
    : public OpSetInvJacSpaceForFaceImpl<2, 1> {
 
  OpBaseDerivativesSetHOInvJacobian(
      boost::shared_ptr<EntitiesFieldData> data_l2,
      boost::shared_ptr<MatrixDouble> inv_jac_ptr);
 
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
 
protected:
  boost::shared_ptr<EntitiesFieldData> dataL2;
};
 
template <int BASE_DIM>
struct OpBaseDerivativesNext;
 
/**
 * @brief Specialisation for calculate derivatives for scalar base functions
 * 
 * @tparam  
 */
template <> struct OpBaseDerivativesNext<1> : public OpBaseDerivativesBase {
 
  OpBaseDerivativesNext(int derivative,
                        boost::shared_ptr<MatrixDouble> base_mass_ptr,
                        boost::shared_ptr<EntitiesFieldData> data_l2,
                        const FieldApproximationBase b, const FieldSpace s,
                        int verb = QUIET, Sev sev = Sev::verbose);
 
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
 
protected:
  int calcBaseDerivative;
  MatrixDouble nF;
 
  template <int BASE_DIM>
  MoFEMErrorCode doWorkImpl(int side, EntityType type,
                            EntitiesFieldData::EntData &data);
 
  template <int BASE_DIM, int SPACE_DIM>
  MoFEMErrorCode setBaseImpl(EntitiesFieldData::EntData &data,
                             EntitiesFieldData::EntData &ent_data);
};
 
/**
 * @brief Specialisation for calculate derivatives for scalar base functions
 * 
 * @tparam  
 */
template <> struct OpBaseDerivativesNext<3> : public OpBaseDerivativesNext<1> {
 
  using OpBaseDerivativesNext<1>::OpBaseDerivativesNext;
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
};
 
} // namespace MoFEM
 
#endif //__BASE_DERIVATIVES_DATA_OPERATORS_HPP__