#include <ReadIDEAS.hpp>

Inheritance diagram for moab::ReadIDEAS:

Collaboration diagram for moab::ReadIDEAS:

Public Member Functions
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)
	Load mesh from a file. More...

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)
	Read tag values from a file. More...

	ReadIDEAS (Interface *impl=NULL)
	Constructor. More...

virtual	~ReadIDEAS ()
	Destructor. More...

Public Member Functions inherited from moab::ReaderIface
virtual	~ReaderIface ()

Static Public Member Functions
static ReaderIface *	factory (Interface *)

Protected Member Functions
ErrorCode	skip_header ()

ErrorCode	create_vertices (EntityHandle &first_vertex, const Tag *file_id_tag)

ErrorCode	create_elements (EntityHandle first_vertex, const Tag *file_id_tag)

Private Attributes
std::ifstream	file

RangeMap< int, EntityHandle >	nodeIdMap

ReadUtilIface *	readMeshIface

Interface *	MBI

Static Private Attributes
static const unsigned	SINGLE_PRECISION_NODES = 15

static const unsigned	DOUBLE_PRECISION_NODES0 = 781

static const unsigned	DOUBLE_PRECISION_NODES1 = 2411

static const unsigned	ELEMENTS0 = 71

static const unsigned	ELEMENTS1 = 780

static const unsigned	ELEMENTS2 = 2412

static const int	ROD0 = 11

static const int	ROD1 = 171

static const int	TRI0 = 41

static const int	TRI1 = 91

static const int	QUAD0 = 44

static const int	QUAD1 = 94

static const int	TET = 111

static const int	WEDGE = 112

static const int	HEX = 115

Detailed Description

Definition at line 24 of file ReadIDEAS.hpp.

Constructor & Destructor Documentation

◆ ReadIDEAS()

moab::ReadIDEAS::ReadIDEAS ( Interface * impl = NULL )

Constructor.

Definition at line 25 of file ReadIDEAS.cpp.

                                       : MBI( impl )
 {
     impl->query_interface( readMeshIface );
 }

References moab::Interface::query_interface(), and readMeshIface.

Referenced by factory().

◆ ~ReadIDEAS()

virtual moab::ReadIDEAS::~ReadIDEAS ( )

inlinevirtual

Destructor.

Definition at line 46 of file ReadIDEAS.hpp.

46 {}

Member Function Documentation

◆ create_elements()

ErrorCode moab::ReadIDEAS::create_elements	(	EntityHandle	first_vertex,
		const Tag *	file_id_tag
	)

protected

