OpSpatialRotationPrec Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpSpatialRotationPrec:

Collaboration diagram for OpSpatialRotationPrec:

Public Types
using	OP = FormsIntegrators< VolUserDataOperator >::Assembly< PETSC >::OpBase

Public Member Functions
	OpSpatialRotationPrec (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_omega_prec)

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

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
const double	alphaOmegaPrec

Detailed Description

Add an intermediate-configuration rotation-gradient block to Schur PC storage.

Definition at line 869 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpSpatialRotationPrec::OP = FormsIntegrators<VolUserDataOperator>::Assembly<PETSC>::OpBase

Definition at line 871 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpSpatialRotationPrec()

OpSpatialRotationPrec::OpSpatialRotationPrec ( const std::string &  field_name, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, const double  alpha_omega_prec  )

inline

Definition at line 873 of file EshelbianOperators.hpp.

      : OP(field_name, field_name, OPROWCOL), dataAtPts(data_ptr),
        alphaOmegaPrec(alpha_omega_prec) {
    sYmm = false;
  }

Member Function Documentation

aSsemble()

MoFEMErrorCode OpSpatialRotationPrec::aSsemble ( EntData &  row_data, EntData &  col_data, const bool    )

overrideprotected

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 3635 of file EshelbianOperators.cpp.

                                                                              {
  MoFEMFunctionBegin;
  // A borrowed pipeline may assemble an exact PETSc tangent for diagnostics.
  PetscBool is_shell = PETSC_FALSE;
  CHKERR PetscObjectTypeCompare(getPetscObject(getTSB()), MATSHELL, &is_shell);
  if (!is_shell)
    MoFEMFunctionReturnHot(0);
  // The form base owns local assembly; only its destination differs.
  CHKERR SchurBackendMatSetValuesPtr::matSetValuesPreconditionedBlockPtr(
      getTSB(), row_data, col_data, locMat, ADD_VALUES);
  MoFEMFunctionReturn(0);
}

iNtegrate()

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

overrideprotected

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 3593 of file EshelbianOperators.cpp.

                                                                   {
  MoFEMFunctionBegin;
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Index<'l', 3> l;
  const auto t_delta = FTensor::Kronecker_Delta<double>();
  const double volume = getVolume();
  auto t_weight = getFTensor0IntegrationWeight();
  auto t_plastic_f = dataAtPts->getFTensorPlasticF(nbIntegrationPts);
  auto t_inv_plastic_f = dataAtPts->getFTensorInvPlasticF(nbIntegrationPts);
  auto t_row_grad = row_data.getFTensor1DiffN<3>();
 
  for (int gg = 0; gg != nbIntegrationPts; ++gg) {
    const double a =
        volume * t_weight * determinantTensor3by3(t_plastic_f) * alphaOmegaPrec;
    int rr = 0;
    for (; rr != nbRows / 3; ++rr) {
      auto t_mat = getLocMat<3>(3 * rr);
      FTensor::Tensor1<double, 3> t_row_intermediate;
      t_row_intermediate(j) = t_row_grad(i) * t_inv_plastic_f(i, j);
      auto t_col_grad = col_data.getFTensor1DiffN<3>(gg, 0);
      for (int cc = 0; cc != nbCols / 3; ++cc) {
        FTensor::Tensor1<double, 3> t_col_intermediate;
        t_col_intermediate(j) = t_col_grad(i) * t_inv_plastic_f(i, j);
        t_mat(k, l) += t_delta(k, l) *
                       (a * (t_row_intermediate(j) * t_col_intermediate(j)));
        ++t_mat;
        ++t_col_grad;
      }
      ++t_row_grad;
    }
    for (; rr != nbRowBaseFunctions; ++rr)
      ++t_row_grad;
    ++t_weight;
    ++t_plastic_f;
    ++t_inv_plastic_f;
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

alphaOmegaPrec

const double OpSpatialRotationPrec::alphaOmegaPrec

private

Definition at line 888 of file EshelbianOperators.hpp.

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpSpatialRotationPrec::dataAtPts

private

Definition at line 887 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp