TS monitore it records temperature at time steps. More...

#include "users_modules/basic_finite_elements/src/ThermalElement.hpp"

Inheritance diagram for ThermalElement::TimeSeriesMonitor:

Collaboration diagram for ThermalElement::TimeSeriesMonitor:

Public Types
using	Ele = ForcesAndSourcesCore

using	VolEle = VolumeElementForcesAndSourcesCore

using	VolOp = VolumeElementForcesAndSourcesCore::UserDataOperator

using	SetPtsData = FieldEvaluatorInterface::SetPtsData

Public Member Functions
	TimeSeriesMonitor (MoFEM::Interface &m_field, const std::string series_name, const std::string temp_name, std::vector< std::array< double, 3 > > eval_points={})

MoFEMErrorCode	postProcess ()

Public Attributes
MoFEM::Interface &	mField

const std::string	seriesName

const std::string	tempName

BitRefLevel	mask

boost::shared_ptr< SetPtsData >	dataFieldEval

boost::shared_ptr< VectorDouble >	tempPtr

std::vector< std::array< double, 3 > >	evalPoints

Detailed Description

TS monitore it records temperature at time steps.

Definition at line 580 of file ThermalElement.hpp.

Member Typedef Documentation

◆ Ele

using ThermalElement::TimeSeriesMonitor::Ele = ForcesAndSourcesCore

Definition at line 587 of file ThermalElement.hpp.

◆ SetPtsData

using ThermalElement::TimeSeriesMonitor::SetPtsData = FieldEvaluatorInterface::SetPtsData

Definition at line 590 of file ThermalElement.hpp.

◆ VolEle

using ThermalElement::TimeSeriesMonitor::VolEle = VolumeElementForcesAndSourcesCore

Definition at line 588 of file ThermalElement.hpp.

◆ VolOp

using ThermalElement::TimeSeriesMonitor::VolOp = VolumeElementForcesAndSourcesCore::UserDataOperator

Definition at line 589 of file ThermalElement.hpp.

Constructor & Destructor Documentation

◆ TimeSeriesMonitor()

ThermalElement::TimeSeriesMonitor::TimeSeriesMonitor	(	MoFEM::Interface &	m_field,
		const std::string	series_name,
		const std::string	temp_name,
		std::vector< std::array< double, 3 > >	eval_points = {} )

inline

Definition at line 597 of file ThermalElement.hpp.

                                                                   {})
        : mField(m_field), seriesName(series_name), tempName(temp_name),
          evalPoints(eval_points),
          dataFieldEval(m_field.getInterface<FieldEvaluatorInterface>()
                            ->getData<VolEle>()) {
      mask.set();
 
      if (!evalPoints.empty()) {
        ierr = m_field.getInterface<FieldEvaluatorInterface>()->buildTree<3>(
            dataFieldEval, "THERMAL_FE");
        CHKERRABORT(PETSC_COMM_WORLD, ierr);
 
        auto no_rule = [](int, int, int) { return -1; };
 
        tempPtr = boost::make_shared<VectorDouble>();
 
        boost::shared_ptr<Ele> vol_ele(dataFieldEval->feMethodPtr.lock());
        vol_ele->getRuleHook = no_rule;
 
        vol_ele->getOpPtrVector().push_back(
            new OpCalculateScalarFieldValues("TEMP", tempPtr));
      }
    }

Member Function Documentation

◆ postProcess()

MoFEMErrorCode ThermalElement::TimeSeriesMonitor::postProcess ( )

