gqt_simple_test.cpp File Reference

#include <iostream>
#include "moab/Interface.hpp"
#include "TestUtil.hpp"
#include "Internals.hpp"
#include "moab/Core.hpp"
#include "moab/GeomTopoTool.hpp"
#include "moab/GeomQueryTool.hpp"

Include dependency graph for gqt_simple_test.cpp:

Go to the source code of this file.

Macros
#define	IS_BUILDING_MB
#define	CHKERR(A)

Functions
void	gqt_load_file ()
void	gqt_load_file_dagmc_build_obb ()
void	gqt_load_file_dagmc_via_moab_build_obb ()
void	gqt_load_file_dagmc_internal_build_obb ()
void	gqt_test_obb_retreval ()
void	gqt_create_impl_compl ()
void	gqt_build_obb ()
void	gqt_num_vols ()
void	gqt_point_in ()
void	gqt_test_obb_retreval_rayfire ()
void	gqt_rayfire ()
void	gqt_closest_to ()
void	gqt_test_boundary ()
void	gqt_point_in_box_1 ()
void	gqt_point_in_box_2 ()
void	gqt_point_in_box_3 ()
void	gqt_point_in_box_4 ()
void	gqt_point_in_box_5 ()
void	cleanup ()
int	main (int, char **)

Variables
GeomTopoTool *	GTT
GeomQueryTool *	GQT
Interface *	MBI
const std::string	input_file = TestDir + "unittest/test_geom.h5m"

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 19 of file gqt_simple_test.cpp.

IS_BUILDING_MB

#define IS_BUILDING_MB

Definition at line 4 of file gqt_simple_test.cpp.

Function Documentation

cleanup()

void cleanup

(

)

Definition at line 338 of file gqt_simple_test.cpp.

{
 // cleanup instances
 delete GQT;
 delete GTT;
 delete MBI;
}

References GQT, GTT, and MBI.

Referenced by main().

gqt_build_obb()

void gqt_build_obb

(

)

Definition at line 126 of file gqt_simple_test.cpp.

{
 ErrorCode rval = GQT->gttool()->construct_obb_trees();CHECK_ERR( rval );
}

References CHECK_ERR, moab::GeomTopoTool::construct_obb_trees(), ErrorCode, GQT, and moab::GeomQueryTool::gttool().

Referenced by main().

gqt_closest_to()

void gqt_closest_to

(

)

Definition at line 242 of file gqt_simple_test.cpp.

{
 const double eps = 1e-6; // epsilon for test, faceting tol?
 
 int vol_idx = 1;
 // note model is cube of side 10, centred at 0,0,0, so ray fire along
 // any unit direction should be exactly 5.0
 double xyz[3] = { -6.0, 0.0, 0.0 };
 double distance; // distance from point to nearest surface
 double expect_distance = 1.0;
 EntityHandle vol_h;
 ErrorCode rval;
 
 // get the volume handle
 vol_h = GQT->gttool()->entity_by_id( 3, vol_idx );
 
 // test closest to location
 rval = GQT->closest_to_location( vol_h, xyz, distance );CHECK_ERR( rval );
 
 // distance should be 1.0 cm
 CHECK_REAL_EQUAL( expect_distance, distance, eps );
}

References CHECK_ERR, CHECK_REAL_EQUAL, moab::GeomQueryTool::closest_to_location(), moab::GeomTopoTool::entity_by_id(), eps, ErrorCode, GQT, and moab::GeomQueryTool::gttool().

Referenced by main().

gqt_create_impl_compl()

void gqt_create_impl_compl

(

)

Definition at line 121 of file gqt_simple_test.cpp.

{
 ErrorCode rval = GQT->gttool()->setup_implicit_complement();CHECK_ERR( rval );
}

References CHECK_ERR, ErrorCode, GQT, moab::GeomQueryTool::gttool(), and moab::GeomTopoTool::setup_implicit_complement().

Referenced by main().

gqt_load_file()

void gqt_load_file

(

)

Definition at line 31 of file gqt_simple_test.cpp.

{
 MBI = new moab::Core();
 ErrorCode rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true );
 GQT = new GeomQueryTool( GTT );
}

