MatOps::OpAxisymmetricLhs< A, DomainEleOp > Struct Template Reference

#include "src/materials/MatAxisymmetric.hpp"

Inheritance diagram for MatOps::OpAxisymmetricLhs< A, DomainEleOp >:

Collaboration diagram for MatOps::OpAxisymmetricLhs< A, DomainEleOp >:

Public Types
using	OpBase = OpAxisymmetricBase< A, DomainEleOp >

Public Member Functions
	OpAxisymmetricLhs (const std::string field_name, boost::shared_ptr< MatrixDouble > full_tangent_ptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data) override

Protected Attributes
boost::shared_ptr< MatrixDouble >	fullTangentPtr

Detailed Description

template<AssemblyType A, typename DomainEleOp>
struct MatOps::OpAxisymmetricLhs< A, DomainEleOp >

Definition at line 177 of file MatAxisymmetric.hpp.

Member Typedef Documentation

OpBase

template<AssemblyType A, typename DomainEleOp >

using MatOps::OpAxisymmetricLhs< A, DomainEleOp >::OpBase = OpAxisymmetricBase<A, DomainEleOp>

Definition at line 178 of file MatAxisymmetric.hpp.

Constructor & Destructor Documentation

OpAxisymmetricLhs()

template<AssemblyType A, typename DomainEleOp >

MatOps::OpAxisymmetricLhs< A, DomainEleOp >::OpAxisymmetricLhs ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  full_tangent_ptr  )

inline

Definition at line 180 of file MatAxisymmetric.hpp.

      : OpBase(field_name, field_name, DomainEleOp::OPROWCOL),
        fullTangentPtr(full_tangent_ptr) {}

Member Function Documentation

iNtegrate()

template<AssemblyType A, typename DomainEleOp >

MoFEMErrorCode MatOps::OpAxisymmetricLhs< A, DomainEleOp >::iNtegrate ( EntData &  row_data, EntData &  col_data  )

inlineoverrideprotected

Definition at line 188 of file MatAxisymmetric.hpp.

                                                                          {
    MoFEMFunctionBegin;
    using DL = DataLayoutTraits<DataLayout::GaussByCoeffs>;
    auto t_D =
        MatrixSizeHelper<GetFTensor4FromMatType<3, 3, 3, 3, -1, DL>, DL>::get(
            *fullTangentPtr, OpBase::nbIntegrationPts)();
    auto t_w = OpBase::getFTensor0IntegrationWeight();
    auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
    auto t_row_base = row_data.getFTensor0N();
    auto t_row_diff = row_data.getFTensor1DiffN<2>();
    FTENSOR_INDEX(3, i);
    FTENSOR_INDEX(3, J);
    FTENSOR_INDEX(3, k);
    FTENSOR_INDEX(3, L);
    FTENSOR_INDEX(2, j);
    FTENSOR_INDEX(2, I);
    FTENSOR_INDEX(2, l);
    FTensor::Number<0> N0;
    FTensor::Number<2> N2;
    FTensor::Tensor2<double, 3, 3> t_hoop;
    t_hoop(i, J) = 0;
    t_hoop(N2, N2) = 1;
    for (int gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
      const double radius = t_coords(0);
      const double alpha = 2 * M_PI * t_w * OpBase::getMeasure();
      int rr = 0;
      for (; rr != OpBase::nbRows / 2; ++rr) {
        auto t_m = OpBase::template getLocMat<2>(2 * rr);
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        auto t_col_diff = col_data.getFTensor1DiffN<2>(gg, 0);
        for (int cc = 0; cc != OpBase::nbCols / 2; ++cc) {
          t_m(j, N0) += alpha * t_col_base * t_row_diff(I) *
                        (t_D(j, I, k, L) * t_hoop(k, L));
          t_m(N0, l) += alpha * t_row_base * t_hoop(i, J) *
                        (t_D(i, J, l, I) * t_col_diff(I));
          t_m(N0, N0) += alpha * t_row_base * t_col_base / radius *
                          t_hoop(i, J) *
                          (t_D(i, J, k, L) * t_hoop(k, L));
          ++t_m;
          ++t_col_base;
          ++t_col_diff;
        }
        ++t_row_base;
        ++t_row_diff;
      }
      for (; rr < OpBase::nbRowBaseFunctions; ++rr) {
        ++t_row_base;
        ++t_row_diff;
      }
      ++t_D;
      ++t_w;
      ++t_coords;
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

fullTangentPtr

template<AssemblyType A, typename DomainEleOp >

boost::shared_ptr<MatrixDouble> MatOps::OpAxisymmetricLhs< A, DomainEleOp >::fullTangentPtr

protected

Definition at line 186 of file MatAxisymmetric.hpp.

---

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

• src/materials/MatAxisymmetric.hpp