#include "tutorials/adv-5/src/SeepageOps.hpp"

Inheritance diagram for SeepageOps::OpDomainRhsHydrostaticStress< DIM >:

Collaboration diagram for SeepageOps::OpDomainRhsHydrostaticStress< DIM >:

Public Member Functions
	OpDomainRhsHydrostaticStress (std::string field_name1, boost::shared_ptr< VectorDouble > h_ptr, double specific_weight_water=9.81)

MoFEMErrorCode	iNtegrate (DataForcesAndSourcesCore::EntData &data)

Private Attributes
double	specificWeightWater

boost::shared_ptr< VectorDouble >	hPtr

Detailed Description

template<int DIM>
struct SeepageOps::OpDomainRhsHydrostaticStress< DIM >

Examples: seepage.cpp.

Definition at line 22 of file SeepageOps.hpp.

Constructor & Destructor Documentation

◆ OpDomainRhsHydrostaticStress()

template<int DIM>

SeepageOps::OpDomainRhsHydrostaticStress< DIM >::OpDomainRhsHydrostaticStress	(	std::string	field_name1,
		boost::shared_ptr< VectorDouble >	h_ptr,
		double	specific_weight_water = 9.81 )

inline

Examples: SeepageOps.hpp.

Definition at line 25 of file SeepageOps.hpp.

28 : AssemblyDomainEleOp(field_name1, field_name1, DomainEleOp::OPROW),

29 hPtr(h_ptr), specificWeightWater(specific_weight_water) {}

SeepageOps::OpDomainRhsHydrostaticStress::hPtr

boost::shared_ptr< VectorDouble > hPtr

Definition SeepageOps.hpp:79

SeepageOps::OpDomainRhsHydrostaticStress::specificWeightWater

double specificWeightWater

Definition SeepageOps.hpp:78

AssemblyDomainEleOp

FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

Definition tensor_divergence_operator.cpp:59

Member Function Documentation

◆ iNtegrate()

template<int DIM>

MoFEMErrorCode SeepageOps::OpDomainRhsHydrostaticStress< DIM >::iNtegrate ( DataForcesAndSourcesCore::EntData & data )

inline

Examples: SeepageOps.hpp.

Definition at line 31 of file SeepageOps.hpp.

                                                                  {
    MoFEMFunctionBegin;
 
    const int nb_dofs = data.getIndices().size();
 
    if (nb_dofs) {
      // locRhs.resize(nb_dofs, false);
      // locRhs.clear();
      auto &nf = AssemblyDomainEleOp::locF;
      // get element area
      const double area = getMeasure();
 
      // get number of integration points
      const int nb_integration_points = getGaussPts().size2();
      // get integration weights
      auto t_w = getFTensor0IntegrationWeight();
 
      // get base function
      auto t_base_diff = data.getFTensor1DiffN<DIM>();
 
      FTensor::Index<'i', DIM> i;
 
      auto t_h = getFTensor0FromVec(*hPtr);
      for (int gg = 0; gg != nb_integration_points; gg++) {
        auto t_nf = getFTensor1FromPtr<DIM>(&nf[0]);
 
        const double a = t_w * area * specificWeightWater * t_h;
 
        for (int rr = 0; rr != nb_dofs / DIM; rr++) {
          t_nf(i) -= t_base_diff(i) * a;
 
          // move to the next base function
          ++t_base_diff; // moves the pointer to the next shape function
          ++t_nf;
        }
 
        // move to the weight of the next integration point
        ++t_w;
        ++t_h;
      }
    }
 
    MoFEMFunctionReturn(0);
  }