Post Process

Collaboration diagram for Post Process:

Files
file	PostProcOnRefMesh.cpp
	Postprocess fields on refined mesh made for 10 Node tets.
file	PostProcOnRefMesh.hpp
	Post-process fields on refined mesh.

Functions
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc (const std::string field_name, Vec v=PETSC_NULL)
	Add operator to post-process L2, H1, Hdiv, Hcurl field value. More...
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc (const std::string field_name, const std::string tag_name, Vec v=PETSC_NULL)
	Add operator to post-process L2 or H1 field value. More...
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, Vec v=PETSC_NULL)
	Add operator to post-process L2 or H1 field gradient. More...
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, const std::string tag_name, Vec v=PETSC_NULL)
	Add operator to post-process L2 or H1 field gradient. More...
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, int space_dim, Vec v=PETSC_NULL)
	Add operator to post-process L2 or H1 field gradient. More...
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::writeFile (const std::string file_name)
	wrote results in (MOAB) format, use "file_name.h5m" More...

Detailed Description

Function Documentation

addFieldValuesGradientPostProc() [1/3]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc	(	const std::string	field_name,
		const std::string	tag_name,
		Vec	v = PETSC_NULL
	)

Add operator to post-process L2 or H1 field gradient.

Parameters

field_name
tag_name	to store results on post-process mesh
v	If vector is given, values from vector are used to set tags on mesh

Definition at line 226 of file PostProcOnRefMesh.hpp.

                                                                    {
    MoFEMFunctionBeginHot;
    ELEMENT::getOpPtrVector().push_back(
        new PostProcCommonOnRefMesh::OpGetFieldGradientValues(
            postProcMesh, mapGaussPts, field_name, tag_name, getCommonData(),
            v));
    MoFEMFunctionReturnHot(0);
  }

addFieldValuesGradientPostProc() [2/3]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc	(	const std::string	field_name,
		int	space_dim,
		Vec	v = PETSC_NULL
	)

Add operator to post-process L2 or H1 field gradient.

Parameters

field_name
space_dim	the dimension of the problem
v	If vector is given, values from vector are used to set tags on mesh

Definition at line 246 of file PostProcOnRefMesh.hpp.

                                                                    {
    MoFEMFunctionBeginHot;
    ELEMENT::getOpPtrVector().push_back(
        new PostProcCommonOnRefMesh::OpGetFieldGradientValues(
            postProcMesh, mapGaussPts, field_name, field_name + "_GRAD",
            getCommonData(), v, space_dim));
    MoFEMFunctionReturnHot(0);
  }

addFieldValuesGradientPostProc() [3/3]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc	(	const std::string	field_name,
		Vec	v = PETSC_NULL
	)

Add operator to post-process L2 or H1 field gradient.

Parameters

field_name
v	If vector is given, values from vector are used to set tags on mesh

Note: Name of the tag to store values on post-process mesh is the same as field name

Definition at line 207 of file PostProcOnRefMesh.hpp.

                                                                    {
    MoFEMFunctionBeginHot;
    ELEMENT::getOpPtrVector().push_back(
        new PostProcCommonOnRefMesh::OpGetFieldGradientValues(
            postProcMesh, mapGaussPts, field_name, field_name + "_GRAD",
            getCommonData(), v));
    MoFEMFunctionReturnHot(0);
  }

addFieldValuesPostProc() [1/2]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc	(	const std::string	field_name,
		const std::string	tag_name,
		Vec	v = PETSC_NULL
	)

Add operator to post-process L2 or H1 field value.

Parameters

field_name
tag_name	to store results on post-process mesh
v	If vector is given, values from vector are used to set tags on mesh

Definition at line 184 of file PostProcOnRefMesh.hpp.

                                                            {
    MoFEMFunctionBeginHot;
    ELEMENT::getOpPtrVector().push_back(
        new PostProcCommonOnRefMesh::OpGetFieldValues(postProcMesh, mapGaussPts,
                                                      field_name, tag_name,
                                                      getCommonData(), v));
    MoFEMFunctionReturnHot(0);
  }

addFieldValuesPostProc() [2/2]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc	(	const std::string	field_name,
		Vec	v = PETSC_NULL
	)

Add operator to post-process L2, H1, Hdiv, Hcurl field value.

Parameters

field_name
v	If vector is given, values from vector are used to set tags on mesh

Note: Name of the tag to store values on post-process mesh is the same as field name

Definition at line 165 of file PostProcOnRefMesh.hpp.

                                                            {
    MoFEMFunctionBeginHot;
    ELEMENT::getOpPtrVector().push_back(
        new PostProcCommonOnRefMesh::OpGetFieldValues(postProcMesh, mapGaussPts,
                                                      field_name, field_name,
                                                      getCommonData(), v));
    MoFEMFunctionReturnHot(0);
  }

writeFile()

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::writeFile

(

const std::string

file_name

)

wrote results in (MOAB) format, use "file_name.h5m"

Parameters

file_name

file name (should always end with .h5m)

Returns

error code

Definition at line 265 of file PostProcOnRefMesh.hpp.

                                                      {
    MoFEMFunctionBegin;
    ParallelComm *pcomm_post_proc_mesh =
        ParallelComm::get_pcomm(&postProcMesh, MYPCOMM_INDEX);
    if (pcomm_post_proc_mesh == NULL)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "ParallelComm not allocated");
    CHKERR postProcMesh.write_file(file_name.c_str(), "MOAB",
                                   "PARALLEL=WRITE_PART");
    MoFEMFunctionReturn(0);
  }