#include <iostream>
#include "moab/Interface.hpp"
#include "TestUtil.hpp"
#include "Internals.hpp"
#include "moab/Core.hpp"
#include "MBTagConventions.hpp"
#include "InitCGMA.hpp"
#include "GeometryQueryTool.hpp"
#include "moab/MeshTopoUtil.hpp"

Include dependency graph for read_cgm_connectivity_test.cpp:

Macros
#define	IS_BUILDING_MB

#define	CHKERR(A)

Functions
void	read_file (Interface moab, const char input_file)

void	cube_verts_connectivity_test ()

void	cube_tris_connectivity_test ()

void	cube_tri_curve_coincidence_test ()

void	cube_edge_adjacencies_test ()

void	cube_tri_vertex_test ()

void	match_tri_edges_w_curve (const Range &tri_edges, const Range &curves)

int	main (int, char **)

Variables
std::string	input_cube = TestDir + "unittest/io/cube.sat"

Macro Definition Documentation

◆ CHKERR

#define CHKERR ( A )

Value:

 do \
 { \
 if( MB_SUCCESS != ( A ) ) \
 { \
 std::cerr << "Failure (error code " << ( A ) << ") at " __FILE__ ":" << __LINE__ << std::endl; \
 return A; \
 } \
 } while( false )

Definition at line 15 of file read_cgm_connectivity_test.cpp.

◆ IS_BUILDING_MB

#define IS_BUILDING_MB

Definition at line 4 of file read_cgm_connectivity_test.cpp.

Function Documentation

◆ cube_edge_adjacencies_test()

void cube_edge_adjacencies_test ( )

Definition at line 176 of file read_cgm_connectivity_test.cpp.

 {
  ErrorCode rval;
  // Open the test file
  Core moab;
  Interface* mb = &moab;
  read_file( mb, input_cube.c_str() );
  
  // Get the curves
  Range curves;
  rval = mb->get_entities_by_type( 0, MBEDGE, curves );CHECK_ERR( rval );
  
  for( Range::const_iterator i = curves.begin(); i != curves.end(); ++i )
  {
  // Get triangle adjacent to each edge
  Range adj_tris;
  rval = mb->get_adjacencies( &( *i ), 1, 2, false, adj_tris );CHECK_ERR( rval );
  
  int num_adj_tris = adj_tris.size();
  // Ensure that no edge is adjacent to more than two triangles
  CHECK( num_adj_tris <= 2 );
  }
 }

References moab::Range::begin(), CHECK, CHECK_ERR, moab::Range::end(), ErrorCode, moab::Core::get_adjacencies(), moab::Core::get_entities_by_type(), input_cube, mb, MBEDGE, read_file(), and moab::Range::size().

◆ cube_tri_curve_coincidence_test()

void cube_tri_curve_coincidence_test ( )

Definition at line 123 of file read_cgm_connectivity_test.cpp.

 {
  
  ErrorCode rval;
  // Open the test file
  Core moab;
  Interface* mb = &moab;
  read_file( mb, input_cube.c_str() );
  
  // Get curves from the mesh
  Range curves;
  rval = mb->get_entities_by_type( 0, MBEDGE, curves );CHECK_ERR( rval );
  curves.print();
  
  // Get triangles from the mesh
  Range tris;
  rval = mb->get_entities_by_type( 0, MBTRI, tris );CHECK_ERR( rval );
  
  for( Range::const_iterator i = tris.begin(); i != tris.end(); ++i )
  {
  // Get the any curve edges that are a part of the triangle
  Range tri_edges;
  rval = mb->get_adjacencies( &( *i ), 1, 1, false, tri_edges );CHECK_ERR( rval );
  // Check that we've retrieved two edges from get_adjacencies
  // For a this file (cube), each triangle should have two curve
  // edges
  int num_of_tri_edges = tri_edges.size();
  CHECK_EQUAL( 2, num_of_tri_edges );
  match_tri_edges_w_curve( tri_edges, curves );CHECK_ERR( rval );
  }
 }

References moab::Range::begin(), CHECK_EQUAL, CHECK_ERR, moab::Range::end(), ErrorCode, moab::Core::get_adjacencies(), moab::Core::get_entities_by_type(), input_cube, match_tri_edges_w_curve(), mb, MBEDGE, MBTRI, moab::Range::print(), read_file(), and moab::Range::size().

Referenced by main().

◆ cube_tri_vertex_test()

void cube_tri_vertex_test ( )

Definition at line 201 of file read_cgm_connectivity_test.cpp.

 {
  ErrorCode rval;
  // Open the test file
  Core moab;
  Interface* mb = &moab;
  read_file( mb, input_cube.c_str() );
  
  // Get all triangles
  Range tris;
  rval = mb->get_entities_by_type( 0, MBTRI, tris );CHECK_ERR( rval );
  
  for( Range::const_iterator i = tris.begin(); i != tris.end(); ++i )
  {
  // Get all triangle vertices
  Range verts;
  rval = mb->get_connectivity( &( *i ), 1, verts );CHECK_ERR( rval );
  // Make sure that each vertex making up
  // the triangle is different
  int number_of_verts = verts.size();
  CHECK( 3 == number_of_verts );
  CHECK( verts[0] != verts[1] );
  CHECK( verts[1] != verts[2] );
  CHECK( verts[2] != verts[0] );
  }
 }

