MixTransport::UnsaturatedFlowElement::preProcessVol Struct Reference

Pre-peprocessing Set head pressure rate and get inital essential boundary conditions. More...

#include "tutorials/cor-0to1/src/UnsaturatedFlow.hpp"

Collaboration diagram for MixTransport::UnsaturatedFlowElement::preProcessVol:

Public Member Functions
	preProcessVol (UnsaturatedFlowElement &ctx, boost::shared_ptr< ForcesAndSourcesCore > &fe_ptr)
MoFEMErrorCode	operator() ()

Public Attributes
UnsaturatedFlowElement &	cTx
boost::shared_ptr< ForcesAndSourcesCore >	fePtr

Detailed Description

Pre-peprocessing Set head pressure rate and get inital essential boundary conditions.

Definition at line 1330 of file UnsaturatedFlow.hpp.

Constructor & Destructor Documentation

preProcessVol()

MixTransport::UnsaturatedFlowElement::preProcessVol::preProcessVol ( UnsaturatedFlowElement &  ctx, boost::shared_ptr< ForcesAndSourcesCore > &  fe_ptr  )

inline

Definition at line 1335 of file UnsaturatedFlow.hpp.

        : cTx(ctx), fePtr(fe_ptr) /*,mArk(mark)*/ {}

Member Function Documentation

operator()()

MoFEMErrorCode MixTransport::UnsaturatedFlowElement::preProcessVol::operator() ( )

inline

Examples

mofem/tutorials/cor-0to1/src/UnsaturatedFlow.hpp.

Definition at line 1340 of file UnsaturatedFlow.hpp.

                                {
      MoFEMFunctionBegin;
      // Update pressure rates
      CHKERR fePtr->mField.getInterface<VecManager>()->setOtherLocalGhostVector(
          fePtr->problemPtr, "VALUES", string("VALUES") + "_t", ROW,
          fePtr->ts_u_t, INSERT_VALUES, SCATTER_REVERSE);
      switch (fePtr->ts_ctx) {
      case TSMethod::CTX_TSSETIFUNCTION:
        CHKERR VecSetOption(fePtr->ts_F, VEC_IGNORE_NEGATIVE_INDICES,
                            PETSC_TRUE);
        if (!cTx.bcIndices.empty()) {
          double scale;
          CHKERR cTx.scaleMethodFlux->getForceScale(fePtr->ts_t, scale);
          if (cTx.bcVecIds.size() != cTx.bcIndices.size()) {
            cTx.bcVecIds.insert(cTx.bcVecIds.begin(), cTx.bcIndices.begin(),
                                cTx.bcIndices.end());
            cTx.bcVecVals.resize(cTx.bcVecIds.size(), false);
            cTx.vecValsOnBc.resize(cTx.bcVecIds.size(), false);
          }
          CHKERR VecGetValues(cTx.D0, cTx.bcVecIds.size(),
                              &*cTx.bcVecIds.begin(), &*cTx.bcVecVals.begin());
          CHKERR VecGetValues(fePtr->ts_u, cTx.bcVecIds.size(),
                              &*cTx.bcVecIds.begin(),
                              &*cTx.vecValsOnBc.begin());
          cTx.bcVecVals *= scale;
          VectorDouble::iterator vit = cTx.bcVecVals.begin();
          const NumeredDofEntity *dof_ptr;
          for (std::vector<int>::iterator it = cTx.bcVecIds.begin();
               it != cTx.bcVecIds.end(); it++, vit++) {
            if (auto dof_ptr =
                    fePtr->problemPtr->getColDofsByPetscGlobalDofIdx(*it)
                        .lock())
              dof_ptr->getFieldData() = *vit;
          }
        } else {
          cTx.bcVecIds.resize(0);
          cTx.bcVecVals.resize(0);
          cTx.vecValsOnBc.resize(0);
        }
        break;
      default:
        // don nothing
        break;
      }
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

cTx

UnsaturatedFlowElement& MixTransport::UnsaturatedFlowElement::preProcessVol::cTx

Examples

mofem/tutorials/cor-0to1/src/UnsaturatedFlow.hpp.

Definition at line 1331 of file UnsaturatedFlow.hpp.

fePtr

boost::shared_ptr<ForcesAndSourcesCore> MixTransport::UnsaturatedFlowElement::preProcessVol::fePtr

Examples

mofem/tutorials/cor-0to1/src/UnsaturatedFlow.hpp.

Definition at line 1332 of file UnsaturatedFlow.hpp.

---

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

• tutorials/cor-0to1/src/UnsaturatedFlow.hpp