Inheritance diagram for OpLhs:

Collaboration diagram for OpLhs:

Public Member Functions
	OpLhs (const std::string field_name, boost::shared_ptr< MatrixDouble > x_ptr, boost::shared_ptr< MatrixDouble > dot_x_ptr)

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

Private Attributes
boost::shared_ptr< MatrixDouble >	xPtr

boost::shared_ptr< MatrixDouble >	xDotPtr

Detailed Description

Examples: mofem/tutorials/vec-4/approx_sphere.cpp, and mofem/users_modules/adolc-plasticity/src/ADOLCPlasticity.hpp.

Definition at line 134 of file approx_sphere.cpp.

Constructor & Destructor Documentation

◆ OpLhs()

OpLhs::OpLhs	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	x_ptr,
		boost::shared_ptr< MatrixDouble >	dot_x_ptr
	)

inline

Definition at line 136 of file approx_sphere.cpp.

      : AssemblyDomainEleOp(field_name, field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        xPtr(x_ptr), xDotPtr(dot_x_ptr) {
    this->sYmm = false;
  }

Member Function Documentation

◆ iNtegrate()

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

inline

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

Definition at line 144 of file approx_sphere.cpp.

                                                               {
    MoFEMFunctionBegin;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_row_base = row_data.getFTensor0N();
 
    auto t_x0 = getFTensor1CoordsAtGaussPts();
    auto t_x = getFTensor1FromMat<3>(*xPtr);
    auto t_dot_x = getFTensor1FromMat<3>(*xDotPtr);
    auto t_normal = getFTensor1NormalsAtGaussPts();
 
    auto get_t_mat = [&](const int rr) {
      return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3>{
          &locMat(rr + 0, 0), &locMat(rr + 0, 1), &locMat(rr + 0, 2),
 
          &locMat(rr + 1, 0), &locMat(rr + 1, 1), &locMat(rr + 1, 2),
 
          &locMat(rr + 2, 0), &locMat(rr + 2, 1), &locMat(rr + 2, 2)};
    };
 
    constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
    const double ts_a = getTSa();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
      FTensor::Tensor1<double, 3> t_n{t_x0(0), t_x0(1), t_x0(2)};
      t_n.normalize();
      FTensor::Tensor2<double, 3, 3> t_P, t_Q;
      t_P(i, j) = t_n(i) * t_n(j);
      t_Q(i, j) = t_kd(i, j) - t_P(i, j);
 
      auto t_J_lhs = lhs_J_dx2(t_x(0), t_x(1), t_x(2));
      double l = std::sqrt(t_normal(i) * t_normal(i));
 
      const double alpha = t_w;
      FTensor::Tensor2<double, 3, 3> t_lhs;
      t_lhs(i, j) =
          (alpha * l) * (t_P(i, k) * t_J_lhs(k, j) + t_Q(i, j) * ts_a);
 
      int rr = 0;
      for (; rr != nbRows / 3; rr++) {
 
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        auto t_mat = get_t_mat(3 * rr);
 
        for (int cc = 0; cc != nbCols / 3; cc++) {
 
          const double rc = t_row_base * t_col_base;
          t_mat(i, j) += rc * t_lhs(i, j);
 
          ++t_col_base;
          ++t_mat;
        }
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr)
        ++t_row_base;
 
      ++t_w;
      ++t_x;
      ++t_x0;
      ++t_normal;
    }
 
    MoFEMFunctionReturn(0);
  }