y = t * t + 1 + u * u z = u return (x, y, z) # surface grid a = D.surface(f, N=50) b = T.splitSize(a, 100) b = X.connectMatch(b, dim=2) t = C.newPyTree(['Surface', 2]) t[2][1][2] += b t = C.initVars(t, 'F', 1.) t = C.initVars(t, 'centers:G', 2.) t[2][1][2][0] = C.addBC2Zone(t[2][1][2][0], 'overlap', 'BCOverlap', 'imin') t = C.fillEmptyBCWith(t, 'wall', 'BCWall', dim=2) b = T.merge(t) t2 = C.newPyTree(['Surface', 2]) t2[2][1][2] += b test.testT(t2, 1) b = T.merge(t, alphaRef=45.) t2 = C.newPyTree(['Surface', 2]) t2[2][1][2] = b test.testT(t2, 2) # # cas d'angles vifs # a1 = G.cart((0, 0, 0), (1, 1, 1), (11, 11, 1)) a3 = G.cart((10, 0, 0), (1, 1, 1), (11, 11, 1)) a2 = T.rotate(a1, (0, 0, 0), (1, 0, 0), 90.) a2 = T.reorder(a2, (-1, 2, 3)) a2[0] = C.getZoneName(a1[0])
# - merge (pyTree) - import Converter.PyTree as C import Generator.PyTree as G import Transform.PyTree as T import Connector.PyTree as X import Geom.PyTree as D def f(t,u): x = t+u y = t*t+1+u*u z = u return (x,y,z) a = D.surface(f) b = T.splitSize(a, 100) b = X.connectMatch(b, dim=2) t = C.newPyTree(['Surface']); t[2][1][2] += b b = T.merge(t) t[2][1][2] = b C.convertPyTree2File(t, "out.cgns")
a = D.surface(f) t = C.newPyTree(['Base', 2]) b = [] i1 = 1 for i in range(10): i2 = (i + 1) * 10 j1 = 1 for j in range(10): j2 = (j + 1) * 10 b.append(T.subzone(a, (i1, j1, 1), (i2, j2, 1))) j1 = j2 i1 = i2 t[2][1][2] += b t = C.initVars(t, 'F', 1.) t = C.initVars(t, 'centers:G', 2.) res = T.merge(t, dir=2) t[2][1][2] = res test.testT(t, 1) res = T.merge(t, dir=1) t[2][1][2] = res test.testT(t, 2) # volume grids distrib = G.cart((0., 0., 0.), (0.1, 1, 1), (10, 1, 1)) a = G.addNormalLayers(a, distrib) i1 = 1 b = [] for i in range(10): i2 = (i + 1) * 10 j1 = 1 for j in range(10):