1 // Copyright (C) 2007-2013 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 #include "MEDMEM_Mesh.hxx"
24 #include "MEDMEM_Group.hxx"
25 #include "MEDMEM_Meshing.hxx"
26 #include "MEDMEM_MedMeshDriver22.hxx"
28 using namespace MED_EN;
30 using namespace MEDMEM;
32 void addMedFacesGroup( MESHING& meshing, int nFaces, int *groupValue,string groupName)
35 int myNumberOfElements[1];
36 medGeometryElement mytypes[1];
39 faces.setName(groupName) ;
40 faces.setMesh(&meshing) ;
41 faces.setEntity(MED_EN::MED_FACE) ;
42 faces.setNumberOfGeometricType(1) ;
43 mytypes[0]=MED_EN::MED_POLYGON;
44 faces.setGeometricType(mytypes);
45 myNumberOfElements[0]=nFaces;
46 faces.setNumberOfElements(myNumberOfElements) ;
47 index[0]=1; index[1]=nFaces;
48 faces.setNumber(index, groupValue) ;
49 meshing.addGroup(faces) ;
52 const int NB_OF_HEXA=1;
53 const int NB_OF_CELL_AXIAL=1;
54 char NomBottom[]="Bottom";
58 MESHING *pmesh=new MESHING;
59 pmesh->setName( "RJHMESH" ) ;
61 double coord3DTab[36]={0.024155000000000003, 0.041837687256826218, -0.30499999999999999, 0.048310000000000013, -1.015761910347357e-17, -0.30499999999999999, 0.096620000000000011, -1.8329792978583059e-18,
62 -0.30499999999999999, 0.12077499999999999, 0.041837687256826225, -0.30499999999999999, 0.096620000000000011, 0.083675374513652451, -0.30499999999999999, 0.048310000000000013, 0.083675374513652465,
63 -0.30499999999999999, 0.024155000000000003, 0.041837687256826218, -0.2863, 0.048310000000000013, -1.015761910347357e-17, -0.2863, 0.096620000000000011, -1.8329792978583059e-18, -0.2863,
64 0.12077499999999999, 0.041837687256826225, -0.2863, 0.096620000000000011, 0.083675374513652451, -0.2863, 0.048310000000000013, 0.083675374513652465, -0.2863};
65 vector<double> coords3D(coord3DTab,coord3DTab+36);
66 unsigned nNodes=coords3D.size()/3;
67 pmesh->setCoordinates(3, nNodes, &(coords3D[0]), "CARTESIAN",MED_EN::MED_FULL_INTERLACE);
68 string coordname[3] = { "x", "y", "z" } ;
69 pmesh->setCoordinatesNames(coordname) ;
70 // Set coordinates units
71 string coordunit[3] = { "cm", "cm", "cm" } ;
72 pmesh->setCoordinatesUnits(coordunit) ;
74 medGeometryElement mytypes3D[1]={MED_POLYHEDRA};
75 medGeometryElement mytypes2D[2]={MED_QUAD4,MED_POLYGON};
76 pmesh->setNumberOfTypes(1, MED_EN::MED_CELL);
77 pmesh->setNumberOfTypes(2,MED_EN::MED_FACE);
78 pmesh->setTypes(mytypes3D,MED_EN::MED_CELL);
79 pmesh->setTypes(mytypes2D,MED_EN::MED_FACE);
81 const int nbPolyH = 1;
82 int polyHInd[2]={1, 9};
83 int polyHFaceConn[43]={8, 9, 10, 11, 12, 7, -1,
92 int conn2DClassicalTab[24]={2, 3, 9, 8, 3, 4, 10, 9, 4, 5, 11, 10, 5, 6, 12, 11, 6, 1, 7, 12, 1, 2, 8, 7};
93 vector<int> conn2DClassical(conn2DClassicalTab,conn2DClassicalTab+24);
95 int conn2DPolyIndTab[3]={1, 7, 13};
96 vector<int > conn2DPolyInd(conn2DPolyIndTab,conn2DPolyIndTab+3);
97 int conn2DPolyTab[12]={1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
99 int nbOfElemsFace[] = { conn2DClassical.size()/4, conn2DPolyInd.size()-1 };
100 pmesh->setNumberOfElements(&nbPolyH,MED_EN::MED_CELL);
101 pmesh->setNumberOfElements(nbOfElemsFace,MED_EN::MED_FACE);
104 pmesh->setConnectivity(MED_EN::MED_CELL, MED_EN::MED_POLYHEDRA, polyHFaceConn,polyHInd);
105 pmesh->setConnectivity(MED_EN::MED_FACE, MED_EN::MED_QUAD4,&conn2DClassical[0]);
106 pmesh->setConnectivity(MED_EN::MED_FACE, MED_EN::MED_POLYGON,&conn2DPolyInd[0],conn2DPolyIndTab);
108 addMedFacesGroup(*pmesh,1,&bottom,NomBottom);
109 int id = pmesh->addDriver(MED_DRIVER, "coeur.med", pmesh->getName());
115 MESH mesh;//(MED_DRIVER,,"RJHMESH");//"/home/geay/test_FGUI.med","FluidMesh_1");//../STEP3/coeur.med","RJHMESH");
116 mesh.setName("RJHMESH");
117 MED_MESH_RDONLY_DRIVER22 *drv=new MED_MESH_RDONLY_DRIVER22("coeur.med",&mesh);
118 drv->desactivateFacesComputation();
119 int newDrv=mesh.addDriver(*drv);
122 cout << " *** " << mesh.getNumberOfTypes(MED_EN::MED_FACE) << endl;
123 cout << "Expected 8 - having " << mesh.getNumberOfElements(MED_FACE,MED_ALL_ELEMENTS) << endl;
124 cout << "Expected 6 - having " << mesh.getNumberOfElements(MED_FACE,MED_QUAD4) << endl;
125 cout << "Expected 1 - having " << mesh.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) << endl;