doc/pyplots/test_prism_quads.py

   1 # Copyright (C) 2009-2023  CEA, EDF
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License, or (at your option) any later version.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19
  20 #
  21
  22 # exemple repris de test_hexa1.cxx
  23
  24 import os
  25 import GEOM
  26 import SALOMEDS
  27 import hexablock
  28
  29 print("test prism...")
  30 doc = hexablock.addDocument()
  31
  32 orig = doc.addVertex(0, 0, 0)
  33 dx = doc.addVector(1, 0, 0)
  34 dy = doc.addVector(0, 1, 0)
  35 dz = doc.addVector(0, 0, 1)
  36
  37 dimx = 11
  38 dimy = 11
  39 dimz = 2
  40 grid = doc.makeCartesian(orig, dx, dy, dz, dimx, dimy, dimz)
  41
  42 mx = dimx/2
  43 my = dimy/2
  44 prems = grid.getQuadIJ(mx, my, dimz)
  45
  46 liste = [prems]
  47 liste.extend([grid.getQuadIJ(nx, my, dimz) for nx in range(dimx) if nx != mx])
  48 liste.extend([grid.getQuadIJ(mx, ny, dimz) for ny in range(dimy) if ny != my])
  49
  50 vec = doc.addVector(1, 1, 1)
  51 elts = doc.prismQuads(liste, vec, 5)
  52
  53 file_name = os.path.join(os.environ['TMP'], 'prisme.vtk')
  54 ####  elts.saveVtk(file_name)
  55
  56 print("...test prism OK")