moab::WriteGMV Class Reference

Output Exodus File for VERDE. More...

#include <WriteGMV.hpp>

Inheritance diagram for moab::WriteGMV:

Collaboration diagram for moab::WriteGMV:

Public Member Functions
	WriteGMV (Interface *impl)
	Constructor. More...
virtual	~WriteGMV ()
	Destructor. More...
ErrorCode	write_file (const char *filename, const bool overwite, const FileOptions &opts, const EntityHandle *output_sets, const int num_output_sets, const std::vector< std::string > &qa_list, const Tag *tag_list=NULL, int num_tags=0, int requested_dimension=3)
	Export mesh to a file. More...
ErrorCode	write_file (const char *file_name, const EntityHandle output_set, const int user_dimension=3, const bool mesh=true, const bool poly_mesh=true)
	writes out a mesh file More...
Public Member Functions inherited from moab::WriterIface
virtual	~WriterIface ()

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

Private Member Functions
ErrorCode	local_write_mesh (const char *file_name, const EntityHandle output_set, const int user_dimension, const bool mesh, const bool poly_mesh)

Private Attributes
Interface *	mbImpl
	interface instance More...
WriteUtilIface *	mWriteIface
Tag	mMaterialSetTag
	Cached tags for reading. Note that all these tags are defined when the core is initialized. More...
Tag	mDirichletSetTag
Tag	mNeumannSetTag
Tag	mHasMidNodesTag
Tag	mGeomDimensionTag
Tag	mGlobalIdTag

Static Private Attributes
static const char *	gmvTypeNames [MBMAXTYPE] = { "", "line", "tri", "quad", "", "tet", "pyramid", "prism", "", "hex", "", "" }

Detailed Description

Output Exodus File for VERDE.

Definition at line 41 of file WriteGMV.hpp.

Constructor & Destructor Documentation

WriteGMV()

moab::WriteGMV::WriteGMV

(

Interface *

impl

)

Constructor.

get and cache predefined tag handles

get and cache predefined tag handles

Definition at line 45 of file WriteGMV.cpp.

 : mbImpl( impl )
{
 assert( impl != NULL );
 
 impl->query_interface( mWriteIface );
 
 // initialize in case tag_get_handle fails below
 mMaterialSetTag = 0;
 mDirichletSetTag = 0;
 mNeumannSetTag = 0;
 mHasMidNodesTag = 0;
 mGeomDimensionTag = 0;
 mGlobalIdTag = 0;
 
 //! get and cache predefined tag handles
 // initialize in case tag_get_handle fails below
 //! get and cache predefined tag handles
 int negone = -1;
 impl->tag_get_handle( MATERIAL_SET_TAG_NAME, 1, MB_TYPE_INTEGER, mMaterialSetTag, MB_TAG_SPARSE | MB_TAG_CREAT,
 &negone );
 
 impl->tag_get_handle( DIRICHLET_SET_TAG_NAME, 1, MB_TYPE_INTEGER, mDirichletSetTag, MB_TAG_SPARSE | MB_TAG_CREAT,
 &negone );
 
 impl->tag_get_handle( NEUMANN_SET_TAG_NAME, 1, MB_TYPE_INTEGER, mNeumannSetTag, MB_TAG_SPARSE | MB_TAG_CREAT,
 &negone );
 
 mGlobalIdTag = impl->globalId_tag();
 
 int dum_val_array[] = { -1, -1, -1, -1 };
 impl->tag_get_handle( HAS_MID_NODES_TAG_NAME, 4, MB_TYPE_INTEGER, mHasMidNodesTag, MB_TAG_SPARSE | MB_TAG_CREAT,
 dum_val_array );
}

References DIRICHLET_SET_TAG_NAME, moab::Interface::globalId_tag(), HAS_MID_NODES_TAG_NAME, MATERIAL_SET_TAG_NAME, MB_TAG_CREAT, MB_TAG_SPARSE, MB_TYPE_INTEGER, mDirichletSetTag, mGeomDimensionTag, mGlobalIdTag, mHasMidNodesTag, mMaterialSetTag, mNeumannSetTag, mWriteIface, NEUMANN_SET_TAG_NAME, moab::Interface::query_interface(), and moab::Interface::tag_get_handle().

Referenced by factory().

~WriteGMV()

moab::WriteGMV::~WriteGMV ( )

virtual

Destructor.

Definition at line 79 of file WriteGMV.cpp.

{
 mbImpl->release_interface( mWriteIface );
}

References mbImpl, mWriteIface, and moab::Interface::release_interface().

Member Function Documentation

factory()

WriterIface * moab::WriteGMV::factory ( Interface *  iface )

static

Definition at line 40 of file WriteGMV.cpp.

{
 return new WriteGMV( iface );
}

References iface, and WriteGMV().

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

local_write_mesh()

ErrorCode moab::WriteGMV::local_write_mesh ( const char *  file_name, const EntityHandle  output_set, const int  user_dimension, const bool  mesh, const bool  poly_mesh  )

private

Definition at line 136 of file WriteGMV.cpp.

{
 std::ofstream ofile;
 ErrorCode result;
 
 if( mesh )
 {
 // need to insert ".gmv"
 std::string tmp_name( file_name );
 tmp_name += ".gmv";
 ofile.open( tmp_name.c_str() );
 }
 else if( poly_mesh )
 {
 // need to insert ".poly.gmv"
 std::string tmp_name( file_name );
 tmp_name += ".poly.gmv";
 ofile.open( tmp_name.c_str() );
 }
 
 ofile << "gmvinput ascii" << std::endl;
 
 // get elements to be output
 Range dum_range, elements, all_verts;
 EntityType otype;
 if( poly_mesh )
 {
 result = mbImpl->get_entities_by_type( output_set, MBPOLYGON, elements, true );
 if( MB_SUCCESS != result ) return result;
 }
 else
 {
 for( otype = CN::TypeDimensionMap[user_dimension].first; otype <= CN::TypeDimensionMap[user_dimension].second;
 otype++ )
 {
 if( otype == MBPOLYGON || otype == MBPOLYHEDRON ) continue;
 dum_range.clear();
 result = mbImpl->get_entities_by_type( output_set, otype, dum_range, true );
 if( MB_SUCCESS != result ) return result;
 
 std::copy( dum_range.begin(), dum_range.end(), range_inserter( elements ) );
 }
 }
 
 // gather the vertices in these elements
 result = mbImpl->get_adjacencies( elements, 0, false, all_verts, Interface::UNION );
 if( MB_SUCCESS != result ) return result;
 
 int num_verts = all_verts.size();
 
 // allocate coordinate arrays and put pointers to them in a list
 double* xcoord = new double[num_verts];
 double* ycoord = new double[num_verts];
 double* zcoord = new double[num_verts];
 std::vector< double* > coord_arrays;
 coord_arrays.push_back( xcoord );
 coord_arrays.push_back( ycoord );
 coord_arrays.push_back( zcoord );
 
 // fill them in, writing id tags at the same time
 result = mWriteIface->get_node_coords( 3, num_verts, all_verts, mGlobalIdTag, 1, coord_arrays );
 if( MB_SUCCESS != result ) return result;
 
 int i, j;
 
 //========================================
 // WRITE COORDINATE DATA TO FILE HERE
 
 ofile << "nodev " << num_verts << std::endl;
 for( i = 0; i < num_verts; i++ )
 ofile << xcoord[i] << " " << ycoord[i] << " " << zcoord[i] << std::endl;
 
 //========================================
 
 delete[] xcoord;
 delete[] ycoord;
 delete[] zcoord;
 
 // iterate over types in selected dimension
 
 std::vector< int > connect;
 std::vector< EntityHandle > connecth;
 
 if( mesh )
 {
 Range sub_range;
 
 ofile << "cells " << elements.size() << std::endl;
 
 for( otype = CN::TypeDimensionMap[user_dimension].first; otype <= CN::TypeDimensionMap[user_dimension].second;
 otype++ )
 {
 
 if( otype == MBPOLYGON || otype == MBPOLYHEDRON ) continue;
 
 // get the first element of this type in the range, and one past the last
 Range::iterator lower =
 Range::lower_bound( elements.begin(), elements.end(), CREATE_HANDLE( otype, MB_START_ID, i ) );
 Range::iterator upper =
 Range::lower_bound( elements.begin(), elements.end(), CREATE_HANDLE( otype + 1, MB_START_ID, i ) );
 
 if( lower == upper ) continue;
 
 // copy these elements into a subrange
 sub_range.clear();
 std::copy( lower, upper, range_inserter( sub_range ) );
 
 // make sure the connectivity array is big enough
 int verts_per = CN::VerticesPerEntity( otype );
 if( connect.size() < verts_per * sub_range.size() ) connect.resize( verts_per * sub_range.size() );
 
 // get the connectivity
 result = mWriteIface->get_element_connect( sub_range.size(), verts_per, mGlobalIdTag, sub_range,
 mGlobalIdTag, 1, &connect[0] );
 if( MB_SUCCESS != result ) return result;
 
 //========================================
 // WRITE CONNECTIVITY DATA TO FILE HERE
 
 for( i = 0; i < (int)sub_range.size(); i++ )
 {
 ofile << gmvTypeNames[otype] << " " << verts_per << std::endl;
 for( j = i * verts_per; j < (int)( i + 1 ) * verts_per; j++ )
 ofile << connect[j] << " ";
 ofile << std::endl;
 }
 
 //========================================
 }
 }
 
 else if( poly_mesh )
 {
 
 // write polygons/hedra, if any
 Range polygons, polyhedra;
 result = mbImpl->get_entities_by_type( output_set, MBPOLYGON, polygons, true );
 if( MB_SUCCESS != result ) return result;
 
 result = mbImpl->get_entities_by_type( output_set, MBPOLYHEDRON, polyhedra, true );
 if( MB_SUCCESS != result ) return result;
 
 if( polygons.size() == 0 ) return result;
 
 // mark polyhedra with global ids
 result = mWriteIface->assign_ids( polyhedra, mGlobalIdTag, 1 );
 if( MB_SUCCESS != result ) return result;
 
 ofile << "faces " << polygons.size() << " " << polyhedra.size() << std::endl;
 
 for( Range::iterator rit = polygons.begin(); rit != polygons.end(); ++rit )
 {
 // get the vertices
 connecth.clear();
 result = mbImpl->get_connectivity( &( *rit ), 1, connecth, true );
 if( MB_SUCCESS != result ) return result;
 
 if( 0 == connecth.size() ) continue;
 
 // get the polyhedra, if any
 if( user_dimension == 3 )
 {
 polyhedra.clear();
 result = mbImpl->get_adjacencies( Range( *rit, *rit ), 3, false, polyhedra );
 if( MB_SUCCESS != result ) return result;
 
 // put them in the connect array
 connecth.push_back( ( polyhedra.size() > 0 ? *polyhedra.begin() : 0 ) );
 connecth.push_back( ( polyhedra.size() > 1 ? *polyhedra.rbegin() : 0 ) );
 }
 
 // replace handles with ids
 connect.resize( connecth.size() + 2 );
 
 // pre-set polyhedra ids in case there aren't any
 connect[connecth.size()] = 0;
 connect[connecth.size() + 1] = 0;
 result =
 mbImpl->tag_get_data( mGlobalIdTag, &connecth[0], connecth.size() - 2 + polyhedra.size(), &connect[0] );
 if( MB_SUCCESS != result ) return result;
 
 // write the data
 ofile << connecth.size() - 2;
 
 for( i = 0; i < (int)connecth.size(); i++ )
 ofile << " " << connect[i];
 
 ofile << std::endl;
 }
 }
 
 ofile << std::endl << "endgmv" << std::endl;
 
 ofile.close();
 
 return MB_SUCCESS;
}

