OpElectrodeCharge< SPACE_DIM > Struct Template Reference

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

Inheritance diagram for OpElectrodeCharge< SPACE_DIM >:

Collaboration diagram for OpElectrodeCharge< SPACE_DIM >:

Public Member Functions
	OpElectrodeCharge (const std::string field_name, boost::shared_ptr< MatrixDouble > d_jump, SmartPetscObj< Vec > alpha_vec, boost::shared_ptr< std::map< int, BlockData >> electrode_block_sets_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	dJumpPtr
boost::shared_ptr< std::map< int, BlockData > >	elecBlockSetsPtr
SmartPetscObj< Vec >	petscVec

Detailed Description

template<int SPACE_DIM>
struct OpElectrodeCharge< SPACE_DIM >

Definition at line 357 of file electrostatics.hpp.

Constructor & Destructor Documentation

OpElectrodeCharge()

template<int SPACE_DIM>

OpElectrodeCharge< SPACE_DIM >::OpElectrodeCharge ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  d_jump, SmartPetscObj< Vec >  alpha_vec, boost::shared_ptr< std::map< int, BlockData >>  electrode_block_sets_ptr  )

inline

Definition at line 358 of file electrostatics.hpp.

      : IntEleOp(field_name, IntEleOp::OPROW, false), dJumpPtr(d_jump),
        petscVec(alpha_vec), elecBlockSetsPtr(electrode_block_sets_ptr) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBVERTEX] = true;
  }

Member Function Documentation

doWork()

template<int SPACE_DIM>

MoFEMErrorCode OpElectrodeCharge< SPACE_DIM >::doWork ( int  side, EntityType  type, EntData &  data  )

inline

Definition at line 367 of file electrostatics.hpp.

                                                                  {
    MoFEMFunctionBegin;
    FTensor::Index<'i', SPACE_DIM> i;
    int index = 0;
    auto fe_ent = getFEEntityHandle();
    auto nb_gauss_pts = getGaussPts().size2();
    double area = getMeasure();
    dJumpPtr->resize(SPACE_DIM, nb_gauss_pts, false);
    auto t_jump = getFTensor1FromMat<SPACE_DIM>(*dJumpPtr);
    auto t_normal = getFTensor1NormalsAtGaussPts();
    auto t_w = getFTensor0IntegrationWeight();
    for (auto &m : *elecBlockSetsPtr) {
      double alphaPart = 0.0;
      if (m.second.electrodeEnts.find(fe_ent) != m.second.electrodeEnts.end()) {
 
        for (int gg = 0; gg != nb_gauss_pts; gg++) {
          FTensor::Tensor1<double, SPACE_DIM> t_r;
          t_r(i) = t_normal(i);
          auto t_r_normalized = t_r.normalize();
          alphaPart += (t_jump(i) * t_r(i)) * t_w * area;
          ++t_jump;
          ++t_normal;
          ++t_w;
        }
        CHKERR ::VecSetValues(petscVec, 1, &index, &alphaPart,
                              ADD_VALUES); // add the alpha value to the vector
      }
      ++index;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dJumpPtr

template<int SPACE_DIM>

boost::shared_ptr<MatrixDouble> OpElectrodeCharge< SPACE_DIM >::dJumpPtr

protected

Definition at line 401 of file electrostatics.hpp.

elecBlockSetsPtr

template<int SPACE_DIM>

boost::shared_ptr<std::map<int, BlockData> > OpElectrodeCharge< SPACE_DIM >::elecBlockSetsPtr

protected

Definition at line 402 of file electrostatics.hpp.

petscVec

template<int SPACE_DIM>

SmartPetscObj<Vec> OpElectrodeCharge< SPACE_DIM >::petscVec

protected

Definition at line 403 of file electrostatics.hpp.

---

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

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