ReadTetGen.hpp

Go to the documentation of this file.

#ifndef READ_TET_GEN_HPP
#define READ_TET_GEN_HPP

#include <iosfwd>
#include "moab/Forward.hpp"
#include "moab/ReaderIface.hpp"
#include <string>

namespace moab {

class ReadUtilIface;

/* TetGen mesh data is typically split into two or three files:
 * name.node : node data
 * name.ele  : tet data
 * name.face : tri data
 * The reader will attempt to guess the correct file names from
 * the single input file name unless explicit file names are
 * specified using file options as described below.
 *
 * File options:
 *   NODE_FILE=filename   : node file name
 *   ELE_FILE[=filename]  : require tet file and optionally specify file name
 *   FACE_FILE[=filename] : reequire tri file and optionally specify file name
 *   EDGE_FILE[=filename] : reequire edge file and optionally specify file name
 *   NODE_ATTR_LIST=name[,name[,...]] : List of tag names in which to store
 *                                     attribute values.  If the same name
 *                                     is repeated multiple times, multiple
 *                                     attribute values will be stored in the
 *                                     same tag as array values in the order
 *                                     they are read from the file.  If a 
 *                                     name is zero-length, the attribute data
 *                                     will be disgarded.  
 */
class ReadTetGen : public ReaderIface
{
   
public:

  static ReaderIface* factory( Interface* );

  ErrorCode load_file( const char* file_name,
                       const EntityHandle* file_set,
                       const FileOptions& opts,
                       const SubsetList* subset_list = 0,
                       const Tag* file_id_tag = 0 );

  ErrorCode read_tag_values( const char* file_name,
                             const char* tag_name,
                             const FileOptions& opts,
                             std::vector<int>& tag_values_out,
                             const SubsetList* subset_list = 0 );
  
  ReadTetGen(Interface* impl = NULL);

  virtual ~ReadTetGen();

private:

  Interface* mbIface;
  ReadUtilIface* readTool;

  ErrorCode open_file( const std::string& input_file_name,
                         const std::string& input_name_base,
                         const std::string& input_name_suffix,
                         const char* file_type_suffix,
                         const char* file_name_option,
                         const FileOptions& opts,
                         std::ifstream& file_stream,
                         bool file_required = false );

  ErrorCode read_line( std::istream& file, std::string& line, int& lineno );

  ErrorCode read_line( std::istream& file, 
                         double* values_out, int num_values,
                         int& lineno );

  ErrorCode parse_attr_list( const std::string& option_str,
                               std::vector<Tag>& tag_list,
                               std::vector<int>& index_list,
                               const char* group_designator = 0 );

  ErrorCode read_node_file( std::istream& file, 
                              const Tag* attr_tag_list,
                              const int* attr_tag_index,
                              int attr_tag_list_len,
                              std::vector<EntityHandle>& nodes );

  ErrorCode read_elem_file( EntityType type,
                              std::istream& file, 
                              const std::vector<EntityHandle>& nodes,
                              Range& elems );
};

} // namespace moab

#endif // defined(READ_TET_GEN_HPP)