-# Copyright (C) 2007-2012 CEA/DEN, EDF R&D
+# Copyright (C) 2007-2024 CEA, EDF
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
-# version 2.1 of the License.
+# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
#
# Author : Anthony Geay (CEA/DEN)
-from MEDCoupling import *
+from medcoupling import *
import math
import os
targetMesh.setName("MyMesh3D");
targetMesh.setDescription("build3DMesh");
targetMesh.allocateCells(12);
- for i in xrange(8):
+ for i in range(8):
targetMesh.insertNextCell(NORM_HEXA8,8,targetConn[8*i:8*(i+1)]);
pass
targetMesh.finishInsertingCells();
targetMesh=MEDCouplingUMesh.New();
targetMesh.setMeshDimension(0);
targetMesh.allocateCells(8);
- targetMesh.setName("Wonderfull 0D mesh");
+ targetMesh.setName("Wonderful 0D mesh");
targetMesh.setDescription("build0DMesh");
targetConn=[]
targetMesh.insertNextCell(NORM_POINT1,1,[0]);
myCoords.setInfoOnComponent(0,"X [m]");
myCoords.setInfoOnComponent(1,"YY [Pm]");
myCoords.setInfoOnComponent(2,"ZZZ [m]");
- targetMesh.checkCoherency();
+ targetMesh.checkConsistencyLight();
return targetMesh;
def buildM1DMesh(self):
- meshM1D=MEDCouplingUMesh.New("wonderfull -1 D mesh",-1);
+ meshM1D=MEDCouplingUMesh.New("wonderful -1 D mesh",-1);
meshM1D.setDescription("buildM1DMesh");
- meshM1D.checkCoherency();
+ meshM1D.checkConsistencyLight();
return meshM1D;
def buildExtrudedMesh(self):
m2D.changeSpaceDimension(3);
m1D=self.build1DMesh();
retu=m2D.buildExtrudedMesh(m1D,0);
- ret=MEDCouplingExtrudedMesh.New(retu,m2D,2);
+ ret=MEDCouplingMappedExtrudedMesh.New(retu,m2D,2);
ret.setName("ExtrudedTestForCorbaTest");
ret.setDescription("buildExtrudedMesh");
return ret;
a2.setInfoOnComponent(0,"SmthZ");
#
targetMesh.setCoordsAt(0,a1);
- targetMesh.setCoordsAt(2,a2);
+ targetMesh.setCoordsAt(1,a2);
#
#
- targetMesh.checkCoherency();
+ targetMesh.checkConsistencyLight();
#
return targetMesh;
+ def buildIMesh(self):
+ targetMesh=MEDCouplingIMesh.New();
+ targetMesh.setTime(2.3,4,5);
+ targetMesh.setTimeUnit("us");
+ targetMesh.setName("Example of IMesh");
+ targetMesh.setDescription("buildIMesh");
+ #
+ targetMesh.setSpaceDimension(3);
+ targetMesh.setNodeStruct([6,7,8]);
+ targetMesh.setOrigin([4.25,3.75,-6.125]);
+ targetMesh.setDXYZ([0.5,0.375,0.75]);
+ targetMesh.setAxisUnit("mm");
+ #
+ return targetMesh
+
+ def buildCLMesh(self):
+ targetMesh=MEDCouplingCurveLinearMesh();
+ targetMesh.setTime(2.3,4,5);
+ targetMesh.setTimeUnit("us");
+ targetMesh.setName("Example of Curve linear mesh");
+ targetMesh.setDescription("buildCLMesh");
+ a1=DataArrayDouble(3*20,1);
+ a1.iota(7.) ; a1.rearrange(3);
+ targetMesh.setCoords(a1);
+ targetMesh.setNodeGridStructure([4,5]);
+ #
+ targetMesh.checkConsistencyLight();
+ #
+ return targetMesh;
+
+ def build1SGTUMesh(self):
+ targetMesh=MEDCoupling1SGTUMesh("Mesh1SGT",NORM_QUAD4)
+ targetMesh.setTime(2.3,44,-55)
+ targetMesh.setTimeUnit("us")
+ targetMesh.setDescription("My Description of 1SGTU");
+ a1=DataArrayDouble([1.,1.,0.,2.,1.,0.,3.,1.,0.,1.,0.,0.,2.,0.,0.,0.,0.,0.,0.,1.,0.,3.,0.,0.,4.,0.,0.,4.,1.,0.],10,3)
+ a1.setInfoOnComponents(["X1 [m]","YY2 [km]","ZZZ3 [km]"])
+ targetMesh.setCoords(a1)
+ a2=DataArrayInt([6,0,3,5,3,0,1,4,1,2,7,4,8,7,2,9],4*4,1)
+ targetMesh.setNodalConnectivity(a2)
+ #
+ targetMesh.checkConsistencyLight();
+ #
+ return targetMesh;
+
+ def build1DGTUMesh(self):
+ targetMesh=MEDCoupling1DGTUMesh("Mesh1DGT",NORM_POLYGON);
+ targetMesh.setTime(2.3,55,-66)
+ targetMesh.setTimeUnit("us")
+ targetMesh.setDescription("My Description of 1DGTU");
+ a1=DataArrayDouble([1.,1.,0.,2.,1.,0.,3.,1.,0.,1.,0.,0.,2.,0.,0.,0.,0.,0.,0.,1.,0.,3.,0.,0.,4.,0.,0.,4.,1.,0.],10,3)
+ a1.setInfoOnComponents(["X1 [m]","YY2 [km]","ZZZ3 [km]"])
+ targetMesh.setCoords(a1)
+ a2=DataArrayInt([6,0,3,5,3,0,1,4,1,2,7,4,8,7,2],15,1)
+ a3=DataArrayInt([0,4,8,12,15],5,1)
+ targetMesh.setNodalConnectivity(a2,a3)
+ #
+ targetMesh.checkConsistencyLight();
+ #
+ return targetMesh;
+
def buildFieldScalarOn2DNT(self):
mesh=self.build2DMesh();
fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
tmp=mesh.getNumberOfCells()*6*[7.]
array.setValues(tmp,mesh.getNumberOfCells(),6);
fieldOnCells.setArray(array)
- fieldOnCells.checkCoherency();
+ fieldOnCells.checkConsistencyLight();
return fieldOnCells;
def buildFieldNodeScalarOn2DNT(self):
tmp=mesh.getNumberOfNodes()*5*[7.1234]
array.setValues(tmp,mesh.getNumberOfNodes(),5);
fieldOnNodes.setArray(array);
- fieldOnNodes.checkCoherency();
+ fieldOnNodes.checkConsistencyLight();
return fieldOnNodes;
def buildFieldScalarOn3DNT(self):
mesh=self.build3DMesh();
fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
- fieldOnCells.setNature(ConservativeVolumic);
+ fieldOnCells.setNature(IntensiveMaximum);
fieldOnCells.setName("toto");
fieldOnCells.setDescription("my wonderful 3D field toto2");
fieldOnCells.setMesh(mesh);
tmp=mesh.getNumberOfCells()*6*[7.]
array.setValues(tmp,mesh.getNumberOfCells(),6);
fieldOnCells.setArray(array);
- fieldOnCells.checkCoherency();
+ fieldOnCells.checkConsistencyLight();
return fieldOnCells;
def buildFieldScalarOn3DSurfWT(self):
array.setInfoOnComponents(["aaa","bbbb","ccccc"])
fieldOnCells.setArray(array);
fieldOnCells.setTime(6.7,1,4);
- fieldOnCells.checkCoherency();
+ fieldOnCells.checkConsistencyLight();
return fieldOnCells;
def buildFieldScalarOn3DSurfCOTI(self):
fieldOnCells.setArray(array);
fieldOnCells.setStartTime(6.7,1,4);
fieldOnCells.setEndTime(7.2,2,8);
- fieldOnCells.checkCoherency();
+ fieldOnCells.checkConsistencyLight();
return fieldOnCells;
def buildFieldScalarOn2DLT(self):
fieldOnCells.setEndArray(array)
fieldOnCells.setStartTime(6.7,25,26);
fieldOnCells.setEndTime(17.2,125,126);
- fieldOnCells.checkCoherency();
+ fieldOnCells.checkConsistencyLight();
return fieldOnCells;
def buildFieldGaussPt2DWT(self):
f.setGaussLocalizationOnType(NORM_QUAD4,_refCoo2,_gsCoo1,_wg1);
array=DataArrayDouble.New();
ptr=18*2*[None]
- for i in xrange(18*2):
+ for i in range(18*2):
ptr[i]=float(i+1);
pass
array.setValues(ptr,18,2);
f.setName("MyFirstFieldOnGaussPoint");
f.setTimeUnit("ms");
f.setDescription("mmmmmmmmmmmm");
- f.checkCoherency();
+ f.checkConsistencyLight();
return f;
def buildFieldGaussPtNE2DWT(self):
f.setDescription("MyDescriptionNE");
array=DataArrayDouble.New();
ptr=18*2*[None]
- for i in xrange(18*2):
+ for i in range(18*2):
ptr[i]=float(i+7)
array.setValues(ptr,18,2);
array.setInfoOnComponent(0,"Power [MW]");
array.setInfoOnComponent(1,"Density [kg/m^3]");
f.setArray(array);
#
- f.checkCoherency();
+ f.checkConsistencyLight();
return f;
def buildFieldVectorOnExtrudedWT(self):
array.setInfoOnComponent(1,"Density [kg/m^3]");
f.setArray(array);
#
- f.checkCoherency();
+ f.checkConsistencyLight();
return f
def buildFieldVectorOnCMeshWT(self):
array.setInfoOnComponent(1,"Density [kg/m^3]");
f.setArray(array);
#
- f.checkCoherency();
+ f.checkConsistencyLight();
return f
def buildFieldTemplateCellOn2D(self):
def buildFieldTemplateNodeOn2D(self):
f1=self.buildFieldNodeScalarOn2DNT();
f2=MEDCouplingFieldTemplate.New(f1);
- f2.setNature(ConservativeVolumic);
+ f2.setNature(IntensiveMaximum);
return f2
def buildFieldTemplateGaussPtOn2D(self):
f1=self.buildFieldGaussPt2DWT();
f2=MEDCouplingFieldTemplate.New(f1);
- f2.setNature(Integral);
+ f2.setNature(ExtensiveMaximum);
return f2
def buildFieldTemplateGaussNEOn2D(self):
f1=self.buildFieldGaussPtNE2DWT();
f2=MEDCouplingFieldTemplate.New(f1);
- f2.setNature(IntegralGlobConstraint);
+ f2.setNature(ExtensiveConservation);
return f2
def buildMultiFields1(self):
d0=DataArrayDouble.New();
d0.setValues(vals0,5,1);
d1=DataArrayDouble.New();
- d1.setValues(vals1,5,1);
+ d1.setValues(vals1[:5],5,1);
d1_1=DataArrayDouble.New();
d1_1.setValues(vals1_1,5,1);
d2=DataArrayDouble.New();
f4.setName("f4");
f4.setTime(2.7,25,26);
ret=MEDCouplingFieldOverTime.New([f0,f1,f2,f3,f4]);
- ret.checkCoherency();
+ ret.checkConsistencyLight();
return ret;
def buildFileNameForIOR(self):
- ret=os.getenv("TMP");
- ret+="/entryPointMEDCouplingCorba.ior";
+ tmpdir=os.getenv("TMP", "/tmp");
+ username="";
+ if os.getenv("USERNAME"):
+ username = os.getenv("USERNAME")+"_";
+ elif os.getenv("USER"):
+ username = os.getenv("USER")+"_";
+ ret=tmpdir+"/"+username+"entryPointMEDCouplingCorba.ior";
return ret;
pass
myCoords=DataArrayDouble.New()
myCoords.setValues(coords,nbOfNodes,3);
mesh.setCoords(myCoords);
- mesh.checkCoherency();
+ mesh.checkConsistencyLight();
myFalseConn=DataArrayInt.New()
myFalseConn.setValues(tab4,6,4)
return mesh
nbOfCells=mesh.getNumberOfCells()
field=MEDCouplingFieldDouble.New(ON_CELLS)
field.setMesh(mesh)
- field.setNature(Integral)
+ field.setNature(ExtensiveMaximum)
myCoords=DataArrayDouble.New()
sampleTab=[]
for i in range(nbOfCells*9):