Definition at line 454 of file ThermalElement.cpp.

                                                            {
  MoFEMFunctionBegin;
 
  CHKERR mField.getInterface<VecManager>()->setGlobalGhostVector(
      problemPtr, ROW, ts_u, INSERT_VALUES, SCATTER_REVERSE);
 
  BitRefLevel proble_bit_level = problemPtr->getBitRefLevel();
 
  SeriesRecorder *recorder_ptr = NULL;
  CHKERR mField.getInterface(recorder_ptr);
  CHKERR recorder_ptr->record_begin(seriesName);
  CHKERR recorder_ptr->record_field(seriesName, tempName, proble_bit_level,
                                    mask);
  CHKERR recorder_ptr->record_end(seriesName, ts_t);
 
  auto post_proc_at_points = [&](std::array<double, 3> point, int num) {
    MoFEMFunctionBegin;
 
    dataFieldEval->setEvalPoints(point.data(), point.size() / 3);
 
    struct OpPrint : public VolOp {
 
      std::array<double, 3> pointCoords;
      int pointNum;
      boost::shared_ptr<VectorDouble> tempPtr;
 
      OpPrint(boost::shared_ptr<VectorDouble> temp_ptr,
              std::array<double, 3> &point_coords, int point_num)
          : VolOp("TEMP", VolOp::OPROW), tempPtr(temp_ptr),
            pointCoords(point_coords), pointNum(point_num) {}
 
      MoFEMErrorCode doWork(int side, EntityType type,
                            EntitiesFieldData::EntData &data) {
        MoFEMFunctionBegin;
        if (type == MBVERTEX) {
          if (getGaussPts().size2()) {
 
            auto t_p = getFTensor0FromVec(*tempPtr);
 
            MOFEM_LOG("THERMALSYNC", Sev::inform)
                << "Pnt: " << std::to_string(pointNum)
                << " Crd: " << getVectorAdaptor(pointCoords.data(), 3)
                << " Tmp: " << t_p;
          }
        }
        MoFEMFunctionReturn(0);
      }
    };
 
    if (auto fe_ptr = dataFieldEval->feMethodPtr.lock()) {
      fe_ptr->getOpPtrVector().push_back(new OpPrint(tempPtr, point, num));
      CHKERR mField.getInterface<FieldEvaluatorInterface>()
          ->evalFEAtThePoint<3>(point.data(), 1e-12, "DMTHERMAL", "THERMAL_FE",
                                dataFieldEval, mField.get_comm_rank(),
                                mField.get_comm_rank(), nullptr, MF_EXIST,
                                QUIET);
      fe_ptr->getOpPtrVector().pop_back();
    }
 
    MoFEMFunctionReturn(0);
  };
 
  if (!evalPoints.empty()) {
    int num = 0;
    for (auto p : evalPoints)
      CHKERR post_proc_at_points(p, num++);
  }
 
  MOFEM_LOG_SYNCHRONISE(mField.get_comm());
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

◆ dataFieldEval

boost::shared_ptr<SetPtsData> ThermalElement::TimeSeriesMonitor::dataFieldEval

Definition at line 592 of file ThermalElement.hpp.

◆ evalPoints

std::vector<std::array<double, 3> > ThermalElement::TimeSeriesMonitor::evalPoints

Definition at line 595 of file ThermalElement.hpp.

◆ mask

BitRefLevel ThermalElement::TimeSeriesMonitor::mask

Definition at line 585 of file ThermalElement.hpp.

◆ mField

MoFEM::Interface& ThermalElement::TimeSeriesMonitor::mField

Definition at line 582 of file ThermalElement.hpp.

◆ seriesName

const std::string ThermalElement::TimeSeriesMonitor::seriesName

Definition at line 583 of file ThermalElement.hpp.

◆ tempName

const std::string ThermalElement::TimeSeriesMonitor::tempName

Definition at line 584 of file ThermalElement.hpp.

◆ tempPtr

boost::shared_ptr<VectorDouble> ThermalElement::TimeSeriesMonitor::tempPtr

Definition at line 593 of file ThermalElement.hpp.

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

users_modules/basic_finite_elements/src/ThermalElement.hpp
users_modules/basic_finite_elements/src/impl/ThermalElement.cpp

Public Types

Public Member Functions

Public Attributes