Mesh Oriented datABase  (version 5.5.1)
An array-based unstructured mesh library
LoadPartial.cpp

Load a part of a file
To run: LoadPartial [file] [tag_name] [val1] [val2] ...

In this example, it is shown how to load only a part of one file; the file must be organized in sets. (cherry-picking only the sets we want) The sets to load are identified by a tag name and the tag values for the sets of interest. This procedure is used when reading in parallel, as each processor will load only its part of the file, identified either by partition or by material/block sets by default, this example will load parallel partition sets with values 1, 2, and 5 from ../MeshFiles/unittest/64bricks_1khex.h5m The example will always write the output to a file name part.h5m

/** @example LoadPartial.cpp \n
* \brief Load a part of a file \n
* <b>To run</b>: LoadPartial [file] [tag_name] [val1] [val2] ...\n
*
* In this example, it is shown how to load only a part of one file; the file must be organized in
* sets. (cherry-picking only the sets we want) The sets to load are identified by a tag name and
* the tag values for the sets of interest. This procedure is used when reading in parallel, as
* each processor will load only its part of the file, identified either by partition or by
* material/block sets by default, this example will load parallel partition sets with values 1, 2,
* and 5 from ../MeshFiles/unittest/64bricks_1khex.h5m The example will always write the output to a
* file name part.h5m
*/
#include <iostream>
#include <vector>
// Include header for MOAB instance and tag conventions for
#include "moab/Core.hpp"
#define NTAGVALS 3
using namespace moab;
using namespace std;
// Function to parse input parameters
ErrorCode get_file_options( int argc, char** argv, string& filename, string& tagName, vector< int >& tagValues )
{
// Get mesh filename
if( argc > 1 )
filename = string( argv[1] );
else
filename = string( MESH_DIR ) + string( "/64bricks_512hex_256part.h5m" );
// Get tag selection options
if( argc > 2 )
tagName = string( argv[2] );
else
tagName = "USERTAG";
if( argc > 3 )
{
tagValues.resize( argc - 3, 0 );
for( int i = 3; i < argc; ++i )
tagValues[i - 3] = atoi( argv[i] );
}
else
{
for( unsigned i = 0; i < tagValues.size(); ++i )
tagValues[i] = 2 * i + 1;
}
if( argc > 1 && argc < 4 ) // print usage
cout << " usage is " << argv[0] << " <file> <tag_name> <value1> <value2> .. \n";
return MB_SUCCESS;
}
int main( int argc, char** argv )
{
// Get MOAB instance
Interface* mb = new( std::nothrow ) Core;
if( NULL == mb ) return 1;
std::string filename, tagname;
vector< int > tagvals( NTAGVALS ); // Allocate for a maximum of 5 tag values
ErrorCode rval = get_file_options( argc, argv, filename, tagname, tagvals );MB_CHK_ERR( rval );
#ifdef MOAB_HAVE_HDF5
// This file is in the mesh files directory
rval = mb->load_file( filename.c_str(), 0, 0, PARALLEL_PARTITION_TAG_NAME, tagvals.data(), (int)tagvals.size() );MB_CHK_SET_ERR( rval, "Failed to read" );
// If HANDLEID tag present, convert to long, and see what we read from file
Tag handleid_tag;
rval = mb->tag_get_handle( "HANDLEID", handleid_tag );
if( MB_SUCCESS == rval )
{
// Convert a few values for a few vertices
Range verts;
rval = mb->get_entities_by_type( 0, MBVERTEX, verts );MB_CHK_SET_ERR( rval, "Failed to get vertices" );
vector< long > valsTag( verts.size() );
rval = mb->tag_get_data( handleid_tag, verts, &valsTag[0] );
if( MB_SUCCESS == rval ) cout << "First 2 long values recovered: " << valsTag[0] << " " << valsTag[1] << "\n";
}
rval = mb->write_file( "part.h5m" );MB_CHK_SET_ERR( rval, "Failed to write partial file" );
cout << "Wrote successfully part.h5m.\n";
#else
std::cout << "Configure MOAB with HDF5 to build and use this example correctly.\n";
#endif
delete mb;
return 0;
}