OpLhsSetInitT_dT< A, I > Struct Template Reference

Lhs for setting initial conditions wrt T. More...

#include "users_modules/multifield-thermoplasticity-private/tutorials/adv-2/src/ThermoElasticOps.hpp"

Inheritance diagram for OpLhsSetInitT_dT< A, I >:

Collaboration diagram for OpLhsSetInitT_dT< A, I >:

Public Member Functions
	OpLhsSetInitT_dT (const std::string field_name, boost::shared_ptr< VectorDouble > T_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Attributes
boost::shared_ptr< VectorDouble >	TPtr

Detailed Description

template<AssemblyType A, IntegrationType I>
struct OpLhsSetInitT_dT< A, I >

Lhs for setting initial conditions wrt T.

Examples

ThermoElasticOps.hpp.

Definition at line 738 of file ThermoElasticOps.hpp.

Constructor & Destructor Documentation

OpLhsSetInitT_dT()

template<AssemblyType A, IntegrationType I>

OpLhsSetInitT_dT< A, I >::OpLhsSetInitT_dT ( const std::string  field_name, boost::shared_ptr< VectorDouble >  T_ptr  )

inline

Definition at line 740 of file ThermoElasticOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        TPtr(T_ptr) {
    sYmm = false;
  }

Member Function Documentation

iNtegrate()

template<AssemblyType A, IntegrationType I>

MoFEMErrorCode OpLhsSetInitT_dT< A, I >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Examples

ThermoElasticOps.hpp.

Definition at line 748 of file ThermoElasticOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    auto t_row_base = row_data.getFTensor0N();
    auto t_row_diff_base = row_data.getFTensor1DiffN<SPACE_DIM>();
 
#ifndef NDEBUG
    if (row_data.getDiffN().size1() != row_data.getN().size1())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 1");
    if (row_data.getDiffN().size2() != row_data.getN().size2() * SPACE_DIM) {
      MOFEM_LOG("SELF", Sev::error)
          << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
      MOFEM_LOG("SELF", Sev::error) << row_data.getN();
      MOFEM_LOG("SELF", Sev::error) << row_data.getDiffN();
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 2");
    }
 
    if (col_data.getDiffN().size1() != col_data.getN().size1())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 1");
    if (col_data.getDiffN().size2() != col_data.getN().size2() * SPACE_DIM) {
      MOFEM_LOG("SELF", Sev::error)
          << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
      MOFEM_LOG("SELF", Sev::error) << col_data.getN();
      MOFEM_LOG("SELF", Sev::error) << col_data.getDiffN();
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 2");
    }
#endif
 
    auto t_T = getFTensor0FromVec(*TPtr);
    // auto t_grad_g = getFTensor1FromMat<SPACE_DIM>(*gradGPtr);
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double alpha = t_w * vol;
 
      int rr = 0;
      for (; rr != nbRows; ++rr) {
 
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        auto t_col_diff_base = col_data.getFTensor1DiffN<SPACE_DIM>(gg, 0);
 
        for (int cc = 0; cc != nbCols; ++cc) {
 
          locMat(rr, cc) += (t_row_base * t_col_base * alpha);
 
          ++t_col_base;
          ++t_col_diff_base;
        }
 
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      ++t_T;
      // ++t_grad_g;
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

TPtr

template<AssemblyType A, IntegrationType I>

boost::shared_ptr<VectorDouble> OpLhsSetInitT_dT< A, I >::TPtr

private

Examples

ThermoElasticOps.hpp.

Definition at line 820 of file ThermoElasticOps.hpp.

---

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

• users_modules/multifield-thermoplasticity-private/tutorials/adv-2/src/ThermoElasticOps.hpp