moab::WriteAns Class Reference

#include <WriteAns.hpp>

Inheritance diagram for moab::WriteAns:

Collaboration diagram for moab::WriteAns:

Classes
struct	DirichletSetData
	struct used to hold data for each nodeset to be output; used by initialize_file to initialize the file header for increased speed More...
struct	MaterialSetData
	struct used to hold data for each block to be output; used by initialize_file to initialize the file header for increased speed More...
class	MeshInfo
	number of dimensions in this file More...
struct	NeumannSetData
	struct used to hold data for each sideset to be output; used by initialize_file to initialize the file header for increased speed More...

Public Member Functions
	WriteAns (Interface *impl)
	Constructor. More...
virtual	~WriteAns ()
	Destructor. More...
ErrorCode	write_file (const char *file_name, const bool overwrite, const FileOptions &opts, const EntityHandle *output_list, const int num_sets, const std::vector< std::string > &qa_list, const Tag *tag_list=NULL, int num_tags=0, int export_dimension=3)
	writes out a file More...
Public Member Functions inherited from moab::WriterIface
virtual	~WriterIface ()

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

Private Member Functions
ErrorCode	write_nodes (const int num_nodes, const Range &nodes, const int dimension, const char *file_name)

Private Attributes
Interface *	mbImpl
	interface instance More...
std::string	fileName
	file name More...
EntityHandle	mCurrentMeshHandle
	Meshset Handle for the mesh that is currently being read. More...
Tag	mMaterialSetTag
	Cached tags for reading. Note that all these tags are defined when the core is initialized. More...
Tag	mDirichletSetTag
Tag	mNeumannSetTag
Tag	mGlobalIdTag
Tag	mMatSetIdTag

Detailed Description

Definition at line 49 of file WriteAns.hpp.

Constructor & Destructor Documentation

WriteAns()

moab::WriteAns::WriteAns

(

Interface *

impl

)

Constructor.

get and cache predefined tag handles

Definition at line 51 of file WriteAns.cpp.

 : mbImpl( impl ), mCurrentMeshHandle( 0 ), mGlobalIdTag( 0 ), mMatSetIdTag( 0 )
{
 assert( impl != NULL );
 
 // impl->query_interface( mWriteIface );
 
 // initialize in case tag_get_handle fails below
 //! get and cache predefined tag handles
 const 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 );
}

References DIRICHLET_SET_TAG_NAME, MATERIAL_SET_TAG_NAME, MB_TAG_CREAT, MB_TAG_SPARSE, MB_TYPE_INTEGER, mDirichletSetTag, mMaterialSetTag, mNeumannSetTag, NEUMANN_SET_TAG_NAME, and moab::Interface::tag_get_handle().

Referenced by factory().

~WriteAns()

moab::WriteAns::~WriteAns ( )

virtual

Destructor.

Definition at line 70 of file WriteAns.cpp.

{
 // mbImpl->release_interface(mWriteIface);
}

Member Function Documentation

factory()

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

static

Definition at line 46 of file WriteAns.cpp.

{
 return new WriteAns( iface );
}

References iface, and WriteAns().

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

write_file()

ErrorCode moab::WriteAns::write_file ( const char *  file_name, const bool  overwrite, const FileOptions &  opts, const EntityHandle *  output_list, const int  num_sets, const std::vector< std::string > &  qa_list, const Tag *  tag_list = NULL, int  num_tags = 0, int  export_dimension = 3  )

virtual

writes out a file

Implements moab::WriterIface.

Definition at line 75 of file WriteAns.cpp.

{
 assert( 0 != mMaterialSetTag && 0 != mNeumannSetTag && 0 != mDirichletSetTag );
 
 ErrorCode result;
 
 // set SOLID45 element type to #60000, hope nobody has already...
 const char* ETSolid45 = "60045";
 const char* ETSolid92 = "60042";
 const char* ETSolid95 = "60095";
 
 // set Material id # to be used as default for all elements
 // will need to be subsequently reassigned inside ANSYS
 // Can, although have not, declare similar defaults for other attributes
 const char* MATDefault = "1";
 
 // create file streams for writing
 std::ofstream node_file;
 std::ofstream elem_file;
 std::ofstream ans_file;
 
 // get base filename from filename.ans
 std::string temp_string;
 std::string base_string;
 base_string.assign( file_name );
 base_string.replace( base_string.find_last_of( ".ans" ) - 3, 4, "" );
 
 // open node file for writing
 temp_string = base_string + ".node";
 node_file.open( temp_string.c_str() );
 node_file.setf( std::ios::scientific, std::ios::floatfield );
 node_file.precision( 13 );
 
 // open elem file for writing
 temp_string = base_string + ".elem";
 elem_file.open( temp_string.c_str() );
 
 // open ans file for writing
 ans_file.open( file_name );
 ans_file << "/prep7" << std::endl;
 
 // gather single output set
 EntityHandle output_set = 0;
 if( ent_handles && num_sets > 0 )
 {
 for( int i = 0; i < num_sets; i++ )
 {
 // from template, maybe can be removed
 result = mbImpl->unite_meshset( output_set, ent_handles[i] );
 if( result != MB_SUCCESS ) return result;
 }
 }
 
 // search for all nodes
 Range node_range;
 result = mbImpl->get_entities_by_type( output_set, MBVERTEX, node_range, true );
 if( result != MB_SUCCESS ) return result;
 
 // Commented out until Seg Fault taken care of in gather_nodes...
 // get any missing nodes which are needed for elements
 // Range all_ent_range,missing_range;
 // result=mbImpl->get_entities_by_handle(output_set,all_ent_range,true);
 // if(result !=MB_SUCCESS) return result;
 // result=mWriteIface->gather_nodes_from_elements(all_ent_range,0,missing_range);
 // node_range.merge(missing_range);
 
 // write the nodes
 double coord[3];
 for( Range::iterator it = node_range.begin(); it != node_range.end(); ++it )
 {
 EntityHandle node_handle = *it;
 
 result = mbImpl->get_coords( &node_handle, 1, coord );
 if( result != MB_SUCCESS ) return result;
 
 node_file.width( 8 );
 node_file << mbImpl->id_from_handle( node_handle );
 node_file.width( 20 );
 node_file << coord[0];
 node_file.width( 20 );
 node_file << coord[1];
 node_file.width( 20 );
 node_file << coord[2] << std::endl;
 }
 
 // update header to load nodes
 ans_file << "nread," << base_string << ",node" << std::endl;
 
 // search for all node sets (Dirichlet Sets)
 Range node_mesh_sets;
 int ns_id;
 result = mbImpl->get_entities_by_type_and_tag( 0, MBENTITYSET, &mDirichletSetTag, NULL, 1, node_mesh_sets );
 if( result != MB_SUCCESS ) return result;
 
 for( Range::iterator ns_it = node_mesh_sets.begin(); ns_it != node_mesh_sets.end(); ++ns_it )
 {
 result = mbImpl->tag_get_data( mDirichletSetTag, &( *ns_it ), 1, &ns_id );
 if( result != MB_SUCCESS ) return result;
 std::vector< EntityHandle > node_vector;
 result = mbImpl->get_entities_by_handle( *ns_it, node_vector, true );
 if( result != MB_SUCCESS ) return result;
 // for every nodeset found, cycle through nodes in set:
 for( std::vector< EntityHandle >::iterator node_it = node_vector.begin(); node_it != node_vector.end();
 ++node_it )
 {
 int ns_node_id = mbImpl->id_from_handle( *node_it );
 if( node_it == node_vector.begin() )
 {
 // select first node in new list
 ans_file << "nsel,s,node,," << std::setw( 8 ) << ns_node_id << std::endl;
 }
 else
 {
 // append node to list
 ans_file << "nsel,a,node,," << std::setw( 8 ) << ns_node_id << std::endl;
 }
 }
 // create NS(#) nodeset
 ans_file << "cm,NS" << ns_id << ",node" << std::endl;
 }
 
 // ANSYS Element format:
 // I, J, K, L, M, N, O, P,etc... MAT, TYPE, REAL, SECNUM, ESYS, IEL
 // I-P are nodes of element
 // MAT = material number
 // TYPE = Element type number
 // REAL = Real constant set number
 // SECNUM = section attribute number
 // ESYS = coordinate system for nodes
 // IEL = element # (unique?)
 // For all nodes past 8, write on second line
 
 // Write all MBTET elements
 Range tet_range;
 result = mbImpl->get_entities_by_type( output_set, MBTET, tet_range, true );
 if( result != MB_SUCCESS ) return result;
 for( Range::iterator elem_it = tet_range.begin(); elem_it != tet_range.end(); ++elem_it )
 {
 EntityHandle elem_handle = *elem_it;
 int elem_id = mbImpl->id_from_handle( elem_handle );
 std::vector< EntityHandle > conn;
 result = mbImpl->get_connectivity( &elem_handle, 1, conn, false );
 if( result != MB_SUCCESS ) return result;
 // make sure 4 or 10 node tet
 if( conn.size() != 4 && conn.size() != 10 )
 {
 std::cout << "Support not added for element type. \n";
 return MB_FAILURE;
 }
 // write information for 4 node tet
 if( conn.size() == 4 )
 {
 elem_file << std::setw( 8 ) << conn[0] << std::setw( 8 ) << conn[1];
 elem_file << std::setw( 8 ) << conn[2] << std::setw( 8 ) << conn[2];
 elem_file << std::setw( 8 ) << conn[3] << std::setw( 8 ) << conn[3];
 elem_file << std::setw( 8 ) << conn[3] << std::setw( 8 ) << conn[3];
 
 elem_file << std::setw( 8 ) << MATDefault << std::setw( 8 ) << ETSolid45;
 elem_file << std::setw( 8 ) << "1" << std::setw( 8 ) << "1";
 elem_file << std::setw( 8 ) << "0" << std::setw( 8 ) << elem_id;
 elem_file << std::endl;
 }
 
 // write information for 10 node tet
 if( conn.size() == 10 )
 {
 elem_file << std::setw( 8 ) << conn[0] << std::setw( 8 ) << conn[1];
 elem_file << std::setw( 8 ) << conn[2] << std::setw( 8 ) << conn[3];
 elem_file << std::setw( 8 ) << conn[4] << std::setw( 8 ) << conn[5];
 elem_file << std::setw( 8 ) << conn[6] << std::setw( 8 ) << conn[7];
 
 elem_file << std::setw( 8 ) << MATDefault << std::setw( 8 ) << ETSolid92;
 elem_file << std::setw( 8 ) << "1" << std::setw( 8 ) << "1";
 elem_file << std::setw( 8 ) << "0" << std::setw( 8 ) << elem_id;
 elem_file << std::endl;
 
 elem_file << std::setw( 8 ) << conn[8] << std::setw( 8 ) << conn[9];
 elem_file << std::endl;
 }
 }
 
 // Write all MBHEX elements
 Range hex_range;
 result = mbImpl->get_entities_by_type( output_set, MBHEX, hex_range, true );
 if( result != MB_SUCCESS ) return result;
 for( Range::iterator elem_it = hex_range.begin(); elem_it != hex_range.end(); ++elem_it )
 {
 EntityHandle elem_handle = *elem_it;
 int elem_id = mbImpl->id_from_handle( elem_handle );
 std::vector< EntityHandle > conn;
 result = mbImpl->get_connectivity( &elem_handle, 1, conn, false );
 if( result != MB_SUCCESS ) return result;
 // make sure supported hex type
 if( conn.size() != 8 && conn.size() != 20 )
 {
 std::cout << "Support not added for element type. \n";
 return MB_FAILURE;
 }
 
 // write information for 8 node hex
 if( conn.size() == 8 )
 {
 elem_file << std::setw( 8 ) << conn[0] << std::setw( 8 ) << conn[1];
 elem_file << std::setw( 8 ) << conn[2] << std::setw( 8 ) << conn[3];
 elem_file << std::setw( 8 ) << conn[4] << std::setw( 8 ) << conn[5];
 elem_file << std::setw( 8 ) << conn[6] << std::setw( 8 ) << conn[7];
 
 elem_file << std::setw( 8 ) << MATDefault << std::setw( 8 ) << ETSolid45;
 elem_file << std::setw( 8 ) << "1" << std::setw( 8 ) << "1";
 elem_file << std::setw( 8 ) << "0" << std::setw( 8 ) << elem_id;
 elem_file << std::endl;
 }
 
 // write information for 20 node hex
 if( conn.size() == 20 )
 {
 
 elem_file << std::setw( 8 ) << conn[4] << std::setw( 8 ) << conn[5];
 elem_file << std::setw( 8 ) << conn[1] << std::setw( 8 ) << conn[0];
 elem_file << std::setw( 8 ) << conn[7] << std::setw( 8 ) << conn[6];
 elem_file << std::setw( 8 ) << conn[2] << std::setw( 8 ) << conn[3];
 
 elem_file << std::setw( 8 ) << MATDefault << std::setw( 8 ) << ETSolid95;
 elem_file << std::setw( 8 ) << "1" << std::setw( 8 ) << "1";
 elem_file << std::setw( 8 ) << "0" << std::setw( 8 ) << elem_id;
 elem_file << std::endl;
 
 elem_file << std::setw( 8 ) << conn[16] << std::setw( 8 ) << conn[13];
 elem_file << std::setw( 8 ) << conn[8] << std::setw( 8 ) << conn[12];
 elem_file << std::setw( 8 ) << conn[18] << std::setw( 8 ) << conn[14];
 elem_file << std::setw( 8 ) << conn[10] << std::setw( 8 ) << conn[15];
 elem_file << std::setw( 8 ) << conn[19] << std::setw( 8 ) << conn[17];
 elem_file << std::setw( 8 ) << conn[9] << std::setw( 8 ) << conn[11];
 elem_file << std::endl;
 }
 }
 // Write all MBPRISM elements
 Range prism_range;
 result = mbImpl->get_entities_by_type( output_set, MBPRISM, prism_range, true );
 if( result != MB_SUCCESS ) return result;
 for( Range::iterator elem_it = prism_range.begin(); elem_it != prism_range.end(); ++elem_it )
 {
 EntityHandle elem_handle = *elem_it;
 int elem_id = mbImpl->id_from_handle( elem_handle );
 std::vector< EntityHandle > conn;
 result = mbImpl->get_connectivity( &elem_handle, 1, conn, false );
 if( result != MB_SUCCESS ) return result;
 // make sure supported prism type
 if( conn.size() != 6 )
 {
 std::cout << "Support not added for element type. \n";
 return MB_FAILURE;
 }
 
 // write information for 6 node prism
 if( conn.size() == 6 )
 {
 elem_file << std::setw( 8 ) << conn[0] << std::setw( 8 ) << conn[3];
 elem_file << std::setw( 8 ) << conn[4] << std::setw( 8 ) << conn[4];
 elem_file << std::setw( 8 ) << conn[1] << std::setw( 8 ) << conn[2];
 elem_file << std::setw( 8 ) << conn[5] << std::setw( 8 ) << conn[5];
 
 elem_file << std::setw( 8 ) << MATDefault << std::setw( 8 ) << ETSolid45;
 elem_file << std::setw( 8 ) << "1" << std::setw( 8 ) << "1";
 elem_file << std::setw( 8 ) << "0" << std::setw( 8 ) << elem_id;
 elem_file << std::endl;
 }
 }
 
 // create element types (for now writes all, even if not used)
 ans_file << "et," << ETSolid45 << ",SOLID45" << std::endl;
 ans_file << "et," << ETSolid92 << ",SOLID92" << std::endl;
 ans_file << "et," << ETSolid95 << ",SOLID95" << std::endl;
 
 // xxx pyramids, other elements later...
 
 // write header to load elements
 ans_file << "eread," << base_string << ",elem" << std::endl;
 
 // search for all side sets (Neumann)
 Range side_mesh_sets;
 int ss_id;
 result = mbImpl->get_entities_by_type_and_tag( 0, MBENTITYSET, &mNeumannSetTag, NULL, 1, side_mesh_sets );
 if( result != MB_SUCCESS ) return result;
 // cycle through all sets found
 for( Range::iterator ss_it = side_mesh_sets.begin(); ss_it != side_mesh_sets.end(); ++ss_it )
 {
 result = mbImpl->tag_get_data( mNeumannSetTag, &( *ss_it ), 1, &ss_id );
 if( result != MB_SUCCESS ) return result;
 std::vector< EntityHandle > elem_vector;
 result = mbImpl->get_entities_by_handle( *ss_it, elem_vector, true );
 if( result != MB_SUCCESS ) return result;
 
 // cycle through elements in current side set
 for( std::vector< EntityHandle >::iterator elem_it = elem_vector.begin(); elem_it != elem_vector.end();
 ++elem_it )
 {
 EntityHandle elem_handle = *elem_it;
 
 // instead of selecting current element in set, select its nodes...
 std::vector< EntityHandle > conn;
 result = mbImpl->get_connectivity( &elem_handle, 1, conn );
 if( result != MB_SUCCESS ) return result;
 if( elem_it == elem_vector.begin() )
 {
 ans_file << "nsel,s,node,," << std::setw( 8 ) << conn[0] << std::endl;
 for( unsigned int i = 1; i < conn.size(); i++ )
 {
 ans_file << "nsel,a,node,," << std::setw( 8 ) << conn[i] << std::endl;
 }
 }
 else
 {
 for( unsigned int i = 0; i < conn.size(); i++ )
 {
 ans_file << "nsel,a,node,," << std::setw( 8 ) << conn[i] << std::endl;
 }
 }
 }
 // create SS(#) node set
 ans_file << "cm,SS" << ss_id << ",node" << std::endl;
 }
 
 // Gather all element blocks
 Range matset;
 int mat_id;
 result = mbImpl->get_entities_by_type_and_tag( 0, MBENTITYSET, &mMaterialSetTag, NULL, 1, matset );
 if( result != MB_SUCCESS ) return result;
 // cycle through all elem blocks
 for( Range::iterator mat_it = matset.begin(); mat_it != matset.end(); ++mat_it )
 {
 EntityHandle matset_handle = *mat_it;
 result = mbImpl->tag_get_data( mMaterialSetTag, &matset_handle, 1, &mat_id );
 if( result != MB_SUCCESS ) return result;
 std::vector< EntityHandle > mat_vector;
 result = mbImpl->get_entities_by_handle( *mat_it, mat_vector, true );
 if( result != MB_SUCCESS ) return result;
 // cycle through elements in current mat set
 for( std::vector< EntityHandle >::iterator elem_it = mat_vector.begin(); elem_it != mat_vector.end();
 ++elem_it )
 {
 EntityHandle elem_handle = *elem_it;
 int elem_id = mbImpl->id_from_handle( elem_handle );
 if( elem_it == mat_vector.begin() )
 {
 ans_file << "esel,s,elem,," << std::setw( 8 ) << elem_id << std::endl;
 }
 else
 {
 ans_file << "esel,a,elem,," << std::setw( 8 ) << elem_id << std::endl;
 }
 }
 // for each matset, write block command
 ans_file << "cm,EB" << mat_id << ",elem" << std::endl;
 }
 
 // close all file streams
 node_file.close();
 elem_file.close();
 ans_file.close();
 
 return MB_SUCCESS;
}

