test_determinism.cpp File Reference

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

Include dependency graph for test_determinism.cpp:

Go to the source code of this file.

Functions
int	main ()

---

Function Documentation

int main

(

)

Definition at line 37 of file test_determinism.cpp.

References moab::Range::begin(), CHECK, CHECK_EQUAL, CHECK_ERR, moab::Range::clear(), moab::Interface::create_element(), moab::Interface::delete_entities(), moab::Range::end(), ErrorCode, filename, moab::Interface::get_adjacencies(), moab::Interface::get_connectivity(), moab::Interface::get_entities_by_dimension(), moab::Range::insert(), moab::Interface::load_file(), MBTET, and moab::Interface::UNION.

{
    Interface* iface0 = new Core;
    Interface* iface1 = new Core;

    std::string filename = TestDir + "unittest/tet_brick.vtk";

    ErrorCode err;
    err = iface0->load_file( filename.c_str() );CHECK_ERR( err );
    err = iface1->load_file( filename.c_str() );CHECK_ERR( err );

    Range tets0, verts0, edges0, edges1, tris0, tris1;
    err = iface0->get_entities_by_dimension( 0, 3, tets0 );CHECK_ERR( err );
    err = iface0->get_adjacencies( tets0, 1, true, edges0, Interface::UNION );CHECK_ERR( err );
    err = iface1->get_adjacencies( tets0, 1, true, edges1, Interface::UNION );CHECK_ERR( err );
    CHECK_EQUAL( edges0, edges1 );
    err = iface0->get_adjacencies( tets0, 2, true, tris0, Interface::UNION );CHECK_ERR( err );
    err = iface1->get_adjacencies( tets0, 2, true, tris1, Interface::UNION );CHECK_ERR( err );
    CHECK_EQUAL( tris0, tris1 );

    std::vector< EntityHandle > conn_squence;

    // At that point iface0 and iface1 should have the same entity handler
    // associated with entities, so we can use one or another, no difference.

    int repeat = 0;
    for( ; repeat != 3; repeat++ )
    {

        std::vector< EntityHandle > conn_seq;

        Range to_delete;
        // Build range of quad to delete form iface0
        int ii = 0;
        for( Range::iterator qit = tets0.begin(); qit != tets0.end(); qit++, ii++ )
        {
            if( ii % 3 )
            {
                to_delete.insert( *qit );
                const EntityHandle* conn;
                int number_nodes = 0;
                iface0->get_connectivity( *qit, conn, number_nodes );
                conn_seq.insert( conn_seq.end(), conn, &conn[number_nodes] );
            }
        }
        // Buidl tange of edges to delete from iface1
        // Create gaps in sequence, to be filled later on.
        for( Range::iterator eit = edges1.begin(); eit != edges1.end(); eit++, ii++ )
        {
            if( ii % 3 )
            {
                to_delete.insert( *eit );
            }
        }
        for( Range::iterator eit = tris1.begin(); eit != tris1.end(); eit++, ii++ )
        {
            if( ii % 3 )
            {
                to_delete.insert( *eit );
            }
        }

        err = iface0->delete_entities( to_delete );
        err = iface1->delete_entities( to_delete );

        for( unsigned qq = 0; qq != conn_seq.size() / 4; qq++ )
        {
            EntityHandle q0, q1;
            err = iface1->create_element( MBTET, &conn_seq[4 * qq], 4, q1 );CHECK_ERR( err );
            err = iface0->create_element( MBTET, &conn_seq[4 * qq], 4, q0 );CHECK_ERR( err );
            CHECK( q0 == q1 );
        }

        tets0.clear();
        err = iface0->get_entities_by_dimension( 0, 3, tets0 );CHECK_ERR( err );
        Range tets1;
        err = iface1->get_entities_by_dimension( 0, 3, tets1 );CHECK_ERR( err );
        CHECK_EQUAL( tets0, tets1 );

        // Finally check adjacency, this finally check if code is deterministic
        for( Range::iterator tit = tets1.begin(); tit != tets1.end(); tit++ )
        {
            std::vector< EntityHandle > adj0, adj1;
            err = iface0->get_adjacencies( &*tit, 1, 1, true, adj0 );CHECK_ERR( err );
            err = iface1->get_adjacencies( &*tit, 1, 1, true, adj1 );CHECK_ERR( err );
            CHECK_EQUAL( adj0, adj1 );
        }
        for( Range::iterator tit = tets1.begin(); tit != tets1.end(); tit++ )
        {
            std::vector< EntityHandle > adj0, adj1;
            err = iface0->get_adjacencies( &*tit, 1, 2, true, adj0 );CHECK_ERR( err );
            err = iface1->get_adjacencies( &*tit, 1, 2, true, adj1 );CHECK_ERR( err );
            CHECK_EQUAL( adj0, adj1 );
        }
    }

    delete iface0;
    delete iface1;

    return 0;
}