1 #include "MEDMEM_Grid.hxx"
3 using namespace MEDMEM;
4 using namespace MED_EN;
6 static void usage(const char * test)
8 cerr << "Usage : " << test
9 << " SpaceDimension nbMaille" << endl << endl
10 << "-> Crée une grille cartesienne en dimension SpaceDimension avec nbMaille suivant chaque direction" << endl;
14 int main (int argc, char ** argv) {
16 /* process the arguments */
20 const int SpaceDimension=atoi(argv[1]);
21 const int nbMaille=atoi(argv[2]);
22 if(SpaceDimension>3 || SpaceDimension<1 || nbMaille<1)
25 // Creation des tableaux necessaires à la construction de GRID
26 std::vector<int> nMaille(SpaceDimension,nbMaille);
27 std::vector<double> Origine(SpaceDimension,0.0);
28 std::vector<double> pas(SpaceDimension,10.0);
29 std::vector<std::vector<double> > XYZ_Array(SpaceDimension);
30 for(int i=0;i!=SpaceDimension;++i)
32 XYZ_Array[i].resize(nMaille[i]+1); // nbre de noeuds = nbre de mailles +1
33 XYZ_Array[i][0]=Origine[i];
34 for(int j=1;j!=XYZ_Array[i].size();++j)
35 XYZ_Array[i][j]=XYZ_Array[i][j-1] + pas[j-1];
38 std::vector<std::string> coord_name(SpaceDimension,"X");
43 std::vector<std::string> coord_unit(SpaceDimension,"cm");
45 // creation du pointeur MESH à partir d'un GRID, test affichage
46 std::auto_ptr<MEDMEM::MESH> Mesh (new MEDMEM::GRID( XYZ_Array, coord_name, coord_unit, MED_CARTESIAN) );
47 const MEDMEM::CONNECTIVITY* conn = Mesh->getConnectivityptr();
48 std::cout << "Affichage du maillage : " << endl << *Mesh << endl;