#include <ReadSms.hpp>

Inheritance diagram for moab::ReadSms:

Collaboration diagram for moab::ReadSms:

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...

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

virtual	~ReadSms ()
	Destructor. More...

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

Static Public Member Functions
static ReaderIface *	factory (Interface *)
	factory method More...

Private Member Functions
ErrorCode	add_entities (EntityHandle start, EntityHandle count, const Tag *file_id_tag)

ErrorCode	load_file_impl (FILE file, const Tag file_id_tag)

ErrorCode	get_set (std::vector< EntityHandle > sets, int set_type, int set_id, Tag set_tag, EntityHandle &this_set, const Tag file_id_tag)

ErrorCode	read_parallel_info (FILE *file_ptr)

Private Attributes
ReadUtilIface *	readMeshIface

Interface *	mdbImpl
	interface instance More...

Tag	globalId

Tag	paramCoords

Tag	geomDimension

int	setId

Detailed Description

Definition at line 15 of file ReadSms.hpp.

Constructor & Destructor Documentation

◆ ReadSms()

ReadSms::ReadSms ( Interface * impl = NULL )

Constructor.

Definition at line 56 of file ReadSms.cpp.

                                   : mdbImpl( impl ), globalId( 0 ), paramCoords( 0 ), geomDimension( 0 ), setId( 0 )
 {
     mdbImpl->query_interface( readMeshIface );
 }

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

Referenced by factory().

◆ ~ReadSms()

ReadSms::~ReadSms ( )

virtual

Destructor.

Definition at line 61 of file ReadSms.cpp.

 {
     if( readMeshIface )
     {
         mdbImpl->release_interface( readMeshIface );
         readMeshIface = 0;
     }
 }

References mdbImpl, readMeshIface, and moab::Interface::release_interface().

Member Function Documentation

◆ add_entities()

ErrorCode ReadSms::add_entities	(	EntityHandle	start,
		EntityHandle	count,
		const Tag *	file_id_tag
	)

private

Definition at line 485 of file ReadSms.cpp.

 {
     if( !count || !file_id_tag ) return MB_FAILURE;
  
     Range range;
     range.insert( start, start + count - 1 );
     return readMeshIface->assign_ids( *file_id_tag, range, 1 );
 }

References moab::ReadUtilIface::assign_ids(), moab::Range::insert(), and readMeshIface.

Referenced by load_file_impl().

◆ factory()

ReaderIface * ReadSms::factory ( Interface * iface )

static

factory method

Definition at line 51 of file ReadSms.cpp.

 {
     return new ReadSms( iface );
 }

References iface, and ReadSms().

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

◆ get_set()

ErrorCode ReadSms::get_set	(	std::vector< EntityHandle > *	sets,
		int	set_type,
		int	set_id,
		Tag	set_tag,
		EntityHandle &	this_set,
		const Tag *	file_id_tag
	)

private

Definition at line 403 of file ReadSms.cpp.

 {
     ErrorCode result = MB_SUCCESS;
  
     if( set_dim < 0 || set_dim > 3 ) return MB_FILE_WRITE_ERROR;
  
     if( (int)sets[set_dim].size() <= set_id || !sets[set_dim][set_id] )
     {
         if( (int)sets[set_dim].size() <= set_id ) sets[set_dim].resize( set_id + 1, 0 );
  
         if( !sets[set_dim][set_id] )
         {
             result = mdbImpl->create_meshset( MESHSET_SET, sets[set_dim][set_id] );
             if( MB_SUCCESS != result ) return result;
             result = mdbImpl->tag_set_data( globalId, &sets[set_dim][set_id], 1, &set_id );
             if( MB_SUCCESS != result ) return result;
             result = mdbImpl->tag_set_data( dim_tag, &sets[set_dim][set_id], 1, &set_dim );
             if( MB_SUCCESS != result ) return result;
  
             if( file_id_tag )
             {
                 result = mdbImpl->tag_set_data( *file_id_tag, &sets[set_dim][set_id], 1, &setId );
                 ++setId;
             }
         }
     }
  
     this_set = sets[set_dim][set_id];
  
     return result;
 }

References moab::Interface::create_meshset(), ErrorCode, globalId, MB_FILE_WRITE_ERROR, MB_SUCCESS, mdbImpl, MESHSET_SET, setId, and moab::Interface::tag_set_data().

Referenced by load_file_impl().

