OpTotalEnergy Struct Reference

#include <tutorials/scl-12/src/electrostatics.hpp>

Inheritance diagram for OpTotalEnergy:

Collaboration diagram for OpTotalEnergy:

Public Member Functions
	OpTotalEnergy (const std::string &field_name, boost::shared_ptr< MatrixDouble > grad_u_negative, boost::shared_ptr< std::map< int, BlockData >> perm_block_sets_ptr, boost::shared_ptr< std::map< int, BlockData >> int_domain_block_sets_ptr, boost::shared_ptr< DataAtIntegrationPts > common_data_ptr, SmartPetscObj< Vec > petsc_vec_energy)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	gradUNegative
boost::shared_ptr< std::map< int, BlockData > >	permBlockSetsPtr
boost::shared_ptr< std::map< int, BlockData > >	intDomainBlocksetsPtr
boost::shared_ptr< DataAtIntegrationPts >	commonDataPtr
SmartPetscObj< Vec >	petscTotalEnergy

Detailed Description

Definition at line 135 of file electrostatics.hpp.

Constructor & Destructor Documentation

OpTotalEnergy()

OpTotalEnergy::OpTotalEnergy ( const std::string &  field_name, boost::shared_ptr< MatrixDouble >  grad_u_negative, boost::shared_ptr< std::map< int, BlockData >>  perm_block_sets_ptr, boost::shared_ptr< std::map< int, BlockData >>  int_domain_block_sets_ptr, boost::shared_ptr< DataAtIntegrationPts >  common_data_ptr, SmartPetscObj< Vec >  petsc_vec_energy  )

inline

Definition at line 136 of file electrostatics.hpp.

      : DomainEleOp(field_name, DomainEleOp::OPROW),
        gradUNegative(grad_u_negative), permBlockSetsPtr(perm_block_sets_ptr),
        intDomainBlocksetsPtr(int_domain_block_sets_ptr),
        commonDataPtr(common_data_ptr), petscTotalEnergy(petsc_vec_energy) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBVERTEX] = true;
  }

Member Function Documentation

doWork()

MoFEMErrorCode OpTotalEnergy::doWork ( int  side, EntityType  type, EntData &  data  )

inline

Definition at line 150 of file electrostatics.hpp.

                                                                  {
    MoFEMFunctionBegin;
 
    auto t_negative_grad_u = getFTensor1FromMat<SPACE_DIM>(*gradUNegative);
    FTensor::Index<'I', SPACE_DIM> I;
    double area = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto nb_gauss_pts = getGaussPts().size2();
    double totalEnergy = 0.0;
    int index = 0;
    // Total Energy Density = 0.5 * E^2 * Perm * Area
    for (const auto &m : *intDomainBlocksetsPtr) {
 
      double blockPermittivity = 0.0;
 
      if (m.second.internalDomainEnts.find(getFEEntityHandle()) !=
          m.second.internalDomainEnts.end()) {
        for (const auto &n : *permBlockSetsPtr) {
          if (n.second.domainEnts.find(getFEEntityHandle()) !=
              n.second.domainEnts.end()) {
            blockPermittivity = n.second.epsPermit;
          }
        }
 
        for (int gg = 0; gg != nb_gauss_pts; gg++) {
          totalEnergy += 0.5 * t_negative_grad_u(I) * t_negative_grad_u(I) *
                         blockPermittivity * t_w * area;
          ++t_negative_grad_u;
          ++t_w;
        }
      }
    }
    CHKERR VecSetValues(petscTotalEnergy, 1, &index, &totalEnergy, ADD_VALUES);
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonDataPtr

boost::shared_ptr<DataAtIntegrationPts> OpTotalEnergy::commonDataPtr

private

Definition at line 191 of file electrostatics.hpp.

gradUNegative

boost::shared_ptr<MatrixDouble> OpTotalEnergy::gradUNegative

private

Definition at line 188 of file electrostatics.hpp.

intDomainBlocksetsPtr

boost::shared_ptr<std::map<int, BlockData> > OpTotalEnergy::intDomainBlocksetsPtr

private

Definition at line 190 of file electrostatics.hpp.

permBlockSetsPtr

boost::shared_ptr<std::map<int, BlockData> > OpTotalEnergy::permBlockSetsPtr

private

Definition at line 189 of file electrostatics.hpp.

petscTotalEnergy

SmartPetscObj<Vec> OpTotalEnergy::petscTotalEnergy

private

Definition at line 192 of file electrostatics.hpp.

---

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

• tutorials/scl-12/src/electrostatics.hpp