#include "users_modules/basic_finite_elements/src/PostProcOnRefMesh.hpp"

Inheritance diagram for PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >:

Collaboration diagram for PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >:

Public Member Functions
	OpGetFieldValuesOnSkinImpl (moab::Interface &post_proc_mesh, std::vector< EntityHandle > &map_gauss_pts, const std::string field_name, const std::string tag_name, boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide > side_fe, const std::string vol_fe_name, boost::shared_ptr< MatrixDouble > mat_ptr, bool save_on_tag)

MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Public Attributes
moab::Interface &	postProcMesh

std::vector< EntityHandle > &	mapGaussPts

boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide >	sideOpFe

const std::string	feVolName

boost::shared_ptr< MatrixDouble >	matPtr

const std::string	tagName

const std::string	fieldName

const bool	saveOnTag

Detailed Description

template<int RANK>
struct PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >

Definition at line 1065 of file PostProcOnRefMesh.hpp.

Constructor & Destructor Documentation

◆ OpGetFieldValuesOnSkinImpl()

template<int RANK>

PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >::OpGetFieldValuesOnSkinImpl	(	moab::Interface &	post_proc_mesh,
		std::vector< EntityHandle > &	map_gauss_pts,
		const std::string	field_name,
		const std::string	tag_name,
		boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide >	side_fe,
		const std::string	vol_fe_name,
		boost::shared_ptr< MatrixDouble >	mat_ptr,
		bool	save_on_tag
	)

inline

Definition at line 1077 of file PostProcOnRefMesh.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              field_name, UserDataOperator::OPCOL),
          postProcMesh(post_proc_mesh), mapGaussPts(map_gauss_pts),
          sideOpFe(side_fe), feVolName(vol_fe_name), matPtr(mat_ptr),
          tagName(tag_name), fieldName(field_name), saveOnTag(save_on_tag) {}

Member Function Documentation

◆ doWork()

template<int RANK>

MoFEMErrorCode PostProcFaceOnRefinedMesh::OpGetFieldValuesOnSkinImpl< RANK >::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

Definition at line 930 of file PostProcOnRefMesh.cpp.

                                                               {
  MoFEMFunctionBegin;
 
  if (type != MBVERTEX)
    MoFEMFunctionReturnHot(0);
 
  CHKERR loopSideVolumes(feVolName, *sideOpFe);
 
  // quit if tag is not needed
  if (!saveOnTag)
    MoFEMFunctionReturnHot(0);
 
  auto &m_field = getPtrFE()->mField;
  auto field_ptr = m_field.get_field_structure(fieldName);
  const int rank = field_ptr->getNbOfCoeffs();
  FieldSpace space = field_ptr->getSpace();
 
  // auto dof_ptr = data.getFieldDofs()[0];
  // int rank = dof_ptr->getNbOfCoeffs();
  // FieldSpace space = dof_ptr->getSpace();
 
  int full_size = rank * RANK;
  if (space == HDIV)
    full_size *= 3;
  // for paraview
  int tag_length = full_size > 3 && full_size < 9 ? 9 : full_size;
  double def_VAL[tag_length];
  bzero(def_VAL, tag_length * sizeof(double));
  Tag th;
  CHKERR postProcMesh.tag_get_handle(tagName.c_str(), tag_length,
                                     MB_TYPE_DOUBLE, th,
                                     MB_TAG_CREAT | MB_TAG_SPARSE, def_VAL);
 
  // zero tags, this for Vertex if H1 and TRI if Hdiv, EDGE for Hcurl
  // no need for L2
  const void *tags_ptr[mapGaussPts.size()];
  int nb_gauss_pts = data.getN().size1();
  if (mapGaussPts.size() != (unsigned int)nb_gauss_pts ||
      nb_gauss_pts != matPtr->size2()) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  switch (space) {
  case H1:
  case L2:
  case HDIV:
 
    CHKERR postProcMesh.tag_get_by_ptr(th, &mapGaussPts[0], mapGaussPts.size(),
                                       tags_ptr);
 
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
      const double *my_ptr2 = static_cast<const double *>(tags_ptr[gg]);
      double *my_ptr = const_cast<double *>(my_ptr2);
      for (int rr = 0; rr != RANK; ++rr) {
        for (int dd = 0; dd != rank; ++dd) {
          my_ptr[rank * rr + dd] = (*matPtr)(rank * rr + dd, gg);
        }
      }
    }
 
    break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
            "field with that space is not implemented");
  }
 
  MoFEMFunctionReturn(0);
}