◆ load_file()

ErrorCode ReadSms::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 79 of file ReadSms.cpp.

 {
     if( subset_list )
     {
         MB_SET_ERR( MB_UNSUPPORTED_OPERATION, "Reading subset of files not supported for Sms" );
     }
  
     setId = 1;
  
     // Open file
     FILE* file_ptr = fopen( filename, "r" );
     if( !file_ptr )
     {
         MB_SET_ERR( MB_FILE_DOES_NOT_EXIST, filename << ": " << strerror( errno ) );
     }
  
     const ErrorCode result = load_file_impl( file_ptr, file_id_tag );
     fclose( file_ptr );
  
     return result;
 }

References ErrorCode, load_file_impl(), MB_FILE_DOES_NOT_EXIST, MB_SET_ERR, MB_UNSUPPORTED_OPERATION, and setId.

◆ load_file_impl()

ErrorCode ReadSms::load_file_impl	(	FILE *	file,
		const Tag *	file_id_tag
	)

private

Definition at line 105 of file ReadSms.cpp.

 {
     bool warned         = false;
     double dum_params[] = { 0.0, 0.0, 0.0 };
  
     globalId = mdbImpl->globalId_tag();
  
     ErrorCode result =
         mdbImpl->tag_get_handle( "PARAMETER_COORDS", 3, MB_TYPE_DOUBLE, paramCoords, MB_TAG_DENSE | MB_TAG_CREAT );
     CHECK( "Failed to create param coords tag." );
  
     int negone = -1;
     result     = mdbImpl->tag_get_handle( GEOM_DIMENSION_TAG_NAME, 1, MB_TYPE_INTEGER, geomDimension,
                                           MB_TAG_SPARSE | MB_TAG_CREAT, &negone );
     CHECK( "Failed to create geom dim tag." );
  
     int n;
     char line[256], all_line[1024];
     int file_type;
  
     if( fgets( all_line, sizeof( all_line ), file_ptr ) == NULL )
     {
         return MB_FAILURE;
     }
     if( sscanf( all_line, "%s %d", line, &file_type ) != 2 )
     {
         return MB_FAILURE;
     }
  
     if( 3 == file_type )
     {
         result = read_parallel_info( file_ptr );
         CHECK( "Failed to read parallel info." );
     }
  
     int nregions, nfaces, nedges, nvertices, npoints;
     n = fscanf( file_ptr, "%d %d %d %d %d", &nregions, &nfaces, &nedges, &nvertices, &npoints );
     CHECKN( 5 );
     if( nregions < 0 || nfaces < 0 || nedges < 0 || nvertices < 0 || npoints < 0 ) return MB_FILE_WRITE_ERROR;
  
     // Create the vertices
     std::vector< double* > coord_arrays;
     EntityHandle vstart = 0;
     result              = readMeshIface->get_node_coords( 3, nvertices, MB_START_ID, vstart, coord_arrays );
     CHECK( "Failed to get node arrays." );
  
     if( file_id_tag )
     {
         result = add_entities( vstart, nvertices, file_id_tag );
         MB_CHK_ERR( result );
     }
  
     EntityHandle this_gent, new_handle;
     std::vector< EntityHandle > gentities[4];
     int gent_id, dum_int;
     int gent_type, num_connections;
  
     for( int i = 0; i < nvertices; i++ )
     {
         n = fscanf( file_ptr, "%d", &gent_id );
         CHECKN( 1 );
         if( !gent_id ) continue;
  
         n = fscanf( file_ptr, "%d %d %lf %lf %lf", &gent_type, &num_connections, coord_arrays[0] + i,
                     coord_arrays[1] + i, coord_arrays[2] + i );
         CHECKN( 5 );
  
         result = get_set( gentities, gent_type, gent_id, geomDimension, this_gent, file_id_tag );
         MB_CHK_ERR( result );
  
         new_handle = vstart + i;
         result     = mdbImpl->add_entities( this_gent, &new_handle, 1 );
         CHECK( "Adding vertex to geom set failed." );
  
         switch( gent_type )
         {
             case 1:
                 n = fscanf( file_ptr, "%le", dum_params );
                 CHECKN( 1 );
                 result = mdbImpl->tag_set_data( paramCoords, &new_handle, 1, dum_params );
                 CHECK( "Failed to set param coords tag for vertex." );
                 break;
             case 2:
                 n = fscanf( file_ptr, "%le %le %d", dum_params, dum_params + 1, &dum_int );
                 CHECKN( 3 );
                 dum_params[2] = dum_int;
                 result        = mdbImpl->tag_set_data( paramCoords, &new_handle, 1, dum_params );
                 CHECK( "Failed to set param coords tag for vertex." );
                 break;
             default:
                 break;
         }
     }  // End of reading vertices
  
     // *******************************
     //  Read Edges
     // *******************************
  
     int vert1, vert2, num_pts;
     std::vector< EntityHandle > everts( 2 );
     EntityHandle estart, *connect;
     result = readMeshIface->get_element_connect( nedges, 2, MBEDGE, 1, estart, connect );
     CHECK( "Failed to create array of edges." );
  
     if( file_id_tag )
     {
         result = add_entities( estart, nedges, file_id_tag );
         if( MB_SUCCESS != result ) return result;
     }
  
     for( int i = 0; i < nedges; i++ )
     {
         n = fscanf( file_ptr, "%d", &gent_id );
         CHECKN( 1 );
         if( !gent_id ) continue;
  
         n = fscanf( file_ptr, "%d %d %d %d %d", &gent_type, &vert1, &vert2, &num_connections, &num_pts );
         CHECKN( 5 );
         if( vert1 < 1 || vert1 > nvertices ) return MB_FILE_WRITE_ERROR;
         if( vert2 < 1 || vert2 > nvertices ) return MB_FILE_WRITE_ERROR;
  
         connect[0] = vstart + vert1 - 1;
         connect[1] = vstart + vert2 - 1;
         if( num_pts > 1 && !warned )
         {
             std::cout << "Warning: num_points > 1 not supported; choosing last one." << std::endl;
             warned = true;
         }
  
         result = get_set( gentities, gent_type, gent_id, geomDimension, this_gent, file_id_tag );
         CHECK( "Problem getting geom set for edge." );
  
         new_handle = estart + i;
         result     = mdbImpl->add_entities( this_gent, &new_handle, 1 );
         CHECK( "Failed to add edge to geom set." );
  
         connect += 2;
  
         for( int j = 0; j < num_pts; j++ )
         {
             switch( gent_type )
             {
                 case 1:
                     n = fscanf( file_ptr, "%le", dum_params );
                     CHECKN( 1 );
                     result = mdbImpl->tag_set_data( paramCoords, &new_handle, 1, dum_params );
                     CHECK( "Failed to set param coords tag for edge." );
                     break;
                 case 2:
                     n = fscanf( file_ptr, "%le %le %d", dum_params, dum_params + 1, &dum_int );
                     CHECKN( 3 );
                     dum_params[2] = dum_int;
                     result        = mdbImpl->tag_set_data( paramCoords, &new_handle, 1, dum_params );
                     CHECK( "Failed to set param coords tag for edge." );
                     break;
                 default:
                     break;
             }
         }
     }  // End of reading edges
  
     // *******************************
     //  Read Faces
     // *******************************
     std::vector< EntityHandle > bound_ents, bound_verts, new_faces;
     int bound_id;
     Range shverts;
     new_faces.resize( nfaces );
     int num_bounding;
  
     for( int i = 0; i < nfaces; i++ )
     {
         n = fscanf( file_ptr, "%d", &gent_id );
         CHECKN( 1 );
         if( !gent_id ) continue;
  
         n = fscanf( file_ptr, "%d %d", &gent_type, &num_bounding );
         CHECKN( 2 );
  
         result = get_set( gentities, gent_type, gent_id, geomDimension, this_gent, file_id_tag );
         CHECK( "Problem getting geom set for face." );
  
         bound_ents.resize( num_bounding + 1 );
         bound_verts.resize( num_bounding );
         for( int j = 0; j < num_bounding; j++ )
         {
             n = fscanf( file_ptr, "%d ", &bound_id );
             CHECKN( 1 );
             if( 0 > bound_id ) bound_id = abs( bound_id );
             assert( 0 < bound_id && bound_id <= nedges );
             if( bound_id < 1 || bound_id > nedges ) return MB_FILE_WRITE_ERROR;
             bound_ents[j] = estart + abs( bound_id ) - 1;
         }
  
         // Convert edge-based model to vertex-based one
         for( int j = 0; j < num_bounding; j++ )
         {
             if( j == num_bounding - 1 ) bound_ents[j + 1] = bound_ents[0];
             result = mdbImpl->get_adjacencies( &bound_ents[j], 2, 0, false, shverts );
             CHECK( "Failed to get vertices bounding edge." );
             assert( shverts.size() == 1 );
             bound_verts[j] = *shverts.begin();
             shverts.clear();
         }
  
         result = mdbImpl->create_element( (EntityType)( MBTRI + num_bounding - 3 ), &bound_verts[0], bound_verts.size(),
                                           new_faces[i] );
         CHECK( "Failed to create edge." );
  
         result = mdbImpl->add_entities( this_gent, &new_faces[i], 1 );
         CHECK( "Failed to add edge to geom set." );
  
         int num_read = fscanf( file_ptr, "%d", &num_pts );
         if( !num_pts || !num_read ) continue;
  
         for( int j = 0; j < num_pts; j++ )
         {
             switch( gent_type )
             {
                 case 1:
                     n = fscanf( file_ptr, "%le", dum_params );
                     CHECKN( 1 );
                     result = mdbImpl->tag_set_data( paramCoords, &new_faces[i], 1, dum_params );
                     CHECK( "Failed to set param coords tag for face." );
                     break;
                 case 2:
                     n = fscanf( file_ptr, "%le %le %d", dum_params, dum_params + 1, &dum_int );
                     CHECKN( 3 );
                     dum_params[2] = dum_int;
                     result        = mdbImpl->tag_set_data( paramCoords, &new_faces[i], 1, dum_params );
                     CHECK( "Failed to set param coords tag for face." );
                     break;
                 default:
                     break;
             }
         }
     }  // End of reading faces
  
     if( file_id_tag )
     {
         result = readMeshIface->assign_ids( *file_id_tag, &new_faces[0], new_faces.size(), 1 );
         if( MB_SUCCESS != result ) return result;
     }
  
     // *******************************
     //  Read Regions
     // *******************************
     int sense[MAX_SUB_ENTITIES];
     bound_verts.resize( MAX_SUB_ENTITIES );
  
     std::vector< EntityHandle > regions;
     if( file_id_tag ) regions.resize( nregions );
     for( int i = 0; i < nregions; i++ )
     {
         n = fscanf( file_ptr, "%d", &gent_id );
         CHECKN( 1 );
         if( !gent_id ) continue;
         result = get_set( gentities, 3, gent_id, geomDimension, this_gent, file_id_tag );
         CHECK( "Couldn't get geom set for region." );
         n = fscanf( file_ptr, "%d", &num_bounding );
         CHECKN( 1 );
         bound_ents.resize( num_bounding );
         for( int j = 0; j < num_bounding; j++ )
         {
             n = fscanf( file_ptr, "%d ", &bound_id );
             CHECKN( 1 );
             assert( abs( bound_id ) < (int)new_faces.size() + 1 && bound_id );
             if( !bound_id || abs( bound_id ) > nfaces ) return MB_FILE_WRITE_ERROR;
             sense[j]      = ( bound_id < 0 ) ? -1 : 1;
             bound_ents[j] = new_faces[abs( bound_id ) - 1];
         }
  
         EntityType etype;
         result = readMeshIface->get_ordered_vertices( &bound_ents[0], sense, num_bounding, 3, &bound_verts[0], etype );
         CHECK( "Failed in get_ordered_vertices." );
  
         // Make the element
         result = mdbImpl->create_element( etype, &bound_verts[0], CN::VerticesPerEntity( etype ), new_handle );
         CHECK( "Failed to create region." );
  
         result = mdbImpl->add_entities( this_gent, &new_handle, 1 );
         CHECK( "Failed to add region to geom set." );
  
         if( file_id_tag ) regions[i] = new_handle;
  
         n = fscanf( file_ptr, "%d ", &dum_int );
         CHECKN( 1 );
     }  // End of reading regions
  
     if( file_id_tag )
     {
         result = readMeshIface->assign_ids( *file_id_tag, &regions[0], regions.size(), 1 );
         if( MB_SUCCESS != result ) return result;
     }
  
     return MB_SUCCESS;
 }

