EshelbianPlasticity::CGGUserPolynomialBase Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianPlasticity.hpp"

Inheritance diagram for EshelbianPlasticity::CGGUserPolynomialBase:

Collaboration diagram for EshelbianPlasticity::CGGUserPolynomialBase:

Public Types
using	CachePhi = boost::tuple< int, int, MatrixDouble >

Public Member Functions
	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 Member Functions
MoFEMErrorCode	getValueHdivForCGGBubble (MatrixDouble &pts)

Private Attributes
MatrixDouble	shapeFun
boost::shared_ptr< CachePhi >	cachePhiPtr

Detailed Description

Definition at line 37 of file EshelbianPlasticity.hpp.

Member Typedef Documentation

CachePhi

using EshelbianPlasticity::CGGUserPolynomialBase::CachePhi = boost::tuple<int, int, MatrixDouble>

Definition at line 39 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

CGGUserPolynomialBase()

EshelbianPlasticity::CGGUserPolynomialBase::CGGUserPolynomialBase

(

boost::shared_ptr< CachePhi >

cache_phi = nullptr

)

Definition at line 2187 of file EshelbianPlasticity.cpp.

    : TetPolynomialBase(), cachePhiPtr(cache_phi_otr) {}

~CGGUserPolynomialBase()

EshelbianPlasticity::CGGUserPolynomialBase::~CGGUserPolynomialBase ( )

default

Member Function Documentation

getValue()

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

Definition at line 2199 of file EshelbianPlasticity.cpp.

                                                                          {
  MoFEMFunctionBeginHot;
 
  this->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");
  }
 
  const auto base = this->cTx->bAse;
  EntitiesFieldData &data = this->cTx->dAta;
 
  switch (this->cTx->sPace) {
  case HDIV:
    CHKERR getValueHdivForCGGBubble(pts);
    break;
  case L2:
    data.dataOnEntities[MBVERTEX][0].getN(base).resize(nb_gauss_pts, 4, false);
    CHKERR Tools::shapeFunMBTET(
        &*data.dataOnEntities[MBVERTEX][0].getN(base).data().begin(),
        &pts(0, 0), &pts(1, 0), &pts(2, 0), nb_gauss_pts);
    data.dataOnEntities[MBVERTEX][0].getDiffN(base).resize(4, 3, false);
    std::copy(Tools::diffShapeFunMBTET.begin(), Tools::diffShapeFunMBTET.end(),
              data.dataOnEntities[MBVERTEX][0].getDiffN(base).data().begin());
    this->cTx->basePolynomialsType0 = Legendre_polynomials;
    CHKERR getValueL2AinsworthBase(pts);
    break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
  }
 
  MoFEMFunctionReturnHot(0);
}

getValueHdivForCGGBubble()

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

private

number of integration points

Definition at line 2242 of file EshelbianPlasticity.cpp.

                                                                 {
  MoFEMFunctionBegin;
 
  // This should be used only in case USER_BASE is selected
  if (this->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 = this->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();
 
  auto check_cache = [this](int order, int nb_gauss_pts) -> bool {
    if (cachePhiPtr) {
      return cachePhiPtr->get<0>() == order &&
             cachePhiPtr->get<1>() == nb_gauss_pts;
    } else {
      return false;
    }
  };
 
  if (check_cache(order, nb_gauss_pts)) {
    auto &mat = cachePhiPtr->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);
    // derivatives 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);
 
    if (cachePhiPtr) {
      cachePhiPtr->get<0>() = order;
      cachePhiPtr->get<1>() = nb_gauss_pts;
      cachePhiPtr->get<2>().resize(phi.size1(), phi.size2(), false);
      noalias(cachePhiPtr->get<2>()) = phi;
    }
  }
 
  MoFEMFunctionReturn(0);
}

query_interface()

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

Definition at line 2191 of file EshelbianPlasticity.cpp.

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

Member Data Documentation

cachePhiPtr

boost::shared_ptr<CachePhi> EshelbianPlasticity::CGGUserPolynomialBase::cachePhiPtr

private

Definition at line 51 of file EshelbianPlasticity.hpp.

shapeFun

MatrixDouble EshelbianPlasticity::CGGUserPolynomialBase::shapeFun

private

Definition at line 50 of file EshelbianPlasticity.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianPlasticity.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp