OpULhs_dH Struct Reference

Inheritance diagram for OpULhs_dH:

Collaboration diagram for OpULhs_dH:

Public Member Functions
	OpULhs_dH (const std::string field_name_row, const std::string field_name_col)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Detailed Description

Examples

shallow_wave.cpp.

Definition at line 267 of file shallow_wave.cpp.

Constructor & Destructor Documentation

OpULhs_dH()

OpULhs_dH::OpULhs_dH ( const std::string field_name_row, const std::string field_name_col )

inline

Examples

shallow_wave.cpp.

Definition at line 269 of file shallow_wave.cpp.

      : AssemblyDomainEleOp(field_name_row, field_name_col,
                            AssemblyDomainEleOp::OPROWCOL) {
    this->sYmm = false;
  }

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpULhs_dH::iNtegrate ( EntitiesFieldData::EntData & row_data, EntitiesFieldData::EntData & col_data )

inline

Examples

shallow_wave.cpp.

Definition at line 275 of file shallow_wave.cpp.

                                                               {
    MoFEMFunctionBegin;
 
    // get element volume
    const double vol = getMeasure();
    // get integration weights
    auto t_w = getFTensor0IntegrationWeight();
    // get base function gradient on rows
    auto t_row_base = row_data.getFTensor0N();
    // normal
    auto t_normal = getFTensor1NormalsAtGaussPts();
 
    auto get_t_vec = [&](const int rr) {
      return FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 3>{
          &locMat(rr + 0, 0),
 
          &locMat(rr + 1, 0),
 
          &locMat(rr + 2, 0)};
    };
 
    // loop over integration points
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      FTensor::Tensor1<double, 3> t_r;
      t_r(i) = t_normal(i);
      t_r.normalize();
 
      constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
      FTensor::Tensor2<double, 3, 3> t_P, t_Q;
      t_P(i, j) = t_r(i) * t_r(j);
      t_Q(i, j) = t_kd(i, j) - t_P(i, j);
 
      const double alpha = t_w * vol;
 
      int rr = 0;
      for (; rr != nbRows / 3; rr++) {
        auto t_vec = get_t_vec(3 * rr);
        auto t_col_diff_base = col_data.getFTensor1DiffN<3>(gg, 0);
        const double a = alpha * g * t_row_base;
 
        for (int cc = 0; cc != nbCols; cc++) {
          t_vec(i) += a * (t_Q(i, m) * t_col_diff_base(m));
          ++t_vec;
          ++t_col_diff_base;
        }
 
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr)
        ++t_row_base;
 
      ++t_w;
      ++t_normal;
    }
 
    MoFEMFunctionReturn(0);
  }

---

The documentation for this struct was generated from the following file:

• tutorials/vec-4/shallow_wave.cpp