References CHECK_ERR, ErrorCode, GQT, GTT, input_file, moab::Interface::load_file(), and MBI.

Referenced by main().

gqt_load_file_dagmc_build_obb()

void gqt_load_file_dagmc_build_obb

(

)

Definition at line 40 of file gqt_simple_test.cpp.

{
 /* 1 - Test with external moab, load file in GTT/GQT*/
 // make new moab core
 ErrorCode rval;
 
 // load a file
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI );
 GQT = new GeomQueryTool( GTT );
 rval = GQT->initialize();CHECK_ERR( rval );
}

References CHECK_ERR, ErrorCode, GQT, GTT, moab::GeomQueryTool::initialize(), input_file, moab::Interface::load_file(), and MBI.

Referenced by main().

gqt_load_file_dagmc_internal_build_obb()

void gqt_load_file_dagmc_internal_build_obb

(

)

Definition at line 72 of file gqt_simple_test.cpp.

{
 /* 3 - Test with internal moab, load file in MOAB*/
 // find geomsets, build implicit complement, and construct obb trees
 ErrorCode rval;
 
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 rval = GQT->gttool()->setup_implicit_complement();CHECK_ERR( rval );
 rval = GQT->gttool()->construct_obb_trees();CHECK_ERR( rval );
}

References CHECK_ERR, moab::GeomTopoTool::construct_obb_trees(), ErrorCode, GQT, GTT, moab::GeomQueryTool::gttool(), input_file, moab::Interface::load_file(), MBI, and moab::GeomTopoTool::setup_implicit_complement().

Referenced by main().

gqt_load_file_dagmc_via_moab_build_obb()

void gqt_load_file_dagmc_via_moab_build_obb

(

)

Definition at line 55 of file gqt_simple_test.cpp.

{
 /* 2 - Test with external moab, load file in MOAB*/
 // build implicit complement and obb trees
 ErrorCode rval;
 
 // load file
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the file
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI );
 GQT = new GeomQueryTool( GTT );
 // create obb trees
 rval = GQT->gttool()->setup_implicit_complement();CHECK_ERR( rval );
 rval = GQT->gttool()->construct_obb_trees();CHECK_ERR( rval );
}

References CHECK_ERR, moab::GeomTopoTool::construct_obb_trees(), ErrorCode, GQT, GTT, moab::GeomQueryTool::gttool(), input_file, moab::Interface::load_file(), MBI, and moab::GeomTopoTool::setup_implicit_complement().

Referenced by main().

gqt_num_vols()

void gqt_num_vols

(

)

Definition at line 131 of file gqt_simple_test.cpp.

{
 int expect_num_vols = 2;
 int num_vols = GQT->gttool()->num_ents_of_dim( 3 );
 CHECK_EQUAL( expect_num_vols, num_vols );
}

References CHECK_EQUAL, GQT, moab::GeomQueryTool::gttool(), and moab::GeomTopoTool::num_ents_of_dim().

Referenced by main().

gqt_point_in()

void gqt_point_in

(

)

Definition at line 138 of file gqt_simple_test.cpp.

{
 int result = 0;
 int expect_result = 1;
 int vol_idx = 1;
 double xyz[3] = { 0.0, 0.0, 0.0 };
 EntityHandle vol_h;
 
 // load file and create instances
 MBI = new moab::Core();
 ErrorCode rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 // initialize for geom queries
 GQT->initialize();
 
 // get volume
 vol_h = GQT->gttool()->entity_by_id( 3, vol_idx );
 // test point in volume
 GQT->point_in_volume( vol_h, xyz, result );CHECK_ERR( rval );
 CHECK_EQUAL( expect_result, result );
}

References CHECK_EQUAL, CHECK_ERR, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, GTT, moab::GeomQueryTool::gttool(), moab::GeomQueryTool::initialize(), input_file, moab::Interface::load_file(), MBI, and moab::GeomQueryTool::point_in_volume().

Referenced by main().

gqt_point_in_box_1()

void gqt_point_in_box_1

(

)

Definition at line 283 of file gqt_simple_test.cpp.

