|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
Inheritance diagram for MeshInterfaceTest: Collaboration diagram for MeshInterfaceTest:Definition at line 56 of file MeshInterfaceTest.cpp.
| MeshInterfaceTest::MeshInterfaceTest | ( | ) | [inline] |
Definition at line 83 of file MeshInterfaceTest.cpp.
: mMesh( 0 ) {}
| MeshInterfaceTest::CPPUNIT_TEST | ( | test_get_geometric_dimension | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST | ( | test_vertices | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST | ( | test_vertex_byte | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST | ( | test_elements | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST | ( | test_elements_get_topology | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST_SUITE | ( | MeshInterfaceTest | ) | [private] |
| MeshInterfaceTest::CPPUNIT_TEST_SUITE_END | ( | ) | [private] |
| void MeshInterfaceTest::setUp | ( | ) | [inline] |
Definition at line 86 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, MBMesquite::MeshImpl::elements_get_attached_vertices(), MBMesquite::MeshImpl::get_all_elements(), mConnectivity, mElements, mErr, mMesh, mOffsets, mVertices, and MBMesquite::MeshImpl::read_vtk().
{
// Read a VTK file -- 1 triangle flanked by 1 quad on each side (1 tri + 3 quads)
mMesh = new MBMesquite::MeshImpl;
mMesh->read_vtk( MESH_FILES_DIR "2D/vtk/mixed/untangled/hybrid_3quad_1tri.vtk", mErr );
CPPUNIT_ASSERT( !mErr );
// Get mesh data
mMesh->get_all_elements( mElements, mErr );
CPPUNIT_ASSERT( !mErr );
mMesh->elements_get_attached_vertices( arrptr( mElements ), mElements.size(), mConnectivity, mOffsets, mErr );
CPPUNIT_ASSERT( !mErr );
// Construct list of vertices w/out duplicates from
// connectivity list.
std::vector< MBMesquite::Mesh::VertexHandle >::iterator new_end;
mVertices = mConnectivity;
std::sort( mVertices.begin(), mVertices.end() );
new_end = std::unique( mVertices.begin(), mVertices.end() );
mVertices.resize( new_end - mVertices.begin() );
}
| void MeshInterfaceTest::tearDown | ( | ) | [inline] |
Definition at line 109 of file MeshInterfaceTest.cpp.
References MBMesquite::MsqError::clear(), mErr, and mMesh.
| void MeshInterfaceTest::test_element_get_attached_vertex_indices | ( | ) | [inline] |
Definition at line 341 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, MBMesquite::MeshImpl::elements_get_topologies(), mConnectivity, mElements, mErr, MBMesquite::MIXED, mMesh, mOffsets, tri_coords(), MBMesquite::TRIANGLE, and MBMesquite::MeshImpl::vertices_get_coordinates().
{
// Find the index of the triangle
MBMesquite::EntityTopology topo = Mesquite::MIXED;
int tri_index = -1;
while( topo != MBMesquite::TRIANGLE )
{
++tri_index;
CPPUNIT_ASSERT( (unsigned)tri_index < mElements.size() );
MBMesquite::Mesh::ElementHandle handle = mElements[tri_index];
mMesh->elements_get_topologies( &handle, &topo, 1, mErr );
CPPUNIT_ASSERT( !mErr );
}
// creates list with correct vertices coordinates for the triangle
std::vector< Mesquite::Vector3D > correct_coords;
correct_coords.push_back( Mesquite::Vector3D( 1., 0., 0. ) );
correct_coords.push_back( Mesquite::Vector3D( 0., 1.732050807, 0. ) );
correct_coords.push_back( Mesquite::Vector3D( -1., 0., 0. ) );
// Creates same list from the mesh implementation
std::vector< MBMesquite::MsqVertex > tri_coords( 3 );
mMesh->vertices_get_coordinates( &mConnectivity[mOffsets[tri_index]], arrptr( tri_coords ), 3, mErr );
CPPUNIT_ASSERT( !mErr );
// Makes sure both list contain the same elements (not necessarily in the same order).
std::vector< Mesquite::Vector3D >::iterator correct_iter;
std::vector< MBMesquite::MsqVertex >::iterator tri_iter;
for( tri_iter = tri_coords.begin(); tri_iter != tri_coords.end(); ++tri_iter )
{
for( correct_iter = correct_coords.begin(); correct_iter != correct_coords.end(); ++correct_iter )
{
if( Mesquite::Vector3D::distance_between( *tri_iter, *correct_iter ) < 10e-4 ) break;
}
// check if a match was found
CPPUNIT_ASSERT( correct_iter != correct_coords.end() );
// remove match from list
correct_coords.erase( correct_iter );
}
CPPUNIT_ASSERT( correct_coords.empty() );
}
| void MeshInterfaceTest::test_elements | ( | ) | [inline] |
Definition at line 268 of file MeshInterfaceTest.cpp.
References CPPUNIT_ASSERT_EQUAL, and mElements.
{
CPPUNIT_ASSERT_EQUAL( 4, (int)mElements.size() );
}
| void MeshInterfaceTest::test_elements_get_attached_vertices | ( | ) | [inline] |
Definition at line 273 of file MeshInterfaceTest.cpp.
References CPPUNIT_ASSERT, mConnectivity, mElements, mOffsets, and mVertices.
{
const size_t nbElem = mElements.size();
// checks we have 3 elements containing 4 vertices
// and 1 element containing 3 vertices.
int n3 = 0;
int n4 = 0;
size_t i;
for( i = 0; i < nbElem; ++i )
{
size_t nve = mOffsets[i + 1] - mOffsets[i];
if( nve == 3 )
++n3;
else if( nve == 4 )
++n4;
else // failure.
CPPUNIT_ASSERT( false );
}
CPPUNIT_ASSERT( n3 == 1 ); // 1 triangle
CPPUNIT_ASSERT( n4 == 3 ); // 3 quads
// Make sure CSR data is valid
std::map< MBMesquite::Mesh::VertexHandle, int > vtx_repeated_occurence;
for( i = 0; i < mVertices.size(); ++i )
vtx_repeated_occurence[mVertices[i]] = 0;
for( i = 0; i < mConnectivity.size(); ++i )
++vtx_repeated_occurence[mConnectivity[i]];
for( i = 0; i < 9; ++i )
{
CPPUNIT_ASSERT( vtx_repeated_occurence[mVertices[i]] <= 3 );
}
// Makes sure CSR offsets are valid
CPPUNIT_ASSERT( mOffsets[0] == 0 );
CPPUNIT_ASSERT( mOffsets[1] >= 3 && mOffsets[1] <= 12 );
CPPUNIT_ASSERT( mOffsets[2] >= 3 && mOffsets[2] <= 12 );
CPPUNIT_ASSERT( mOffsets[3] >= 3 && mOffsets[3] <= 12 );
CPPUNIT_ASSERT( mOffsets[4] == 15 );
}
| void MeshInterfaceTest::test_elements_get_topology | ( | ) | [inline] |
Definition at line 314 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, CPPUNIT_ASSERT_EQUAL, MBMesquite::MeshImpl::elements_get_topologies(), mElements, mErr, mMesh, MBMesquite::QUADRILATERAL, and MBMesquite::TRIANGLE.
{
const size_t nbElem = mElements.size();
int nb_quads = 0;
int nb_tri = 0;
MBMesquite::EntityTopology* topos = new MBMesquite::EntityTopology[nbElem];
mMesh->elements_get_topologies( arrptr( mElements ), topos, nbElem, mErr );
CPPUNIT_ASSERT( !mErr );
for( size_t i = 0; i < nbElem; ++i )
{
switch( topos[i] )
{
case MBMesquite::TRIANGLE:
++nb_tri;
break;
case MBMesquite::QUADRILATERAL:
++nb_quads;
break;
default:
CPPUNIT_FAIL( "Topology should be quad or Hex only." );
}
}
CPPUNIT_ASSERT_EQUAL( 1, nb_tri );
CPPUNIT_ASSERT_EQUAL( 3, nb_quads );
delete[] topos;
}
| void MeshInterfaceTest::test_get_geometric_dimension | ( | ) | [inline] |
Definition at line 118 of file MeshInterfaceTest.cpp.
References CPPUNIT_ASSERT_EQUAL, MBMesquite::MeshImpl::get_geometric_dimension(), mErr, and mMesh.
{
int d = mMesh->get_geometric_dimension( mErr );
CPPUNIT_ASSERT_EQUAL( d, 3 );
}
| void MeshInterfaceTest::test_vertex_byte | ( | ) | [inline] |
Definition at line 183 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, mErr, mMesh, mVertices, MBMesquite::MeshImpl::vertex_get_byte(), MBMesquite::MeshImpl::vertex_set_byte(), MBMesquite::MeshImpl::vertices_get_byte(), and MBMesquite::MeshImpl::vertices_set_byte().
{
size_t nbVert = mVertices.size();
size_t i;
unsigned char* bytes = new unsigned char[nbVert];
mMesh->vertices_get_byte( arrptr( mVertices ), bytes, nbVert, mErr );
CPPUNIT_ASSERT( !mErr );
// Asserts all vertex bytes are initialised to 0.
for( i = 0; i < nbVert; ++i )
CPPUNIT_ASSERT( bytes[i] == 0 );
// Test various vertex byte read / write routines.
bytes[3] |= 4;
mMesh->vertices_set_byte( arrptr( mVertices ), bytes, nbVert, mErr );
CPPUNIT_ASSERT( !mErr );
mMesh->vertex_set_byte( mVertices[5], 8, mErr );
CPPUNIT_ASSERT( !mErr );
unsigned char byte;
mMesh->vertex_get_byte( mVertices[3], &byte, mErr );
CPPUNIT_ASSERT( !mErr );
CPPUNIT_ASSERT( byte == 4 );
mMesh->vertices_get_byte( arrptr( mVertices ), bytes, nbVert, mErr );
CPPUNIT_ASSERT( !mErr );
for( i = 0; i < nbVert; ++i )
{
if( i == 3 )
CPPUNIT_ASSERT( bytes[i] == 4 );
else if( i == 5 )
CPPUNIT_ASSERT( bytes[i] == 8 );
else
CPPUNIT_ASSERT( bytes[i] == 0 );
}
delete[] bytes;
}
| void MeshInterfaceTest::test_vertex_get_attached_elements | ( | ) | [inline] |
Definition at line 220 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, CPPUNIT_ASSERT_EQUAL, mErr, mMesh, mVertices, and MBMesquite::MeshImpl::vertices_get_attached_elements().
{
size_t i;
const size_t nbVert = mVertices.size();
std::vector< MBMesquite::Mesh::ElementHandle > elements;
std::vector< size_t > offsets;
mMesh->vertices_get_attached_elements( arrptr( mVertices ), mVertices.size(), elements, offsets, mErr );
CPPUNIT_ASSERT( !mErr );
CPPUNIT_ASSERT_EQUAL( offsets.size(), mVertices.size() + 1 );
// checks we have 6 vertices contained in 1 element only
// and 3 vertices contained in 3 elements.
int n1 = 0;
int n3 = 0;
for( i = 1; i <= mVertices.size(); ++i )
{
const size_t nev = offsets[i] - offsets[i - 1];
if( nev == 1 )
++n1;
else if( nev == 3 )
++n3;
else // failure.
CPPUNIT_ASSERT( false );
}
CPPUNIT_ASSERT( n1 == 6 );
CPPUNIT_ASSERT( n3 == 3 );
// gets the index of a vertex in a corner
int one_corner_vertex_index = 0;
for( i = 0; i < nbVert; ++i )
{
const size_t nev = offsets[i + 1] - offsets[i];
if( 1 == nev ) break;
}
CPPUNIT_ASSERT( i < nbVert );
one_corner_vertex_index = i;
// retrieve the attached element.
// This is a poor test. We allow an element handle of zero,
// and just testing that the function returned something isn't
// that useful. - J.K.
// MBMesquite::Mesh::ElementHandle elem=0;
// mMesh->vertex_get_attached_elements(mVertices[one_corner_vertex_index], &elem, 1, mErr);
// CPPUNIT_ASSERT(!mErr);
// CPPUNIT_ASSERT(elem!=0);
}
| void MeshInterfaceTest::test_vertex_is_fixed | ( | ) | [inline] |
Definition at line 169 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, CPPUNIT_ASSERT_EQUAL, fixed, mErr, mMesh, mVertices, and MBMesquite::MeshImpl::vertices_get_fixed_flag().
{
size_t nbVert = mVertices.size();
bool correct_fixed[9] = { false, false, false, true, true, true, true, true, true };
std::vector< bool > fixed;
mMesh->vertices_get_fixed_flag( arrptr( mVertices ), fixed, 9, mErr );
CPPUNIT_ASSERT( !mErr );
CPPUNIT_ASSERT_EQUAL( (size_t)8, fixed.size() );
for( size_t i = 0; i < nbVert; ++i )
{
CPPUNIT_ASSERT_EQUAL( correct_fixed[i], (bool)fixed[i] );
}
}
| void MeshInterfaceTest::test_vertices | ( | ) | [inline] |
Definition at line 124 of file MeshInterfaceTest.cpp.
References MBMesquite::arrptr(), CPPUNIT_ASSERT, CPPUNIT_ASSERT_DOUBLES_EQUAL, CPPUNIT_ASSERT_EQUAL, mErr, mMesh, mVertices, MBMesquite::Vector3D::set(), MBMesquite::MeshImpl::vertex_set_coordinates(), and MBMesquite::MeshImpl::vertices_get_coordinates().
{
size_t nbVert = mVertices.size();
CPPUNIT_ASSERT_EQUAL( 9, (int)nbVert );
MBMesquite::MsqVertex correct_coords[9], coords[9];
correct_coords[0].set( 1, 0, 0 );
correct_coords[1].set( 0, 1.732, 0 );
correct_coords[2].set( -1, 0, 0 );
correct_coords[3].set( -1, -2, 0 );
correct_coords[4].set( 1, -2, 0 );
correct_coords[5].set( 2.732, 1, 0 );
correct_coords[6].set( 1.732, 2.732, 0 );
correct_coords[7].set( -1.732, 2.732, 0 );
correct_coords[8].set( -2.732, 1, 0 );
mMesh->vertices_get_coordinates( arrptr( mVertices ), coords, nbVert, mErr );
CPPUNIT_ASSERT( !mErr );
for( size_t i = 0; i < nbVert; ++i )
{
for( int j = 0; j < 3; ++j )
CPPUNIT_ASSERT_DOUBLES_EQUAL( coords[i][j], correct_coords[i][j], .01 );
}
coords[3].set( 2., 3., 4. );
mMesh->vertex_set_coordinates( mVertices[3], coords[3], mErr );
CPPUNIT_ASSERT( !mErr );
MBMesquite::MsqVertex coords_2;
mMesh->vertices_get_coordinates( &mVertices[3], &coords_2, 1, mErr );
CPPUNIT_ASSERT( !mErr );
for( int j = 0; j < 3; ++j )
CPPUNIT_ASSERT_DOUBLES_EQUAL( coords[3][j], coords_2[j], 1e-6 );
}
std::vector< MBMesquite::Mesh::VertexHandle > MeshInterfaceTest::mConnectivity [private] |
Definition at line 76 of file MeshInterfaceTest.cpp.
Referenced by setUp(), test_element_get_attached_vertex_indices(), and test_elements_get_attached_vertices().
std::vector< MBMesquite::Mesh::ElementHandle > MeshInterfaceTest::mElements [private] |
Definition at line 78 of file MeshInterfaceTest.cpp.
Referenced by setUp(), test_element_get_attached_vertex_indices(), test_elements(), test_elements_get_attached_vertices(), and test_elements_get_topology().
MBMesquite::MsqError MeshInterfaceTest::mErr [private] |
Definition at line 80 of file MeshInterfaceTest.cpp.
Referenced by setUp(), tearDown(), test_element_get_attached_vertex_indices(), test_elements_get_topology(), test_get_geometric_dimension(), test_vertex_byte(), test_vertex_get_attached_elements(), test_vertex_is_fixed(), and test_vertices().
MBMesquite::MeshImpl* MeshInterfaceTest::mMesh [private] |
Definition at line 75 of file MeshInterfaceTest.cpp.
Referenced by setUp(), tearDown(), test_element_get_attached_vertex_indices(), test_elements_get_topology(), test_get_geometric_dimension(), test_vertex_byte(), test_vertex_get_attached_elements(), test_vertex_is_fixed(), and test_vertices().
std::vector< size_t > MeshInterfaceTest::mOffsets [private] |
Definition at line 79 of file MeshInterfaceTest.cpp.
Referenced by setUp(), test_element_get_attached_vertex_indices(), and test_elements_get_attached_vertices().
std::vector< MBMesquite::Mesh::VertexHandle > MeshInterfaceTest::mVertices [private] |
Definition at line 77 of file MeshInterfaceTest.cpp.
Referenced by setUp(), test_elements_get_attached_vertices(), test_vertex_byte(), test_vertex_get_attached_elements(), test_vertex_is_fixed(), and test_vertices().
1.7.6.1