h5regression.cpp

Go to the documentation of this file.

#include "moab/Core.hpp"
#include "TestUtil.hpp"
#include "moab/Range.hpp"
#include "moab/ReadUtilIface.hpp"
#include "WriteHDF5.hpp"
#include "moab/FileOptions.hpp"

#ifdef MOAB_HAVE_MPI
#include "moab_mpi.h"
#endif

#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstdlib>
#include <cmath>

using namespace moab;

const char filename[] = "bad.h5m";

void test_write_invalid_elem();
void test_write_read_many_tags();

int main( int argc, char* argv[] )
{
#ifdef MOAB_HAVE_MPI
    int fail = MPI_Init( &argc, &argv );
    if( fail ) return fail;
#else
    argv[0] = argv[argc - argc];  // warning in serial
#endif

    int exitval = 0;
    exitval += RUN_TEST( test_write_invalid_elem );
    exitval += RUN_TEST( test_write_read_many_tags );

#ifdef MOAB_HAVE_MPI
    fail = MPI_Finalize();
    if( fail ) return fail;
#endif

    return exitval;
}

void test_write_invalid_elem()
{
    Core mbcore;
    Interface& moab         = mbcore;
    ReadUtilIface* readtool = 0;
    ErrorCode rval;

    rval = moab.query_interface( readtool );CHECK_ERR( rval );
    CHECK( readtool != 0 );

    // create two nodes
    EntityHandle first_node;
    std::vector< double* > coords;
    rval = readtool->get_node_coords( 3, 2, 1, first_node, coords );CHECK_ERR( rval );
    // Cppcheck warning (false positive): variable coords is assigned a value that is never used
    coords[0][0] = coords[0][1] = 0.0;
    coords[1][0] = coords[1][1] = 0.0;
    coords[2][0] = coords[2][1] = 0.0;

    // create a triangle with an invalid node handle for its
    // third vertex
    EntityHandle tri;
    EntityHandle* conn = 0;
    rval               = readtool->get_element_connect( 1, 3, MBTRI, 1, tri, conn );CHECK_ERR( rval );
    conn[0] = first_node;      // valid
    conn[1] = first_node + 1;  // valid
    conn[2] = first_node + 2;  // invalid

    // try to write the file (should fail)
    WriteHDF5 writer( &moab );
    FileOptions opts( 0 );
    rval = writer.write_file( filename, true, opts, 0, 0, std::vector< std::string >() );
    CHECK( MB_SUCCESS != rval );
}

void test_write_read_many_tags()
{
    const int N = 200;
    Core mbcore;
    Interface& mb = mbcore;
    ErrorCode rval;

    double coords[3] = { 0, 0, 0 };
    EntityHandle node;
    rval = mb.create_vertex( coords, node );CHECK_ERR( rval );

    // create a lot of tags
    std::vector< Tag > tags;
    for( int i = 0; i < N; ++i )
    {
        Tag t;
        std::ostringstream name( "IntTag" );
        name << i;
        rval = mb.tag_get_handle( name.str().c_str(), 1, MB_TYPE_INTEGER, t,
                                  ( i % 2 ? MB_TAG_SPARSE : MB_TAG_DENSE ) | MB_TAG_EXCL, &i );CHECK_ERR( rval );
        tags.push_back( t );
    }

    // write the file
    rval = mb.write_file( filename, "MOAB" );CHECK_ERR( rval );

    // clear moab instance
    rval = mb.delete_mesh();CHECK_ERR( rval );
    for( int i = 0; i < N; ++i )
    {
        rval = mb.tag_delete( tags[i] );CHECK_ERR( rval );
    }

    // read the file
    rval = mb.load_file( filename );CHECK_ERR( rval );
    remove( filename );

    // check that we have the expected tags
    for( int i = 0; i < N; ++i )
    {
        Tag t;
        std::ostringstream name( "IntTag" );
        name << i;
        rval = mb.tag_get_handle( name.str().c_str(), 1, MB_TYPE_INTEGER, t );CHECK_ERR( rval );

        TagType storage;
        rval = mb.tag_get_type( t, storage );CHECK_ERR( rval );
        CHECK_EQUAL( ( i % 2 ) ? MB_TAG_SPARSE : MB_TAG_DENSE, storage );

        DataType type;
        rval = mb.tag_get_data_type( t, type );CHECK_ERR( rval );
        CHECK_EQUAL( MB_TYPE_INTEGER, type );

        int size;
        rval = mb.tag_get_length( t, size );CHECK_ERR( rval );
        CHECK_EQUAL( 1, size );

        int def;
        rval = mb.tag_get_default_value( t, &def );CHECK_ERR( rval );
        CHECK_EQUAL( i, def );
    }
}