ReadUtilIface.hpp

Go to the documentation of this file.

#ifndef MOAB_READ_UTIL_IFACE_HPP
#define MOAB_READ_UTIL_IFACE_HPP

#include <vector>
#include <string>
#include "moab/Types.hpp"
#include "moab/Compiler.hpp"

namespace moab {

class Range;

class MB_DLL_EXPORT ReadUtilIface
{
public:

  ReadUtilIface(){}

  virtual ~ReadUtilIface(){}

  virtual ErrorCode get_node_coords(
    const int num_arrays,
    const int num_nodes, 
    const int preferred_start_id,
    EntityHandle& actual_start_handle, 
    std::vector<double*>& arrays,
    const int sequence_size = -1
    ) = 0;

  virtual ErrorCode get_element_connect(
    const int num_elements, 
    const int verts_per_element,
    const EntityType mdb_type,
    const int preferred_start_id, 
    EntityHandle& actual_start_handle, 
    EntityHandle*& array,
    int sequence_size = -1
    ) = 0;

  virtual ErrorCode gather_related_ents(Range &partition,
                                        Range &related_ents,
                                        EntityHandle *file_set = NULL) = 0;

  virtual ErrorCode create_entity_sets(
    EntityID num_sets,
    const unsigned* set_flags,
    EntityID preffered_start_id,
    EntityHandle& actual_start_handle
  ) = 0;

  virtual ErrorCode update_adjacencies(
    const EntityHandle start_handle,
    const int number_elements,
    const int number_vertices_per_element,
    const EntityHandle* conn_array
    ) = 0;

  static inline 
  void reorder( const int* order, EntityHandle* conn, 
                int num_elem, int node_per_elem );

  virtual ErrorCode report_error( const std::string& error ) = 0;

  virtual ErrorCode report_error( const char* error, ... ) MB_PRINTF(1) = 0;

  virtual ErrorCode get_ordered_vertices(EntityHandle *bound_ents, 
                                           int *sense, 
                                           int num_bound,
                                           int dim,
                                           EntityHandle *bound_verts, 
                                           EntityType &etype) = 0;

  virtual ErrorCode assign_ids( Tag id_tag, const Range& ents, 
                                  int start = 0 ) = 0;

  virtual ErrorCode assign_ids( Tag id_tag, const EntityHandle* ents, 
                                  size_t num_ents, int start = 0 ) = 0;

  virtual ErrorCode create_gather_set(EntityHandle& gather_set) = 0;

  virtual ErrorCode get_gather_set(EntityHandle& gather_set) = 0;
};

inline 
void ReadUtilIface::reorder( const int* order, EntityHandle* conn, 
                               int num_elem, int node_per_elem )
{
  std::vector<EntityHandle> elem(node_per_elem);
  EntityHandle* const end = conn + num_elem*node_per_elem;
  while (conn != end) {
    std::copy( conn, conn+node_per_elem, elem.begin() );
    for (int j = 0; j < node_per_elem; ++j) 
      conn[order[j]] = elem[j];
    conn += node_per_elem;
  }
}

} // namespace moab 

#endif