#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 1063 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 1075 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);
 }