{
 int vol_id = 1;
 double origin[3] = { 0.0, 0.0, 0.0 };
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_id );
 int result;
 ErrorCode rval = GQT->point_in_box( vol_h, origin, result );
 CHECK_EQUAL( rval, moab::MB_SUCCESS );
 CHECK_EQUAL( result, 1 ); // inside
}

References CHECK_EQUAL, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), MB_SUCCESS, origin, and moab::GeomQueryTool::point_in_box().

gqt_point_in_box_2()

void gqt_point_in_box_2

(

)

Definition at line 294 of file gqt_simple_test.cpp.

{
 int vol_id = 1;
 double origin[3] = { 5.0, 0.0, 0.0 };
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_id );
 int result;
 ErrorCode rval = GQT->point_in_box( vol_h, origin, result );
 CHECK_EQUAL( rval, moab::MB_SUCCESS );
 CHECK_EQUAL( result, 1 ); // inside
}

References CHECK_EQUAL, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), MB_SUCCESS, origin, and moab::GeomQueryTool::point_in_box().

gqt_point_in_box_3()

void gqt_point_in_box_3

(

)

Definition at line 305 of file gqt_simple_test.cpp.

{
 int vol_id = 1;
 double origin[3] = { 5.1, 0.0, 0.0 };
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_id );
 int result;
 ErrorCode rval = GQT->point_in_box( vol_h, origin, result );
 CHECK_EQUAL( rval, moab::MB_SUCCESS );
 CHECK_EQUAL( result, 0 ); // outside
}

References CHECK_EQUAL, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), MB_SUCCESS, origin, and moab::GeomQueryTool::point_in_box().

gqt_point_in_box_4()

void gqt_point_in_box_4

(

)

Definition at line 316 of file gqt_simple_test.cpp.

{
 int vol_id = 1;
 double origin[3] = { 5.1, 5.1, 5.1 };
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_id );
 int result;
 ErrorCode rval = GQT->point_in_box( vol_h, origin, result );
 CHECK_EQUAL( rval, moab::MB_SUCCESS );
 CHECK_EQUAL( result, 0 ); // outside
}

References CHECK_EQUAL, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), MB_SUCCESS, origin, and moab::GeomQueryTool::point_in_box().

gqt_point_in_box_5()

void gqt_point_in_box_5

(

)

Definition at line 327 of file gqt_simple_test.cpp.

{
 int vol_id = 1;
 double origin[3] = { 5.0 - 1.e-308, 5.0 - 1.e-308, 5.0 - 1.e-308 };
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_id );
 int result;
 ErrorCode rval = GQT->point_in_box( vol_h, origin, result );
 CHECK_EQUAL( rval, moab::MB_SUCCESS );
 CHECK_EQUAL( result, 1 ); // outside
}

References CHECK_EQUAL, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), MB_SUCCESS, origin, and moab::GeomQueryTool::point_in_box().

gqt_rayfire()

void gqt_rayfire

(

)

Definition at line 215 of file gqt_simple_test.cpp.

{
 const double eps = 1e-6; // epsilon for test, faceting tol?
 
 int vol_idx = 1;
 // note model is cube of side 10, centred at 0,0,0, so ray fire along
 // any unit direction should be exactly 5.0
 double xyz[3] = { 0.0, 0.0, 0.0 };
 double dir[3] = { 0.0, 0.0, 1.0 };
 EntityHandle next_surf;
 double next_surf_dist;
 double expect_next_surf_dist = 5.0;
 EntityHandle vol_h;
 ErrorCode rval;
 
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 GQT->initialize();CHECK_ERR( rval );
 
 vol_h = GQT->gttool()->entity_by_id( 3, vol_idx );
 rval = GQT->ray_fire( vol_h, xyz, dir, next_surf, next_surf_dist );CHECK_ERR( rval );
 CHECK_REAL_EQUAL( expect_next_surf_dist, next_surf_dist, eps );
}

References CHECK_ERR, CHECK_REAL_EQUAL, moab::GeomTopoTool::entity_by_id(), eps, ErrorCode, GQT, GTT, moab::GeomQueryTool::gttool(), moab::GeomQueryTool::initialize(), input_file, moab::Interface::load_file(), MBI, and moab::GeomQueryTool::ray_fire().

Referenced by main().

gqt_test_boundary()

