X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FwhichSide.py;h=c88a9097d03f6c13660f3064087e88ca89877022;hp=1cf3856028e8ad88fb99237f76da4ed19feedafb;hb=HEAD;hpb=385d4cede5f752d0eec26c306f3b5e14511e2a3d diff --git a/src/Tools/blocFissure/gmu/whichSide.py b/src/Tools/blocFissure/gmu/whichSide.py index 1cf385602..e17ad331f 100644 --- a/src/Tools/blocFissure/gmu/whichSide.py +++ b/src/Tools/blocFissure/gmu/whichSide.py @@ -1,5 +1,5 @@ # -*- coding: utf-8 -*- -# Copyright (C) 2014-2019 CEA/DEN, EDF R&D +# Copyright (C) 2014-2024 EDF # # This library is free software; you can redistribute it and/or # modify it under the terms of the GNU Lesser General Public @@ -17,13 +17,11 @@ # # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # +"""Calcul de la position d'une shape par rapport à une face (dessus, dessous, sur la surface même)""" import logging from .geomsmesh import geompy -# ----------------------------------------------------------------------------- -# --- calcul de la position d'une shape par rapport à une face (dessus, dessous, sur la surface même) - def whichSide(face, obj, tol = 1.e-3): """ Position d'un objet par rapport a une face (non nécessairement plane). @@ -40,6 +38,7 @@ def whichSide(face, obj, tol = 1.e-3): nbFaces = geompy.NbShapes(obj, geompy.ShapeType["FACE"]) # --- attention ! pour une seule face presente, renvoie 2 logging.debug(" nbFaces %s", nbFaces) vertices = geompy.ExtractShapes(obj, geompy.ShapeType["VERTEX"], False) + if nbEdges > 0 and nbFaces == 0: # --- edges if nbEdges <= 2: point = geompy.MakeVertexOnCurve(obj, 0.5) @@ -49,17 +48,20 @@ def whichSide(face, obj, tol = 1.e-3): for anEdge in edges: point = geompy.MakeVertexOnCurve(anEdge, 0.5) vertices.append(point) + elif nbFaces >0: # --- faces if nbFaces <=2: point = geompy.MakeVertexOnSurface(obj, 0.5, 0.5) vertices.append(point) - if nbFaces > 2: + else: faces = geompy.ExtractShapes(obj, geompy.ShapeType["FACE"], False) for aFace in faces: point = geompy.MakeVertexOnSurface(aFace, 0.5, 0.5) vertices.append(point) + else: # --- vertices vertices = [obj] + for vertex in vertices: distance = geompy.MinDistance(vertex, face) logging.debug(" distance %s", distance) @@ -74,5 +76,5 @@ def whichSide(face, obj, tol = 1.e-3): side = -1 break logging.debug(" side %s", side) - return side + return side