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: mofem/users_modules/basic_finite_elements/atom_tests/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: mofem/users_modules/basic_finite_elements/atom_tests/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, -1, CoeffsByGauss>(*rhoGradAtGaussPtsPtr);
 
    FTensor::Index<'i', 3> i;
    auto set_density = [&](auto t_coords) {
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        t_rho = 1 + t_coords(i) * t_coords(i); // 1+x^2+y^2+z^2
        t_grad_rho(i) = 2 * t_coords(i);
 
        ++t_rho;
        ++t_coords;
        ++t_grad_rho;
      }
    };
 
    if constexpr (ALE) {
      set_density(getFTensor1FromMat<3>(*matCoordsPtr));
    } else {
      MatrixDouble coords =
          trans(getCoordsAtGaussPts()); // because the size is (nb_gg,3)
      set_density(getFTensor1FromMat<3, -1, CoeffsByGauss>(coords));
    }
 
    MoFEMFunctionReturn(0);
  }