References moab::Interface::add_entities(), add_entities(), moab::ReadUtilIface::assign_ids(), moab::Range::begin(), CHECK, CHECKN, moab::Range::clear(), moab::Interface::create_element(), ErrorCode, GEOM_DIMENSION_TAG_NAME, geomDimension, moab::Interface::get_adjacencies(), moab::ReadUtilIface::get_element_connect(), moab::ReadUtilIface::get_node_coords(), moab::ReadUtilIface::get_ordered_vertices(), get_set(), globalId, moab::Interface::globalId_tag(), moab::MAX_SUB_ENTITIES, MB_CHK_ERR, MB_FILE_WRITE_ERROR, MB_START_ID, MB_SUCCESS, MB_TAG_CREAT, MB_TAG_DENSE, MB_TAG_SPARSE, MB_TYPE_DOUBLE, MB_TYPE_INTEGER, MBEDGE, MBTRI, mdbImpl, paramCoords, read_parallel_info(), readMeshIface, moab::Range::size(), moab::Interface::tag_get_handle(), moab::Interface::tag_set_data(), and moab::CN::VerticesPerEntity().

Referenced by load_file().

◆ read_parallel_info()

ErrorCode ReadSms::read_parallel_info ( FILE * file_ptr )

private

Definition at line 440 of file ReadSms.cpp.

 {
     // ErrorCode result;
  
     // Read partition info
     int nparts, part_id, num_ifaces, num_corner_ents;
     int num_read = fscanf( file_ptr, "%d %d %d %d", &nparts, &part_id, &num_ifaces, &num_corner_ents );
     if( !num_read ) return MB_FAILURE;
  
     // Read interfaces
     int iface_id, iface_dim, iface_own, num_iface_corners;
     // EntityHandle this_iface;
     std::vector< int >* iface_corners = NULL;
     for( int i = 0; i < num_ifaces; i++ )
     {
         num_read = fscanf( file_ptr, "%d %d %d %d", &iface_id, &iface_dim, &iface_own, &num_iface_corners );
         if( !num_read ) return MB_FAILURE;
  
         // result = get_set(sets, iface_dim, iface_id, dim_tag, iface_own, this_iface);
         // CHECK("Failed to make iface set.");
  
         // Read the corner ids and store them on the set for now
         iface_corners = new std::vector< int >( num_iface_corners );
         for( int j = 0; j < num_iface_corners; j++ )
         {
             num_read = fscanf( file_ptr, "%d", &( *iface_corners )[j] );
             if( !num_read )
             {
                 delete iface_corners;
                 return MB_FAILURE;
             }
         }
  
         // result = tag_set_data(ifaceCornerTag, &this_iface, 1,
         //&iface_corners);
         // CHECK("Failed to set iface corner tag.");
  
         delete iface_corners;
         iface_corners = NULL;
     }
  
     // Interface data has been read
     return MB_SUCCESS;
 }

