definition of volume element More...

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

Inheritance diagram for ConvectiveMassElement::MyVolumeFE:

Collaboration diagram for ConvectiveMassElement::MyVolumeFE:

Public Member Functions
	MyVolumeFE (MoFEM::Interface &m_field)

int	getRule (int order)
	it is used to calculate nb. of Gauss integration points More...

MoFEMErrorCode	preProcess ()

MoFEMErrorCode	postProcess ()

Public Attributes
Mat	A

Vec	F

SmartPetscObj< Vec >	V

double	eNergy

Detailed Description

definition of volume element

Examples: nonlinear_dynamics.cpp.

Definition at line 31 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

◆ MyVolumeFE()

ConvectiveMassElement::MyVolumeFE::MyVolumeFE ( MoFEM::Interface & m_field )

Definition at line 22 of file ConvectiveMassElement.cpp.

    : VolumeElementForcesAndSourcesCore(m_field), A(PETSC_NULL), F(PETSC_NULL) {
  meshPositionsFieldName = "NoNE";
 
  auto create_vec = [&]() {
    constexpr int ghosts[] = {0};
    if (mField.get_comm_rank() == 0) {
      return createVectorMPI(mField.get_comm(), 1, 1);
    } else {
      return createVectorMPI(mField.get_comm(), 0, 1);
    }
  };
 
  V = create_vec();
}

Member Function Documentation

◆ getRule()

int ConvectiveMassElement::MyVolumeFE::getRule ( int order )

it is used to calculate nb. of Gauss integration points

for more details pleas look Reference:

Albert Nijenhuis, Herbert Wilf, Combinatorial Algorithms for Computers and Calculators, Second Edition, Academic Press, 1978, ISBN: 0-12-519260-6, LC: QA164.N54.

More details about algorithm http://people.sc.fsu.edu/~jburkardt/cpp_src/gm_rule/gm_rule.html

Definition at line 38 of file ConvectiveMassElement.cpp.

38{ return 2 * order; };

order

constexpr int order

Definition: dg_projection.cpp:18

◆ postProcess()

MoFEMErrorCode ConvectiveMassElement::MyVolumeFE::postProcess ( )

Definition at line 56 of file ConvectiveMassElement.cpp.

                                                            {
  MoFEMFunctionBeginHot;
 
  CHKERR VolumeElementForcesAndSourcesCore::postProcess();
 
  const double *array;
  switch (ts_ctx) {
  case CTX_TSNONE:
    CHKERR VecAssemblyBegin(V);
    CHKERR VecAssemblyEnd(V);
    CHKERR VecSum(V, &eNergy);
    break;
  default:
    break;
  }
 
  MoFEMFunctionReturnHot(0);
}

◆ preProcess()

MoFEMErrorCode ConvectiveMassElement::MyVolumeFE::preProcess ( )

Definition at line 40 of file ConvectiveMassElement.cpp.

                                                           {
  MoFEMFunctionBeginHot;
 
  CHKERR VolumeElementForcesAndSourcesCore::preProcess();
 
  switch (ts_ctx) {
  case CTX_TSNONE:
    CHKERR VecZeroEntries(V);
    break;
  default:
    break;
  }
 
  MoFEMFunctionReturnHot(0);
}