EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM > Struct Template Reference

Rhs for testing EP mapping with initial conditions. More...

#include "users_modules/multifield-thermoplasticity-private/src/EdgeFlippingOps.hpp"

Inheritance diagram for EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >:

Collaboration diagram for EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >:

Public Member Functions
	OpRhsSetInitEP (const std::string field_name, boost::shared_ptr< MatrixDouble > dot_EP_ptr, boost::shared_ptr< MatrixDouble > EP_ptr, boost::shared_ptr< MatrixDouble > grad_EP_ptr, boost::shared_ptr< double > initial_EP_ptr, boost::shared_ptr< double > peak_EP_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	dotEPPtr
boost::shared_ptr< MatrixDouble >	EPPtr
boost::shared_ptr< MatrixDouble >	gradEPPtr
boost::shared_ptr< MatrixDouble >	gradQPtr
boost::shared_ptr< double >	initialEPPtr
boost::shared_ptr< double >	peakEPPtr

Detailed Description

template<AssemblyType A, IntegrationType I, int DIM>
struct EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >

Rhs for testing EP mapping with initial conditions.

Examples

thermoplastic.cpp.

Definition at line 521 of file EdgeFlippingOps.hpp.

Constructor & Destructor Documentation

OpRhsSetInitEP()

template<AssemblyType A, IntegrationType I, int DIM>

EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::OpRhsSetInitEP ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  dot_EP_ptr, boost::shared_ptr< MatrixDouble >  EP_ptr, boost::shared_ptr< MatrixDouble >  grad_EP_ptr, boost::shared_ptr< double >  initial_EP_ptr, boost::shared_ptr< double >  peak_EP_ptr  )

inline

Definition at line 523 of file EdgeFlippingOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        dotEPPtr(dot_EP_ptr), EPPtr(EP_ptr), gradEPPtr(grad_EP_ptr),
        initialEPPtr(initial_EP_ptr), peakEPPtr(peak_EP_ptr) {}

Member Function Documentation

iNtegrate()

template<AssemblyType A, IntegrationType I, int DIM>

MoFEMErrorCode EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::iNtegrate ( EntitiesFieldData::EntData &  data )

inline

Examples

EdgeFlippingOps.hpp.

Definition at line 533 of file EdgeFlippingOps.hpp.

                                                           {
    MoFEMFunctionBegin;
 
    FTENSOR_INDEX(DIM, i);
    FTENSOR_INDEX(DIM, j);
    constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
    FTENSOR_INDEX(size_symm, L);
 
    const double vol = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    auto t_base = data.getFTensor0N();
    auto t_diff_base = data.getFTensor1DiffN<SPACE_DIM>();
 
#ifndef NDEBUG
    if (data.getDiffN().size1() != data.getN().size1())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 1");
    if (data.getDiffN().size2() != data.getN().size2() * SPACE_DIM) {
      MOFEM_LOG("SELF", Sev::error)
          << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
      MOFEM_LOG("SELF", Sev::error) << data.getN();
      MOFEM_LOG("SELF", Sev::error) << data.getDiffN();
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 2");
    }
#endif
 
    auto t_EP = getFTensor2SymmetricFromMat<SPACE_DIM>(*EPPtr);
    // auto t_grad_temp = getFTensor1FromMat<SPACE_DIM>(*gradTPtr);
 
    const double alpha = t_w * vol;
 
    FTensor::Tensor2_symmetric<double, DIM> t_set_EP;
 
    auto t_L = PlasticOps::symm_L_tensor(FTensor::Number<DIM>());
    auto &nf = AssemblyDomainEleOp::locF;
 
    for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
      t_set_EP(0, 0) = init_T(*initialEPPtr, *peakEPPtr, t_coords(0),
                              t_coords(1), t_coords(2));
      t_set_EP(1, 1) = t_set_EP(0, 0) - 0.01;
      t_set_EP(0, 1) = t_set_EP(0, 0) - 0.02;
 
      // if (DIM == 3) {
      //   t_set_EP(2, 2) = t_set_EP(0, 0) - 3.0;
      //   t_set_EP(0, 2) = t_set_EP(0, 0) - 4.0;
      //   t_set_EP(1, 2) = t_set_EP(0, 0) - 5.0;
      // }
 
      const double alpha = t_w * vol;
 
      FTensor::Tensor1<double, size_symm> t_set_EP_L;
      t_set_EP_L(L) = alpha * t_set_EP(i, j) * t_L(i, j, L);
 
      auto t_nf = getFTensor1FromArray<size_symm, size_symm>(nf);
      int bb = 0;
      for (; bb != AssemblyDomainEleOp::nbRows / size_symm; ++bb) {
        t_nf(L) -= t_base * (t_set_EP_L(L));
        ++t_base;
        ++t_diff_base;
        ++t_nf;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb) {
        ++t_base;
        ++t_diff_base;
      }
 
      ++t_EP;
      // ++t_grad_g;
 
      ++t_coords;
      ++t_w;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dotEPPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<MatrixDouble> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::dotEPPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 612 of file EdgeFlippingOps.hpp.

EPPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<MatrixDouble> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::EPPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 613 of file EdgeFlippingOps.hpp.

gradEPPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<MatrixDouble> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::gradEPPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 614 of file EdgeFlippingOps.hpp.

gradQPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<MatrixDouble> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::gradQPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 615 of file EdgeFlippingOps.hpp.

initialEPPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<double> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::initialEPPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 616 of file EdgeFlippingOps.hpp.

peakEPPtr

template<AssemblyType A, IntegrationType I, int DIM>

boost::shared_ptr<double> EdgeFlippingOps::OpRhsSetInitEP< A, I, DIM >::peakEPPtr

private

Examples

EdgeFlippingOps.hpp.

Definition at line 617 of file EdgeFlippingOps.hpp.

---

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

• users_modules/multifield-thermoplasticity-private/src/EdgeFlippingOps.hpp