MoFEM::FieldEntity_change_order Struct Reference

structure to change FieldEntity order More...

#include <src/multi_indices/FieldEntsMultiIndices.hpp>

Collaboration diagram for MoFEM::FieldEntity_change_order:

Public Member Functions
	FieldEntity_change_order (const ApproximationOrder order, const bool reduce_tag_size=false)
void	operator() (boost::shared_ptr< FieldEntity > &e)
void	operator() (FieldEntity *e)

Private Attributes
const ApproximationOrder	order
const bool	reduceTagSize
std::vector< FieldData >	data

Detailed Description

structure to change FieldEntity order

Definition at line 374 of file FieldEntsMultiIndices.hpp.

Constructor & Destructor Documentation

FieldEntity_change_order()

MoFEM::FieldEntity_change_order::FieldEntity_change_order ( const ApproximationOrder  order, const bool  reduce_tag_size = false  )

inline

Definition at line 376 of file FieldEntsMultiIndices.hpp.

      : order(order), reduceTagSize(reduce_tag_size) {}

Member Function Documentation

operator()() [1/2]

void MoFEM::FieldEntity_change_order::operator() ( boost::shared_ptr< FieldEntity > &  e )

inline

Definition at line 379 of file FieldEntsMultiIndices.hpp.

                                                          {
    (*this)(e.get());
  }

operator()() [2/2]

void MoFEM::FieldEntity_change_order::operator()

(

FieldEntity *

e

)

Definition at line 63 of file FieldEntsMultiIndices.cpp.

                                                        {
 
  moab::Interface &moab = e->sPtr->getBasicDataPtr()->moab;
  const EntityHandle ent = e->getEnt();
  *const_cast<ApproximationOrder *>(e->getMaxOrderPtr()) = order;
  std::size_t nb_dofs = e->getOrderNbDofs(order) * e->getNbOfCoeffs();
 
  double *tag_field_data = nullptr;
  int tag_field_data_size = 0;
 
  auto set_verts = [&]() {
    if (e->getFieldRawPtr()->tagFieldDataVertsType == MB_TAG_SPARSE) {
      // Get pointer and size of field values tag
      rval = moab.tag_get_by_ptr(e->getFieldRawPtr()->th_FieldDataVerts, &ent,
                                 1, (const void **)&tag_field_data,
                                 &tag_field_data_size);
      if (nb_dofs) {
        if (nb_dofs != tag_field_data_size) {
          rval = moab.tag_set_data(e->getFieldRawPtr()->th_FieldDataVerts, &ent,
                                   1, &*data.begin());
          MOAB_THROW(rval);
        }
      } else if (rval == MB_SUCCESS) {
        rval = moab.tag_delete_data(e->getFieldRawPtr()->th_FieldDataVerts,
                                    &ent, 1);
        MOAB_THROW(rval);
      }
    } else {
      rval = moab.tag_get_by_ptr(e->getFieldRawPtr()->th_FieldDataVerts, &ent,
                                 1, (const void **)&tag_field_data);
      MOAB_THROW(rval);
      rval = moab.tag_set_data(e->getFieldRawPtr()->th_FieldDataVerts, &ent, 1,
                               tag_field_data);
      MOAB_THROW(rval);
    }
  };
 
  auto set_default = [&]() {
    if (reduceTagSize || nb_dofs) {
 
      // Get pointer and size of field values tag
      rval = moab.tag_get_by_ptr(e->getFieldRawPtr()->th_FieldData, &ent, 1,
                                 (const void **)&tag_field_data,
                                 &tag_field_data_size);
 
      if ((reduceTagSize && nb_dofs != tag_field_data_size) ||
          nb_dofs > tag_field_data_size) {
 
        // Tag exist and are some data on it
        if (rval == MB_SUCCESS) {
          // Size of tag is different than new size, so copy data to new
          // container
          data.resize(tag_field_data_size);
          FieldData *ptr_begin = static_cast<FieldData *>(tag_field_data);
          FieldData *ptr_end =
              static_cast<FieldData *>(tag_field_data) + tag_field_data_size;
          std::copy(ptr_begin, ptr_end, data.begin());
        }
 
        if (rval != MB_SUCCESS || nb_dofs) {
 
          // Set field dof data
          data.resize(nb_dofs, 0);
          int tag_size[1];
          tag_size[0] = data.size();
          void const *tag_data[] = {&data[0]};
          rval = moab.tag_set_by_ptr(e->getFieldRawPtr()->th_FieldData, &ent, 1,
                                     tag_data, tag_size);
          MOAB_THROW(rval);
          rval = moab.tag_get_by_ptr(e->getFieldRawPtr()->th_FieldData, &ent, 1,
                                     (const void **)&tag_field_data,
                                     &tag_field_data_size);
          MOAB_THROW(rval);
 
        } else {
 
          rval =
              moab.tag_delete_data(e->getFieldRawPtr()->th_FieldData, &ent, 1);
          MOAB_THROW(rval);
        }
      }
    }
  };
 
  switch (e->getEntType()) {
  case MBVERTEX:
    set_verts();
    break;
  default:
    set_default();
  }
 
  if (nb_dofs)
    const_cast<double *&>(*(e->getEntFieldDataPtr())) = tag_field_data;
  else
    const_cast<double *&>(*(e->getEntFieldDataPtr())) = nullptr;
}

Member Data Documentation

data

std::vector<FieldData> MoFEM::FieldEntity_change_order::data

private

Definition at line 387 of file FieldEntsMultiIndices.hpp.

order

const ApproximationOrder MoFEM::FieldEntity_change_order::order

private

Definition at line 385 of file FieldEntsMultiIndices.hpp.

reduceTagSize

const bool MoFEM::FieldEntity_change_order::reduceTagSize

private

Definition at line 386 of file FieldEntsMultiIndices.hpp.

---

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

• src/multi_indices/FieldEntsMultiIndices.hpp
• src/multi_indices/impl/FieldEntsMultiIndices.cpp