References moab::WriteUtilIface::assign_ids(), moab::Range::begin(), moab::Range::clear(), moab::CREATE_HANDLE(), moab::Range::end(), ErrorCode, moab::GeomUtil::first(), moab::Interface::get_adjacencies(), moab::Interface::get_connectivity(), moab::WriteUtilIface::get_element_connect(), moab::Interface::get_entities_by_type(), moab::WriteUtilIface::get_node_coords(), gmvTypeNames, moab::Range::lower_bound(), MB_START_ID, MB_SUCCESS, mbImpl, MBPOLYGON, MBPOLYHEDRON, mGlobalIdTag, mWriteIface, moab::Range::rbegin(), moab::Range::size(), moab::Interface::tag_get_data(), moab::CN::TypeDimensionMap, moab::Interface::UNION, and moab::CN::VerticesPerEntity().

Referenced by write_file().

write_file() [1/2]

ErrorCode moab::WriteGMV::write_file	(	const char *	file_name,
		const EntityHandle	output_set,
		const int	user_dimension = 3,
		const bool	mesh = true,
		const bool	poly_mesh = true
	)

writes out a mesh file

Definition at line 84 of file WriteGMV.cpp.

{
 // general function for writing a mesh
 
 ErrorCode result = MB_SUCCESS;
 
 // initialize file
 
 if( mesh )
 {
 result = local_write_mesh( file_name, output_set, user_dimension, true, false );
 if( MB_SUCCESS != result ) return result;
 }
 
 if( poly_mesh )
 {
 result = local_write_mesh( file_name, output_set, user_dimension, false, true );
 if( MB_SUCCESS != result ) return result;
 }
 
 return result;
}

