2 .. _python_testMEDCoupling2D_solution:
9 import medcoupling as mc
15 print("1 ********************")
18 # Coordinates of central polygon
19 X = [1.,0.5,-0.5,-1.,-0.5,0.5]
20 Y = [0.,sqrt(3.)/2.,sqrt(3.)/2.,0.,-sqrt(3.)/2.,-sqrt(3.)/2.]
24 coordinates.append(0.)
25 coordinates.append(0.)
29 coordinates.append(X[i])
30 coordinates.append(Y[i])
32 # Coordinates of second couron
44 coordinates.append(xtmp+X[(i-1)%6])
45 coordinates.append(ytmp+Y[(i-1)%6])
46 coordinates.append(xtmp+X[i%6])
47 coordinates.append(ytmp+Y[i%6])
48 coordinates.append(xtmp+X[(i+1)%6])
49 coordinates.append(ytmp+Y[(i+1)%6])
51 print("2 ********************")
53 mesh=mc.MEDCoupplingUMesh.New()
54 mesh.setMeshDimension(2)
55 mesh.allocateCells(numberOfCells)
56 mesh.setName("MaFleur")
58 myCoords=mc.DataArrayDouble.New()
59 myCoords.setValues(coordinates,numberOfNodes,2)
60 mesh.setCoords(myCoords)
62 print("3 ********************")
63 # Connectivity of triangular meshing
66 connectivity.append(0)
67 connectivity.append(i%6+1)
68 connectivity.append((i+1)%6+1)
70 mesh.insertNextCell(mc.NORM_TRI3,3,connectivity[3*i:3*(i+1)])
73 print("4 ********************")
74 # Connectivity of hexagons
78 connectivity.append(start)
79 connectivity.append(start+2*(i+3))
80 connectivity.append(start+2*(i+3)+1)
81 connectivity.append(start+2*(i+3)+2)
83 connectivity.append(7)
85 connectivity.append(start+2*(i+3)+3)
86 connectivity.append((i+1)%6+1)
88 mesh.insertNextCell(mc.NORM_POLYGON,6,connectivity[6*i:6*(i+1)])
91 print("5 ********************")
92 mesh.checkConsistencyLight()
94 medFileName = "MEDCouppling_Fleur.med"
95 mc.WriteUMesh(medFileName,mesh,True)