FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts Struct Reference

Evaluate density at integration points. More...

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

Inheritance diagram for FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts:

Collaboration diagram for FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts:

Public Member Functions
	OpMWLSRhoAtGaussPts (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

Evaluate density at integration points.

Definition at line 330 of file MWLS.hpp.

Constructor & Destructor Documentation

OpMWLSRhoAtGaussPts()

FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts::OpMWLSRhoAtGaussPts ( 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 336 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::OpMWLSRhoAtGaussPts::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 1383 of file MWLS.cpp.

                                                                      {
  MoFEMFunctionBegin;
  if (testing)
    cerr << "Element " << getFEEntityHandle() << endl;
 
  FTensor::Index<'i', 3> i;
 
  const int nb_gauss_pts = data.getN().size1();
  const bool use_global_base =
      mwlsApprox->getUseGlobalBaseAtMaterialReferenceConfiguration();
  Tag th;
  CHKERR mwlsApprox->mwlsMoab.tag_get_handle(rhoTagName.c_str(), th);
 
  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);
    FTensor::Tensor1<double, 3> t_approx_base_pos;
    t_approx_base_pos(i) = t_approx_base_point(i);
 
    if (testing) {
      cerr << "material_positions " << t_mwls_material_positions << endl;
      cerr << "approx_base_point " << t_approx_base_point << endl;
    }
    CHKERR mwlsApprox->getInfluenceNodes(&t_approx_base_pos(0));
    CHKERR mwlsApprox->calculateApproxCoeff(true, use_global_base);
    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(ii + 3 * jj, gg) = diff_diff_vals(jj, ii);
      }
    }
 
    ++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::OpMWLSRhoAtGaussPts::calculate2ndDerivative

Definition at line 334 of file MWLS.hpp.

calculateDerivative

const bool FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts::calculateDerivative

Definition at line 333 of file MWLS.hpp.

rhoTagName

std::string FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts::rhoTagName

Definition at line 331 of file MWLS.hpp.

testing

const bool FractureMechanics::MWLSApprox::OpMWLSRhoAtGaussPts::testing

Definition at line 335 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