#include <users_modules/basic_finite_elements/mix_transport/src/UnsaturatedFlow.hpp>

Inheritance diagram for MixTransport::UnsaturatedFlowElement::OpIntegrateFluxes:

Collaboration diagram for MixTransport::UnsaturatedFlowElement::OpIntegrateFluxes:

Public Member Functions
	OpIntegrateFluxes (UnsaturatedFlowElement &ctx, const std::string fluxes_name)
	Constructor. More...

MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)
	Integrate boundary condition. More...

Public Attributes
UnsaturatedFlowElement &	cTx

FTensor::Index< 'i', 3 >	i

Detailed Description

Definition at line 901 of file UnsaturatedFlow.hpp.

Constructor & Destructor Documentation

◆ OpIntegrateFluxes()

MixTransport::UnsaturatedFlowElement::OpIntegrateFluxes::OpIntegrateFluxes	(	UnsaturatedFlowElement &	ctx,
		const std::string	fluxes_name
	)

Constructor.

Examples: UnsaturatedFlow.hpp.

Definition at line 908 of file UnsaturatedFlow.hpp.

         : MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator(
               fluxes_name, UserDataOperator::OPROW),
           cTx(ctx) {}

Member Function Documentation

◆ doWork()

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

Integrate boundary condition.

Parameters

side	local index of entity
type	type of entity
data	data on entity

Returns: error code

Examples: UnsaturatedFlow.hpp.

Definition at line 923 of file UnsaturatedFlow.hpp.

                                                                  {
       MoFEMFunctionBegin;
       int nb_dofs = data.getFieldData().size();
       if (nb_dofs == 0)
         MoFEMFunctionReturnHot(0);
       // Get base function
       auto t_n_hdiv = data.getFTensor1N<3>();
       // get normal of face
       auto t_normal = getFTensor1Normal();
       // Integration weight
       auto t_w = getFTensor0IntegrationWeight();
       double flux_on_entity = 0;
       int nb_gauss_pts = data.getN().size1();
       for (int gg = 0; gg < nb_gauss_pts; gg++) {
         auto t_data = data.getFTensor0FieldData();
         for (int rr = 0; rr != nb_dofs; rr++) {
           flux_on_entity -= (0.5 * t_data * t_w) * (t_n_hdiv(i) * t_normal(i));
           ++t_n_hdiv;
           ++t_data;
         }
         ++t_w;
       }
       CHKERR VecSetValue(cTx.ghostFlux, 0, flux_on_entity, ADD_VALUES);
       MoFEMFunctionReturn(0);
     }