Definition at line 210 of file ReadIDEAS.cpp.

 {
     char line1[10000], line2[10000];
     int il1, il2;
     char *ctmp1, *ctmp2;
     std::string s1, s2;
     ErrorCode rval;
     EntityHandle handle;
  
     Tag mat_tag, phys_tag, id_tag;
     rval = MBI->tag_get_handle( MAT_PROP_TABLE_TAG, 1, MB_TYPE_INTEGER, mat_tag, MB_TAG_DENSE | MB_TAG_CREAT );
     if( MB_SUCCESS != rval && MB_ALREADY_ALLOCATED != rval ) return rval;
     rval = MBI->tag_get_handle( PHYS_PROP_TABLE_TAG, 1, MB_TYPE_INTEGER, phys_tag, MB_TAG_DENSE | MB_TAG_CREAT );
     if( MB_SUCCESS != rval && MB_ALREADY_ALLOCATED != rval ) return rval;
     id_tag = MBI->globalId_tag();
  
     for( ;; )
     {
         if( !file.getline( line1, 10000 ) || !file.getline( line2, 10000 ) ) return MB_FAILURE;
  
         // Check if we are at the end of the block
         il1 = std::strtol( line1, &ctmp1, 10 );
         il2 = std::strtol( line2, &ctmp2, 10 );
         if( ( il1 == -1 ) && ( il2 == -1 ) )
         {
             s1 = ctmp1;
             s2 = ctmp2;
             if( ( s1.empty() ) && ( s2.empty() ) ) return MB_SUCCESS;
         }
  
         // The first line describes attributes of the element other than connectivity.
         const int element_id = strtol( line1 + 1, &ctmp1, 10 );
         const int ideas_type = strtol( line1 + 11, &ctmp1, 10 );
         const int phys_table = strtol( line1 + 21, &ctmp1, 10 );
         const int mat_table  = strtol( line1 + 31, &ctmp1, 10 );
  
         // Determine the element type.
         EntityType mb_type;
         if( TRI0 == ideas_type || TRI1 == ideas_type )
             mb_type = MBTRI;
         else if( QUAD0 == ideas_type || QUAD1 == ideas_type )
             mb_type = MBQUAD;
         else if( TET == ideas_type )
             mb_type = MBTET;
         else if( HEX == ideas_type )
             mb_type = MBHEX;
         else if( WEDGE == ideas_type )
             mb_type = MBPRISM;
         else
         {
             std::cout << "IDEAS element type not yet added to MOAB reader." << std::endl;
             return MB_NOT_IMPLEMENTED;
         }
  
         // Get the connectivity out of the 2nd line
         std::stringstream ss( line2 );
         const int n_conn = CN::VerticesPerEntity( mb_type );
         EntityHandle conn[CN::MAX_NODES_PER_ELEMENT];
         EntityHandle vert;
         for( int i = 0; i < n_conn; ++i )
         {
             ss >> vert;
             conn[i] = vstart + vert - 1;
         }
  
         // Make the element. According to the Gmsh 2.2.3 source code, the IDEAS
         // canonical numbering is the same as MBCN.
         MB_CHK_SET_ERR( MBI->create_element( mb_type, conn, n_conn, handle ),
                         "can't create elements of type " << mb_type );
  
         // If the phys set does not already exist, create it.
         Range phys_sets;
         EntityHandle phys_set;
         const void* const phys_set_id_val[] = { &phys_table };
         MB_CHK_SET_ERR( MBI->get_entities_by_type_and_tag( 0, MBENTITYSET, &phys_tag, phys_set_id_val, 1, phys_sets ),
                         "can't get phys sets" );
         if( phys_sets.empty() )
         {
             MB_CHK_SET_ERR( MBI->create_meshset( MESHSET_SET, phys_set ), "can't create phys set" );
             MB_CHK_SET_ERR( MBI->tag_set_data( phys_tag, &phys_set, 1, &phys_table ), "can't set tag to phys set" );
         }
         else if( 1 == phys_sets.size() )
         {
             phys_set = phys_sets.front();
         }
         else
         {
             return MB_MULTIPLE_ENTITIES_FOUND;
         }
         MB_CHK_SET_ERR( MBI->add_entities( phys_set, &handle, 1 ), "can't add entities to phys set" );
  
         // If the material set does not already exist, create it.
         Range mat_sets;
         EntityHandle mat_set;
         const void* const mat_set_id_val[] = { &mat_table };
         rval = MBI->get_entities_by_type_and_tag( 0, MBENTITYSET, &mat_tag, mat_set_id_val, 1, mat_sets );
         if( MB_SUCCESS != rval ) return rval;
         if( mat_sets.empty() )
         {
             rval = MBI->create_meshset( MESHSET_SET, mat_set );
             if( MB_SUCCESS != rval ) return rval;
             rval = MBI->tag_set_data( mat_tag, &mat_set, 1, &mat_table );
             if( MB_SUCCESS != rval ) return rval;
         }
         else if( 1 == mat_sets.size() )
         {
             mat_set = mat_sets.front();
         }
         else
         {
             return MB_MULTIPLE_ENTITIES_FOUND;
         }
         rval = MBI->add_entities( mat_set, &handle, 1 );
         if( MB_SUCCESS != rval ) return rval;
  
         // Tag the element with its id
         MB_CHK_SET_ERR( MBI->tag_set_data( id_tag, &handle, 1, &element_id ), "Failed to assign IDs" );
         if( file_id_tag )
         {
             MB_CHK_SET_ERR( MBI->tag_set_data( *file_id_tag, &handle, 1, &element_id ), "Failed to assign file IDs" );
         }
     }
 }

