OpEnergyDensity Struct Reference

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

Inheritance diagram for OpEnergyDensity:

Collaboration diagram for OpEnergyDensity:

Public Member Functions
	OpEnergyDensity (const std::string &field_name, boost::shared_ptr< MatrixDouble > grad_u_negative, boost::shared_ptr< VectorDouble > energy_densiity_ptr, boost::shared_ptr< std::map< int, BlockData >> perm_block_sets_ptr, boost::shared_ptr< DataAtIntegrationPts > common_data_ptr)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Attributes
boost::shared_ptr< MatrixDouble >	gradUNegative
boost::shared_ptr< VectorDouble >	energyDensity
boost::shared_ptr< std::map< int, BlockData > >	permBlockSetsPtr
boost::shared_ptr< DataAtIntegrationPts >	commonDataPtr

Detailed Description

Definition at line 195 of file electrostatics.hpp.

Constructor & Destructor Documentation

OpEnergyDensity()

OpEnergyDensity::OpEnergyDensity ( const std::string &  field_name, boost::shared_ptr< MatrixDouble >  grad_u_negative, boost::shared_ptr< VectorDouble >  energy_densiity_ptr, boost::shared_ptr< std::map< int, BlockData >>  perm_block_sets_ptr, boost::shared_ptr< DataAtIntegrationPts >  common_data_ptr  )

inline

Definition at line 197 of file electrostatics.hpp.

      : DomainEleOp(field_name, field_name, OPROWCOL, false),
        energyDensity(energy_densiity_ptr), gradUNegative(grad_u_negative),
        permBlockSetsPtr(perm_block_sets_ptr), commonDataPtr(common_data_ptr) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBVERTEX] = true;
  }

Member Function Documentation

doWork()

MoFEMErrorCode OpEnergyDensity::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Definition at line 210 of file electrostatics.hpp.

                                                            {
    MoFEMFunctionBegin;
 
    size_t nb_gauss_pts = getGaussPts().size2();
    energyDensity->resize(nb_gauss_pts, false);
    energyDensity->clear();
    auto t_energy_density = getFTensor0FromVec(*energyDensity);
    auto t_negative_grad_u = getFTensor1FromMat<SPACE_DIM>(*gradUNegative);
    FTensor::Index<'I', SPACE_DIM> I;
    double blockPermittivity = 0.0;
    for (const auto &n : *permBlockSetsPtr) {
      if (n.second.domainEnts.find(getFEEntityHandle()) !=
          n.second.domainEnts.end()) {
        blockPermittivity = n.second.epsPermit;
      }
    }
    // E = 0.5 * E^2 * Perm
    for (int gg = 0; gg != nb_gauss_pts; gg++) {
      t_energy_density =
          0.5 * t_negative_grad_u(I) * t_negative_grad_u(I) * blockPermittivity;
      ++t_negative_grad_u;
      ++t_energy_density;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonDataPtr

boost::shared_ptr<DataAtIntegrationPts> OpEnergyDensity::commonDataPtr

private

Definition at line 244 of file electrostatics.hpp.

energyDensity

boost::shared_ptr<VectorDouble> OpEnergyDensity::energyDensity

private

Definition at line 242 of file electrostatics.hpp.

gradUNegative

boost::shared_ptr<MatrixDouble> OpEnergyDensity::gradUNegative

private

Definition at line 241 of file electrostatics.hpp.

permBlockSetsPtr

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

private

Definition at line 243 of file electrostatics.hpp.

---

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

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