CGGUserPolynomialBase.hpp

Go to the documentation of this file.

/**
 * @brief CGG User Polynomial Base
 *
 */
 
struct CGGUserPolynomialBase : public TetPolynomialBase {
 
  using CachePhi = boost::tuple<int, int, MatrixDouble, MatrixDouble>;
 
  CGGUserPolynomialBase(boost::shared_ptr<CachePhi> cache_phi = nullptr);
  ~CGGUserPolynomialBase() = default;
 
  MoFEMErrorCode query_interface(boost::typeindex::type_index type_index,
                                 BaseFunctionUnknownInterface **iface) const;
  MoFEMErrorCode getValue(MatrixDouble &pts,
                          boost::shared_ptr<BaseFunctionCtx> ctx_ptr);
 
private:
  MatrixDouble shapeFun;
  boost::shared_ptr<CachePhi> cachePhiPtr;
 
  MoFEMErrorCode getValueHdivForCGGBubble(MatrixDouble &pts);
};