Inheritance diagram for OpCalculateDeformationGradient< DIM_0, DIM_1 >:

Collaboration diagram for OpCalculateDeformationGradient< DIM_0, DIM_1 >:

Public Member Functions
	OpCalculateDeformationGradient (boost::shared_ptr< MatrixDouble > def_grad_ptr, boost::shared_ptr< MatrixDouble > grad_tensor_ptr)

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

Private Attributes
boost::shared_ptr< MatrixDouble >	gradTensorPtr

boost::shared_ptr< MatrixDouble >	defGradPtr

Detailed Description

template<int DIM_0, int DIM_1>
struct OpCalculateDeformationGradient< DIM_0, DIM_1 >

Examples: dynamic_first_order_con_law.cpp.

Definition at line 338 of file dynamic_first_order_con_law.cpp.

Constructor & Destructor Documentation

◆ OpCalculateDeformationGradient()

template<int DIM_0, int DIM_1>

OpCalculateDeformationGradient< DIM_0, DIM_1 >::OpCalculateDeformationGradient	(	boost::shared_ptr< MatrixDouble >	def_grad_ptr,
		boost::shared_ptr< MatrixDouble >	grad_tensor_ptr )

inline

Examples: dynamic_first_order_con_law.cpp.

Definition at line 340 of file dynamic_first_order_con_law.cpp.

343 : ForcesAndSourcesCore::UserDataOperator(NOSPACE, OPLAST),

344 defGradPtr(def_grad_ptr), gradTensorPtr(grad_tensor_ptr) {}

NOSPACE

@ NOSPACE

Definition definitions.h:83

OpCalculateDeformationGradient::defGradPtr

boost::shared_ptr< MatrixDouble > defGradPtr

Definition dynamic_first_order_con_law.cpp:379

OpCalculateDeformationGradient::gradTensorPtr

boost::shared_ptr< MatrixDouble > gradTensorPtr

Definition dynamic_first_order_con_law.cpp:378

Member Function Documentation

◆ doWork()

template<int DIM_0, int DIM_1>

MoFEMErrorCode OpCalculateDeformationGradient< DIM_0, DIM_1 >::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data )

inline

Examples: dynamic_first_order_con_law.cpp.

Definition at line 346 of file dynamic_first_order_con_law.cpp.

                                                               {
    MoFEMFunctionBegin;
    // Define Indicies
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
 
    // Define Kronecker Delta
    constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
 
    // Number of Gauss points
    const size_t nb_gauss_pts = getGaussPts().size2();
 
    // Resize Piola
    defGradPtr->resize(DIM_0 * DIM_1, nb_gauss_pts, false);
    defGradPtr->clear();
 
    // Extract matrix from data matrix
    auto t_F = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*defGradPtr);
    auto t_H = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*gradTensorPtr);
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
      t_F(i, j) = t_H(i, j) + t_kd(i, j);
 
      ++t_F;
      ++t_H;
    }
 
    MoFEMFunctionReturn(0);
  }