References moab::Interface::add_entities(), moab::Interface::create_element(), moab::Interface::create_meshset(), moab::Range::empty(), ErrorCode, file, moab::Range::front(), moab::Interface::get_entities_by_type_and_tag(), moab::Interface::globalId_tag(), HEX, MAT_PROP_TABLE_TAG, moab::CN::MAX_NODES_PER_ELEMENT, MB_ALREADY_ALLOCATED, MB_CHK_SET_ERR, MB_MULTIPLE_ENTITIES_FOUND, MB_NOT_IMPLEMENTED, MB_SUCCESS, MB_TAG_CREAT, MB_TAG_DENSE, MB_TYPE_INTEGER, MBENTITYSET, MBHEX, MBI, MBPRISM, MBQUAD, MBTET, MBTRI, MESHSET_SET, PHYS_PROP_TABLE_TAG, QUAD0, QUAD1, moab::Range::size(), moab::Interface::tag_get_handle(), moab::Interface::tag_set_data(), TET, TRI0, TRI1, moab::CN::VerticesPerEntity(), and WEDGE.

Referenced by load_file().

◆ create_vertices()

ErrorCode moab::ReadIDEAS::create_vertices	(	EntityHandle &	first_vertex,
		const Tag *	file_id_tag
	)

protected

Definition at line 129 of file ReadIDEAS.cpp.

 {
     // Read two lines: first has some data, second has coordinates
     char line1[10000], line2[10000];
     int il1, il2;
     char *ctmp1, *ctmp2;
     std::string s1, s2;
  
     ErrorCode rval;
  
     std::streampos top_of_block = file.tellg();
     unsigned int num_verts      = 0;
  
     for( ;; )
     {
         // Read both lines
         if( !file.getline( line1, 10000 ) ) return MB_FAILURE;
         if( !file.getline( line2, 10000 ) ) return MB_FAILURE;
  
         // Check if we are at the end of the block
         il1 = std::strtol( line1, &ctmp1, 10 );
         il2 = std::strtol( line2, &ctmp2, 10 );
         if( ( il1 == -1 ) && ( il2 == -1 ) )
         {
             s1 = ctmp1;
             s2 = ctmp2;
             if( ( s1.empty() ) && ( s2.empty() ) ) break;
         }
         num_verts++;
     }
  
     file.seekg( top_of_block );
  
     std::vector< double* > arrays;
     rval = readMeshIface->get_node_coords( 3, num_verts, 0, first_vertex, arrays );
     if( MB_SUCCESS != rval ) return rval;
  
     Range verts;
     verts.insert( first_vertex, first_vertex + num_verts - 1 );
  
     double* x = arrays[0];
     double* y = arrays[1];
     double* z = arrays[2];
  
     // For now, assume ids are sequential and begin with 1
     Tag id_tag                  = MBI->globalId_tag();
     const int beginning_node_id = 1;
     int node_id                 = beginning_node_id;
  
     for( unsigned int i = 0; i < num_verts; i++ )
     {
         if( !file.getline( line1, 10000 ) ) return MB_FAILURE;
         if( !file.getline( line2, 10000 ) ) return MB_FAILURE;
  
         // Get the id out of the 1st line. Check the assumption that node ids are
         // sequential and begin with 1.
         if( node_id != std::strtol( line1, &ctmp1, 10 ) )
             MB_SET_ERR( MB_FAILURE, "node_id " << node_id << " line2:" << line2 << " ctmp1:" << ctmp1 );
         else
             ++node_id;
  
         // Get the doubles out of the 2nd line
         x[i] = std::strtod( line2, &ctmp2 );
         y[i] = std::strtod( ctmp2 + 1, &ctmp2 );
         z[i] = std::strtod( ctmp2 + 1, NULL );
     }
  
     if( !file.getline( line1, 10000 ) ) MB_SET_ERR( MB_FAILURE, " expect more lines" );
     if( !file.getline( line2, 10000 ) ) MB_SET_ERR( MB_FAILURE, " expect more lines 2" );
  
     // Tag the nodes with ids
     MB_CHK_SET_ERR( readMeshIface->assign_ids( id_tag, verts, beginning_node_id ), "Failed to assign IDs" );
     if( file_id_tag )
     {
         MB_CHK_SET_ERR( readMeshIface->assign_ids( *file_id_tag, verts, beginning_node_id ),
                         "Failed to assign file IDs" );
     }
  
     return MB_SUCCESS;
 }

References moab::ReadUtilIface::assign_ids(), ErrorCode, file, moab::ReadUtilIface::get_node_coords(), moab::Interface::globalId_tag(), moab::Range::insert(), MB_CHK_SET_ERR, MB_SET_ERR, MB_SUCCESS, MBI, and readMeshIface.

Referenced by load_file().

◆ factory()

ReaderIface * moab::ReadIDEAS::factory ( Interface * iface )

