|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
#include "TestUtil.hpp"#include "moab/Core.hpp"#include "moab/ParallelComm.hpp"#include "moab/ProgOptions.hpp"#include "MBParallelConventions.h"#include "moab/ReadUtilIface.hpp"#include "MBTagConventions.hpp"#include <sstream>
Include dependency graph for gcrm_par.cpp:Go to the source code of this file.
Functions | |
| void | test_read_onevar_trivial () |
| void | test_read_mesh_parallel_trivial () |
| void | test_gather_onevar_on_rank0 () |
| void | test_gather_onevar_on_rank1 () |
| void | test_multiple_loads_of_same_file () |
| void | read_one_cell_var (bool rcbzoltan) |
| void | read_mesh_parallel (bool rcbzoltan) |
| void | gather_one_cell_var (int gather_set_rank) |
| void | multiple_loads_of_same_file () |
| int | main (int argc, char *argv[]) |
| void | test_read_onevar_rcbzoltan () |
| void | test_read_mesh_parallel_rcbzoltan () |
Variables | |
| std::string | example = "unittest/io/gcrm_r3.nc" |
| std::string | read_options |
| const double | eps = 1e-6 |
| const int | layers = 3 |
| void gather_one_cell_var | ( | int | gather_set_rank | ) |
Definition at line 470 of file gcrm_par.cpp.
References CHECK_EQUAL, CHECK_ERR, CHECK_REAL_EQUAL, moab::Interface::create_meshset(), eps, ErrorCode, example, moab::ParallelComm::filter_pstatus(), moab::ParallelComm::gather_data(), moab::Interface::get_entities_by_type(), moab::ReadUtilIface::get_gather_set(), moab::ParallelComm::get_pcomm(), moab::Interface::globalId_tag(), layers, moab::Interface::load_file(), mb, MB_TAG_DENSE, MB_TYPE_DOUBLE, MBPOLYGON, MESHSET_SET, moab::ParallelComm::proc_config(), moab::ProcConfig::proc_rank(), moab::ProcConfig::proc_size(), moab::Range::psize(), PSTATUS_NOT, PSTATUS_NOT_OWNED, moab::Interface::query_interface(), rank, read_options, moab::Range::size(), moab::Interface::tag_get_data(), and moab::Interface::tag_get_handle().
Referenced by test_gather_onevar_on_rank0(), and test_gather_onevar_on_rank1().
{
Core moab;
Interface& mb = moab;
EntityHandle file_set;
ErrorCode rval = mb.create_meshset( MESHSET_SET, file_set );CHECK_ERR( rval );
read_options = "PARALLEL=READ_PART;PARTITION_METHOD=TRIVIAL;PARALLEL_RESOLVE_SHARED_ENTS";
std::ostringstream gather_set_option;
gather_set_option << ";GATHER_SET=" << gather_set_rank;
read_options += gather_set_option.str();
rval = mb.load_file( example.c_str(), &file_set, read_options.c_str() );CHECK_ERR( rval );
ParallelComm* pcomm = ParallelComm::get_pcomm( &mb, 0 );
int procs = pcomm->proc_config().proc_size();
int rank = pcomm->proc_config().proc_rank();
// Make sure gather_set_rank is valid
if( gather_set_rank < 0 || gather_set_rank >= procs ) return;
Range cells, cells_owned;
rval = mb.get_entities_by_type( file_set, MBPOLYGON, cells );CHECK_ERR( rval );
// Get local owned cells
rval = pcomm->filter_pstatus( cells, PSTATUS_NOT_OWNED, PSTATUS_NOT, -1, &cells_owned );CHECK_ERR( rval );
EntityHandle gather_set = 0;
if( gather_set_rank == rank )
{
// Get gather set
ReadUtilIface* readUtilIface;
mb.query_interface( readUtilIface );
rval = readUtilIface->get_gather_set( gather_set );CHECK_ERR( rval );
assert( gather_set != 0 );
}
Tag vorticity_tag0, gid_tag;
rval = mb.tag_get_handle( "vorticity0", layers, MB_TYPE_DOUBLE, vorticity_tag0, MB_TAG_DENSE );CHECK_ERR( rval );
gid_tag = mb.globalId_tag();
pcomm->gather_data( cells_owned, vorticity_tag0, gid_tag, gather_set, gather_set_rank );
if( gather_set_rank == rank )
{
// Get gather set cells
Range gather_set_cells;
rval = mb.get_entities_by_type( gather_set, MBPOLYGON, gather_set_cells );CHECK_ERR( rval );
CHECK_EQUAL( (size_t)642, gather_set_cells.size() );
CHECK_EQUAL( (size_t)1, gather_set_cells.psize() );
// Check vorticity0 tag values on 4 gather set cells: first cell, two median cells, and last
// cell
EntityHandle cell_ents[] = { gather_set_cells[0], gather_set_cells[320], gather_set_cells[321],
gather_set_cells[641] };
double vorticity0_val[4 * layers];
rval = mb.tag_get_data( vorticity_tag0, &cell_ents[0], 4, vorticity0_val );CHECK_ERR( rval );
// Only check first two layers
// Layer 0
CHECK_REAL_EQUAL( 3.629994, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( 0.131688, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[2 * layers], eps );
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[3 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.629944, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.131686, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[2 * layers + 1], eps );
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[3 * layers + 1], eps );
}
}
| int main | ( | int | argc, |
| char * | argv[] | ||
| ) |
Definition at line 40 of file gcrm_par.cpp.
References RUN_TEST, test_gather_onevar_on_rank0(), test_gather_onevar_on_rank1(), test_multiple_loads_of_same_file(), test_read_mesh_parallel_rcbzoltan(), test_read_mesh_parallel_trivial(), test_read_onevar_rcbzoltan(), and test_read_onevar_trivial().
{
MPI_Init( &argc, &argv );
int result = 0;
result += RUN_TEST( test_read_onevar_trivial );
#if defined( MOAB_HAVE_MPI ) && defined( MOAB_HAVE_ZOLTAN )
result += RUN_TEST( test_read_onevar_rcbzoltan );
#endif
result += RUN_TEST( test_read_mesh_parallel_trivial );
#if defined( MOAB_HAVE_MPI ) && defined( MOAB_HAVE_ZOLTAN )
result += RUN_TEST( test_read_mesh_parallel_rcbzoltan );
#endif
result += RUN_TEST( test_gather_onevar_on_rank0 );
result += RUN_TEST( test_gather_onevar_on_rank1 );
result += RUN_TEST( test_multiple_loads_of_same_file );
MPI_Finalize();
return result;
}
| void multiple_loads_of_same_file | ( | ) |
Definition at line 544 of file gcrm_par.cpp.
References CHECK_EQUAL, CHECK_ERR, CHECK_REAL_EQUAL, moab::Interface::create_meshset(), eps, ErrorCode, example, moab::Interface::get_entities_by_type(), moab::ParallelComm::get_pcomm(), layers, moab::Interface::load_file(), mb, MB_TYPE_DOUBLE, MBEDGE, MBPOLYGON, MBVERTEX, MESHSET_SET, moab::Range::psize(), rank, read_options, moab::Range::size(), moab::Interface::tag_get_data(), and moab::Interface::tag_get_handle().
Referenced by test_multiple_loads_of_same_file(), and test_multiple_loads_of_same_file_no_mixed_elements().
{
Core moab;
Interface& mb = moab;
// Need a file set for nomesh to work right
EntityHandle file_set;
ErrorCode rval;
rval = mb.create_meshset( MESHSET_SET, file_set );CHECK_ERR( rval );
// Read first only header information, no mesh, no variable
read_options = "PARALLEL=READ_PART;PARTITION;NOMESH;VARIABLE=;PARTITION_METHOD=TRIVIAL";
rval = mb.load_file( example.c_str(), &file_set, read_options.c_str() );CHECK_ERR( rval );
// Create mesh, no variable
read_options = "PARALLEL=READ_PART;PARTITION;PARALLEL_RESOLVE_SHARED_ENTS;PARTITION_METHOD="
"TRIVIAL;VARIABLE=";
rval = mb.load_file( example.c_str(), &file_set, read_options.c_str() );CHECK_ERR( rval );
// Read variable vorticity at timestep 0, no mesh
read_options = "PARALLEL=READ_PART;PARTITION;PARTITION_METHOD=TRIVIAL;NOMESH;VARIABLE="
"vorticity;TIMESTEP=0";
rval = mb.load_file( example.c_str(), &file_set, read_options.c_str() );CHECK_ERR( rval );
Range local_verts;
rval = mb.get_entities_by_type( file_set, MBVERTEX, local_verts );CHECK_ERR( rval );
Range local_edges;
rval = mb.get_entities_by_type( file_set, MBEDGE, local_edges );CHECK_ERR( rval );
Range local_cells;
rval = mb.get_entities_by_type( file_set, MBPOLYGON, local_cells );CHECK_ERR( rval );
// No mixed elements
CHECK_EQUAL( (size_t)1, local_cells.psize() );
ParallelComm* pcomm = ParallelComm::get_pcomm( &mb, 0 );
int procs = pcomm->proc_config().proc_size();
int rank = pcomm->proc_config().proc_rank();
// Make check runs this test on two processors
if( 2 == procs )
{
CHECK_EQUAL( (size_t)321, local_cells.size() );
// Check tag for cell variable vorticity at timestep 0
Tag vorticity_tag0;
rval = mb.tag_get_handle( "vorticity0", layers, MB_TYPE_DOUBLE, vorticity_tag0 );CHECK_ERR( rval );
// Get vorticity0 tag values on 3 local cells
double vorticity0_val[3 * layers];
EntityHandle cell_ents[] = { local_cells[0], local_cells[160], local_cells[320] };
rval = mb.tag_get_data( vorticity_tag0, cell_ents, 3, vorticity0_val );CHECK_ERR( rval );
if( 0 == rank )
{
CHECK_EQUAL( (size_t)687, local_verts.size() );
CHECK_EQUAL( (size_t)1007, local_edges.size() );
// Layer 0
CHECK_REAL_EQUAL( 3.629994, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.708188, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( 0.131688, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.629944, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.708164, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.131686, vorticity0_val[2 * layers + 1], eps );
}
else if( 1 == rank )
{
CHECK_EQUAL( (size_t)688, local_verts.size() );
CHECK_EQUAL( (size_t)1008, local_edges.size() );
// Layer 0
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( 2.434397, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( 2.434363, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[2 * layers + 1], eps );
}
}
}
| void read_mesh_parallel | ( | bool | rcbzoltan | ) |
Definition at line 288 of file gcrm_par.cpp.
References moab::ParallelComm::check_all_shared_handles(), CHECK_EQUAL, CHECK_ERR, ErrorCode, example, moab::ParallelComm::filter_pstatus(), moab::Interface::get_entities_by_type(), moab::ParallelComm::get_pcomm(), moab::Interface::load_file(), mb, MBEDGE, MBPOLYGON, MBVERTEX, moab::ProcConfig::proc_comm(), moab::ParallelComm::proc_config(), moab::ProcConfig::proc_rank(), moab::ProcConfig::proc_size(), moab::Range::psize(), PSTATUS_NOT, PSTATUS_NOT_OWNED, rank, read_options, moab::Range::size(), and moab::Interface::write_file().
Referenced by test_read_mesh_parallel_rcbzoltan(), test_read_mesh_parallel_rcbzoltan_no_mixed_elements(), test_read_mesh_parallel_trivial(), and test_read_mesh_parallel_trivial_no_mixed_elements().
{
Core moab;
Interface& mb = moab;
read_options = "PARALLEL=READ_PART;PARTITION_METHOD=TRIVIAL;PARALLEL_RESOLVE_SHARED_ENTS;VARIABLE=";
if( rcbzoltan )
read_options = "PARALLEL=READ_PART;PARTITION_METHOD=RCBZOLTAN;PARALLEL_RESOLVE_SHARED_ENTS;VARIABLE=";
ErrorCode rval = mb.load_file( example.c_str(), NULL, read_options.c_str() );CHECK_ERR( rval );
ParallelComm* pcomm = ParallelComm::get_pcomm( &mb, 0 );
int procs = pcomm->proc_config().proc_size();
int rank = pcomm->proc_config().proc_rank();
rval = pcomm->check_all_shared_handles();CHECK_ERR( rval );
// Get local vertices
Range local_verts;
rval = mb.get_entities_by_type( 0, MBVERTEX, local_verts );CHECK_ERR( rval );
int verts_num = local_verts.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 688, verts_num );
else if( 1 == rank )
CHECK_EQUAL( 687, verts_num ); // Not owned vertices included
}
else
{
if( 0 == rank )
CHECK_EQUAL( 687, verts_num );
else if( 1 == rank )
CHECK_EQUAL( 688, verts_num ); // Not owned vertices included
}
}
rval = pcomm->filter_pstatus( local_verts, PSTATUS_NOT_OWNED, PSTATUS_NOT );CHECK_ERR( rval );
verts_num = local_verts.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 688, verts_num );
else if( 1 == rank )
CHECK_EQUAL( 592, verts_num ); // Not owned vertices excluded
}
else
{
if( 0 == rank )
CHECK_EQUAL( 687, verts_num );
else if( 1 == rank )
CHECK_EQUAL( 593, verts_num ); // Not owned vertices excluded
}
}
// Get local edges
Range local_edges;
rval = mb.get_entities_by_type( 0, MBEDGE, local_edges );CHECK_ERR( rval );
int edges_num = local_edges.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 1008, edges_num );
else if( 1 == rank )
CHECK_EQUAL( 1007, edges_num ); // Not owned edges included
}
else
{
if( 0 == rank )
CHECK_EQUAL( 1007, edges_num );
else if( 1 == rank )
CHECK_EQUAL( 1008, edges_num ); // Not owned edges included
}
}
rval = pcomm->filter_pstatus( local_edges, PSTATUS_NOT_OWNED, PSTATUS_NOT );CHECK_ERR( rval );
edges_num = local_edges.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 1008, edges_num );
else if( 1 == rank )
CHECK_EQUAL( 912, edges_num ); // Not owned edges excluded
}
else
{
if( 0 == rank )
CHECK_EQUAL( 1007, edges_num );
else if( 1 == rank )
CHECK_EQUAL( 913, edges_num ); // Not owned edges excluded
}
}
// Get local cells
Range local_cells;
rval = mb.get_entities_by_type( 0, MBPOLYGON, local_cells );CHECK_ERR( rval );
// No mixed elements
CHECK_EQUAL( (size_t)1, local_cells.psize() );
int cells_num = local_cells.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 321, cells_num );
else
CHECK_EQUAL( 321, cells_num );
}
else
CHECK_EQUAL( 321, cells_num );
}
rval = pcomm->filter_pstatus( local_cells, PSTATUS_NOT_OWNED, PSTATUS_NOT );CHECK_ERR( rval );
cells_num = local_cells.size();
if( 2 == procs )
{
if( rcbzoltan )
{
if( 0 == rank )
CHECK_EQUAL( 321, cells_num );
else
CHECK_EQUAL( 321, cells_num );
}
else
CHECK_EQUAL( 321, cells_num );
}
std::cout << "proc: " << rank << " verts:" << verts_num << "\n";
int total_verts_num;
MPI_Reduce( &verts_num, &total_verts_num, 1, MPI_INT, MPI_SUM, 0, pcomm->proc_config().proc_comm() );
if( 0 == rank )
{
std::cout << "total vertices: " << total_verts_num << "\n";
CHECK_EQUAL( 1280, total_verts_num );
}
std::cout << "proc: " << rank << " edges:" << edges_num << "\n";
int total_edges_num;
MPI_Reduce( &edges_num, &total_edges_num, 1, MPI_INT, MPI_SUM, 0, pcomm->proc_config().proc_comm() );
if( 0 == rank )
{
std::cout << "total edges: " << total_edges_num << "\n";
CHECK_EQUAL( 1920, total_edges_num );
}
std::cout << "proc: " << rank << " cells:" << cells_num << "\n";
int total_cells_num;
MPI_Reduce( &cells_num, &total_cells_num, 1, MPI_INT, MPI_SUM, 0, pcomm->proc_config().proc_comm() );
if( 0 == rank )
{
std::cout << "total cells: " << total_cells_num << "\n";
CHECK_EQUAL( 642, total_cells_num );
}
#ifdef MOAB_HAVE_HDF5_PARALLEL
std::string write_options( "PARALLEL=WRITE_PART;" );
std::string output_file = "test_gcrm";
if( rcbzoltan ) output_file += "_rcbzoltan";
output_file += ".h5m";
mb.write_file( output_file.c_str(), NULL, write_options.c_str() );
#endif
}
| void read_one_cell_var | ( | bool | rcbzoltan | ) |
Definition at line 100 of file gcrm_par.cpp.
References CHECK_EQUAL, CHECK_ERR, CHECK_REAL_EQUAL, eps, ErrorCode, example, moab::Interface::get_entities_by_type(), moab::ParallelComm::get_pcomm(), moab::Interface::globalId_tag(), layers, moab::Interface::load_file(), mb, MB_TYPE_DOUBLE, MBEDGE, MBPOLYGON, moab::Range::psize(), rank, read_options, moab::Range::size(), moab::Interface::tag_get_data(), and moab::Interface::tag_get_handle().
Referenced by test_read_onevar_rcbzoltan(), test_read_onevar_rcbzoltan_no_mixed_elements(), test_read_onevar_trivial(), and test_read_onevar_trivial_no_mixed_elements().
{
Core moab;
Interface& mb = moab;
read_options = "PARALLEL=READ_PART;PARTITION_METHOD=TRIVIAL;NO_EDGES;VARIABLE=vorticity";
if( rcbzoltan )
read_options = "PARALLEL=READ_PART;PARTITION_METHOD=RCBZOLTAN;NO_EDGES;VARIABLE=vorticity;DEBUG_IO=1";
ErrorCode rval = mb.load_file( example.c_str(), NULL, read_options.c_str() );CHECK_ERR( rval );
// Get local edges
Range local_edges;
rval = mb.get_entities_by_type( 0, MBEDGE, local_edges );CHECK_ERR( rval );
CHECK_EQUAL( (size_t)0, local_edges.size() );
// Get local cells
Range local_cells;
rval = mb.get_entities_by_type( 0, MBPOLYGON, local_cells );CHECK_ERR( rval );
// No mixed elements
CHECK_EQUAL( (size_t)1, local_cells.psize() );
Tag gid_tag = mb.globalId_tag();
std::vector< int > gids( local_cells.size() );
rval = mb.tag_get_data( gid_tag, local_cells, &gids[0] );
Range local_cell_gids;
std::copy( gids.rbegin(), gids.rend(), range_inserter( local_cell_gids ) );
ParallelComm* pcomm = ParallelComm::get_pcomm( &mb, 0 );
int procs = pcomm->proc_config().proc_size();
int rank = pcomm->proc_config().proc_rank();
// Make check runs this test on two processors
if( 2 == procs )
{
// Check tag for cell variable vorticity at timestep 0
Tag vorticity_tag0;
rval = mb.tag_get_handle( "vorticity0", layers, MB_TYPE_DOUBLE, vorticity_tag0 );CHECK_ERR( rval );
// Check tag for cell variable vorticity at timestep 1
Tag vorticity_tag1;
rval = mb.tag_get_handle( "vorticity1", layers, MB_TYPE_DOUBLE, vorticity_tag1 );CHECK_ERR( rval );
// Get vorticity0 and vorticity1 tag values on 3 local cells
double vorticity0_val[3 * layers];
double vorticity1_val[3 * layers];
if( rcbzoltan )
{
CHECK_EQUAL( (size_t)14, local_cell_gids.psize() );
if( 0 == rank )
{
CHECK_EQUAL( (size_t)321, local_cells.size() );
CHECK_EQUAL( (size_t)321, local_cell_gids.size() );
CHECK_EQUAL( 3, (int)local_cell_gids[0] );
CHECK_EQUAL( 162, (int)local_cell_gids[159] );
CHECK_EQUAL( 640, (int)local_cell_gids[318] );
EntityHandle cell_ents[] = { local_cells[0], local_cells[159], local_cells[318] };
rval = mb.tag_get_data( vorticity_tag0, cell_ents, 3, vorticity0_val );CHECK_ERR( rval );
// Timestep 0
// Layer 0
CHECK_REAL_EQUAL( -0.725999, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.814997, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( 0.131688, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( -0.725989, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.814972, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.131686, vorticity0_val[2 * layers + 1], eps );
rval = mb.tag_get_data( vorticity_tag1, cell_ents, 3, vorticity1_val );CHECK_ERR( rval );
// Timestep 1
// Layer 0
CHECK_REAL_EQUAL( -0.706871, vorticity1_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.767178, vorticity1_val[1 * layers], eps );
CHECK_REAL_EQUAL( 0.128218, vorticity1_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( -0.706861, vorticity1_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.767153, vorticity1_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.128216, vorticity1_val[2 * layers + 1], eps );
}
else if( 1 == rank )
{
CHECK_EQUAL( (size_t)321, local_cells.size() );
CHECK_EQUAL( (size_t)321, local_cell_gids.size() );
CHECK_EQUAL( 1, (int)local_cell_gids[0] );
CHECK_EQUAL( 366, (int)local_cell_gids[161] );
CHECK_EQUAL( 1, (int)local_cell_gids[322] );
EntityHandle cell_ents[] = { local_cells[0], local_cells[161], local_cells[322] };
rval = mb.tag_get_data( vorticity_tag0, cell_ents, 3, vorticity0_val );CHECK_ERR( rval );
// Timestep 0
// Layer 0
CHECK_REAL_EQUAL( 3.629994, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.512985, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( 3.629994, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.629944, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.512964, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 3.629944, vorticity0_val[2 * layers + 1], eps );
rval = mb.tag_get_data( vorticity_tag1, cell_ents, 3, vorticity1_val );CHECK_ERR( rval );
// Timestep 1
// Layer 0
CHECK_REAL_EQUAL( 3.534355, vorticity1_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.473122, vorticity1_val[1 * layers], eps );
CHECK_REAL_EQUAL( 3.534355, vorticity1_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.534306, vorticity1_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.473102, vorticity1_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 3.534306, vorticity1_val[2 * layers + 1], eps );
}
}
else
{
CHECK_EQUAL( (size_t)321, local_cells.size() );
CHECK_EQUAL( (size_t)321, local_cell_gids.size() );
CHECK_EQUAL( (size_t)1, local_cell_gids.psize() );
EntityHandle cell_ents[] = { local_cells[0], local_cells[160], local_cells[320] };
rval = mb.tag_get_data( vorticity_tag0, cell_ents, 3, vorticity0_val );CHECK_ERR( rval );
rval = mb.tag_get_data( vorticity_tag1, cell_ents, 3, vorticity1_val );CHECK_ERR( rval );
if( 0 == rank )
{
CHECK_EQUAL( 1, (int)local_cell_gids[0] );
CHECK_EQUAL( 161, (int)local_cell_gids[160] );
CHECK_EQUAL( 321, (int)local_cell_gids[320] );
// Timestep 0
// Layer 0
CHECK_REAL_EQUAL( 3.629994, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.708188, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( 0.131688, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.629944, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.708164, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.131686, vorticity0_val[2 * layers + 1], eps );
// Timestep 1
// Layer 0
CHECK_REAL_EQUAL( 3.534355, vorticity1_val[0 * layers], eps );
CHECK_REAL_EQUAL( -1.663182, vorticity1_val[1 * layers], eps );
CHECK_REAL_EQUAL( 0.128218, vorticity1_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( 3.534306, vorticity1_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( -1.663160, vorticity1_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( 0.128216, vorticity1_val[2 * layers + 1], eps );
}
else if( 1 == rank )
{
CHECK_EQUAL( 322, (int)local_cell_gids[0] );
CHECK_EQUAL( 482, (int)local_cell_gids[160] );
CHECK_EQUAL( 642, (int)local_cell_gids[320] );
// Timestep 0
// Layer 0
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[0 * layers], eps );
CHECK_REAL_EQUAL( 2.434397, vorticity0_val[1 * layers], eps );
CHECK_REAL_EQUAL( -0.554888, vorticity0_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( 2.434363, vorticity0_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( -0.554881, vorticity0_val[2 * layers + 1], eps );
// Timestep 1
// Layer 0
CHECK_REAL_EQUAL( -0.540269, vorticity1_val[0 * layers], eps );
CHECK_REAL_EQUAL( 2.370258, vorticity1_val[1 * layers], eps );
CHECK_REAL_EQUAL( -0.540269, vorticity1_val[2 * layers], eps );
// Layer 1
CHECK_REAL_EQUAL( -0.540262, vorticity1_val[0 * layers + 1], eps );
CHECK_REAL_EQUAL( 2.370226, vorticity1_val[1 * layers + 1], eps );
CHECK_REAL_EQUAL( -0.540262, vorticity1_val[2 * layers + 1], eps );
}
}
}
}
| void test_gather_onevar_on_rank0 | ( | ) |
Definition at line 84 of file gcrm_par.cpp.
References gather_one_cell_var().
Referenced by main().
{
gather_one_cell_var( 0 );
}
| void test_gather_onevar_on_rank1 | ( | ) |
Definition at line 89 of file gcrm_par.cpp.
References gather_one_cell_var().
Referenced by main().
{
gather_one_cell_var( 1 );
}
| void test_multiple_loads_of_same_file | ( | ) |
Definition at line 94 of file gcrm_par.cpp.
References multiple_loads_of_same_file().
Referenced by main().
{
multiple_loads_of_same_file();
}
| void test_read_mesh_parallel_rcbzoltan | ( | ) |
Definition at line 79 of file gcrm_par.cpp.
References read_mesh_parallel().
Referenced by main().
{
read_mesh_parallel( true );
}
| void test_read_mesh_parallel_trivial | ( | ) |
Definition at line 74 of file gcrm_par.cpp.
References read_mesh_parallel().
Referenced by main().
{
read_mesh_parallel( false );
}
| void test_read_onevar_rcbzoltan | ( | ) |
Definition at line 69 of file gcrm_par.cpp.
References read_one_cell_var().
Referenced by main().
{
read_one_cell_var( true );
}
| void test_read_onevar_trivial | ( | ) |
Definition at line 64 of file gcrm_par.cpp.
References read_one_cell_var().
Referenced by main().
{
read_one_cell_var( false );
}
| const double eps = 1e-6 |
Definition at line 37 of file gcrm_par.cpp.
| std::string example = "unittest/io/gcrm_r3.nc" |
Definition at line 13 of file gcrm_par.cpp.
| const int layers = 3 |
Definition at line 38 of file gcrm_par.cpp.
| std::string read_options |
Definition at line 36 of file gcrm_par.cpp.
Referenced by adj_perf(), ahf_test(), closedsurface_uref_hirec_convergence_study(), gather_one_cell_var(), iMOAB_LoadMesh(), load_meshset_hirec(), main(), multiple_loads_of_same_file(), perf_inmesh(), read_mesh_parallel(), read_one_cell_var(), test_closedsurface_mesh(), test_mesh(), and test_new_pcomm_instance().
1.7.6.1