Class used to calculate base functions at integration points. More...

Inheritance diagram for SomeUserPolynomialBase:

Collaboration diagram for SomeUserPolynomialBase:

Public Member Functions
	SomeUserPolynomialBase ()=default

	~SomeUserPolynomialBase ()=default

MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const
	Return interface to this class when one ask for for tetrahedron, otherisw return interface class for generic class. More...

MoFEMErrorCode	getValue (MatrixDouble &pts, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
	Calculate base functions at intergeneration points. More...

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 MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const =0

virtual	~BaseFunctionUnknownInterface ()=default

Public Member Functions inherited from MoFEM::UnknownInterface
virtual MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const =0

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

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

Class used to calculate base functions at integration points.

Examples: forces_and_sources_testing_users_base.cpp.

Definition at line 33 of file forces_and_sources_testing_users_base.cpp.

Constructor & Destructor Documentation

◆ SomeUserPolynomialBase()

SomeUserPolynomialBase::SomeUserPolynomialBase ( )

default

Examples: forces_and_sources_testing_users_base.cpp.

◆ ~SomeUserPolynomialBase()

SomeUserPolynomialBase::~SomeUserPolynomialBase ( )

default

Examples: forces_and_sources_testing_users_base.cpp.

Member Function Documentation

◆ getValue()

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

Calculate base functions at intergeneration points.

Parameters

pts
ctx_ptr

Returns: MoFEMErrorCode

Examples: forces_and_sources_testing_users_base.cpp.

Definition at line 59 of file forces_and_sources_testing_users_base.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 SomeUserPolynomialBase::getValueHdivForCGGBubble ( MatrixDouble & pts )

private

Examples: forces_and_sources_testing_users_base.cpp.

Definition at line 92 of file forces_and_sources_testing_users_base.cpp.

                                                             {
    MoFEMFunctionBegin;
 
    const FieldApproximationBase base = cTx->bAse;
    // 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");
    }
 
    // This is example, simply use Demkowicz HDiv base to generate base
    // functions
 
    EntitiesFieldData &data = cTx->dAta;
    int nb_gauss_pts = pts.size2();
 
    // 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);
 
    int volume_order = data.dataOnEntities[MBTET][0].getOrder();
 
    int p_f[4];
    double *phi_f[4];
    double *diff_phi_f[4];
 
    // Calculate base function on tet faces
    for (int ff = 0; ff != 4; ff++) {
      int face_order = data.dataOnEntities[MBTRI][ff].getOrder();
      int order = volume_order > face_order ? volume_order : face_order;
      data.dataOnEntities[MBTRI][ff].getN(base).resize(
          nb_gauss_pts, 3 * NBFACETRI_DEMKOWICZ_HDIV(order), false);
      data.dataOnEntities[MBTRI][ff].getDiffN(base).resize(
          nb_gauss_pts, 9 * NBFACETRI_DEMKOWICZ_HDIV(order), false);
      p_f[ff] = order;
      phi_f[ff] = &*data.dataOnEntities[MBTRI][ff].getN(base).data().begin();
      diff_phi_f[ff] =
          &*data.dataOnEntities[MBTRI][ff].getDiffN(base).data().begin();
      if (NBFACETRI_DEMKOWICZ_HDIV(order) == 0)
        continue;
      CHKERR Hdiv_Demkowicz_Face_MBTET_ON_FACE(
          &data.facesNodes(ff, 0), order, &*shapeFun.data().begin(),
          diff_shape_fun, phi_f[ff], diff_phi_f[ff], nb_gauss_pts, 4);
    }
 
    // Calculate base functions in tet interior
    if (NBVOLUMETET_DEMKOWICZ_HDIV(volume_order) > 0) {
      data.dataOnEntities[MBTET][0].getN(base).resize(
          nb_gauss_pts, 3 * NBVOLUMETET_DEMKOWICZ_HDIV(volume_order), false);
      data.dataOnEntities[MBTET][0].getDiffN(base).resize(
          nb_gauss_pts, 9 * NBVOLUMETET_DEMKOWICZ_HDIV(volume_order), false);
      double *phi_v = &*data.dataOnEntities[MBTET][0].getN(base).data().begin();
      double *diff_phi_v =
          &*data.dataOnEntities[MBTET][0].getDiffN(base).data().begin();
      CHKERR Hdiv_Demkowicz_Interior_MBTET(
          volume_order, &*shapeFun.data().begin(), diff_shape_fun, p_f, phi_f,
          diff_phi_f, phi_v, diff_phi_v, nb_gauss_pts);
    }
 
    // Set size of face base correctly
    for (int ff = 0; ff != 4; ff++) {
      int face_order = data.dataOnEntities[MBTRI][ff].getOrder();
      data.dataOnEntities[MBTRI][ff].getN(base).resize(
          nb_gauss_pts, 3 * NBFACETRI_DEMKOWICZ_HDIV(face_order), true);
      data.dataOnEntities[MBTRI][ff].getDiffN(base).resize(
          nb_gauss_pts, 9 * NBFACETRI_DEMKOWICZ_HDIV(face_order), true);
    }
 
    MoFEMFunctionReturn(0);
  }

◆ query_interface()

MoFEMErrorCode SomeUserPolynomialBase::query_interface	(	boost::typeindex::type_index	type_index,
		UnknownInterface **	iface
	)		const

virtual

Return interface to this class when one ask for for tetrahedron, otherisw return interface class for generic class.

Parameters

iface interface class

Returns: MoFEMErrorCode

Implements MoFEM::BaseFunctionUnknownInterface.

Examples: forces_and_sources_testing_users_base.cpp.

Definition at line 45 of file forces_and_sources_testing_users_base.cpp.

                                                                 {
    MoFEMFunctionBegin;
    *iface = const_cast<SomeUserPolynomialBase *>(this);
    MoFEMFunctionReturn(0);
  }