Post Process

Collaboration diagram for Post Process:

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

Classes
struct	PostProcCommonOnRefMesh
	Set of operators and data structures used for post-processing. More...
struct	PostProcCommonOnRefMesh::CommonData
struct	PostProcCommonOnRefMesh::CommonDataForVolume
struct	PostProcCommonOnRefMesh::OpGetFieldValues
	operator to post-process (save gradients on refined post-processing mesh) field gradient More...
struct	PostProcCommonOnRefMesh::OpGetFieldGradientValues
	operator to post-process (save gradients on refined post-processing mesh) field gradient More...
struct	PostProcTemplateOnRefineMesh< ELEMENT >
	Generic post-processing class. More...
struct	PostProcVolumeOnRefinedMesh
	Post processing. More...
struct	PostProcFatPrismOnRefinedMesh
	Postprocess on prism. More...
struct	PostProcFatPrismOnRefinedMesh::PointsMap3D
struct	PostProcFaceOnRefinedMesh
	Postprocess on face. More...
struct	PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >
struct	PostProcEdgeOnRefinedMesh
	Postprocess on edge. More...

Functions
MoFEMErrorCode	MoFEM::PostProcBrokenMeshInMoabBase< E >::writeFile (const std::string file_name)
	wrote results in (MOAB) format, use "file_name.h5m"
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc (const std::string field_name, Vec v=PETSC_NULLPTR)
	Add operator to post-process L2, H1, Hdiv, Hcurl field value.
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesPostProc (const std::string field_name, const std::string tag_name, Vec v=PETSC_NULLPTR)
	Add operator to post-process L2 or H1 field value.
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, Vec v=PETSC_NULLPTR)
	Add operator to post-process L2 or H1 field gradient.
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, const std::string tag_name, Vec v=PETSC_NULLPTR)
	Add operator to post-process L2 or H1 field gradient.
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::addFieldValuesGradientPostProc (const std::string field_name, int space_dim, Vec v=PETSC_NULLPTR)
	Add operator to post-process L2 or H1 field gradient.
MoFEMErrorCode	PostProcTemplateOnRefineMesh< ELEMENT >::writeFile (const std::string file_name, const char *file_type="MOAB", const char *file_options="PARALLEL=WRITE_PART")
	wrote results in (MOAB) format, use "file_name.h5m"

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_NULLPTR  )

inline

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

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 214 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_NULLPTR  )

inline

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

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 234 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_NULLPTR  )

inline

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

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

Examples

mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp.

Definition at line 195 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_NULLPTR  )

inline

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

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 172 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_NULLPTR  )

inline

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

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

Examples

mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp.

Definition at line 153 of file PostProcOnRefMesh.hpp.

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

writeFile() [1/2]

template<typename E >

MoFEMErrorCode MoFEM::PostProcBrokenMeshInMoabBase< E >::writeFile

(

const std::string

file_name

)

#include <src/post_proc/PostProcBrokenMeshInMoabBase.hpp>

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 668 of file PostProcBrokenMeshInMoabBase.hpp.

                                                                    {
  MoFEMFunctionBegin;
  auto pcomm_post_proc_mesh = getPostProcMeshPcommPtr();
  if (!pcomm_post_proc_mesh)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "ParallelComm not allocated");
 
  auto &post_proc_mesh = getPostProcMesh();
  Tag part_tag = pcomm_post_proc_mesh->part_tag();
  Range tagged_sets;
  CHK_MOAB_THROW(post_proc_mesh.get_entities_by_type_and_tag(
                     0, MBENTITYSET, &part_tag, NULL, 1, tagged_sets,
                     moab::Interface::UNION),
                 "get partition sets");
  if (tagged_sets.empty()) {
    EntityHandle part_set;
    CHK_MOAB_THROW(post_proc_mesh.create_meshset(MESHSET_SET, part_set),
                   "create partition meshset");
 
    Range ents;
    CHK_MOAB_THROW(post_proc_mesh.get_entities_by_handle(0, ents, false),
                   "get post-process entities");
    ents = subtract(ents, ents.subset_by_type(MBENTITYSET));
    CHK_MOAB_THROW(post_proc_mesh.add_entities(part_set, ents),
                   "add entities to partition meshset");
 
    int part = pcomm_post_proc_mesh->rank();
    CHK_MOAB_THROW(post_proc_mesh.tag_set_data(part_tag, &part_set, 1, &part),
                   "tag partition meshset");
  }
 
  CHKERR getPostProcMesh().write_file(file_name.c_str(), "MOAB",
                                      "PARALLEL=WRITE_PART");
  MoFEMFunctionReturn(0);
};

writeFile() [2/2]

template<class ELEMENT >

MoFEMErrorCode PostProcTemplateOnRefineMesh< ELEMENT >::writeFile ( const std::string  file_name, const char *  file_type = "MOAB", const char *  file_options = "PARALLEL=WRITE_PART"  )

inline

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

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

Parameters

file_name

file name (should always end with .h5m)

Returns

error code

Examples

mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp.

Definition at line 253 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(), file_type, file_options);
    MoFEMFunctionReturn(0);
  }