void gqt_test_boundary

(

)

Definition at line 265 of file gqt_simple_test.cpp.

{
 int vol_idx = 1;
 EntityHandle vol_h = GQT->gttool()->entity_by_id( 3, vol_idx );
 int surf_idx = 1;
 EntityHandle surf_h = GQT->gttool()->entity_by_id( 2, surf_idx );
 
 double xyz[3] = { 0.0, 0.0, 5.0 };
 double dir[3] = { 0.0, 0.0, 1.0 };
 int result;
 int expect_result = 0;
 
 // test vol boundary
 ErrorCode rval = GQT->test_volume_boundary( vol_h, surf_h, xyz, dir, result );CHECK_ERR( rval );
 // check ray leaving volume
 CHECK_EQUAL( expect_result, result );
}

References CHECK_EQUAL, CHECK_ERR, moab::GeomTopoTool::entity_by_id(), ErrorCode, GQT, moab::GeomQueryTool::gttool(), and moab::GeomQueryTool::test_volume_boundary().

Referenced by main().

gqt_test_obb_retreval()

void gqt_test_obb_retreval

(

)

Definition at line 87 of file gqt_simple_test.cpp.

{
 // make new dagmc
 std::cout << "test_obb_retreval" << std::endl;
 
 ErrorCode rval;
 // load a file
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 GQT->initialize();CHECK_ERR( rval );
 
 // write to file
 rval = MBI->write_file( "fcad" );CHECK_ERR( rval );
 
 // cleanup instances
 delete GQT;
 delete GTT;
 delete MBI;
 
 // re-load written file
 MBI = new moab::Core();
 rval = MBI->load_file( "fcad" ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 rval = GQT->initialize();CHECK_ERR( rval );
 
 // remove file
 remove( "fcad" );
}

References CHECK_ERR, ErrorCode, GQT, GTT, moab::GeomQueryTool::initialize(), input_file, moab::Interface::load_file(), MBI, and moab::Interface::write_file().

Referenced by main().

gqt_test_obb_retreval_rayfire()

void gqt_test_obb_retreval_rayfire

(

)

Definition at line 162 of file gqt_simple_test.cpp.

{
 // make new dagmc
 std::cout << "test_obb_retreval and ray_fire" << std::endl;
 
 ErrorCode rval;
 
 // load file and initialize
 MBI = new moab::Core();
 rval = MBI->load_file( input_file.c_str() ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 GQT->initialize();CHECK_ERR( rval );
 
 // write to file
 rval = MBI->write_file( "fcad" );CHECK_ERR( rval );
 
 // delete instances
 delete GTT;
 delete GQT;
 delete MBI;
 
 // load the file that was written and initialize
 MBI = new moab::Core();
 rval = MBI->load_file( "fcad" ); // open the
 CHECK_ERR( rval );
 GTT = new GeomTopoTool( MBI, true, 0 );
 GQT = new GeomQueryTool( GTT );
 rval = GQT->initialize();CHECK_ERR( rval );
 
 // remove generated test file
 remove( "fcad" );
 
 // now perform full ray fire
 double eps = 1.e-6;
 int vol_idx = 1;
 // note model is cube of side 10, centred at 0,0,0, so ray fire along
 // any unit direction should be exactly 5.0
 double xyz[3] = { 0.0, 0.0, 0.0 };
 double dir[3] = { 0.0, 0.0, 1.0 };
 EntityHandle next_surf;
 double next_surf_dist;
 double expect_next_surf_dist = 5.0;
 EntityHandle vol_h;
 
 // get volume handle to fire on
 vol_h = GQT->gttool()->entity_by_id( 3, vol_idx );
 // test ray fire
 rval = GQT->ray_fire( vol_h, xyz, dir, next_surf, next_surf_dist );CHECK_ERR( rval );
 CHECK_REAL_EQUAL( expect_next_surf_dist, next_surf_dist, eps );
}

References CHECK_ERR, CHECK_REAL_EQUAL, moab::GeomTopoTool::entity_by_id(), eps, ErrorCode, GQT, GTT, moab::GeomQueryTool::gttool(), moab::GeomQueryTool::initialize(), input_file, moab::Interface::load_file(), MBI, moab::GeomQueryTool::ray_fire(), and moab::Interface::write_file().

Referenced by main().

main()

int main	(	int	,
		char **
	)

Definition at line 346 of file gqt_simple_test.cpp.

{
 int result = 0;
 
 result += RUN_TEST( gqt_load_file ); // test ray fire
 result += RUN_TEST( gqt_build_obb ); // build the obb
 result += RUN_TEST( gqt_create_impl_compl ); // build the obb
 result += RUN_TEST( gqt_num_vols ); // make sure the num of vols correct
 cleanup();
 
 result += RUN_TEST( gqt_load_file_dagmc_build_obb ); //
 cleanup();
 
 result += RUN_TEST( gqt_load_file_dagmc_via_moab_build_obb ); //
 cleanup();
 
 result += RUN_TEST( gqt_load_file_dagmc_internal_build_obb ); //
 cleanup();
 
 result += RUN_TEST( gqt_test_obb_retreval ); // check that we are retreving loaded obbs
 cleanup();
 
 result += RUN_TEST( gqt_test_obb_retreval_rayfire ); // check that we can ray fire on loaded obbs
 cleanup();
 
 result += RUN_TEST( gqt_point_in ); // check entity by point
 cleanup();
 
 result += RUN_TEST( gqt_rayfire ); // ensure ray fire distance is correct
 result += RUN_TEST( gqt_closest_to ); // check the distance to surface nearest point
 result += RUN_TEST( gqt_test_boundary ); // check particle entering leaving
 cleanup();
 
 return result;
}

References cleanup(), gqt_build_obb(), gqt_closest_to(), gqt_create_impl_compl(), gqt_load_file(), gqt_load_file_dagmc_build_obb(), gqt_load_file_dagmc_internal_build_obb(), gqt_load_file_dagmc_via_moab_build_obb(), gqt_num_vols(), gqt_point_in(), gqt_rayfire(), gqt_test_boundary(), gqt_test_obb_retreval(), gqt_test_obb_retreval_rayfire(), and RUN_TEST.

Variable Documentation

GQT

GeomQueryTool* GQT

Definition at line 16 of file gqt_simple_test.cpp.

Referenced by cleanup(), gqt_build_obb(), gqt_closest_to(), gqt_create_impl_compl(), gqt_load_file(), gqt_load_file_dagmc_build_obb(), gqt_load_file_dagmc_internal_build_obb(), gqt_load_file_dagmc_via_moab_build_obb(), gqt_num_vols(), gqt_point_in(), gqt_point_in_box_1(), gqt_point_in_box_2(), gqt_point_in_box_3(), gqt_point_in_box_4(), gqt_point_in_box_5(), gqt_rayfire(), gqt_test_boundary(), gqt_test_obb_retreval(), and gqt_test_obb_retreval_rayfire().

GTT

GeomTopoTool* GTT

Definition at line 15 of file gqt_simple_test.cpp.

Referenced by cleanup(), gqt_load_file(), gqt_load_file_dagmc_build_obb(), gqt_load_file_dagmc_internal_build_obb(), gqt_load_file_dagmc_via_moab_build_obb(), gqt_point_in(), gqt_rayfire(), gqt_test_obb_retreval(), and gqt_test_obb_retreval_rayfire().

input_file

const std::string input_file = TestDir + "unittest/test_geom.h5m"

Definition at line 29 of file gqt_simple_test.cpp.

Referenced by gqt_load_file(), gqt_load_file_dagmc_build_obb(), gqt_load_file_dagmc_internal_build_obb(), gqt_load_file_dagmc_via_moab_build_obb(), gqt_point_in(), gqt_rayfire(), gqt_test_obb_retreval(), and gqt_test_obb_retreval_rayfire().

MBI

Interface* MBI

Definition at line 17 of file gqt_simple_test.cpp.

Referenced by cleanup(), gqt_load_file(), gqt_load_file_dagmc_build_obb(), gqt_load_file_dagmc_internal_build_obb(), gqt_load_file_dagmc_via_moab_build_obb(), gqt_point_in(), gqt_rayfire(), gqt_test_obb_retreval(), and gqt_test_obb_retreval_rayfire().