SplitVertices.cpp

Go to the documentation of this file.

#include "SplitVertices.hpp"
#include "RefinerTagManager.hpp"

#include "MBParallelConventions.h"

namespace moab {

SplitVerticesBase::SplitVerticesBase( RefinerTagManager* tag_mgr )
{
  this->tag_manager = tag_mgr;
  this->mesh_out = tag_mgr->get_output_mesh();
}

SplitVerticesBase::~SplitVerticesBase()
{
}

EntitySource::EntitySource( int nc, RefinerTagManager* tag_mgr )
{
  this->tag_manager = tag_mgr;
  this->mesh_out = tag_mgr->get_output_mesh();
  this->num_corners = nc;
}

EntitySource::~EntitySource()
{
}

bool EntitySource::create_element(
  EntityType etyp, int nconn, const EntityHandle* elem_verts, EntityHandle& elem_handle,
  std::map<ProcessSet,int>& proc_partition_counts )
{
  // Get the global IDs of the input vertices
  //int stat;
  proc_partition_counts[this->tag_manager->get_element_procs()]++;
  if ( this->mesh_out->create_element( etyp, elem_verts, nconn, elem_handle ) != MB_SUCCESS )
    {
    return false;
    }
  this->push_back( EntitySourceRecord( this->num_corners, elem_handle, this->tag_manager->get_element_procs() ) );
  this->tag_manager->set_sharing( elem_handle, this->tag_manager->get_element_procs() );
  return true;
}

void EntitySource::assign_global_ids( std::map<ProcessSet,int>& gids )
{
  std::vector<EntityHandle> adjacencies;
  adjacencies.resize( this->num_corners );
  std::vector<EntitySourceRecord>::iterator it;
  for ( it = this->begin(); it != this->end(); ++ it )
    {
    int num_nodes;
    const EntityHandle* conn;
    this->mesh_out->get_connectivity( it->handle, conn, num_nodes );
    this->tag_manager->get_output_gids( this->num_corners, conn, it->ids );
    std::sort( it->ids.begin(), it->ids.end() );
    }
  std::sort( this->begin(), this->end() );
  for ( it = this->begin(); it != this->end(); ++ it )
    {
    int gid = gids[it->process_set] ++;
    this->tag_manager->set_gid( it->handle, gid );
#ifdef MB_DEBUG
    std::cout << "Assigning entity: " << it->handle << " GID: " << gid << "\n";
#endif // MB_DEBUG
    }
}

} // namespace moab