ThermalElement::OpGetFieldAtGaussPts< OP > Struct Template Reference

operator to calculate temperature and rate of temperature at Gauss points More...

#include <users_modules/basic_finite_elements/src/ThermalElement.hpp>

Inheritance diagram for ThermalElement::OpGetFieldAtGaussPts< OP >:

Collaboration diagram for ThermalElement::OpGetFieldAtGaussPts< OP >:

Public Member Functions
	OpGetFieldAtGaussPts (const std::string field_name, VectorDouble &field_at_gauss_pts)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	operator calculating temperature and rate of temperature More...

Public Attributes
VectorDouble &	fieldAtGaussPts

Detailed Description

template<typename OP>
struct ThermalElement::OpGetFieldAtGaussPts< OP >

operator to calculate temperature and rate of temperature at Gauss points

Definition at line 177 of file ThermalElement.hpp.

Constructor & Destructor Documentation

OpGetFieldAtGaussPts()

template<typename OP >

ThermalElement::OpGetFieldAtGaussPts< OP >::OpGetFieldAtGaussPts ( const std::string  field_name, VectorDouble &  field_at_gauss_pts  )

inline

Definition at line 180 of file ThermalElement.hpp.

        : OP::UserDataOperator(field_name, OP::UserDataOperator::OPROW),
          fieldAtGaussPts(field_at_gauss_pts) {}

Member Function Documentation

doWork()

template<typename OP >

MoFEMErrorCode ThermalElement::OpGetFieldAtGaussPts< OP >::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

operator calculating temperature and rate of temperature

temperature temperature or rate of temperature is calculated multiplying shape functions by degrees of freedom

Definition at line 190 of file ThermalElement.hpp.

                                                          {
      MoFEMFunctionBeginHot;
      try {
 
        if (data.getFieldData().size() == 0)
          MoFEMFunctionReturnHot(0);
        int nb_dofs = data.getFieldData().size();
        int nb_gauss_pts = data.getN().size1();
 
        // initialize
        fieldAtGaussPts.resize(nb_gauss_pts);
        if (type == MBVERTEX) {
          // loop over shape functions on entities always start from
          // vertices, so if nodal shape functions are processed, vector of
          // field values is zero at initialization
          std::fill(fieldAtGaussPts.begin(), fieldAtGaussPts.end(), 0);
        }
 
        for (int gg = 0; gg < nb_gauss_pts; gg++) {
          fieldAtGaussPts[gg] +=
              inner_prod(data.getN(gg, nb_dofs), data.getFieldData());
        }
 
      } catch (const std::exception &ex) {
        std::ostringstream ss;
        ss << "throw in method: " << ex.what() << std::endl;
        SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, ss.str().c_str());
      }
 
      MoFEMFunctionReturnHot(0);
    }

Member Data Documentation

fieldAtGaussPts

template<typename OP >

VectorDouble& ThermalElement::OpGetFieldAtGaussPts< OP >::fieldAtGaussPts

Definition at line 179 of file ThermalElement.hpp.

---

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

• users_modules/basic_finite_elements/src/ThermalElement.hpp