tests/test_enforced_vertex.py

   1 # -*- coding: utf-8 -*-
   2
   3 import salome
   4 import math
   5
   6 import GEOM
   7 from salome.geom import geomBuilder
   8 geompy = geomBuilder.New(salome.myStudy)
   9
  10 r = 10
  11 dist = 10
  12 dist_coin = 10.1
  13
  14 p1 = geompy.MakeVertex(0., 0., 0.)
  15 p2 = geompy.MakeVertex(100., 100., 100.)
  16 box = geompy.MakeBoxTwoPnt(p1, p2)
  17 geompy.addToStudy(box, "box")
  18
  19 p3 = geompy.MakeVertex(dist_coin, 0, dist_coin)
  20 geompy.addToStudy(p3, "p3")
  21
  22 left = geompy.GetFaceNearPoint(box, p3)
  23 geompy.addToStudyInFather(box, left, "left")
  24
  25 # Mesh
  26 # ====
  27
  28 import SMESH
  29 from salome.smesh import smeshBuilder
  30 smesh = smeshBuilder.New(salome.myStudy)
  31
  32 Mesh = smesh.Mesh(box, "Mesh")
  33
  34
  35 algo2d = Mesh.Triangle(algo=smeshBuilder.MG_CADSurf)
  36 algo2d.SetGeometricMesh( 1 )
  37 algo2d.SetAngleMesh( 4 )
  38 algo2d.SetPhySize( 8 )
  39
  40 algo2d.SetEnforcedVertex(left, dist_coin, 0, dist_coin)
  41
  42 Mesh.Compute()
  43
  44 id_node = Mesh.FindNodeClosestTo(dist_coin, 0, dist_coin)
  45
  46 x, y, z = Mesh.GetNodeXYZ(id_node)
  47
  48 assert("%.2f, %.2f, %.2f"%(x, y, z) == "%.2f, %.2f, %.2f"%(dist_coin, 0, dist_coin))
  49
  50
  51 salome.sg.updateObjBrowser(True)
  52