EshelbianPlasticity::CGGUserPolynomialBase Struct Reference

Inheritance diagram for EshelbianPlasticity::CGGUserPolynomialBase:

Collaboration diagram for EshelbianPlasticity::CGGUserPolynomialBase:

Public Member Functions
	CGGUserPolynomialBase ()
	~CGGUserPolynomialBase ()
MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, BaseFunctionUnknownInterface **iface) const
MoFEMErrorCode	getValue (MatrixDouble &pts, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
Public Member Functions inherited from MoFEM::BaseFunction
MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, MoFEM::UnknownInterface **iface) const
virtual MoFEMErrorCode	getValue (MatrixDouble &pts, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
virtual MoFEMErrorCode	getValue (MatrixDouble &pts_x, MatrixDouble &pts_t, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
Public Member Functions inherited from MoFEM::BaseFunctionUnknownInterface
virtual	~BaseFunctionUnknownInterface ()=default
Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (bool error_if_registration_failed=true)
	Register interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface refernce to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface. More...
template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface. More...
virtual	~UnknownInterface ()=default

Private Member Functions
MoFEMErrorCode	getValueHdivForCGGBubble (MatrixDouble &pts)

Private Attributes
EntPolynomialBaseCtx *	cTx
MatrixDouble	shapeFun
boost::tuple< int, int, MatrixDouble >	cachePhi = {0, 0, MatrixDouble()}

Additional Inherited Members
Static Public Member Functions inherited from MoFEM::UnknownInterface
static MoFEMErrorCode	getLibVersion (Version &version)
	Get library version. More...
static MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version)
	Get database major version. More...
static MoFEMErrorCode	setFileVersion (moab::Interface &moab, Version version=Version(MoFEM_VERSION_MAJOR, MoFEM_VERSION_MINOR, MoFEM_VERSION_BUILD))
	Get database major version. More...
static MoFEMErrorCode	getInterfaceVersion (Version &version)
	Get database major version. More...

Detailed Description

Definition at line 752 of file EshelbianPlasticity.cpp.

Constructor & Destructor Documentation

CGGUserPolynomialBase()

EshelbianPlasticity::CGGUserPolynomialBase::CGGUserPolynomialBase ( )

inline

Definition at line 754 of file EshelbianPlasticity.cpp.

{}

~CGGUserPolynomialBase()

EshelbianPlasticity::CGGUserPolynomialBase::~CGGUserPolynomialBase ( )

inline

Definition at line 755 of file EshelbianPlasticity.cpp.

{}

Member Function Documentation

getValue()

MoFEMErrorCode EshelbianPlasticity::CGGUserPolynomialBase::getValue ( MatrixDouble &  pts, boost::shared_ptr< BaseFunctionCtx >  ctx_ptr  )

inline

Definition at line 763 of file EshelbianPlasticity.cpp.

                                                                    {
    MoFEMFunctionBeginHot;
 
    cTx = ctx_ptr->getInterface<EntPolynomialBaseCtx>();
 
    int nb_gauss_pts = pts.size2();
    if (!nb_gauss_pts) {
      MoFEMFunctionReturnHot(0);
    }
 
    if (pts.size1() < 3) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Wrong dimension of pts, should be at least 3 rows with "
              "coordinates");
    }
 
    switch (cTx->sPace) {
    case HDIV:
      CHKERR getValueHdivForCGGBubble(pts);
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
    }
 
    MoFEMFunctionReturnHot(0);
  }

getValueHdivForCGGBubble()

MoFEMErrorCode EshelbianPlasticity::CGGUserPolynomialBase::getValueHdivForCGGBubble ( MatrixDouble &  pts )

inlineprivate

number of integration points

Definition at line 797 of file EshelbianPlasticity.cpp.

                                                             {
    MoFEMFunctionBegin;
 
    // This should be used only in case USER_BASE is selected
    if (cTx->bAse != USER_BASE) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Wrong base, should be USER_BASE");
    }
    // get access to data structures on element
    EntitiesFieldData &data = cTx->dAta;
    // Get approximation order on element. Note that bubble functions are only
    // on tetrahedron.
    const int order = data.dataOnEntities[MBTET][0].getOrder();
    /// number of integration points
    const int nb_gauss_pts = pts.size2();
 
    if(cachePhi.get<0>() == order && cachePhi.get<1>() == nb_gauss_pts) {
      auto &mat = cachePhi.get<2>();
      auto &phi = data.dataOnEntities[MBTET][0].getN(USER_BASE);
      phi.resize(mat.size1(), mat.size2(), false);
      noalias(phi) = mat;
    } else {
      // calculate shape functions, i.e. barycentric coordinates
      shapeFun.resize(nb_gauss_pts, 4, false);
      CHKERR ShapeMBTET(&*shapeFun.data().begin(), &pts(0, 0), &pts(1, 0),
                        &pts(2, 0), nb_gauss_pts);
      // direvatives of shape functions
      double diff_shape_fun[12];
      CHKERR ShapeDiffMBTET(diff_shape_fun);
 
      const int nb_base_functions = NBVOLUMETET_CCG_BUBBLE(order);
      // get base functions and set size
      MatrixDouble &phi = data.dataOnEntities[MBTET][0].getN(USER_BASE);
      phi.resize(nb_gauss_pts, 9 * nb_base_functions, false);
      // finally calculate base functions
      FTensor::Tensor2<FTensor::PackPtr<double *, 9>, 3, 3> t_phi(
          &phi(0, 0), &phi(0, 1), &phi(0, 2),
 
          &phi(0, 3), &phi(0, 4), &phi(0, 5),
 
          &phi(0, 6), &phi(0, 7), &phi(0, 8));
      CHKERR CGG_BubbleBase_MBTET(order, &shapeFun(0, 0), diff_shape_fun, t_phi,
                                  nb_gauss_pts);
 
      cachePhi.get<0>() = order;
      cachePhi.get<1>() = nb_gauss_pts;
      cachePhi.get<2>().resize(phi.size1(), phi.size2(), false);
      noalias(cachePhi.get<2>()) = phi;
    }
 
    MoFEMFunctionReturn(0);
  }

query_interface()

MoFEMErrorCode EshelbianPlasticity::CGGUserPolynomialBase::query_interface ( boost::typeindex::type_index  type_index, BaseFunctionUnknownInterface **  iface  ) const

inline

Definition at line 757 of file EshelbianPlasticity.cpp.

                                                                             {
    *iface = const_cast<CGGUserPolynomialBase *>(this);
    return 0;
  }

Member Data Documentation

cachePhi

boost::tuple<int, int, MatrixDouble> EshelbianPlasticity::CGGUserPolynomialBase::cachePhi = {0, 0, MatrixDouble()}

private

Definition at line 795 of file EshelbianPlasticity.cpp.

cTx

EntPolynomialBaseCtx* EshelbianPlasticity::CGGUserPolynomialBase::cTx

private

Definition at line 792 of file EshelbianPlasticity.cpp.

shapeFun

MatrixDouble EshelbianPlasticity::CGGUserPolynomialBase::shapeFun

private

Definition at line 794 of file EshelbianPlasticity.cpp.

---

The documentation for this struct was generated from the following file:

• users_modules/eshelbian_plasticit/src/impl/EshelbianPlasticity.cpp