References moab::Range::begin(), moab::Range::end(), ErrorCode, moab::Interface::get_connectivity(), moab::Interface::get_coords(), moab::Interface::get_entities_by_handle(), moab::Interface::get_entities_by_type(), moab::Interface::get_entities_by_type_and_tag(), moab::Interface::id_from_handle(), MB_SUCCESS, MBENTITYSET, MBHEX, mbImpl, MBPRISM, MBTET, MBVERTEX, mDirichletSetTag, mMaterialSetTag, mNeumannSetTag, moab::Interface::tag_get_data(), and moab::Interface::unite_meshset().

write_nodes()

ErrorCode moab::WriteAns::write_nodes ( const int  num_nodes, const Range &  nodes, const int  dimension, const char *  file_name  )

private

Member Data Documentation

fileName

std::string moab::WriteAns::fileName

private

file name

Definition at line 137 of file WriteAns.hpp.

mbImpl

Interface* moab::WriteAns::mbImpl

private

interface instance

Definition at line 133 of file WriteAns.hpp.

Referenced by write_file().

mCurrentMeshHandle

EntityHandle moab::WriteAns::mCurrentMeshHandle

private

Meshset Handle for the mesh that is currently being read.

Definition at line 140 of file WriteAns.hpp.

mDirichletSetTag

Tag moab::WriteAns::mDirichletSetTag

private

Definition at line 145 of file WriteAns.hpp.

Referenced by write_file(), and WriteAns().

mGlobalIdTag

Tag moab::WriteAns::mGlobalIdTag

private

Definition at line 147 of file WriteAns.hpp.

mMaterialSetTag

Tag moab::WriteAns::mMaterialSetTag

private

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

Definition at line 144 of file WriteAns.hpp.

Referenced by write_file(), and WriteAns().

mMatSetIdTag

Tag moab::WriteAns::mMatSetIdTag

private

Definition at line 148 of file WriteAns.hpp.

mNeumannSetTag

Tag moab::WriteAns::mNeumannSetTag

private

Definition at line 146 of file WriteAns.hpp.

Referenced by write_file(), and WriteAns().

---

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

• WriteAns.hpp
• WriteAns.cpp