DGProjection.hpp

Go to the documentation of this file.

/*/**
 * @file DGProjection.hpp
 * @author your name (you@domain.com)
 * @brief
 * @version 0.1
 * @date 2023-08-23
 *
 * @copyright Copyright (c) 2023
 *
 */
 
#ifndef DGPROJECTION_HPP
#define DGPROJECTION_HPP
 
namespace MoFEM {
 
template <int DIM>
MoFEMErrorCode
setDGSetIntegrationPoints(ForcesAndSourcesCore::GaussHookFun &set_hook,
                          ForcesAndSourcesCore::RuleHookFun rule,
                          boost::shared_ptr<Range> edges_ptr);
 
template <>
MoFEMErrorCode
setDGSetIntegrationPoints<2>(ForcesAndSourcesCore::GaussHookFun &set_hook,
                             ForcesAndSourcesCore::RuleHookFun rule,
                             boost::shared_ptr<Range> edges_ptr);
 
struct OpDGProjectionMassMatrix : public OpBaseDerivativesBase {
  OpDGProjectionMassMatrix(int order, boost::shared_ptr<MatrixDouble> 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);
 
private:
  int baseOrder;
};
 
/**
 * @brief Evaluate field for given DG projection coefficients.
 * 
 */
struct OpDGProjectionCoefficients : public OpBaseDerivativesBase {
 
  OpDGProjectionCoefficients(
      boost::shared_ptr<MatrixDouble> data_ptr,
      boost::shared_ptr<MatrixDouble> coeffs_ptr,
      boost::shared_ptr<MatrixDouble> mass_ptr,
      boost::shared_ptr<EntitiesFieldData> data_l2,
      const FieldApproximationBase b, const FieldSpace s,
      const LogManager::SeverityLevel sev = Sev::noisy,
      DataLayout input_layout = DataLayout::GaussByCoeffs);
 
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
 
protected:
  boost::shared_ptr<MatrixDouble> dataPtr;
  boost::shared_ptr<MatrixDouble> coeffsPtr;
  DataLayout inputLayout;
};
 
/**
 * @brief Evaluate right hand side for given data coefficients.
 * 
 */
struct OpDGProjectionEvaluation : public OpDGProjectionCoefficients {
 
  using OutputLayout = DataLayout;
 
  OpDGProjectionEvaluation(boost::shared_ptr<MatrixDouble> data_ptr,
                           boost::shared_ptr<MatrixDouble> coeffs_ptr,
                           boost::shared_ptr<EntitiesFieldData> data_l2,
                           const FieldApproximationBase b, const FieldSpace s,
                           const LogManager::SeverityLevel sev = Sev::noisy,
                           DataLayout output_layout = DataLayout::GaussByCoeffs);
 
  MoFEMErrorCode doWork(int side, EntityType type,
                        EntitiesFieldData::EntData &data);
 
private:
  DataLayout outputLayout;
};
}; // namespace MoFEM
 
#endif // DGPROJECTION_HPP