References moab::Range::begin(), CHECK, CHECK_ERR, moab::Range::end(), ErrorCode, moab::Core::get_connectivity(), moab::Core::get_entities_by_type(), input_cube, mb, MBTRI, read_file(), and moab::Range::size().

Referenced by main().

◆ cube_tris_connectivity_test()

void cube_tris_connectivity_test ( )

Definition at line 92 of file read_cgm_connectivity_test.cpp.

 {
  ErrorCode rval;
  // Open the test file
  Core moab;
  Interface* mb = &moab;
  read_file( mb, input_cube.c_str() );
  
  // Get triangles from the mesh
  Range tris;
  rval = mb->get_entities_by_type( 0, MBTRI, tris );CHECK_ERR( rval );
  
  int expected_num_of_adj_tris = 3;
  
  for( Range::const_iterator i = tris.begin() + 1; i != tris.end(); ++i )
  {
  Range adj_tris;
  moab::MeshTopoUtil mu( mb );
  // Use Triangle edges to get all adjacent triangles
  rval = mu.get_bridge_adjacencies( *i, 1, 2, adj_tris );CHECK_ERR( rval );
  CHECK_EQUAL( expected_num_of_adj_tris, (int)adj_tris.size() );
  
  // Check that the entities we found from bridge_adjacencies
  // are triangles
  Range adj_tri_test = adj_tris.subset_by_type( MBTRI );
  CHECK_EQUAL( (int)adj_tris.size(), (int)adj_tri_test.size() );
  }
 }

References moab::Range::begin(), CHECK_EQUAL, CHECK_ERR, moab::Range::end(), ErrorCode, moab::MeshTopoUtil::get_bridge_adjacencies(), moab::Core::get_entities_by_type(), input_cube, mb, MBTRI, read_file(), moab::Range::size(), and moab::Range::subset_by_type().

Referenced by main().

◆ cube_verts_connectivity_test()

void cube_verts_connectivity_test ( )

Definition at line 66 of file read_cgm_connectivity_test.cpp.

 {
  
  ErrorCode rval;
  // Open the test file
  Core moab;
  Interface* mb = &moab;
  read_file( mb, input_cube.c_str() );
  
  // Get all vertex handles from the mesh
  Range verts;
  rval = mb->get_entities_by_type( 0, MBVERTEX, verts );CHECK_ERR( rval );
  
  // Check that each vertex connects to at least 4 and no more than 6 triangles
  for( Range::const_iterator i = verts.begin(); i != verts.end(); ++i )
  {
  std::vector< EntityHandle > adj_tris;
  rval = mb->get_adjacencies( &( *i ), 1, 2, false, adj_tris );CHECK_ERR( rval );
  
  int adj_size = adj_tris.size();
  CHECK( adj_size >= 4 && adj_size <= 6 );
  }
 }

References moab::Range::begin(), CHECK, CHECK_ERR, moab::Range::end(), ErrorCode, moab::Core::get_adjacencies(), moab::Core::get_entities_by_type(), input_cube, mb, MBVERTEX, and read_file().

Referenced by main().

◆ main()

int main	(	int	,
		char **
	)

Definition at line 44 of file read_cgm_connectivity_test.cpp.

 {
  int result = 0;
  
  result += RUN_TEST( cube_verts_connectivity_test );
  result += RUN_TEST( cube_tris_connectivity_test );
  result += RUN_TEST( cube_tri_curve_coincidence_test );
  result += RUN_TEST( cube_tri_vertex_test );
  
  return result;
 }

References cube_tri_curve_coincidence_test(), cube_tri_vertex_test(), cube_tris_connectivity_test(), cube_verts_connectivity_test(), and RUN_TEST.

◆ match_tri_edges_w_curve()

void match_tri_edges_w_curve	(	const Range &	tri_edges,
		const Range &	curves
	)

Definition at line 155 of file read_cgm_connectivity_test.cpp.

 {
  int match_counter = 0;
  int num_of_tri_edges = tri_edges.size();
  CHECK( num_of_tri_edges );
  for( Range::const_iterator i = tri_edges.begin(); i != tri_edges.end(); ++i )
  {
  for( Range::const_iterator j = curves.begin(); j != curves.end(); ++j )
  {
  // If the edge handle matches a curve handle, increment the number
  // matches
  if( *i == *j ) match_counter++;
  }
  }
  // Make sure that each edge returned from triangle edges
  // has been matched to a curve
  CHECK_EQUAL( num_of_tri_edges, match_counter );
 }

References moab::Range::begin(), CHECK, CHECK_EQUAL, moab::Range::end(), and moab::Range::size().

Referenced by cube_tri_curve_coincidence_test().

◆ read_file()

void read_file	(	Interface *	moab,
		const char *	input_file
	)

Definition at line 56 of file read_cgm_connectivity_test.cpp.

 {
  InitCGMA::initialize_cgma();
  GeometryQueryTool::instance()->delete_geometry();
  
  ErrorCode rval = moab->load_file( input_file );CHECK_ERR( rval );
 }