|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
#include <cstdio>#include <cmath>#include "moab/Core.hpp"#include "moab/Range.hpp"#include "MBTagConventions.hpp"#include "moab/CN.hpp"#include <iostream>#include <sstream>
Include dependency graph for h5file_test.cpp:Go to the source code of this file.
Defines | |
| #define | filename "h5test.h5m" |
| #define | tagname "foo/\\/\\" |
| #define | bitname "bar\n" |
| #define | intname "int tag" |
| #define | dblname " dbl " |
| #define | handlename "hanlde" |
Functions | |
| void | create () |
| bool | compare () |
| void | moab_error (const char *function) |
| int | main (int argc, char *argv[]) |
| EntityHandle | vtx (double x, double y, double z) |
| EntityHandle | pent (EntityHandle *vtx_list, int i1, int i2, int i3, int i4, int i5) |
| EntityHandle | make_set (unsigned int options, EntityHandle *entities, size_t num_entities, bool reverse, int id) |
| bool | compare_conn (std::vector< EntityHandle > &conn1, std::vector< EntityHandle > &conn2) |
| bool | compare_sets (int id, const char *tag_name=0) |
| bool | compare_tags (EntityHandle dod[]) |
Variables | |
| const int | FACE_SET_ID = 1101 |
| const int | VERTEX_SET_ID = 1102 |
| const int | REGION_SET_ID = 1103 |
| const int | EMPTY_SET_ID = 1100 |
| const int | SET_SET_ID = 1105 |
| Interface * | iface |
| #define bitname "bar\n" |
Definition at line 34 of file h5file_test.cpp.
| #define dblname " dbl " |
Definition at line 36 of file h5file_test.cpp.
Referenced by compare_tags(), and create().
| #define filename "h5test.h5m" |
MOAB, a Mesh-Oriented datABase, is a software component for creating, storing and accessing finite element mesh data.
Copyright 2004 Sandia Corporation. Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains certain rights in this software.
This library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version.
Definition at line 24 of file h5file_test.cpp.
Referenced by do_file(), do_ray_fire_test(), iMesh_load(), iMesh_save(), moab::ReadCGNS::load_file(), moab::ReadTemplate::load_file(), moab::ReadTetGen::load_file(), main(), mesh_size_test(), moab::ReadTetGen::open_file(), parse_options(), read_write(), read_write_file(), save_entset_test(), signal_handler(), split_test_across(), test_correct_ghost(), test_eul_check_across_files(), test_eul_check_append(), test_eul_check_T(), test_eul_check_timestep(), test_fv_check_T(), test_gcrm_check_vars(), test_homme_check_T(), test_mpas_check_vars(), test_native_read(), test_read_and_ghost_after(), test_read_with_ghost(), test_read_with_ghost_no_augment(), test_read_with_thin_ghost_layer(), test_sets_fileids(), test_tags_retrieval(), test_var_length_parallel(), test_write_invalid_elem(), test_write_read_many_tags(), and verdict_test1().
| #define handlename "hanlde" |
Definition at line 37 of file h5file_test.cpp.
Referenced by compare_tags(), and create().
| #define intname "int tag" |
Definition at line 35 of file h5file_test.cpp.
Referenced by compare_tags(), and create().
| #define tagname "foo/\\/\\" |
Definition at line 33 of file h5file_test.cpp.
Referenced by compare(), compare_conn(), create(), create_mesh(), moab::AdaptiveKDTree::init(), main(), mb_common_tag_test(), test_delete_tag_data(), test_delete_type_tag(), test_var_length_parallel(), and test_write_dense_tags().
| bool compare | ( | ) |
Definition at line 714 of file h5file_test.cpp.
References moab::Range::begin(), bitname, moab::Range::clear(), compare_conn(), compare_sets(), compare_tags(), EMPTY_SET_ID, FACE_SET_ID, moab::Interface::get_connectivity(), moab::Interface::get_entities_by_type_and_tag(), GLOBAL_ID_TAG_NAME, iface, MB_SUCCESS, MB_TYPE_BIT, MBHEX, MBPOLYHEDRON, moab_error(), REGION_SET_ID, SET_SET_ID, size, moab::Range::size(), moab::Interface::tag_get_handle(), tagname, and VERTEX_SET_ID.
Referenced by moab::Bvh_tree< _Entity_handles, _Box, _Moab, _Parametrizer >::choose_best_split(), moab::find_map_values(), moab::find_tag_values(), main(), and moab::ReadRTT::split_string().
{
Range range;
EntityHandle hex[2];
EntityHandle dod[2];
Tag elemtag;
Range::iterator iter;
// get tag
if( MB_SUCCESS != iface->tag_get_handle( bitname, 2, MB_TYPE_BIT, elemtag ) ) moab_error( "tag_get_handle" );
// get two hexes
char two = '\002';
const void* tarray[] = { &two };
if( MB_SUCCESS != iface->get_entities_by_type_and_tag( 0, MBHEX, &elemtag, tarray, 1, range ) )
moab_error( "get_entities_by_type_and_tag" );
if( range.size() != 2 )
{
fprintf( stderr, "Expected 2 Hexes. Got %lu\n", (unsigned long)range.size() );
exit( 1 );
}
iter = range.begin();
hex[0] = *iter;
hex[1] = *++iter;
// get two polyhedra
range.clear();
char one = '\001';
const void* oarray[] = { &one };
if( MB_SUCCESS != iface->get_entities_by_type_and_tag( 0, MBPOLYHEDRON, &elemtag, oarray, 1, range ) )
moab_error( "get_entities_by_type_and_tag" );
if( range.size() != 2 )
{
fprintf( stderr, "Expected 2 Polyhedra. Got %lu\n", (unsigned long)range.size() );
exit( 1 );
}
iter = range.begin();
dod[0] = *iter;
dod[1] = *++iter;
// compare hexes
std::vector< EntityHandle > conn[2];
if( MB_SUCCESS != iface->get_connectivity( hex, 1, conn[0] ) ||
MB_SUCCESS != iface->get_connectivity( hex + 1, 1, conn[1] ) )
moab_error( "get_connectivity" );
if( !compare_conn( conn[0], conn[1] ) ) return false;
// compare polyhedra
std::vector< EntityHandle > face[2];
conn[0].clear();
conn[1].clear();
if( MB_SUCCESS != iface->get_connectivity( dod, 1, conn[0], false ) ||
MB_SUCCESS != iface->get_connectivity( dod + 1, 1, conn[1], false ) )
moab_error( "get_connectivity" );
if( conn[0].size() != 12 || conn[1].size() != 12 )
{
fprintf( stderr,
"Expected two dodecahedrons. Got polyhedrons with "
"%lu and %lu faces respectively.\n",
(unsigned long)conn[0].size(), (unsigned long)conn[1].size() );
return false;
}
for( int i = 0; i < 12; ++i )
{
face[0].clear();
face[1].clear();
if( MB_SUCCESS != iface->get_connectivity( &conn[0][i], 1, face[0], false ) ||
MB_SUCCESS != iface->get_connectivity( &conn[1][i], 1, face[1], false ) )
moab_error( "get_connectivity" );
if( !compare_conn( face[0], face[1] ) ) return false;
}
// compare sets
if( !compare_sets( VERTEX_SET_ID, tagname ) || !compare_sets( FACE_SET_ID ) || !compare_sets( REGION_SET_ID ) ||
!compare_sets( EMPTY_SET_ID ) || !compare_sets( SET_SET_ID, GLOBAL_ID_TAG_NAME ) )
return false;
// check tags
if( !compare_tags( dod ) ) return false;
return true;
}
| bool compare_conn | ( | std::vector< EntityHandle > & | conn1, |
| std::vector< EntityHandle > & | conn2 | ||
| ) |
Definition at line 367 of file h5file_test.cpp.
References moab::Interface::get_coords(), iface, MB_SUCCESS, MB_TYPE_INTEGER, moab_error(), size, moab::Interface::tag_get_data(), moab::Interface::tag_get_handle(), and tagname.
Referenced by compare().
{
unsigned i;
if( conn1.size() != conn2.size() || conn1.size() == 0 )
{
fprintf( stderr, "Error comparing connectivity: sizes %lu and %lu\n", (unsigned long)conn1.size(),
(unsigned long)conn2.size() );
return false;
}
std::vector< double > coords[2];
coords[0].resize( 3 * conn1.size() );
coords[1].resize( 3 * conn2.size() );
if( MB_SUCCESS != iface->get_coords( &conn1[0], conn1.size(), &coords[0][0] ) ||
MB_SUCCESS != iface->get_coords( &conn2[0], conn2.size(), &coords[1][0] ) ||
coords[0].size() != coords[1].size() )
moab_error( "get_coords" );
std::vector< double >::iterator citer1 = coords[0].begin(), citer2 = coords[1].begin();
for( i = 0; i < conn1.size(); i++ )
{
double x1 = *( citer1++ ), y1 = *( citer1++ ), z1 = *( citer1++ );
double x2 = *( citer2++ ), y2 = *( citer2++ ), z2 = *( citer2++ );
if( x1 != x2 || y1 != y2 || z1 != z2 )
{
fprintf( stderr, "Vertex coords don't match: ( %f, %f, %f ) and ( %f, %f, %f )\n", x1, y1, z1, x2, y2, z2 );
return false;
}
}
std::vector< int > tags[2];
tags[0].resize( conn1.size() );
tags[1].resize( conn2.size() );
Tag tag;
if( MB_SUCCESS != iface->tag_get_handle( tagname, 1, MB_TYPE_INTEGER, tag ) ) moab_error( "tag_get_handle" );
if( MB_SUCCESS != iface->tag_get_data( tag, &conn1[0], conn1.size(), &tags[0][0] ) ||
MB_SUCCESS != iface->tag_get_data( tag, &conn2[0], conn2.size(), &tags[1][0] ) )
moab_error( "tag_get_data" );
std::vector< int >::iterator titer1 = tags[0].begin(), titer2 = tags[1].begin();
for( i = 0; i < conn1.size(); i++ )
{
int t1 = *( titer1++ );
int t2 = *( titer2++ );
if( t1 != t2 )
{
fprintf( stderr, "Vertex tags don't match: %d != %d\n", t1, t2 );
return false;
}
}
return true;
}
| bool compare_sets | ( | int | id, |
| const char * | tag_name = 0 |
||
| ) |
Definition at line 427 of file h5file_test.cpp.
References moab::Range::begin(), moab::CN::EntityTypeName(), ErrorCode, moab::Interface::get_child_meshsets(), moab::Interface::get_entities_by_handle(), moab::Interface::get_entities_by_type_and_tag(), moab::Interface::get_meshset_options(), moab::Interface::get_parent_meshsets(), moab::Interface::globalId_tag(), id_tag, iface, MB_ENTITY_NOT_FOUND, MB_SUCCESS, MB_TAG_ANY, MB_TAG_NOT_FOUND, MB_TYPE_INTEGER, MB_TYPE_OPAQUE, MBENTITYSET, MBMAXTYPE, MESHSET_SET, MESHSET_TRACK_OWNER, moab_error(), moab::Range::size(), moab::Interface::tag_get_data(), moab::Interface::tag_get_data_type(), moab::Interface::tag_get_handle(), and moab::Interface::type_from_handle().
Referenced by compare().
{
bool ok;
ErrorCode rval;
// get sets
Tag id_tag = iface->globalId_tag();
Range range;
const void* tag_data[] = { &id };
rval = iface->get_entities_by_type_and_tag( 0, MBENTITYSET, &id_tag, tag_data, 1, range );
if( MB_ENTITY_NOT_FOUND == rval || range.size() != 2 )
{
fprintf( stderr, "Could not find set with id %d pair in file\n", id );
return false;
}
else if( MB_SUCCESS != rval )
moab_error( "get_entities_by_type_and_tag" );
EntityHandle set1 = *range.begin();
EntityHandle set2 = *++range.begin();
// Compare set descriptions
unsigned opt1 = 0, opt2 = 0;
rval = iface->get_meshset_options( set1, opt1 );
if( MB_SUCCESS != rval ) moab_error( "get_meshset_options for set1 failed" );
rval = iface->get_meshset_options( set2, opt2 );
if( MB_SUCCESS != rval ) moab_error( "get_meshset_options for set2 failed" );
if( opt1 != opt2 )
{
fprintf( stderr,
"Sets with id %d do not have matching options.\n"
"Set 1: track_owner=%s set=%s ordered=%s\n"
"Set 2: track_owner=%s set=%s ordered=%s\n",
id, opt1 & MESHSET_TRACK_OWNER ? "yes" : "no", opt1 & MESHSET_SET ? "yes" : "no",
opt1 & MESHSET_ORDERED ? "yes" : "no", opt2 & MESHSET_TRACK_OWNER ? "yes" : "no",
opt2 & MESHSET_SET ? "yes" : "no", opt2 & MESHSET_ORDERED ? "yes" : "no" );
return false;
}
// Compare set contents
// First check if same number of entities.
// Then select from three possible methods to compare set contents
// o If caller gave us a tag to use, compare entities by tag value
// o If set is ordered, compare entity types in order
// o Otherwise compare counts of entity types in sets
std::vector< EntityHandle > list1, list2;
rval = iface->get_entities_by_handle( set1, list1 );
if( MB_SUCCESS != rval ) moab_error( "get_entities_by_handle for set1 failed" );
rval = iface->get_entities_by_handle( set2, list2 );
if( MB_SUCCESS != rval ) moab_error( "get_entities_by_handle for set2 failed" );
if( list1.size() != list2.size() )
{
fprintf( stderr,
"Sets with id %d do not have the same number of entities.\n"
"Set 1 : %u Set 2 : %u\n",
id, (unsigned)list1.size(), (unsigned)list2.size() );
return false;
}
if( tag_name ) // compare contents using tag value
{
Tag tag;
if( MB_SUCCESS != iface->tag_get_handle( tag_name, 0, MB_TYPE_OPAQUE, tag, MB_TAG_ANY ) )
{
fprintf( stderr, "Could not find tag \"%s\" in file.\n", tag_name );
return false;
}
// Make sure tag is integer type
DataType type;
if( MB_SUCCESS != iface->tag_get_data_type( tag, type ) ) moab_error( "tag_get_data_type" );
if( MB_TYPE_INTEGER != type && MB_TYPE_OPAQUE != type ) moab_error( "compare_sets" );
std::vector< int > data1( list1.size() ), data2( list2.size() );
if( MB_SUCCESS != iface->tag_get_data( tag, &list1[0], list1.size(), &data1[0] ) ||
MB_SUCCESS != iface->tag_get_data( tag, &list2[0], list2.size(), &data2[0] ) )
moab_error( "tag_get_data" );
if( !( opt1 & MESHSET_ORDERED ) )
{
std::sort( data1.begin(), data1.end() );
std::sort( data2.begin(), data2.end() );
}
for( unsigned i = 0; i < data1.size(); ++i )
if( data1[i] != data2[i] )
{
fprintf( stderr, "Entities in sets with id %d do not have matching tag values.\n", id );
return false;
}
}
else if( opt1 & MESHSET_ORDERED ) // compare types of each entity in set
{
ok = true;
for( unsigned i = 0; i < list1.size(); ++i )
if( iface->type_from_handle( list1[i] ) != iface->type_from_handle( list2[i] ) )
{
fprintf( stderr,
"Entities at position %u in ordered sets with id %d\n"
"have different types.\n",
i, id );
ok = false;
}
if( !ok ) return false;
}
else // compare count of entity types in each set
{
unsigned counts1[MBMAXTYPE], counts2[MBMAXTYPE];
memset( counts1, 0, MBMAXTYPE * sizeof( unsigned ) );
memset( counts2, 0, MBMAXTYPE * sizeof( unsigned ) );
for( unsigned i = 0; i < list1.size(); ++i )
{
counts1[iface->type_from_handle( list1[i] )]++;
counts2[iface->type_from_handle( list2[i] )]++;
}
ok = true;
for( int j = 0; j < MBMAXTYPE; ++j )
if( counts1[j] != counts2[j] )
{
fprintf( stderr, "Sets with id %d have differing numbers of %s: %u and %u\n", id,
CN::EntityTypeName( (EntityType)j ), counts1[j], counts2[j] );
ok = false;
}
if( !ok ) return false;
}
// Compare set parent/child links using global id
ok = true;
const char* words[] = { "children", "parents" };
std::vector< EntityHandle > adj1, adj2;
for( unsigned two = 0; two < 2; ++two )
{
adj1.clear();
adj2.clear();
if( two )
{
if( MB_SUCCESS != iface->get_parent_meshsets( set1, adj1 ) ||
MB_SUCCESS != iface->get_parent_meshsets( set2, adj2 ) )
moab_error( "get_parent_meshsets" );
}
else
{
if( MB_SUCCESS != iface->get_child_meshsets( set1, adj1 ) ||
MB_SUCCESS != iface->get_child_meshsets( set2, adj2 ) )
moab_error( "get_child_meshsets" );
}
if( adj1.size() != adj2.size() )
{
fprintf( stderr, "Sets with id %d have different number of %s: %u and %u\n", id, words[two],
(unsigned)adj1.size(), (unsigned)adj2.size() );
ok = false;
continue;
}
std::vector< int > ids1( adj1.size() ), ids2( adj2.size() );
unsigned i;
for( i = 0; i < adj1.size(); ++i )
{
rval = iface->tag_get_data( id_tag, &adj1[i], 1, &ids1[i] );
if( MB_TAG_NOT_FOUND == rval )
ids1[i] = 0;
else if( MB_SUCCESS != rval )
moab_error( "tag_get_data" );
rval = iface->tag_get_data( id_tag, &adj2[i], 1, &ids2[i] );
if( MB_TAG_NOT_FOUND == rval )
ids2[i] = 0;
else if( MB_SUCCESS != rval )
moab_error( "tag_get_data" );
}
std::sort( ids1.begin(), ids1.end() );
std::sort( ids2.begin(), ids2.end() );
for( i = 0; i < ids1.size(); ++i )
{
if( ids1[i] != ids2[i] )
{
fprintf( stderr, "Sets with id %d have non-matching %s\n", id, words[two] );
ok = false;
break;
}
}
}
return ok;
}
| bool compare_tags | ( | EntityHandle | dod[] | ) |
Definition at line 620 of file h5file_test.cpp.
References dblname, moab::Interface::get_connectivity(), handlename, iface, intname, MB_ENTITY_NOT_FOUND, MB_SUCCESS, MB_TAG_DENSE, MB_TAG_SPARSE, MB_TAG_STORE, MB_TYPE_DOUBLE, MB_TYPE_HANDLE, MB_TYPE_INTEGER, MBHEX, MBPOLYHEDRON, moab_error(), moab::Interface::tag_get_data(), moab::Interface::tag_get_default_value(), moab::Interface::tag_get_handle(), and moab::Interface::type_from_handle().
Referenced by compare(), test_global_opaque_common(), test_global_value_common(), test_huge_var_length(), test_var_length_big_data(), test_var_length_data_common(), test_var_length_handle_tag(), and test_var_length_no_data().
{
Tag tag;
// Get integer tag handle and characterstics
if( MB_SUCCESS != iface->tag_get_handle( intname, 2, MB_TYPE_INTEGER, tag, MB_TAG_SPARSE | MB_TAG_STORE ) )
moab_error( "tag_get_handle(intname)" );
// integer tag should not have a default value
int idata[4];
if( MB_ENTITY_NOT_FOUND != iface->tag_get_default_value( tag, idata ) )
{
fprintf( stderr, "tag_get_default_value() for nonexistant default integer tag value did "
"not fail correctly.\n" );
return false;
}
// check data for integer tag on both dodecahedrons
if( MB_SUCCESS != iface->tag_get_data( tag, dod, 2, idata ) ) moab_error( "tag_get_data(itag)" );
if( idata[0] != (int)0xDEADBEEF || idata[1] != (int)0xDEFACED || idata[2] != idata[0] || idata[3] != idata[1] )
{
fprintf( stderr, "Incorrect values for integer tag data.\n" );
return false;
}
// Get double tag handle and characterstics
if( MB_SUCCESS != iface->tag_get_handle( dblname, 1, MB_TYPE_DOUBLE, tag, MB_TAG_DENSE | MB_TAG_STORE ) )
moab_error( "tag_get_handle(dblname)" );
// check default value of double tag
double ddata[2];
if( MB_SUCCESS != iface->tag_get_default_value( tag, ddata ) ) moab_error( "tag_get_default_value" );
if( ddata[0] != 3.14159 )
{
fprintf( stderr, "incorrect default value for double tag.\n" );
return false;
}
// check data for double tag on both dodecahedrons
if( MB_SUCCESS != iface->tag_get_data( tag, dod, 2, ddata ) ) moab_error( "tag_get_data()" );
if( ddata[0] != .333 || ddata[0] != ddata[1] )
{
fprintf( stderr, "Incorrect values for double tag data.\n" );
return false;
}
// Get handle tag handle and characterstics
if( MB_SUCCESS != iface->tag_get_handle( handlename, 3, MB_TYPE_HANDLE, tag, MB_TAG_SPARSE | MB_TAG_STORE ) )
moab_error( "tag_get_handle(handlename)" );
// check default value of handle tag
// default value will not change after tag is created. As we
// do multiple iterations of save/restore, we shouldn't expect
// the dodecahedron handles to point to the most recent two.
EntityHandle hdata[6];
if( MB_SUCCESS != iface->tag_get_default_value( tag, hdata ) ) moab_error( "tag_get_default_value" );
if( iface->type_from_handle( hdata[0] ) != MBHEX || hdata[2] != 0 ||
iface->type_from_handle( hdata[1] ) != MBPOLYHEDRON )
{
fprintf( stderr, "incorrect default value for handle tag '%s'\n", handlename );
return false;
}
// check global value for handle tag
// global value should be changed each time a new global is read,
// so expect one of the two dodecahedrons in the last slot.
const EntityHandle root = 0;
if( MB_SUCCESS != iface->tag_get_data( tag, &root, 1, hdata ) ) moab_error( "tag_get_data(mesh)" );
if( iface->type_from_handle( hdata[1] ) != MBHEX || hdata[0] != 0 || ( hdata[2] != dod[0] && hdata[2] != dod[1] ) )
{
fprintf( stderr, "incorrect global/mesh value for handle tag.\n" );
return false;
}
// check data on each dodecahedron
const EntityHandle* conn;
int len;
if( MB_SUCCESS != iface->tag_get_data( tag, dod, 2, hdata ) ) moab_error( "tag_get_data()" );
if( MB_SUCCESS != iface->get_connectivity( dod[0], conn, len ) ) moab_error( "get_connectivity" );
if( memcmp( conn, hdata, 3 * sizeof( EntityHandle ) ) )
{
fprintf( stderr, "Incorrect values for handle tag data.\n" );
return false;
}
if( MB_SUCCESS != iface->get_connectivity( dod[1], conn, len ) ) moab_error( "get_connectivity" );
if( memcmp( conn, hdata + 3, 3 * sizeof( EntityHandle ) ) )
{
fprintf( stderr, "Incorrect values for handle tag data.\n" );
return false;
}
return true;
}
| void create | ( | ) |
Definition at line 217 of file h5file_test.cpp.
References moab::Interface::add_child_meshset(), moab::Interface::add_parent_child(), moab::Interface::add_parent_meshset(), bitname, moab::Interface::create_element(), dblname, EMPTY_SET_ID, FACE_SET_ID, handlename, iface, intname, make_set(), MB_SUCCESS, MB_TAG_DENSE, MB_TAG_EXCL, MB_TAG_SPARSE, MB_TYPE_BIT, MB_TYPE_DOUBLE, MB_TYPE_HANDLE, MB_TYPE_INTEGER, MBHEX, MBPOLYHEDRON, MESHSET_SET, MESHSET_TRACK_OWNER, moab_error(), pent(), REGION_SET_ID, SET_SET_ID, moab::Interface::tag_get_handle(), moab::Interface::tag_set_data(), tagname, VERTEX_SET_ID, and vtx().
Referenced by main(), and moab::WriteNC::write_file().
{
// Create dodecahedron
// radius
const double r = 50.;
// center
const double x = 0., y = 0., z = 0.;
// length of edge of inscribed cube
const double cube = r * ( 2.0 / sqrt( 3.0 ) );
// length of dodecahedron edge
const double edge = cube * ( 2.0 / ( 1.0 + sqrt( 5.0 ) ) );
// distance of projection of a dodecahedron vertex to
// closest edge of inscribed cube
const double p = ( cube - edge ) / 2.0;
// distance of projection of a dodecahedron vertex to
// closest face of inscribed cube
const double d = sqrt( edge * edge - cube * cube / 4 - p );
// coordinate values
const double c = cube / 2.0;
const double a = c + d;
const double b = edge / 2.0;
// list of vertex handles
EntityHandle vertices[20];
// list of pentagon handles
EntityHandle faces[12];
// Dodecahedron handle
EntityHandle dodec;
// Inscribed Hex handle
EntityHandle hex;
// Create vertices if inscribed cube
vertices[0] = vtx( x - c, y + c, z + c );
vertices[1] = vtx( x + c, y + c, z + c );
vertices[2] = vtx( x + c, y - c, z + c );
vertices[3] = vtx( x - c, y - c, z + c );
vertices[4] = vtx( x - c, y + c, z - c );
vertices[5] = vtx( x + c, y + c, z - c );
vertices[6] = vtx( x + c, y - c, z - c );
vertices[7] = vtx( x - c, y - c, z - c );
// Create inscribed hex
if( MB_SUCCESS != iface->create_element( MBHEX, vertices, 8, hex ) ) moab_error( "create_element" );
// Create vertices, 2 "above" each face of inscribed cube
// +z face
vertices[8] = vtx( x - b, y, z + a );
vertices[9] = vtx( x + b, y, z + a );
// +x face
vertices[10] = vtx( x + a, y + b, z );
vertices[11] = vtx( x + a, y - b, z );
// -z face
vertices[12] = vtx( x - b, y, z - a );
vertices[13] = vtx( x + b, y, z - a );
// -x face
vertices[14] = vtx( x - a, y + b, z );
vertices[15] = vtx( x - a, y - b, z );
// +y face
vertices[16] = vtx( x, y + a, z + b );
vertices[17] = vtx( x, y + a, z - b );
// -y face
vertices[18] = vtx( x, y - a, z + b );
vertices[19] = vtx( x, y - a, z - b );
// Create petagons
faces[0] = pent( vertices, 0, 8, 9, 1, 16 );
faces[1] = pent( vertices, 3, 18, 2, 9, 8 );
faces[2] = pent( vertices, 2, 11, 10, 1, 9 );
faces[3] = pent( vertices, 2, 18, 19, 6, 11 );
faces[4] = pent( vertices, 1, 10, 5, 17, 16 );
faces[5] = pent( vertices, 5, 10, 11, 6, 13 );
faces[6] = pent( vertices, 4, 17, 5, 13, 12 );
faces[7] = pent( vertices, 7, 12, 13, 6, 19 );
faces[8] = pent( vertices, 4, 12, 7, 15, 14 );
faces[9] = pent( vertices, 0, 16, 17, 4, 14 );
faces[10] = pent( vertices, 3, 15, 7, 19, 18 );
faces[11] = pent( vertices, 0, 14, 15, 3, 8 );
// Create dodecahedron
if( MB_SUCCESS != iface->create_element( MBPOLYHEDRON, faces, 12, dodec ) ) moab_error( "create_element" );
// Create a dense tag
int zero = 0;
Tag tag;
if( MB_SUCCESS != iface->tag_get_handle( tagname, 1, MB_TYPE_INTEGER, tag, MB_TAG_DENSE | MB_TAG_EXCL, &zero ) )
moab_error( "tag_get_handle" );
// Put dense tag on all vertices.
for( int i = 0; i < 20; ++i )
if( MB_SUCCESS != iface->tag_set_data( tag, vertices + i, 1, &i ) ) moab_error( "tag_set_data" );
// Create bit tag
Tag tag2;
if( MB_SUCCESS != iface->tag_get_handle( bitname, 2, MB_TYPE_BIT, tag2, MB_TAG_EXCL ) )
moab_error( "tag_get_handle" );
// Set tag to 0 on Hex
char two = '\002';
if( MB_SUCCESS != iface->tag_set_data( tag2, &hex, 1, &two ) ) moab_error( "tag_set_data" );
// set tag to 1 on dodecahedron
char one = '\001';
if( MB_SUCCESS != iface->tag_set_data( tag2, &dodec, 1, &one ) ) moab_error( "tag_set_data" );
// Create an integer array tag and set some values on the dodecahedron
Tag itag;
if( MB_SUCCESS != iface->tag_get_handle( intname, 2, MB_TYPE_INTEGER, itag, MB_TAG_SPARSE | MB_TAG_EXCL ) )
moab_error( "tag_get_handle(MB_TYPE_INT)" );
int idata[] = { static_cast< int >( 0xDEADBEEF ), static_cast< int >( 0xDEFACED ) };
if( MB_SUCCESS != iface->tag_set_data( itag, &dodec, 1, idata ) ) moab_error( "tag_set_data(itag)" );
// Create a double tag with a non-zero default value, and set on dodecahedron
Tag dtag;
double ddef = 3.14159;
if( MB_SUCCESS != iface->tag_get_handle( dblname, 1, MB_TYPE_DOUBLE, dtag, MB_TAG_DENSE | MB_TAG_EXCL, &ddef ) )
moab_error( "tag_get_handle(dtag)" );
double dval = 0.333;
if( MB_SUCCESS != iface->tag_set_data( dtag, &dodec, 1, &dval ) ) moab_error( "tag_set_data(dtag)" );
// Create a tag containing entity handles, with default values
Tag htag;
EntityHandle hdef[] = { hex, dodec, 0 };
if( MB_SUCCESS != iface->tag_get_handle( handlename, 3, MB_TYPE_HANDLE, htag, MB_TAG_SPARSE | MB_TAG_EXCL, hdef ) )
moab_error( "tag_get_handle(htag)" );
// Set global (mesh) value for tag
EntityHandle hgbl[] = { 0, hex, dodec };
const EntityHandle root = 0;
if( MB_SUCCESS != iface->tag_set_data( htag, &root, 1, hgbl ) ) moab_error( "tag_set_data(hgbl)" );
// Store first three entiries of dodec connectivity on dodec
EntityHandle hval[] = { faces[0], faces[1], faces[2] };
if( MB_SUCCESS != iface->tag_set_data( htag, &dodec, 1, hval ) ) moab_error( "tag_set_data(hgbl)" );
// create some sets
EntityHandle face_set, vertex_set, region_set, empty_set;
EntityHandle regions[] = { dodec, hex };
const unsigned empty_flags = MESHSET_ORDERED | MESHSET_TRACK_OWNER;
face_set = make_set( MESHSET_SET, faces, 12, false, FACE_SET_ID );
vertex_set = make_set( MESHSET_ORDERED, vertices, 20, true, VERTEX_SET_ID );
region_set = make_set( MESHSET_SET, regions, 2, false, REGION_SET_ID );
empty_set = make_set( empty_flags, 0, 0, true, EMPTY_SET_ID );
EntityHandle sets[] = { face_set, vertex_set, region_set, empty_set };
make_set( MESHSET_ORDERED, sets, 4, false, SET_SET_ID );
// create some set parent-child links
if( MB_SUCCESS != iface->add_parent_child( face_set, vertex_set ) ) moab_error( "add_parent_child" );
if( MB_SUCCESS != iface->add_child_meshset( region_set, face_set ) ) moab_error( "add_child_meshset" );
if( MB_SUCCESS != iface->add_parent_meshset( vertex_set, region_set ) ) moab_error( "add_parent_meshet" );
}
| int main | ( | int | argc, |
| char * | argv[] | ||
| ) |
Definition at line 53 of file h5file_test.cpp.
References compare(), create(), moab::Interface::delete_mesh(), ErrorCode, filename, moab::Interface::get_last_error(), iface, moab::Interface::load_file(), MB_SUCCESS, moab_error(), and moab::Interface::write_file().
{
ErrorCode rval;
std::string msg;
std::string read_opt;
std::string write_opt( "DEBUG_BINIO" );
for( int i = 1; i < argc; ++i )
{
long val;
char* endptr;
if( argv[i][0] == '-' && ( argv[i][1] == 'r' || argv[i][1] == 'w' ) && i + 1 < argc &&
( val = strtol( argv[i + 1], &endptr, 0 ) ) >= 0 && !*endptr )
{
std::string& s = argv[i][1] == 'r' ? read_opt : write_opt;
std::ostringstream str;
str << "DEBUG_IO=" << val;
if( s.empty() )
s = str.str();
else
{
s += ';';
s += str.str();
}
++i;
}
else
{
std::cerr << "Usage: " << argv[0] << " [-r <n>] [-w <n>]" << std::endl;
return 1;
}
}
iface = new Core();
// create a dodecahedron and inscribed hex
fprintf( stderr, "creating... " );
create();
// write out the dodecahedron
fprintf( stderr, "writing... " );
rval = iface->write_file( filename, 0, write_opt.c_str() );
if( MB_SUCCESS != rval )
{
fprintf( stderr, "Failed to write \"%s\"\n", filename );
if( MB_SUCCESS == iface->get_last_error( msg ) ) fprintf( stderr, "%s\n", msg.c_str() );
delete iface;
return 1;
}
// Read back in as a copy of the original
fprintf( stderr, "reading... " );
rval = iface->load_file( filename, 0, read_opt.c_str() );
if( MB_SUCCESS != rval )
{
fprintf( stderr, "Failed to read \"%s\"\n", filename );
if( MB_SUCCESS == iface->get_last_error( msg ) ) fprintf( stderr, "%s\n", msg.c_str() );
delete iface;
return 1;
}
// Compare the two.
fprintf( stderr, "comparing... " );
if( !compare() )
{
fprintf( stderr, "Comparison failed.\n" );
delete iface;
return 1;
}
fprintf( stderr, "success!\n" );
// Write both the original and copy to a file
fprintf( stderr, "writing... " );
rval = iface->write_file( filename, 0, write_opt.c_str() );
if( MB_SUCCESS != rval )
{
fprintf( stderr, "Failed to write \"%s\"\n", filename );
if( MB_SUCCESS == iface->get_last_error( msg ) ) fprintf( stderr, "%s\n", msg.c_str() );
delete iface;
return 1;
}
// Delete the mesh
fprintf( stderr, "clearing db... " );
rval = iface->delete_mesh();
if( MB_SUCCESS != rval ) moab_error( "delete_mesh" );
// Read the two dodecahedrons from the file
fprintf( stderr, "reading... " );
rval = iface->load_file( filename, 0, read_opt.c_str() );
if( MB_SUCCESS != rval )
{
fprintf( stderr, "Failed to read \"%s\"\n", filename );
if( MB_SUCCESS == iface->get_last_error( msg ) ) fprintf( stderr, "%s\n", msg.c_str() );
delete iface;
return 1;
}
// Compare them
fprintf( stderr, "comparing... " );
if( !compare() )
{
fprintf( stderr, "Comparison failed.\n" );
delete iface;
return 1;
}
fprintf( stderr, "success!\n" );
// Delete the mesh
fprintf( stderr, "cleaning up... " );
rval = iface->delete_mesh();
if( MB_SUCCESS != rval ) moab_error( "delete_mesh" );
// Clean up the file.
remove( filename );
fprintf( stderr, "done.\n" );
delete iface;
return 0;
}
| EntityHandle make_set | ( | unsigned int | options, |
| EntityHandle * | entities, | ||
| size_t | num_entities, | ||
| bool | reverse, | ||
| int | id | ||
| ) |
Definition at line 196 of file h5file_test.cpp.
References moab::Interface::add_entities(), moab::Interface::create_meshset(), moab::Interface::globalId_tag(), id_tag, iface, MB_SUCCESS, moab_error(), and moab::Interface::tag_set_data().
Referenced by create().
{
EntityHandle handle;
if( MB_SUCCESS != iface->create_meshset( options, handle ) ) moab_error( "create_meshset" );
if( reverse )
{
for( int i = (int)num_entities - 1; i >= 0; --i )
if( MB_SUCCESS != iface->add_entities( handle, entities + i, 1 ) ) moab_error( "add_entities" );
}
else
{
if( MB_SUCCESS != iface->add_entities( handle, entities, num_entities ) ) moab_error( "add_entities" );
}
Tag id_tag = iface->globalId_tag();
if( MB_SUCCESS != iface->tag_set_data( id_tag, &handle, 1, &id ) ) moab_error( "tag_set_data" );
return handle;
}
| void moab_error | ( | const char * | function | ) |
Definition at line 800 of file h5file_test.cpp.
References moab::Interface::get_last_error(), iface, and MB_SUCCESS.
Referenced by compare(), compare_conn(), compare_sets(), compare_tags(), create(), main(), make_set(), pent(), and vtx().
{
std::string msg;
fprintf( stderr, "%s() failed.\n", str );
if( MB_SUCCESS == iface->get_last_error( msg ) ) fprintf( stderr, "%s\n", msg.c_str() );
delete iface;
exit( 1 );
}
| EntityHandle pent | ( | EntityHandle * | vtx_list, |
| int | i1, | ||
| int | i2, | ||
| int | i3, | ||
| int | i4, | ||
| int | i5 | ||
| ) |
Definition at line 182 of file h5file_test.cpp.
References moab::Interface::create_element(), iface, MB_SUCCESS, MBPOLYGON, and moab_error().
Referenced by create().
{
const EntityHandle conn[5] = { vtx_list[i1], vtx_list[i2], vtx_list[i3], vtx_list[i4], vtx_list[i5] };
EntityHandle result;
if( MB_SUCCESS != iface->create_element( MBPOLYGON, conn, 5, result ) ) moab_error( "create_element" );
return result;
}
| EntityHandle vtx | ( | double | x, |
| double | y, | ||
| double | z | ||
| ) |
Definition at line 174 of file h5file_test.cpp.
References moab::Interface::create_vertex(), iface, MB_SUCCESS, and moab_error().
Referenced by check_one_layer(), compute_dual_mesh(), create(), moab::ReadCGM::create_surface_facets(), moab::BSPTreePoly::is_valid(), mb_skin_adjacent_surf_patches(), PartMap::part_from_coords(), moab::ReadVtk::read_vertices(), moab::BSPTreePoly::VertexUse::set_vertex(), test_entity_copy_parts(), and test_ghost_elements().
{
const double p[3] = { x, y, z };
EntityHandle result;
if( MB_SUCCESS != iface->create_vertex( p, result ) ) moab_error( "create_vertex" );
return result;
}
| const int EMPTY_SET_ID = 1100 |
Definition at line 42 of file h5file_test.cpp.
| const int FACE_SET_ID = 1101 |
Definition at line 39 of file h5file_test.cpp.
Definition at line 45 of file h5file_test.cpp.
| const int REGION_SET_ID = 1103 |
Definition at line 41 of file h5file_test.cpp.
| const int SET_SET_ID = 1105 |
Definition at line 43 of file h5file_test.cpp.
| const int VERTEX_SET_ID = 1102 |
Definition at line 40 of file h5file_test.cpp.
1.7.6.1