static

Definition at line 20 of file ReadIDEAS.cpp.

 {
     return new ReadIDEAS( iface );
 }

References iface, and ReadIDEAS().

Referenced by moab::ReaderWriterSet::ReaderWriterSet().

◆ load_file()

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

virtual

Load mesh from a file.

Method all readers must provide to import a mesh.

Parameters

file_name	The file to read.
file_set	Optional pointer to entity set representing file. If this is not NULL, reader may optionally tag the pointed-to set with format-specific meta-data.
subset_list	An optional struct pointer specifying the tags identifying entity sets to be read.
file_id_tag	If specified, reader should store for each entity it reads, a unique integer ID for this tag.

Author: Jason Kraftcheck

Implements moab::ReaderIface.

Definition at line 39 of file ReadIDEAS.cpp.

 {
     if( subset_list )
     {
         MB_SET_ERR( MB_UNSUPPORTED_OPERATION, "Reading subset of files not supported for IDEAS" );
     }
  
     file.open( fname );
     if( !file.good() )
     {
         MB_SET_ERR( MB_FILE_DOES_NOT_EXIST, "Failed to open file: " << fname );
     }
  
     ErrorCode rval;
  
     char line[10000];
     file.getline( line, 10000 );
     char* liter = line;
     while( *liter && isspace( *liter ) )
         ++liter;
     if( *liter != '-' ) return MB_FAILURE;
     ++liter;
     if( *liter != '1' ) return MB_FAILURE;
     while( *++liter )
         if( !isspace( *liter ) ) return MB_FAILURE;
  
     EntityHandle first_vertex = 0;
     while( !file.eof() )
     {
         file.getline( line, 10000 );
         unsigned int header_id = (unsigned int)strtol( line, NULL, 10 );
  
         // Create vertices
         if( DOUBLE_PRECISION_NODES0 == header_id || DOUBLE_PRECISION_NODES1 == header_id )
         {
             if( first_vertex )  // multiple vertex blocks?
                 return MB_FAILURE;
             MB_CHK_SET_ERR( create_vertices( first_vertex, file_id_tag ), "Failed to read vertices" );
         }
         // Create elements
         else if( ELEMENTS0 == header_id || ELEMENTS1 == header_id || ELEMENTS2 == header_id )
         {
             if( !first_vertex )  // Need to read vertices first
                 return MB_FAILURE;
             MB_CHK_SET_ERR( create_elements( first_vertex, file_id_tag ), "Failed to read elements" );
         }
         // Skip everything else
         else
         {
             rval = skip_header();
             if( MB_SUCCESS != rval ) return MB_FAILURE;
         }
     }
  
     file.close();
     return MB_SUCCESS;
 }

References create_elements(), create_vertices(), DOUBLE_PRECISION_NODES0, DOUBLE_PRECISION_NODES1, ELEMENTS0, ELEMENTS1, ELEMENTS2, ErrorCode, file, MB_CHK_SET_ERR, MB_FILE_DOES_NOT_EXIST, MB_SET_ERR, MB_SUCCESS, MB_UNSUPPORTED_OPERATION, and skip_header().

◆ read_tag_values()

