adj_time.cpp

Go to the documentation of this file.

#include "moab/Core.hpp"
#include "moab/Range.hpp"
#include "TestUtil.hpp"
#include <iostream>
#include <cassert>
#include <ctime>

using namespace moab;

static const std::string default_input_file = TestDir + "unittest/64bricks_512hex.h5m";

int main()
{
    ErrorCode rval;
    Core moab;
    Interface& mb = moab;

    // load test file
    rval = mb.load_file( default_input_file.c_str() );
    if( MB_SUCCESS != rval )
    {
        std::cerr << default_input_file << ": failed to load file." << std::endl;
        return 1;
    }

    // get all region elements
    Range vols;
    rval = mb.get_entities_by_dimension( 0, 3, vols );
    if( MB_SUCCESS != rval ) return 2;
    if( vols.empty() ) return 1;

    // create internal face elements
    Range faces;
    rval = mb.get_adjacencies( vols, 2, true, faces, Interface::UNION );
    if( MB_SUCCESS != rval ) return 2;
    assert( faces.size() > vols.size() );

    // time query of all adjacent volumes
    std::vector< EntityHandle > adj;
    clock_t t_0 = clock();
    for( Range::iterator i = faces.begin(); i != faces.end(); ++i )
    {
        adj.clear();
        rval = mb.get_adjacencies( &*i, 1, 3, false, adj );
        if( MB_SUCCESS != rval ) return 2;
        assert( adj.size() == 1 || adj.size() == 2 );
    }
    clock_t t_up = clock() - t_0;
    std::cout << "Querying of volumes for " << faces.size() << " faces: " << t_up / (double)CLOCKS_PER_SEC << " seconds"
              << std::endl;

    // time downward adjacency query from volumes to faces
    t_0 = clock();
    for( Range::iterator i = vols.begin(); i != vols.end(); ++i )
    {
        adj.clear();
        rval = mb.get_adjacencies( &*i, 1, 1, false, adj );
        if( MB_SUCCESS != rval ) return 2;
        assert( adj.size() > 3 );
    }
    clock_t t_down = clock() - t_0;
    std::cout << "Querying of faces for " << vols.size() << " volumes: " << t_down / (double)CLOCKS_PER_SEC
              << " seconds" << std::endl;

    return 0;
}