OpElasticTools::FePrePostProcess Struct Reference

Not used at this stage. Could be used to do some calculations, before assembly of local elements. More...

#include <users_modules/multifield_plasticity/src/ElasticOperators.hpp>

Inheritance diagram for OpElasticTools::FePrePostProcess:

Collaboration diagram for OpElasticTools::FePrePostProcess:

Public Member Functions
	FePrePostProcess (MoFEM::Interface &m_field, boost::ptr_vector< DataFromMove > &bc_data, boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	preProcess ()
MoFEMErrorCode	postProcess ()

Public Attributes
MoFEM::Interface &	mField
boost::ptr_vector< DataFromMove > &	bcData
boost::shared_ptr< MethodForForceScaling >	methodsOpForMove
boost::shared_ptr< MethodForForceScaling >	methodsOpForRollersDisp
boost::shared_ptr< CommonData >	commonDataPtr
bool	runOnceFlag

Detailed Description

Not used at this stage. Could be used to do some calculations, before assembly of local elements.

Definition at line 1220 of file ElasticOperators.hpp.

Constructor & Destructor Documentation

FePrePostProcess()

OpElasticTools::FePrePostProcess::FePrePostProcess ( MoFEM::Interface &  m_field, boost::ptr_vector< DataFromMove > &  bc_data, boost::shared_ptr< CommonData >  common_data_ptr  )

inline

Definition at line 1235 of file ElasticOperators.hpp.

      : mField(m_field), bcData(bc_data), commonDataPtr(common_data_ptr) {}

Member Function Documentation

postProcess()

MoFEMErrorCode OpElasticTools::FePrePostProcess::postProcess ( )

inline

Definition at line 1317 of file ElasticOperators.hpp.

                               {
    MoFEMFunctionBeginHot;
    // if (ts_ctx == CTX_TSSETIJACOBIAN)
    //   CHKERR MatView(ts_B, PETSC_VIEWER_STDOUT_SELF);
    auto assemble = [](Mat a) {
      MoFEMFunctionBeginHot;
      if (a) {
        CHKERR MatAssemblyBegin(a, MAT_FINAL_ASSEMBLY);
        CHKERR MatAssemblyEnd(a, MAT_FINAL_ASSEMBLY);
      }
      MoFEMFunctionReturnHot(0);
    };
 
    if (ts_ctx == CTX_TSSETIJACOBIAN) {
 
      if (commonDataPtr->SEpEp)
        CHKERR assemble(commonDataPtr->SEpEp);
      if (commonDataPtr->STauTau)
        CHKERR assemble(commonDataPtr->STauTau);
 
      // if (true)
      //   CHKERR MatView(ts_B, PETSC_VIEWER_DRAW_WORLD);
      // if (commonDataPtr->SEpEp)
      //   CHKERR MatView(commonDataPtr->SEpEp, PETSC_VIEWER_DRAW_WORLD);
      // if (commonDataPtr->STauTau)
      //   CHKERR MatView(commonDataPtr->STauTau, PETSC_VIEWER_DRAW_WORLD);
 
    }
 
    MoFEMFunctionReturnHot(0);
  }

preProcess()

MoFEMErrorCode OpElasticTools::FePrePostProcess::preProcess ( )

inline

Definition at line 1240 of file ElasticOperators.hpp.

                              {
    MoFEMFunctionBegin;
    runOnceFlag = true;
    CHKERR VecSetOption(ts_F, VEC_IGNORE_NEGATIVE_INDICES, PETSC_TRUE);
 
    switch (ts_ctx) {
    case CTX_TSSETIJACOBIAN: {
      snes_ctx = CTX_SNESSETJACOBIAN;
      snes_B = ts_B;
      // CHKERR MatSetOption(ts_B, MAT_NO_OFF_PROC_ZERO_ROWS, PETSC_TRUE);
      // CHKERR MatSetOption(ts_B, MAT_KEEP_NONZERO_PATTERN, PETSC_TRUE);
 
      if (commonDataPtr->SEpEp)
        CHKERR MatZeroEntries(commonDataPtr->SEpEp);
      if (commonDataPtr->STauTau)
        CHKERR MatZeroEntries(commonDataPtr->STauTau);
 
      break;
    }
    case CTX_TSSETIFUNCTION: {
      snes_ctx = CTX_SNESSETFUNCTION;
      snes_f = ts_F;
      break;
    }
    default:
      break;
    }
 
    auto time = ts_t;
 
    // forceDataScaled = forceData;
    // pressureDataScaled = pressureData;
    if (methodsOpForMove)
      for (auto &bdata : bcData) {
        bdata.scaledValues = bdata.values;
        CHKERR MethodForForceScaling::applyScale(this, methodsOpForMove,
                                                 bdata.scaledValues);
      }
 
    //FIXME: obsolete, for debugging get rid of this flag
    if (time <= (*cache).rollers_stop_time) {
 
      if (methodsOpForRollersDisp)
        for (auto &roller : (*cache).rollerDataVec) {
          roller.BodyDispScaled = roller.rollerDisp;
          CHKERR MethodForForceScaling::applyScale(
              this, methodsOpForRollersDisp, roller.BodyDispScaled);
        }
 
      for (auto &roller : (*cache).rollerDataVec) {
        if (roller.methodOpForRollerPosition) {
          roller.BodyDispScaled.clear();
          CHKERR MethodForForceScaling::applyScale(
              this, roller.methodOpForRollerPosition, roller.BodyDispScaled);
        }
        if (roller.methodOpForRollerDirection) {
          roller.BodyDirectionScaled.clear();
          CHKERR MethodForForceScaling::applyScale(
              this, roller.methodOpForRollerDirection,
              roller.BodyDirectionScaled);
        }
      }
    }
 
    // if (snes_ctx == CTX_SNESNONE) {
    //   if (!dofsIndices.empty()) {
    //     CHKERR VecSetValues(snes_x, dofsIndices.size(),
    //     &*dofsIndices.begin(),
    //                         &*dofsValues.begin(), INSERT_VALUES);
    //   }
    //   CHKERR VecAssemblyBegin(snes_x);
    //   CHKERR VecAssemblyEnd(snes_x);
    // }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

bcData

boost::ptr_vector<DataFromMove>& OpElasticTools::FePrePostProcess::bcData

Definition at line 1223 of file ElasticOperators.hpp.

commonDataPtr

boost::shared_ptr<CommonData> OpElasticTools::FePrePostProcess::commonDataPtr

Definition at line 1228 of file ElasticOperators.hpp.

methodsOpForMove

boost::shared_ptr<MethodForForceScaling> OpElasticTools::FePrePostProcess::methodsOpForMove

Definition at line 1225 of file ElasticOperators.hpp.

methodsOpForRollersDisp

boost::shared_ptr<MethodForForceScaling> OpElasticTools::FePrePostProcess::methodsOpForRollersDisp

Definition at line 1226 of file ElasticOperators.hpp.

mField

MoFEM::Interface& OpElasticTools::FePrePostProcess::mField

Definition at line 1222 of file ElasticOperators.hpp.

runOnceFlag

bool OpElasticTools::FePrePostProcess::runOnceFlag

Definition at line 1233 of file ElasticOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/ElasticOperators.hpp