From f4910d91d61fb1552e8711c8cee02f56e76a36ae Mon Sep 17 00:00:00 2001 From: ageay Date: Fri, 5 Apr 2013 10:00:32 +0000 Subject: [PATCH] MEDMEM suppression --- src/RENUMBER/renumbering.cxx | 22 +- src/RENUMBER/testRenumbering.py | 439 ++++++++++++++------------------ 2 files changed, 205 insertions(+), 256 deletions(-) diff --git a/src/RENUMBER/renumbering.cxx b/src/RENUMBER/renumbering.cxx index c9670f191..c073f1feb 100644 --- a/src/RENUMBER/renumbering.cxx +++ b/src/RENUMBER/renumbering.cxx @@ -101,10 +101,24 @@ int main(int argc, char** argv) // Fields cout << "Reordering fields and writing : " << flush; MEDFileFields *fs=fd->getFields(); - MEDCouplingAutoRefCountObjectPtr o2n=DataArrayInt::New(); - o2n->useArray(&perm[0],false,CPP_DEALLOC,perm.size(),1); - fs->renumberEntitiesLyingOnMesh(meshname.c_str(),code,code,o2n); - fs->write(filename_out.c_str(),0); + if(fs) + { + for(int i=0;igetNumberOfFields();i++) + { + MEDFileFieldMultiTS *fmts=fs->getFieldAtPos(i); + if(fmts->getMeshName()==meshname) + { + for(int j=0;jgetNumberOfTS();j++) + { + MEDFileField1TS *f1ts=fmts->getTimeStepAtPos(j); + DataArrayDouble *arr=f1ts->getUndergroundDataArray(); + arr->renumberInPlace(&perm[0]); + } + } + } + fs->write(filename_out.c_str(),0); + //fs->renumberEntitiesLyingOnMesh(meshname.c_str(),code,code,o2n); bugged + } t_field=clock(); cout << " : " << (t_field-t_family)/(double) CLOCKS_PER_SEC << "s" << endl << flush; return 0; diff --git a/src/RENUMBER/testRenumbering.py b/src/RENUMBER/testRenumbering.py index 782af1d3f..745720213 100755 --- a/src/RENUMBER/testRenumbering.py +++ b/src/RENUMBER/testRenumbering.py @@ -19,257 +19,192 @@ # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # -from libMEDMEM_Swig import * +from MEDLoader import * import os import sys - -## *************** -## *** TEST 2D *** -## *************** -srcdir = os.getenv("srcdir") -med_root = os.getenv("MED_ROOT_DIR") -if srcdir: - # make check is being performed - dir_renumber="./renumber" - dir_mesh = os.path.join( srcdir, "../../resources") -elif med_root: - # hope renumber has been already installed - dir_renumber=os.path.join( med_root, "bin/salome/renumber") - dir_mesh = os.path.join( med_root, "share/salome/resources/med") -else: - # initial version - dir_renumber="../../../MED_INSTALL/bin/salome/renumber" - dir_mesh="../../resources" - -filename="Test2D.med" -meshname="Mesh_1" - -print "TEST 2D Boost" -method="BOOST" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) - -field_ini=[2,3,12,13,14,15,4,5,6,7,8,9,16,17,0,1,10,11] -s = m.getSupportOnAll(MED_CELL) -f = FIELDDOUBLE(s,2) -id=f.addDriver(MED_DRIVER,dir_mesh+"/out_"+filename,"Test field") -f.read(id); -field=True -for i in range(9): - field=field&(f.getValueIJ(i+1,1)==field_ini[i*2]) - field=field&(f.getValueIJ(i+1,2)==field_ini[i*2+1]) -f.rmDriver(id) - -nbcell2dboost=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[2,6,13,11,11,13,14,12,6,5,15,13,12,14,10,4,13,15,16,14,5,1,7,15,14,16,9,10,15,7,8,16,16,8,3,9] -connectivite_index=[1,5,9,13,17,21,25,29,33,37] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_QUAD4) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn2dboost=(list(conn)==connectivite) # convert numpy.ndarray to list -conn_index2dboost=(list(conn_index)==connectivite_index) -Boost2D=conn2dboost and conn_index2dboost and (nbcell2dboost==9) and field -os.remove(dir_mesh+"/out_"+filename) - - -print "TEST 2D Metis" -method="METIS" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell2dmetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[12,14,10,4,2,6,13,11,11,13,14,12,16,8,3,9,5,1,7,15,15,7,8,16,14,16,9,10,6,5,15,13,13,15,16,14] -connectivite_index=[1,5,9,13,17,21,25,29,33,37] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_QUAD4) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn2dmetis=(list(conn)==connectivite) -conn_index2dmetis=(list(conn_index)==connectivite_index) -Metis2D=conn2dmetis and conn_index2dmetis and (nbcell2dmetis==9) -os.remove(dir_mesh+"/out_"+filename) - -## *** Avec polygone *** - -filename="Test2Dpoly.med" -meshname="Mesh_1" - - -print "TEST 2D Boost with polygons" -method="BOOST" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell2dpolyboost=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[2,5,9,10,11,10,9,12,5,6,8,9,4,11,12,16,12,9,8,13,6,1,7,8,16,12,13,15,13,8,7,14,15,13,14,3] -connectivite_index=[1,5,9,13,17,21,25,29,33,37] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYGON) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL) -conn2dpolyboost=(list(conn)==connectivite) -conn_index2dpolyboost=(list(conn_index)==connectivite_index) -PolyBoost2D=conn2dpolyboost and conn_index2dpolyboost and (nbcell2dpolyboost==9) -os.remove(dir_mesh+"/out_"+filename) - -print "TEST 2D Metis with polygons" -method="METIS" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell2dpolymetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[6,1,7,8,2,5,9,10,5,6,8,9,15,13,14,3,4,11,12,16,16,12,13,15,11,10,9,12,12,9,8,13,13,8,7,14] -connectivite_index=[1,5,9,13,17,21,25,29,33,37] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYGON) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL) -conn2dpolymetis=(list(conn)==connectivite) -conn_index2dpolymetis=(list(conn_index)==connectivite_index) -PolyMetis2D=conn2dpolymetis and conn_index2dpolymetis and (nbcell2dpolymetis==9) -os.remove(dir_mesh+"/out_"+filename) - - -## *************** -## *** TEST 3D *** -## *************** - - -filename="Test3D.med" -meshname="Mesh_1" - - -print "TEST 3D Boost" -method="BOOST" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell3dboost=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[23,13,5,18,27,22,14,26,17,6,13,23,25,16,22,27,27,22,14,26,24,15,7,20,9,23,18,1,21,27,26,10,25,16,22,27,19,8,15,24,2,17,23,9,12,25,27,21,21,27,26,10,11,24,20,3,12,25,27,21,4,19,24,11] -connectivite_index=[1,9,17,25,33,41,49,57,65] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_HEXA8) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn3dboost=(list(conn)==connectivite) -conn_index3dboost=(list(conn_index)==connectivite_index) -Boost3D=conn3dboost and conn_index3dboost and (nbcell3dboost==8) -os.remove(dir_mesh+"/out_"+filename) - - -print "TEST 3D Metis" -method="METIS" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell3dmetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[12,25,27,21,4,19,24,11,27,22,14,26,24,15,7,20,17,6,13,23,25,16,22,27,9,23,18,1,21,27,26,10,23,13,5,18,27,22,14,26,25,16,22,27,19,8,15,24,2,17,23,9,12,25,27,21,21,27,26,10,11,24,20,3] -connectivite_index=[1,9,17,25,33,41,49,57,65] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_HEXA8) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn3dmetis=(list(conn)==connectivite) -conn_index3dmetis=(list(conn_index)==connectivite_index) -Metis3D=conn3dmetis&conn_index3dmetis&(nbcell3dmetis==8) -os.remove(dir_mesh+"/out_"+filename) - - -## *** Avec polyedres *** - -## 23,13,5,18,27,26,14,22,23,27,22,13,13,22,14,5,5,14,26,18,18,26,27,23, -## 21,27,26,10,11,3,20,24,21,11,24,27,27,24,20,26,26,20,3,10,10,3,11,21, -## 12,25,27,21,4,11,24,19,12,4,19,25,25,19,24,27,27,24,11,21,21,11,4,12, -## 9,23,18,1,21,10,26,27,9,21,27,23,23,27,26,18,18,26,10,1,1,10,21,9, -## 2,17,23,9,12,21,27,25,2,12,25,17,17,25,27,23,23,27,21,9,9,21,12,2, -## 25,16,22,27,19,24,15,8,25,19,8,16,16,8,15,22,22,15,24,27,27,24,19,25, -## 17,6,13,23,25,27,22,16,17,25,16,6,6,16,22,13,13,22,27,23,23,27,25,17, -## 27,22,14,26,24,20,7,15,27,24,15,22,22,15,7,14,14,7,20,26,26,20,24,27, - - -filename="Test3Dpoly.med" -meshname="Mesh_1" - - -print "TEST 3D Boost with polyhedra" -method="BOOST" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell3dpolyboost=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[23,13,5,18,-1,27,26,14,22,-1,23,27,22,13,-1,13,22,14,5,-1,5,14,26,18,-1,18,26,27,23, - 17,6,13,23,-1,25,27,22,16,-1,17,25,16,6,-1,6,16,22,13,-1,13,22,27,23,-1,23,27,25,17, - 27,22,14,26,-1,24,20,7,15,-1,27,24,15,22,-1,22,15,7,14,-1,14,7,20,26,-1,26,20,24,27, - 9,23,18,1,-1,21,10,26,27,-1,9,21,27,23,-1,23,27,26,18,-1,18,26,10,1,-1,1,10,21,9, - 25,16,22,27,-1,19,24,15,8,-1,25,19,8,16,-1,16,8,15,22,-1,22,15,24,27,-1,27,24,19,25, - 2,17,23,9,-1,12,21,27,25,-1,2,12,25,17,-1,17,25,27,23,-1,23,27,21,9,-1,9,21,12,2, - 21,27,26,10,-1,11,3,20,24,-1,21,11,24,27,-1,27,24,20,26,-1,26,20,3,10,-1,10,3,11,21, - 12,25,27,21,-1,4,11,24,19,-1,12,4,19,25,-1,25,19,24,27,-1,27,24,11,21,-1,21,11,4,12] -connectivite_index=[1, 30, 59, 88, 117, 146, 175, 204, 233] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYHEDRA) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn3dpolyboost=(connectivite==list(conn)) -conn_index3dpolyboost=(connectivite_index==list(conn_index)) -PolyBoost3D=(conn3dpolyboost and conn_index3dpolyboost and (nbcell3dpolyboost==8)) -os.remove(dir_mesh+"/out_"+filename) - - -print "TEST 3D Metis with polyhedra" -method="METIS" -string_to_execute="'"+dir_renumber+" "+dir_mesh+"/"+filename+" "+meshname+" "+method+" "+dir_mesh+"/out_"+filename+"'" -eval("os.system("+string_to_execute+")") -m = MESH(MED_DRIVER,dir_mesh+"/out_"+filename,meshname) -nbcell3dpolymetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) -connectivite=[12,25,27,21,-1,4,11,24,19,-1,12,4,19,25,-1,25,19,24,27,-1,27,24,11,21,-1,21,11,4,12, - 27,22,14,26,-1,24,20,7,15,-1,27,24,15,22,-1,22,15,7,14,-1,14,7,20,26,-1,26,20,24,27, - 17,6,13,23,-1,25,27,22,16,-1,17,25,16,6,-1,6,16,22,13,-1,13,22,27,23,-1,23,27,25,17, - 9,23,18,1,-1,21,10,26,27,-1,9,21,27,23,-1,23,27,26,18,-1,18,26,10,1,-1,1,10,21,9, - 23,13,5,18,-1,27,26,14,22,-1,23,27,22,13,-1,13,22,14,5,-1,5,14,26,18,-1,18,26,27,23, - 25,16,22,27,-1,19,24,15,8,-1,25,19,8,16,-1,16,8,15,22,-1,22,15,24,27,-1,27,24,19,25, - 2,17,23,9,-1,12,21,27,25,-1,2,12,25,17,-1,17,25,27,23,-1,23,27,21,9,-1,9,21,12,2, - 21,27,26,10,-1,11,3,20,24,-1,21,11,24,27,-1,27,24,20,26,-1,26,20,3,10,-1,10,3,11,21] -connectivite_index=[1, 30, 59, 88, 117, 146, 175, 204, 233] -conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYHEDRA) -conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); -conn3dpolymetis=(list(conn)==connectivite) -conn_index3dpolymetis=(list(conn_index)==connectivite_index) -PolyMetis3D=(conn3dpolymetis and conn_index3dpolymetis and (nbcell3dpolymetis==8)) -os.remove(dir_mesh+"/out_"+filename) - - - - -print "" -if Boost2D: - print "Boost 2D ok" -else: - print "ERROR Boost 2D" -if Metis2D: - print "Metis 2D ok" -else: - print "ERROR Metis 2D" -if PolyBoost2D: - print "Poly Boost 2D ok" -else: - print "ERROR Poly Boost 2D" -if PolyMetis2D: - print "Poly Metis 2D ok" -else: - print "ERROR Poly Metis 2D" -if Boost3D: - print "Boost 3D ok" -else: - print "ERROR Boost 3D" -if Metis3D: - print "Metis 3D ok" -else: - print "ERROR Metis 3D" -if PolyBoost3D: - print "Poly Boost 3D ok" -else: - print "ERROR Poly Boost 3D" -if PolyMetis3D: - print "Poly Metis 3D ok" -else: - print "ERROR Poly Metis 3D" - - -print "" -if Boost2D&Metis2D&PolyBoost2D&PolyMetis2D&Boost3D&Metis3D&PolyBoost3D&PolyMetis3D: - print "Every mesh correctly renumbered" - sys.exit() -else: - print "Error" - sys.exit("Error in the renumbering test") - - +import unittest + +class RenumberingTest(unittest.TestCase): + def testBoost2D(self): + filename="Test2D.med" + meshname="Mesh_1" + method="BOOST" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + mm=MEDFileMesh.New(self.dir_mesh+"/out_"+filename,meshname) + m=mm.getMeshAtLevel(0) + ff=MEDFileField1TS(self.dir_mesh+"/out_"+filename,"Test field") + field_ini=DataArrayDouble([(2,3),(12,13),(14,15),(4,5),(6,7),(8,9),(16,17),(0,1),(10,11)]) + ff.getFieldOnMeshAtLevel(ON_CELLS,0,mm) + f=ff.getFieldOnMeshAtLevel(ON_CELLS,0,mm) + field=f.getArray().isEqual(field_ini,1e-15) + connectivite=[4,1,5,12,10,4,10,12,13,11,4,5,4,14,12,4,11,13,9,3,4,12,14,15,13,4,4,0,6,14,4,13,15,8,9,4,14,6,7,15,4,15,7,2,8] + connectivite_index=[0,5,10,15,20,25,30,35,40,45] + Boost2D=m.getNodalConnectivity().getValues()==connectivite and m.getNodalConnectivityIndex().getValues()==connectivite_index and field + self.assertTrue(Boost2D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def tessMetis2D(self):#not activated yet + filename="Test2D.med" + meshname="Mesh_1" + method="METIS" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + m = MESH(MED_DRIVER,self.dir_mesh+"/out_"+filename,meshname) + nbcell2dmetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) + connectivite=[12,14,10,4,2,6,13,11,11,13,14,12,16,8,3,9,5,1,7,15,15,7,8,16,14,16,9,10,6,5,15,13,13,15,16,14] + connectivite_index=[1,5,9,13,17,21,25,29,33,37] + conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_QUAD4) + conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); + conn2dmetis=(list(conn)==connectivite) + conn_index2dmetis=(list(conn_index)==connectivite_index) + Metis2D=conn2dmetis and conn_index2dmetis and (nbcell2dmetis==9) + self.assertTrue(Metis2D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def testBoost2DPolygon(self): + filename="Test2Dpoly.med" + meshname="Mesh_1" + method="BOOST" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + mm=MEDFileMesh.New(self.dir_mesh+"/out_"+filename,meshname) + m=mm.getMeshAtLevel(0) + nbcell2dpolyboost=m.getNumberOfCells() + connectivite=[5,1,4,8,9,5,10,9,8,11,5,4,5,7,8,5,3,10,11,15,5,11,8,7,12,5,5,0,6,7,5,15,11,12,14,5,12,7,6,13,5,14,12,13,2] + connectivite_index=[0,5,10,15,20,25,30,35,40,45] + conn=m.getNodalConnectivity().getValues() + conn_index=m.getNodalConnectivityIndex().getValues() + conn2dpolyboost=(list(conn)==connectivite) + conn_index2dpolyboost=(list(conn_index)==connectivite_index) + PolyBoost2D=conn2dpolyboost and conn_index2dpolyboost and (nbcell2dpolyboost==9) + self.assertTrue(PolyBoost2D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def tessMetis2DPolygon(self):#not activated yet + filename="Test2Dpoly.med" + meshname="Mesh_1" + method="METIS" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + m = MESH(MED_DRIVER,self.dir_mesh+"/out_"+filename,meshname) + nbcell2dpolymetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) + connectivite=[6,1,7,8,2,5,9,10,5,6,8,9,15,13,14,3,4,11,12,16,16,12,13,15,11,10,9,12,12,9,8,13,13,8,7,14] + connectivite_index=[1,5,9,13,17,21,25,29,33,37] + conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYGON) + conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL) + conn2dpolymetis=(list(conn)==connectivite) + conn_index2dpolymetis=(list(conn_index)==connectivite_index) + PolyMetis2D=conn2dpolymetis and conn_index2dpolymetis and (nbcell2dpolymetis==9) + self.assertTrue(PolyMetis2D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def testBoost3D(self): + filename="Test3D.med" + meshname="Mesh_1" + method="BOOST" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + mm=MEDFileMesh.New(self.dir_mesh+"/out_"+filename,meshname) + m=mm.getMeshAtLevel(0) + nbcell3dboost=m.getNumberOfCells() + connectivite=[18,22,12,4,17,26,21,13,25,18,16,5,12,22,24,15,21,26,18,26,21,13,25,23,14,6,19,18,8,22,17,0,20,26,25,9,18,24,15,21,26,18,7,14,23,18,1,16,22,8,11,24,26,20,18,20,26,25,9,10,23,19,2,18,11,24,26,20,3,18,23,10] + connectivite_index=[0,9,18,27,36,45,54,63,72] + conn=m.getNodalConnectivity().getValues() + conn_index=m.getNodalConnectivityIndex().getValues() + conn3dboost=(list(conn)==connectivite) + conn_index3dboost=(list(conn_index)==connectivite_index) + Boost3D=conn3dboost and conn_index3dboost and (nbcell3dboost==8) + self.assertTrue(Boost3D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def tessMetis3D(self):#not activated yet + filename="Test3D.med" + meshname="Mesh_1" + method="METIS" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + m = MESH(MED_DRIVER,self.dir_mesh+"/out_"+filename,meshname) + nbcell3dmetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) + connectivite=[12,25,27,21,4,19,24,11,27,22,14,26,24,15,7,20,17,6,13,23,25,16,22,27,9,23,18,1,21,27,26,10,23,13,5,18,27,22,14,26,25,16,22,27,19,8,15,24,2,17,23,9,12,25,27,21,21,27,26,10,11,24,20,3] + connectivite_index=[1,9,17,25,33,41,49,57,65] + conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_HEXA8) + conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); + conn3dmetis=(list(conn)==connectivite) + conn_index3dmetis=(list(conn_index)==connectivite_index) + Metis3D=conn3dmetis&conn_index3dmetis&(nbcell3dmetis==8) + self.assertTrue(Metis3D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def testBoost3DPoly(self): + filename="Test3Dpoly.med" + meshname="Mesh_1" + method="BOOST" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + mm=MEDFileMesh.New(self.dir_mesh+"/out_"+filename,meshname) + m=mm.getMeshAtLevel(0) + nbcell3dpolyboost=m.getNumberOfCells() + connectivite=[31,22,12,4,17,-1,26,25,13,21,-1,22,26,21,12,-1,12,21,13,4,-1,4,13,25,17,-1,17,25,26,22,31,16,5,12,22,-1,24,26,21,15,-1,16,24,15,5,-1,5,15,21,12,-1,12,21,26,22,-1,22,26,24,16,31,26,21,13,25,-1,23,19,6,14,-1,26,23,14,21,-1,21,14,6,13,-1,13,6,19,25,-1,25,19,23,26,31,8,22,17,0,-1,20,9,25,26,-1,8,20,26,22,-1,22,26,25,17,-1,17,25,9,0,-1,0,9,20,8,31,24,15,21,26,-1,18,23,14,7,-1,24,18,7,15,-1,15,7,14,21,-1,21,14,23,26,-1,26,23,18,24,31,1,16,22,8,-1,11,20,26,24,-1,1,11,24,16,-1,16,24,26,22,-1,22,26,20,8,-1,8,20,11,1,31,20,26,25,9,-1,10,2,19,23,-1,20,10,23,26,-1,26,23,19,25,-1,25,19,2,9,-1,9,2,10,20,31,11,24,26,20,-1,3,10,23,18,-1,11,3,18,24,-1,24,18,23,26,-1,26,23,10,20,-1,20,10,3,11] + connectivite_index=[0,30,60,90,120,150,180,210,240] + conn=m.getNodalConnectivity().getValues() + conn_index=m.getNodalConnectivityIndex().getValues() + conn3dpolyboost=(connectivite==list(conn)) + conn_index3dpolyboost=(connectivite_index==list(conn_index)) + PolyBoost3D=(conn3dpolyboost and conn_index3dpolyboost and (nbcell3dpolyboost==8)) + self.assertTrue(PolyBoost3D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def tessBoost3DPoly(self):#not activated yet + filename="Test3Dpoly.med" + meshname="Mesh_1" + method="METIS" + string_to_execute=self.dir_renumber+" "+self.dir_mesh+"/"+filename+" "+meshname+" "+method+" "+self.dir_mesh+"/out_"+filename + os.system(string_to_execute) + m = MESH(MED_DRIVER,self.dir_mesh+"/out_"+filename,meshname) + nbcell3dpolymetis=m.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS) + connectivite=[12,25,27,21,-1,4,11,24,19,-1,12,4,19,25,-1,25,19,24,27,-1,27,24,11,21,-1,21,11,4,12, + 27,22,14,26,-1,24,20,7,15,-1,27,24,15,22,-1,22,15,7,14,-1,14,7,20,26,-1,26,20,24,27, + 17,6,13,23,-1,25,27,22,16,-1,17,25,16,6,-1,6,16,22,13,-1,13,22,27,23,-1,23,27,25,17, + 9,23,18,1,-1,21,10,26,27,-1,9,21,27,23,-1,23,27,26,18,-1,18,26,10,1,-1,1,10,21,9, + 23,13,5,18,-1,27,26,14,22,-1,23,27,22,13,-1,13,22,14,5,-1,5,14,26,18,-1,18,26,27,23, + 25,16,22,27,-1,19,24,15,8,-1,25,19,8,16,-1,16,8,15,22,-1,22,15,24,27,-1,27,24,19,25, + 2,17,23,9,-1,12,21,27,25,-1,2,12,25,17,-1,17,25,27,23,-1,23,27,21,9,-1,9,21,12,2, + 21,27,26,10,-1,11,3,20,24,-1,21,11,24,27,-1,27,24,20,26,-1,26,20,3,10,-1,10,3,11,21] + connectivite_index=[1, 30, 59, 88, 117, 146, 175, 204, 233] + conn=m.getConnectivity(MED_NODAL,MED_CELL,MED_POLYHEDRA) + conn_index=m.getConnectivityIndex(MED_NODAL,MED_CELL); + conn3dpolymetis=(list(conn)==connectivite) + conn_index3dpolymetis=(list(conn_index)==connectivite_index) + PolyMetis3D=(conn3dpolymetis and conn_index3dpolymetis and (nbcell3dpolymetis==8)) + self.assertTrue(PolyMetis3D) + os.remove(self.dir_mesh+"/out_"+filename) + pass + + def setUp(self): + srcdir = os.getenv("srcdir") + med_root = os.getenv("MED_ROOT_DIR") + if srcdir: + # make check is being performed + self.dir_renumber="./renumber" + self.dir_mesh = os.path.join( srcdir, "../../resources") + elif med_root: + # hope renumber has been already installed + self.dir_renumber=os.path.join( med_root, "bin/salome/renumber") + self.dir_mesh = os.path.join( med_root, "share/salome/resources/med") + else: + # initial version + self.dir_renumber="../../../MED_INSTALL/bin/salome/renumber" + self.dir_mesh="../../resources" + pass + pass + pass + +unittest.main() -- 2.39.2