PoissonExample::AuxFunctions Struct Reference

#include <users_modules/basic_finite_elements/poisson/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 26 of file AuxPoissonFunctions.hpp.

Constructor & Destructor Documentation

AuxFunctions()

PoissonExample::AuxFunctions::AuxFunctions

(

const MoFEM::Interface &

m_field

)

Examples

AuxPoissonFunctions.hpp.

Definition at line 28 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

Assemble error vector.

Examples

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

Definition at line 55 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

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 42 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

)

Print error.

Examples

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

Definition at line 72 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

)

Test error.

Examples

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

Definition at line 85 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 101 of file AuxPoissonFunctions.hpp.

rAnk

const int PoissonExample::AuxFunctions::rAnk

private

Examples

AuxPoissonFunctions.hpp.

Definition at line 102 of file AuxPoissonFunctions.hpp.

---

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

• users_modules/basic_finite_elements/poisson/src/AuxPoissonFunctions.hpp