PoissonExample::AuxFunctions Struct Reference

#include <tutorials/cor-2to5/src/AuxPoissonFunctions.hpp>

Collaboration diagram for PoissonExample::AuxFunctions:

Public Member Functions
	AuxFunctions (const MoFEM::Interface &m_field)
MoFEMErrorCode	createGhostVec (Vec *ghost_vec) const
MoFEMErrorCode	assembleGhostVector (Vec ghost_vec) const
	Assemble error vector. More...
MoFEMErrorCode	printError (Vec ghost_vec)
	Print error. More...
MoFEMErrorCode	testError (Vec ghost_vec)
	Test error. More...

Private Attributes
MPI_Comm	cOmm
const int	rAnk

Detailed Description

Examples

analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, and analytical_poisson_field_split.cpp.

Definition at line 14 of file AuxPoissonFunctions.hpp.

Constructor & Destructor Documentation

AuxFunctions()

PoissonExample::AuxFunctions::AuxFunctions ( const MoFEM::Interface &  m_field )

inline

Examples

AuxPoissonFunctions.hpp.

Definition at line 16 of file AuxPoissonFunctions.hpp.

      : cOmm(m_field.get_comm()), rAnk(m_field.get_comm_rank()) {}

Member Function Documentation

assembleGhostVector()

MoFEMErrorCode PoissonExample::AuxFunctions::assembleGhostVector ( Vec  ghost_vec ) const

inline

Assemble error vector.

Examples

analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, and AuxPoissonFunctions.hpp.

Definition at line 43 of file AuxPoissonFunctions.hpp.

                                                          {
 
    MoFEMFunctionBegin;
    CHKERR VecAssemblyBegin(ghost_vec);
    CHKERR VecAssemblyEnd(ghost_vec);
    // accumulate errors from processors
    CHKERR VecGhostUpdateBegin(ghost_vec, ADD_VALUES, SCATTER_REVERSE);
    CHKERR VecGhostUpdateEnd(ghost_vec, ADD_VALUES, SCATTER_REVERSE);
    // scatter errors to all processors
    CHKERR VecGhostUpdateBegin(ghost_vec, INSERT_VALUES, SCATTER_FORWARD);
    CHKERR VecGhostUpdateEnd(ghost_vec, INSERT_VALUES, SCATTER_FORWARD);
    MoFEMFunctionReturn(0);
  }

createGhostVec()

MoFEMErrorCode PoissonExample::AuxFunctions::createGhostVec ( Vec *  ghost_vec ) const

inline

Create ghost vector to assemble errors from all element on distributed mesh. Ghost vector has size 1, where one element is owned by processor 0, other processor have one ghost element of zero element at processor 0.

[createGhostVec description]

Parameters

ghost_vec

pointer to created ghost vector

Returns

error code

Examples

analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, and AuxPoissonFunctions.hpp.

Definition at line 30 of file AuxPoissonFunctions.hpp.

                                                      {
 
    MoFEMFunctionBegin;
    int ghosts[] = {0};
    int nb_locals = rAnk == 0 ? 1 : 0;
    int nb_ghosts = rAnk > 0 ? 1 : 0;
    CHKERR VecCreateGhost(cOmm, nb_locals, 1, nb_ghosts, ghosts, ghost_vec);
    MoFEMFunctionReturn(0);
  }

printError()

MoFEMErrorCode PoissonExample::AuxFunctions::printError ( Vec  ghost_vec )

inline

Print error.

Examples

analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, and AuxPoissonFunctions.hpp.

Definition at line 60 of file AuxPoissonFunctions.hpp.

                                           {
 
    MoFEMFunctionBegin;
    double *e;
    CHKERR VecGetArray(ghost_vec, &e);
    CHKERR PetscPrintf(cOmm, "Approximation error %4.3e\n", sqrt(e[0]));
    CHKERR VecRestoreArray(ghost_vec, &e);
    MoFEMFunctionReturn(0);
  }

testError()

MoFEMErrorCode PoissonExample::AuxFunctions::testError ( Vec  ghost_vec )

inline

Test error.

Examples

analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, and AuxPoissonFunctions.hpp.

Definition at line 73 of file AuxPoissonFunctions.hpp.

                                          {
 
    MoFEMFunctionBegin;
    double *e;
    CHKERR VecGetArray(ghost_vec, &e);
    // Check if error is zero, otherwise throw error
    const double eps = 1e-8;
    if ((sqrt(e[0]) > eps) || (!boost::math::isnormal(e[0]))) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
              "Test failed, error too big");
    }
    CHKERR VecRestoreArray(ghost_vec, &e);
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

cOmm

MPI_Comm PoissonExample::AuxFunctions::cOmm

private

Examples

AuxPoissonFunctions.hpp.

Definition at line 89 of file AuxPoissonFunctions.hpp.

rAnk

const int PoissonExample::AuxFunctions::rAnk

private

Examples

AuxPoissonFunctions.hpp.

Definition at line 90 of file AuxPoissonFunctions.hpp.

---

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

• tutorials/cor-2to5/src/AuxPoissonFunctions.hpp