References ErrorCode, local_write_mesh(), and MB_SUCCESS.

write_file() [2/2]

ErrorCode moab::WriteGMV::write_file ( const char *  file_name, const bool  overwrite, const FileOptions &  opts, const EntityHandle *  meshset_list, const int  num_sets, const std::vector< std::string > &  qa_records, const Tag *  tag_list = NULL, int  num_tags = 0, int  requested_output_dimension = 3  )

virtual

Export mesh to a file.

Method all writers must provide to export a mesh.

Parameters

file_name	The name of the file to create.
overwrite	If false, reader should fail if the file already exists.
meshset_list	A list of meshsets to export, or NULL if the whole mesh is to be exported.
num_sets	The length of meshset_list or zero if the whole mesh is to be exported.
qa_records	File history metadata
tag_list	Array of handles for tags to write. If null, write all tags. If non-NULL but num_tags is zero, write no tags.
requseted_output_dimension	The geometric dimension of the output mesh (coord values per vertex.) If zero, the dimension of the mesh as returned from Interface should be used.

Author

Jason Kraftcheck

Implements moab::WriterIface.

Definition at line 111 of file WriteGMV.cpp.

{
 EntityHandle output_set = 0;
 if( output_sets && num_output_sets > 0 )
 {
 if( num_output_sets > 1 ) return MB_FAILURE;
 output_set = output_sets[0];
 }
 
 if( dimension == 0 )
 {
 mbImpl->get_dimension( dimension );
 }
 
 return write_file( filename, output_set, dimension, true, true );
}

References moab::Interface::get_dimension(), and mbImpl.

Member Data Documentation

gmvTypeNames

const char * moab::WriteGMV::gmvTypeNames = { "", "line", "tri", "quad", "", "tet", "pyramid", "prism", "", "hex", "", "" }

staticprivate

Definition at line 85 of file WriteGMV.hpp.

Referenced by local_write_mesh().

mbImpl

Interface* moab::WriteGMV::mbImpl

private

interface instance

Definition at line 73 of file WriteGMV.hpp.

Referenced by local_write_mesh(), write_file(), and ~WriteGMV().

mDirichletSetTag

Tag moab::WriteGMV::mDirichletSetTag

private

Definition at line 79 of file WriteGMV.hpp.

Referenced by WriteGMV().

mGeomDimensionTag

Tag moab::WriteGMV::mGeomDimensionTag

private

Definition at line 82 of file WriteGMV.hpp.

Referenced by WriteGMV().

mGlobalIdTag

Tag moab::WriteGMV::mGlobalIdTag

private

Definition at line 83 of file WriteGMV.hpp.

Referenced by local_write_mesh(), and WriteGMV().

mHasMidNodesTag

Tag moab::WriteGMV::mHasMidNodesTag

private

Definition at line 81 of file WriteGMV.hpp.

Referenced by WriteGMV().

mMaterialSetTag

Tag moab::WriteGMV::mMaterialSetTag

private

Cached tags for reading. Note that all these tags are defined when the core is initialized.

Definition at line 78 of file WriteGMV.hpp.

Referenced by WriteGMV().

mNeumannSetTag

Tag moab::WriteGMV::mNeumannSetTag

private

Definition at line 80 of file WriteGMV.hpp.

Referenced by WriteGMV().

mWriteIface

WriteUtilIface* moab::WriteGMV::mWriteIface

private

Definition at line 74 of file WriteGMV.hpp.

Referenced by local_write_mesh(), WriteGMV(), and ~WriteGMV().

---

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

• WriteGMV.hpp
• WriteGMV.cpp