Inheritance diagram for OpGetDensityField< ALE >:

Collaboration diagram for OpGetDensityField< ALE >:

Public Member Functions
	OpGetDensityField (const std::string row_field, boost::shared_ptr< MatrixDouble > mat_coords_ptr, boost::shared_ptr< VectorDouble > density_at_pts, boost::shared_ptr< MatrixDouble > rho_grad_at_gauss_pts_ptr)

MoFEMErrorCode	doWork (int row_side, EntityType row_type, HookeElement::EntData &row_data)

Public Attributes
boost::shared_ptr< MatrixDouble >	matCoordsPtr

boost::shared_ptr< VectorDouble >	rhoAtGaussPtsPtr

boost::shared_ptr< MatrixDouble >	rhoGradAtGaussPtsPtr

Detailed Description

template<bool ALE>
struct OpGetDensityField< ALE >

Examples: testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 19 of file testing_jacobian_of_hook_scaled_with_density_element.cpp.

Constructor & Destructor Documentation

◆ OpGetDensityField()

template<bool ALE>

OpGetDensityField< ALE >::OpGetDensityField	(	const std::string	row_field,
		boost::shared_ptr< MatrixDouble >	mat_coords_ptr,
		boost::shared_ptr< VectorDouble >	density_at_pts,
		boost::shared_ptr< MatrixDouble >	rho_grad_at_gauss_pts_ptr
	)

inline

Definition at line 24 of file testing_jacobian_of_hook_scaled_with_density_element.cpp.

       : HookeElement::VolUserDataOperator(row_field, OPROW),
         matCoordsPtr(mat_coords_ptr), rhoAtGaussPtsPtr(density_at_pts),
         rhoGradAtGaussPtsPtr(rho_grad_at_gauss_pts_ptr) {}

Member Function Documentation

◆ doWork()

template<bool ALE>

MoFEMErrorCode OpGetDensityField< ALE >::doWork	(	int	row_side,
		EntityType	row_type,
		HookeElement::EntData &	row_data
	)

inline

Definition at line 32 of file testing_jacobian_of_hook_scaled_with_density_element.cpp.

                                                        {
     MoFEMFunctionBegin;
     if (row_type != MBVERTEX)
       MoFEMFunctionReturnHot(0);
     // get number of integration points
     const int nb_integration_pts = getGaussPts().size2();
     rhoAtGaussPtsPtr->resize(nb_integration_pts, false);
     rhoAtGaussPtsPtr->clear();
     rhoGradAtGaussPtsPtr->resize(3, nb_integration_pts, false);
     rhoGradAtGaussPtsPtr->clear();
  
     auto t_rho = getFTensor0FromVec(*rhoAtGaussPtsPtr);
     auto t_grad_rho = getFTensor1FromMat<3>(*rhoGradAtGaussPtsPtr);
  
     MatrixDouble coords;
     if (ALE)
       coords = *matCoordsPtr;
     else
       coords = trans(getCoordsAtGaussPts()); // because the size is (nb_gg,3)
  
     FTensor::Index<'i', 3> i;
     auto t_coords = getFTensor1FromMat<3>(coords);
  
     for (int gg = 0; gg != nb_integration_pts; ++gg) {
       t_rho = 1 + t_coords(i) * t_coords(i); // density is equal to
                                              // 1+x^2+y^2+z^2 (x,y,z coordines)
       t_grad_rho(i) = 2 * t_coords(i);
  
       ++t_rho;
       ++t_coords;
       ++t_grad_rho;
     }
  
     MoFEMFunctionReturn(0);
   }