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; };

◆ 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);
 }