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

Examples: testing_jacobian_of_hook_scaled_with_density_element.cpp.

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);
  }