ErrorCode moab::ReadIDEAS::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`
	)

virtual

Read tag values from a file.

Read the list if all integer tag values from the file for a tag that is a single integer value per entity.

Parameters

file_name	The file to read.
tag_name	The tag for which to read values
tag_values_out	Output: The list of tag values.
subset_list	An array of tag name and value sets specifying the subset of the file to read. If multiple tags are specified, the sets that match all tags (intersection) should be read.
subset_list_length	The length of the 'subset_list' array.

Implements moab::ReaderIface.

Definition at line 30 of file ReadIDEAS.cpp.

 {
     return MB_NOT_IMPLEMENTED;
 }

References MB_NOT_IMPLEMENTED.

◆ skip_header()

ErrorCode moab::ReadIDEAS::skip_header ( )

protected

Definition at line 101 of file ReadIDEAS.cpp.

 {
     // Go until finding a pair of -1 lines
     char* ctmp;
     char line[10000];
     std::string s;
  
     int end_of_block = 0;
  
     long int il;
  
     while( file.getline( line, 10000 ) )
     {
         il = std::strtol( line, &ctmp, 10 );
         if( il == -1 )
         {
             s = ctmp;
             if( s.empty() ) end_of_block++;
         }
         else
             end_of_block = 0;
  
         if( end_of_block >= 2 ) return MB_SUCCESS;
     }
  
     return MB_SUCCESS;
 }

References file, and MB_SUCCESS.

Referenced by load_file().

Member Data Documentation

◆ DOUBLE_PRECISION_NODES0

const unsigned moab::ReadIDEAS::DOUBLE_PRECISION_NODES0 = 781

staticprivate

Definition at line 82 of file ReadIDEAS.hpp.

Referenced by load_file().

◆ DOUBLE_PRECISION_NODES1

const unsigned moab::ReadIDEAS::DOUBLE_PRECISION_NODES1 = 2411

staticprivate

Definition at line 83 of file ReadIDEAS.hpp.

Referenced by load_file().

◆ ELEMENTS0

const unsigned moab::ReadIDEAS::ELEMENTS0 = 71

staticprivate

Definition at line 91 of file ReadIDEAS.hpp.

Referenced by load_file().

◆ ELEMENTS1

const unsigned moab::ReadIDEAS::ELEMENTS1 = 780

staticprivate

Definition at line 92 of file ReadIDEAS.hpp.

Referenced by load_file().

◆ ELEMENTS2

const unsigned moab::ReadIDEAS::ELEMENTS2 = 2412

staticprivate

Definition at line 93 of file ReadIDEAS.hpp.

Referenced by load_file().

◆ file

std::ifstream moab::ReadIDEAS::file

private

Definition at line 54 of file ReadIDEAS.hpp.

Referenced by create_elements(), create_vertices(), load_file(), and skip_header().

◆ HEX

const int moab::ReadIDEAS::HEX = 115

staticprivate

Definition at line 107 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ MBI

Interface* moab::ReadIDEAS::MBI

private

Definition at line 61 of file ReadIDEAS.hpp.

Referenced by create_elements(), and create_vertices().

◆ nodeIdMap

RangeMap< int, EntityHandle > moab::ReadIDEAS::nodeIdMap

private

Definition at line 55 of file ReadIDEAS.hpp.

◆ QUAD0

const int moab::ReadIDEAS::QUAD0 = 44

staticprivate

Definition at line 103 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ QUAD1

const int moab::ReadIDEAS::QUAD1 = 94

staticprivate

Definition at line 104 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ readMeshIface

ReadUtilIface* moab::ReadIDEAS::readMeshIface

private

Definition at line 58 of file ReadIDEAS.hpp.

Referenced by create_vertices(), and ReadIDEAS().

◆ ROD0

const int moab::ReadIDEAS::ROD0 = 11

staticprivate

Definition at line 99 of file ReadIDEAS.hpp.

◆ ROD1

const int moab::ReadIDEAS::ROD1 = 171

staticprivate

Definition at line 100 of file ReadIDEAS.hpp.

◆ SINGLE_PRECISION_NODES

const unsigned moab::ReadIDEAS::SINGLE_PRECISION_NODES = 15

staticprivate

Definition at line 76 of file ReadIDEAS.hpp.

◆ TET

const int moab::ReadIDEAS::TET = 111

staticprivate

Definition at line 105 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ TRI0

const int moab::ReadIDEAS::TRI0 = 41

staticprivate

Definition at line 101 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ TRI1

const int moab::ReadIDEAS::TRI1 = 91

staticprivate

Definition at line 102 of file ReadIDEAS.hpp.

Referenced by create_elements().

◆ WEDGE

const int moab::ReadIDEAS::WEDGE = 112

staticprivate

Definition at line 106 of file ReadIDEAS.hpp.

Referenced by create_elements().

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ReadIDEAS()

◆ ~ReadIDEAS()

Member Function Documentation

◆ create_elements()

◆ create_vertices()

◆ factory()

◆ load_file()

◆ read_tag_values()

◆ skip_header()

Member Data Documentation

◆ DOUBLE_PRECISION_NODES0

◆ DOUBLE_PRECISION_NODES1

◆ ELEMENTS0

◆ ELEMENTS1

◆ ELEMENTS2

◆ file

◆ HEX

◆ MBI

◆ nodeIdMap

◆ QUAD0

◆ QUAD1

◆ readMeshIface

◆ ROD0

◆ ROD1

◆ SINGLE_PRECISION_NODES

◆ TET

◆ TRI0

◆ TRI1

◆ WEDGE