Author: tautges Date: 2009-10-29 15:17:58 -0500 (Thu, 29 Oct 2009) New Revision: 3251 Added: MOAB/trunk/parallel/uber_parallel_test.cpp Log: Adding the file I'd fogotten to commit, thanks Jed. Added: MOAB/trunk/parallel/uber_parallel_test.cpp =================================================================== --- MOAB/trunk/parallel/uber_parallel_test.cpp (rev 0) +++ MOAB/trunk/parallel/uber_parallel_test.cpp 2009-10-29 20:17:58 UTC (rev 3251) @@ -0,0 +1,139 @@ +#include "MBParallelComm.hpp" +#include "MBParallelConventions.h" +#include "ReadParallel.hpp" +#include "FileOptions.hpp" +#include "MBTagConventions.hpp" +#include "MBCore.hpp" +#include "mpi.h" +#include <iostream> +#include <algorithm> +#include <sstream> +#include <assert.h> +#if !defined(_MSC_VER) && !defined(__MINGW32__) +#include <unistd.h> +#endif + + +#define STRINGIFY_(X) #X +#define STRINGIFY(X) STRINGIFY_(X) + + +#define CHKERR(a) do { \ + MBErrorCode val = (a); \ + if (MB_SUCCESS != val) { \ + std::cerr << "Error code " << val << " at " << __FILE__ << ":" << __LINE__ << std::endl;\ + return val; \ + } \ +} while (false) + +#define PCHECK(A) if (is_any_proc_error(!(A))) return report_error(__FILE__,__LINE__) + +MBErrorCode report_error( const char* file, int line ) +{ + std::cerr << "Failure at " << file << ':' << line << std::endl; + return MB_FAILURE; +} + +MBErrorCode test_read(const char *filename, const char *option); + +#define RUN_TEST(A, B, C) run_test( &A, #A, B, C) + +int is_any_proc_error( int is_my_error ) +{ + int result = 0; + int err = MPI_Allreduce( &is_my_error, &result, 1, MPI_INT, MPI_MAX, MPI_COMM_WORLD ); + return err || result; +} + +int run_test( MBErrorCode (*func)(const char*, const char*), + const char* func_name, + const char* file_name, + const char *option) +{ + MBErrorCode result = (*func)(file_name, option); + int is_err = is_any_proc_error( (MB_SUCCESS != result) ); + int rank; + MPI_Comm_rank( MPI_COMM_WORLD, &rank ); + if (rank == 0) { + if (is_err) + std::cout << func_name << " : FAILED!!" << std::endl; + else + std::cout << func_name << " : success" << std::endl; + } + + return is_err; +} + +int main( int argc, char* argv[] ) +{ + int rank, size; + MPI_Init(&argc, &argv); + MPI_Comm_rank( MPI_COMM_WORLD, &rank ); + MPI_Comm_size( MPI_COMM_WORLD, &size ); + int num_errors = 0; + + const char* filename, *option; + if (1 < argc) + filename = argv[1]; + else { +#ifdef SRCDIR + filename = STRINGIFY(SRCDIR) "/../test/64bricks_512hex.h5m"; +#else + filename = "../test/64bricks_512hex.h5m"; +#endif + } + + //=========== read_delete, geom_dimension, resolve_shared + option = "PARALLEL=READ_DELETE;PARTITION=GEOM_DIMENSION;PARTITION_VAL=3;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== read_delete, material_set, resolve_shared + option = "PARALLEL=READ_DELETE;PARTITION=MATERIAL_SET;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== bcast_delete, geom_dimension, resolve_shared + option = "PARALLEL=READ_DELETE;PARTITION=GEOM_DIMENSION;PARTITION_VAL=3;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== bcast_delete, material_set, resolve_shared + option = "PARALLEL=READ_DELETE;PARTITION=MATERIAL_SET;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== read_delete, geom_dimension, resolve_shared, exch ghost + option = "PARALLEL=READ_DELETE;PARTITION=GEOM_DIMENSION;PARTITION_VAL=3;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;PARALLEL_GHOSTS=3.0.1;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== read_delete, material_set, resolve_shared, exch ghost + option = "PARALLEL=READ_DELETE;PARTITION=MATERIAL_SET;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;PARALLEL_GHOSTS=3.0.1;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== bcast_delete, geom_dimension, resolve_shared, exch ghost + option = "PARALLEL=READ_DELETE;PARTITION=GEOM_DIMENSION;PARTITION_VAL=3;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;PARALLEL_GHOSTS=3.0.1;"; + num_errors += RUN_TEST( test_read, filename, option ); + + //=========== bcast_delete, material_set, resolve_shared, exch ghost + option = "PARALLEL=READ_DELETE;PARTITION=MATERIAL_SET;PARTITION_DISTRIBUTE;PARALLEL_RESOLVE_SHARED_ENTS;PARALLEL_GHOSTS=3.0.1;"; + num_errors += RUN_TEST( test_read, filename, option ); + + MPI_Finalize(); + + return num_errors; +} + +MBErrorCode test_read(const char *filename, const char *option) +{ + MBCore mb_instance; + MBInterface& moab = mb_instance; + MBErrorCode rval; + MBEntityHandle set; + + rval = moab.load_file( filename, set, option); + CHKERR(rval); + + MBParallelComm* pcomm = MBParallelComm::get_pcomm(&moab, 0); + + rval = pcomm->check_all_shared_handles(); + CHKERR(rval); + + return MB_SUCCESS; +}