References MB_SUCCESS, and nparts.

Referenced by load_file_impl().

◆ read_tag_values()

ErrorCode ReadSms::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 70 of file ReadSms.cpp.

 {
     return MB_NOT_IMPLEMENTED;
 }

References MB_NOT_IMPLEMENTED.

Member Data Documentation

◆ geomDimension

Tag moab::ReadSms::geomDimension

private

Definition at line 59 of file ReadSms.hpp.

Referenced by load_file_impl().

◆ globalId

Tag moab::ReadSms::globalId

private

Definition at line 59 of file ReadSms.hpp.

Referenced by get_set(), and load_file_impl().

◆ mdbImpl

Interface* moab::ReadSms::mdbImpl

private

interface instance

Definition at line 57 of file ReadSms.hpp.

Referenced by get_set(), load_file_impl(), ReadSms(), and ~ReadSms().

◆ paramCoords

Tag moab::ReadSms::paramCoords

private

Definition at line 59 of file ReadSms.hpp.

Referenced by load_file_impl().

◆ readMeshIface

ReadUtilIface* moab::ReadSms::readMeshIface

private

Definition at line 54 of file ReadSms.hpp.

Referenced by add_entities(), load_file_impl(), ReadSms(), and ~ReadSms().

◆ setId

int moab::ReadSms::setId

private

Definition at line 61 of file ReadSms.hpp.

Referenced by get_set(), and load_file().

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ReadSms()

◆ ~ReadSms()

Member Function Documentation

◆ add_entities()

◆ factory()

◆ get_set()

◆ load_file()

◆ load_file_impl()

◆ read_parallel_info()

◆ read_tag_values()

Member Data Documentation

◆ geomDimension

◆ globalId

◆ mdbImpl

◆ paramCoords

◆ readMeshIface

◆ setId