#include <users_modules/fracture_mechanics/src/MWLS.hpp>

Inheritance diagram for FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs:

Collaboration diagram for FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs:

Public Member Functions
	OpMWLSRhoAtGaussUsingPrecalulatedCoeffs (boost::shared_ptr< MatrixDouble > mat_pos_at_pts_ptr, boost::shared_ptr< MatrixDouble > mat_grad_pos_at_pts_ptr, boost::shared_ptr< CrackFrontElement > fe_singular_ptr, boost::shared_ptr< MWLSApprox > mwls_approx, const std::string rho_tag_name, const bool calculate_derivative, const bool calculate_2nd_derivative, bool testing=false)

Public Member Functions inherited from FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >
	OpMWLSBase (boost::shared_ptr< MatrixDouble > mat_pos_at_pts_ptr, boost::shared_ptr< MatrixDouble > mat_grad_pos_at_pts_ptr, boost::shared_ptr< CrackFrontElement > fe_singular_ptr, boost::shared_ptr< MWLSApprox > mwls_approx, bool testing=false)

virtual	~OpMWLSBase ()=default

MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
std::string	rhoTagName

const bool	calculateDerivative

const bool	calculate2ndDerivative

const bool	testing

Public Attributes inherited from FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >
boost::shared_ptr< MatrixDouble >	matPosAtPtsPtr

boost::shared_ptr< MatrixDouble >	matGradPosAtPtsPtr

boost::weak_ptr< CrackFrontElement >	feSingularPtr

boost::shared_ptr< MWLSApprox >	mwlsApprox

Protected Member Functions
MoFEMErrorCode	doMWLSWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)
	Do specific work in operator. More...

Detailed Description

Definition at line 302 of file MWLS.hpp.

Constructor & Destructor Documentation

◆ OpMWLSRhoAtGaussUsingPrecalulatedCoeffs()

FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs	(	boost::shared_ptr< MatrixDouble >	mat_pos_at_pts_ptr,
		boost::shared_ptr< MatrixDouble >	mat_grad_pos_at_pts_ptr,
		boost::shared_ptr< CrackFrontElement >	fe_singular_ptr,
		boost::shared_ptr< MWLSApprox >	mwls_approx,
		const std::string	rho_tag_name,
		const bool	calculate_derivative,
		const bool	calculate_2nd_derivative,
		bool	testing = `false`
	)

inline

Definition at line 309 of file MWLS.hpp.

         : OpMWLSBase<VolumeElementForcesAndSourcesCore>(
               mat_pos_at_pts_ptr, mat_grad_pos_at_pts_ptr, fe_singular_ptr,
               mwls_approx, testing),
           rhoTagName(rho_tag_name), calculateDerivative(calculate_derivative),
           calculate2ndDerivative(calculate_2nd_derivative), testing(testing) {}

Member Function Documentation

◆ doMWLSWork()

MoFEMErrorCode FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::doMWLSWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

protectedvirtual

Do specific work in operator.

Is protected to prevent direct call. This can be only called by doWork method.

Parameters

side
type
data

Returns: MoFEMErrorCode

Implements FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >.

Definition at line 1293 of file MWLS.cpp.

                                                                       {
   MoFEMFunctionBegin;
  
   if (testing)
     cerr << "Element " << getFEEntityHandle() << endl;
   const bool use_global_base =
       mwlsApprox->getUseGlobalBaseAtMaterialReferenceConfiguration();
   if (use_global_base)
     SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
             "Coefficients can be precalculated only for local base in "
             "reference configuration");
  
   FTensor::Index<'i', 3> i;
  
   Tag th;
   CHKERR mwlsApprox->mwlsMoab.tag_get_handle(rhoTagName.c_str(), th);
  
   const int nb_gauss_pts = data.getN().size1();
   auto &inv_AB_map = mwlsApprox->invABMap.at(getFEEntityHandle());
   auto &influence_nodes_map =
       mwlsApprox->influenceNodesMap.at(getFEEntityHandle());
   auto &dm_nodes_map = mwlsApprox->dmNodesMap.at(this->getFEEntityHandle());
  
   VectorDouble &rho = *mwlsApprox->rhoAtGaussPts;
   rho.resize(nb_gauss_pts, false);
   MatrixDouble &diff_rho = *mwlsApprox->diffRhoAtGaussPts;
   diff_rho.resize(3, nb_gauss_pts, false);
   MatrixDouble &diff_diff_rho = *mwlsApprox->diffDiffRhoAtGaussPts;
   if (calculate2ndDerivative)
     diff_diff_rho.resize(9, nb_gauss_pts, false);
  
   auto t_mwls_material_positions =
       getFTensor1FromMat<3>(mwlsApprox->mwlsMaterialPositions);
   auto t_approx_base_point = getFTensor1FromMat<3>(mwlsApprox->approxBasePoint);
  
   for (int gg = 0; gg != nb_gauss_pts; ++gg) {
  
     FTensor::Tensor1<double, 3> t_mat_pos;
     t_mat_pos(i) = t_mwls_material_positions(i);
  
     if (testing) {
       cerr << "material_positions " << t_mwls_material_positions << endl;
     }
  
     for (int dd : {0, 1, 2})
       mwlsApprox->approxPointCoords[dd] = t_approx_base_point(dd);
  
     mwlsApprox->influenceNodes = influence_nodes_map[gg];
     mwlsApprox->invAB = inv_AB_map[gg];
     mwlsApprox->shortenDm = dm_nodes_map[gg];
  
     CHKERR mwlsApprox->evaluateApproxFun(&t_mat_pos(0));
     CHKERR mwlsApprox->evaluateFirstDiffApproxFun(&t_mat_pos(0),
                                                   use_global_base);
     if (calculate2ndDerivative)
       CHKERR mwlsApprox->evaluateSecondDiffApproxFun(&t_mat_pos(0),
                                                      use_global_base);
  
     CHKERR mwlsApprox->getTagData(th);
     const auto &vals = mwlsApprox->getDataApprox();
  
     rho(gg) = vals[0];
     if (calculateDerivative) {
       const auto &diff_vals = mwlsApprox->getDiffDataApprox();
       for (int ii = 0; ii != 3; ++ii) {
         diff_rho(ii, gg) = diff_vals(ii, 0);
       }
     }
  
     if (calculate2ndDerivative) {
       const auto &diff_diff_vals = mwlsApprox->getDiffDiffDataApprox();
       for (int ii = 0; ii != 3; ++ii) {
         for (int jj = 0; jj != 3; ++jj)
           diff_diff_rho(jj * 3 + ii, gg) = diff_diff_vals(jj * 3 + ii, 0);
       }
     }
  
     ++t_mwls_material_positions;
     ++t_approx_base_point;
   }
  
   if (testing) {
     cerr << "rho " << rho << endl;
     cerr << endl;
   }
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ calculate2ndDerivative

const bool FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::calculate2ndDerivative

Definition at line 307 of file MWLS.hpp.

◆ calculateDerivative

const bool FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::calculateDerivative

Definition at line 306 of file MWLS.hpp.

◆ rhoTagName

std::string FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::rhoTagName

Definition at line 305 of file MWLS.hpp.

◆ testing

const bool FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussUsingPrecalulatedCoeffs::testing

Definition at line 308 of file MWLS.hpp.

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

users_modules/fracture_mechanics/src/MWLS.hpp
users_modules/fracture_mechanics/src/impl/MWLS.cpp

Public Member Functions

Public Attributes

Protected Member Functions