Inheritance diagram for OpURhs:

Collaboration diagram for OpURhs:

Public Member Functions
	OpURhs (const std::string field_name, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > u_dot_ptr, boost::shared_ptr< MatrixDouble > grad_u_ptr, boost::shared_ptr< MatrixDouble > grad_h_ptr)

MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< MatrixDouble >	uPtr

boost::shared_ptr< MatrixDouble >	uDotPtr

boost::shared_ptr< MatrixDouble >	uGradPtr

boost::shared_ptr< MatrixDouble >	hGradPtr

Detailed Description

Examples: mofem/tutorials/vec-4/shallow_wave.cpp.

Definition at line 82 of file shallow_wave.cpp.

Constructor & Destructor Documentation

◆ OpURhs()

OpURhs::OpURhs	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	u_ptr,
		boost::shared_ptr< MatrixDouble >	u_dot_ptr,
		boost::shared_ptr< MatrixDouble >	grad_u_ptr,
		boost::shared_ptr< MatrixDouble >	grad_h_ptr
	)

inline

Definition at line 84 of file shallow_wave.cpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        uPtr(u_ptr), uDotPtr(u_dot_ptr), uGradPtr(grad_u_ptr),
        hGradPtr(grad_h_ptr) {}

Member Function Documentation

◆ iNtegrate()

MoFEMErrorCode OpURhs::iNtegrate ( EntitiesFieldData::EntData & row_data )

inline

Examples: mofem/tutorials/vec-4/shallow_wave.cpp.

Definition at line 92 of file shallow_wave.cpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_row_base = row_data.getFTensor0N();
    auto t_row_diff_base = row_data.getFTensor1DiffN<3>();
    auto t_dot_u = getFTensor1FromMat<3>(*uDotPtr);
    auto t_u = getFTensor1FromMat<3>(*uPtr);
    auto t_grad_u = getFTensor2FromMat<3, 3>(*uGradPtr);
    auto t_grad_h = getFTensor1FromMat<3>(*hGradPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
    auto t_normal = getFTensor1NormalsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double alpha = t_w * vol;
      auto t_nf = getFTensor1FromArray<3, 3>(locF);
 
      FTensor::Tensor1<double, 3> t_rhs;
      FTensor::Tensor2<double, 3, 3> t_rhs_grad;
 
      const auto a = std::sqrt(t_coords(i) * t_coords(i));
      const auto sin_fi = t_coords(2) / a;
      const auto f = 2 * omega * sin_fi;
 
      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);
 
      FTensor::Tensor2<double, 3, 3> t_A;
      t_A(i, m) = levi_civita(i, j, m) * t_r(j);
 
      t_rhs(m) = t_Q(m, i) * (t_dot_u(i) + t_grad_u(i, j) * t_u(j) +
                              f * t_A(i, j) * t_u(j) + g * t_grad_h(i));
      t_rhs_grad(m, j) = t_Q(m, i) * (mu * t_grad_u(i, j));
 
      t_rhs(m) += t_P(m, j) * t_u(j);
 
      int rr = 0;
      for (; rr != nbRows / 3; ++rr) {
        t_nf(i) += alpha * t_row_base * t_rhs(i);
        t_nf(i) += alpha * t_row_diff_base(j) * t_rhs_grad(i, j);
        ++t_row_base;
        ++t_row_diff_base;
        ++t_nf;
      }
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      ++t_w;
      ++t_u;
      ++t_dot_u;
      ++t_grad_u;
      ++t_grad_h;
      ++t_coords;
      ++t_normal;
    }
 
    MoFEMFunctionReturn(0);
  }