|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
#include <iostream>#include "moab/Interface.hpp"#include "TestUtil.hpp"#include "Internals.hpp"#include "moab/Core.hpp"#include "DagMC.hpp"
Include dependency graph for dagmc_simple_test.cpp:Go to the source code of this file.
Defines | |
| #define | IS_BUILDING_MB |
| #define | CHKERR(A) |
Functions | |
| void | dagmc_load_file () |
| void | dagmc_load_file_dagmc () |
| void | dagmc_load_file_dagmc_via_moab () |
| void | dagmc_load_file_dagmc_internal () |
| void | dagmc_load_file_dagmc_build_obb () |
| void | dagmc_load_file_dagmc_via_moab_build_obb () |
| void | dagmc_load_file_dagmc_internal_build_obb () |
| void | dagmc_test_obb_retreval () |
| void | dagmc_build_obb () |
| void | dagmc_num_vols () |
| void | dagmc_entity_handle () |
| void | dagmc_point_in () |
| void | dagmc_test_obb_retreval_rayfire () |
| void | dagmc_rayfire () |
| void | dagmc_closest_to () |
| void | dagmc_test_boundary () |
| int | main (int, char **) |
Variables | |
| DagMC * | DAG |
| const std::string | input_file = "unittest/test_geom.h5m" |
| #define CHKERR | ( | A | ) |
do \ { \ if( MB_SUCCESS != ( A ) ) \ { \ std::cerr << "Failure (error code " << ( A ) << ") at " __FILE__ ":" << __LINE__ << std::endl; \ return A; \ } \ } while( false )
Definition at line 18 of file dagmc_simple_test.cpp.
| #define IS_BUILDING_MB |
Definition at line 4 of file dagmc_simple_test.cpp.
| void dagmc_build_obb | ( | ) |
| void dagmc_closest_to | ( | ) |
Definition at line 252 of file dagmc_simple_test.cpp.
References CHECK_ERR, CHECK_REAL_EQUAL, DAG, eps, and ErrorCode.
Referenced by main().
{
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 = DAG->entity_by_index( 3, vol_idx );
ErrorCode rval = DAG->closest_to_location( vol_h, xyz, distance );CHECK_ERR( rval );
// distance should be 1.0 cm
CHECK_REAL_EQUAL( expect_distance, distance, eps );
}
| void dagmc_entity_handle | ( | ) |
Definition at line 168 of file dagmc_simple_test.cpp.
{
/*
int num_vols = DAG->num_entities(3);
EntityHandle vol_h;
for (int i = 0; i < num_vols; i++)
vol_h = DAG->entity_by_index(3, i);
*/
// EntityHandle expect_vol_h = 12682136550675316765;
// CHECK_EQUAL(expect_vol_h, vol_h);
}
| void dagmc_load_file | ( | ) |
Definition at line 30 of file dagmc_simple_test.cpp.
References CHECK_ERR, DAG, ErrorCode, and input_file.
Referenced by main().
{
ErrorCode rval = DAG->load_file( input_file.c_str() ); // open the Dag file
CHECK_ERR( rval );
}
| void dagmc_load_file_dagmc | ( | ) |
Definition at line 36 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, and input_file.
Referenced by main().
{
/* 1 - Test with external moab, load file in DAGMC*/
// make new moab core
Core* mbi = new moab::Core();
// make new dagmc into that moab
DagMC* dagmc = new moab::DagMC( mbi );
ErrorCode rval;
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
// delete dagmc
delete dagmc;
delete mbi;
}
| void dagmc_load_file_dagmc_build_obb | ( | ) |
Definition at line 81 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, and input_file.
Referenced by main().
{
/* 1 - Test with external moab, load file in DAGMC*/
// make new moab core
ErrorCode rval;
moab::Core* mbi = new moab::Core();
// make new dagmc into that moab
DagMC* dagmc = new moab::DagMC( mbi );
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// delete dagmc
delete dagmc;
delete mbi;
}
| void dagmc_load_file_dagmc_internal | ( | ) |
Definition at line 69 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, and input_file.
Referenced by main().
{
/* 3 - Test with internal moab, load file in DAG*/
// make new dagmc into that moab
ErrorCode rval;
moab::DagMC* dagmc = new moab::DagMC();
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
delete dagmc;
}
Definition at line 116 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, and input_file.
Referenced by main().
{
/* 3 - Test with internal moab, load file in DAG*/
// make new dagmc into that moab
ErrorCode rval;
moab::DagMC* dagmc = new moab::DagMC();
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
delete dagmc;
}
| void dagmc_load_file_dagmc_via_moab | ( | ) |
Definition at line 53 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, input_file, and moab::Core::load_file().
Referenced by main().
{
/* 2 - Test with external moab, load file in MOAB*/
// load the file into moab rather than dagmc
ErrorCode rval;
moab::Core* mbi = new moab::Core();
rval = mbi->load_file( input_file.c_str() );CHECK_ERR( rval );
moab::DagMC* dagmc = new moab::DagMC( mbi );
rval = dagmc->load_existing_contents();CHECK_ERR( rval );
// delete dagmc;
delete dagmc;
delete mbi;
}
Definition at line 99 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, input_file, and moab::Core::load_file().
Referenced by main().
{
/* 2 - Test with external moab, load file in MOAB*/
// load the file into moab rather than dagmc
ErrorCode rval;
moab::Core* mbi = new moab::Core();
rval = mbi->load_file( input_file.c_str() );CHECK_ERR( rval );
moab::DagMC* dagmc = new moab::DagMC( mbi );
rval = dagmc->load_existing_contents();CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// delete dagmc;
delete dagmc;
delete mbi;
}
| void dagmc_num_vols | ( | ) |
Definition at line 161 of file dagmc_simple_test.cpp.
References CHECK_EQUAL, and DAG.
Referenced by main().
{
int expect_num_vols = 2;
int num_vols = DAG->num_entities( 3 );
CHECK_EQUAL( expect_num_vols, num_vols );
}
| void dagmc_point_in | ( | ) |
Definition at line 180 of file dagmc_simple_test.cpp.
References CHECK_EQUAL, CHECK_ERR, DAG, and ErrorCode.
{
int result = 0;
int expect_result = 1;
int vol_idx = 1;
double xyz[3] = { 0.0, 0.0, 0.0 };
EntityHandle vol_h = DAG->entity_by_index( 3, vol_idx );
ErrorCode rval = DAG->point_in_volume( vol_h, xyz, result );CHECK_ERR( rval );
CHECK_EQUAL( expect_result, result );
}
| void dagmc_rayfire | ( | ) |
Definition at line 234 of file dagmc_simple_test.cpp.
References CHECK_ERR, CHECK_REAL_EQUAL, DAG, eps, and ErrorCode.
Referenced by main().
{
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 = DAG->entity_by_index( 3, vol_idx );
ErrorCode rval = DAG->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 );
}
| void dagmc_test_boundary | ( | ) |
Definition at line 269 of file dagmc_simple_test.cpp.
References CHECK_EQUAL, CHECK_ERR, DAG, and ErrorCode.
Referenced by main().
{
int vol_idx = 1;
EntityHandle vol_h = DAG->entity_by_index( 3, vol_idx );
int surf_idx = 1;
EntityHandle surf_h = DAG->entity_by_index( 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;
ErrorCode rval = DAG->test_volume_boundary( vol_h, surf_h, xyz, dir, result );CHECK_ERR( rval );
// check ray leaving volume
CHECK_EQUAL( expect_result, result );
}
| void dagmc_test_obb_retreval | ( | ) |
Definition at line 129 of file dagmc_simple_test.cpp.
References CHECK_ERR, ErrorCode, and input_file.
Referenced by main().
{
// make new dagmc
std::cout << "test_obb_retreval" << std::endl;
DagMC* dagmc = new moab::DagMC();
ErrorCode rval;
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// write the file
rval = dagmc->write_mesh( "fcad", 4 );
// now remove the dagmc instance a
delete dagmc;
dagmc = new moab::DagMC();
rval = dagmc->load_file( "fcad" );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// delete the fcad file
remove( "fcad" );
delete dagmc;
}
| void dagmc_test_obb_retreval_rayfire | ( | ) |
Definition at line 191 of file dagmc_simple_test.cpp.
References CHECK_ERR, CHECK_REAL_EQUAL, DAG, eps, ErrorCode, and input_file.
Referenced by main().
{
// make new dagmc
std::cout << "test_obb_retreval and ray_fire" << std::endl;
DagMC* dagmc = new moab::DagMC();
ErrorCode rval;
// load a file
rval = dagmc->load_file( input_file.c_str() );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// write the file
rval = dagmc->write_mesh( "fcad", 4 );
// now remove the dagmc instance a
delete dagmc;
// now create new DAGMC
dagmc = new moab::DagMC();
rval = dagmc->load_file( "fcad" );CHECK_ERR( rval );
rval = dagmc->init_OBBTree();CHECK_ERR( rval );
// delete the fcad 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 = DAG->entity_by_index( 3, vol_idx );
rval = DAG->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 );
delete dagmc;
}
| int main | ( | int | , |
| char ** | |||
| ) |
Definition at line 286 of file dagmc_simple_test.cpp.
References DAG, dagmc_build_obb(), dagmc_closest_to(), dagmc_load_file(), dagmc_load_file_dagmc(), dagmc_load_file_dagmc_build_obb(), dagmc_load_file_dagmc_internal(), dagmc_load_file_dagmc_internal_build_obb(), dagmc_load_file_dagmc_via_moab(), dagmc_load_file_dagmc_via_moab_build_obb(), dagmc_num_vols(), dagmc_point_in(), dagmc_rayfire(), dagmc_test_boundary(), dagmc_test_obb_retreval(), dagmc_test_obb_retreval_rayfire(), and RUN_TEST.
{
int result = 0;
DAG = new moab::DagMC();
result += RUN_TEST( dagmc_load_file ); // test ray fire
result += RUN_TEST( dagmc_build_obb ); // build the obb
result += RUN_TEST( dagmc_num_vols ); // make sure the num of vols correct
result += RUN_TEST( dagmc_load_file_dagmc ); //
result += RUN_TEST( dagmc_load_file_dagmc_via_moab ); //
result += RUN_TEST( dagmc_load_file_dagmc_internal ); //
result += RUN_TEST( dagmc_load_file_dagmc_build_obb ); //
result += RUN_TEST( dagmc_load_file_dagmc_via_moab_build_obb ); //
result += RUN_TEST( dagmc_load_file_dagmc_internal_build_obb ); //
result += RUN_TEST( dagmc_test_obb_retreval ); // check that we are retreving loaded obbs
result += RUN_TEST( dagmc_test_obb_retreval_rayfire ); // check that we can ray fire on loaded obbs
result += RUN_TEST( dagmc_point_in ); // check entity by point
result += RUN_TEST( dagmc_rayfire ); // ensure ray fire distance is correct
result += RUN_TEST( dagmc_closest_to ); // check the distance to surface nearest point
result += RUN_TEST( dagmc_test_boundary ); // check particle entering leaving
delete DAG;
return result;
}
| DagMC* DAG |
Definition at line 16 of file dagmc_simple_test.cpp.
| const std::string input_file = "unittest/test_geom.h5m" |
Definition at line 28 of file dagmc_simple_test.cpp.
1.7.6.1