1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142 | #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<--- Same expression on both sides of '-'. [+]Finding the same expression on both sides of an operator is suspicious and might indicate a cut and paste or logic error. Please examine this code carefully to determine if it is correct.
#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 );<--- Clarify calculation precedence for '%' and '?'. [+]Suspicious calculation. Please use parentheses to clarify the code. The code ''a%b?c:d'' should be written as either ''(a%b)?c:d'' or ''a%(b?c:d)''.
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 );
}
}
|