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
|
// -------------------------------------------------------------------
// MAdLib - Copyright (C) 2008-2009 Universite catholique de Louvain
//
// See the Copyright.txt and License.txt files for license information.
// You should have received a copy of these files along with MAdLib.
// If not, see <http://www.madlib.be/license/>
//
// Please report all bugs and problems to <contrib@madlib.be>
//
// Authors: Gaetan Compere, Jean-Francois Remacle
// -------------------------------------------------------------------
#include "MAdLib.h"
#include <iostream>
using std::cout;
#include <sys/resource.h>
#include <string>
using std::string;
#include <stdlib.h>
using namespace MAd;
#ifdef PARALLEL
#include "mpi.h"
#endif
// ----------------------------------------------------------------------
int main(int argc, char* argv[]) {
#ifdef PARALLEL
MPI_Init(&argc, &argv);
#endif
// Check input
// ------------
if ( argc != 3 && argc != 2 ) {
printf("Error: usage: \'executable mshFile [geoFile]\'\n");
exit(0);
}
string meshFile = argv[1];
// Build tools
// ------------
printf ("Checking mesh %s...\n\n",meshFile.c_str());
// --- Reading model ---
pGModel model = 0;
GM_create(&model,"theModel");
if ( argc == 3 ) {
string geoFile = argv[2];
GM_read(model, geoFile.c_str());
}
else {
GM_readFromMSH(model, meshFile.c_str());
}
pMesh mesh = M_new(model);
M_load(mesh,meshFile.c_str());
MeshStatus status = VALID;
if ( !checkMesh(mesh, CHECK_ALL, 1, std::cout, &status) )
{
printf("\nThe mesh is invalid, status = %d\n\n",status);
}
else
{
printf("\nThe mesh is valid, status = %d\n\n",status);
}
#ifdef PARALLEL
MPI_Finalize();
#endif
}
// ----------------------------------------------------------------------
|