X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FwhichSideMulti.py;h=74cb172182de373c5b686b305f7679312c5aa788;hp=db013217cefb26f878d57bb469a6170ac9a33cc8;hb=93738e779d1a2bd6628e3f545a748ec2a930c695;hpb=385d4cede5f752d0eec26c306f3b5e14511e2a3d diff --git a/src/Tools/blocFissure/gmu/whichSideMulti.py b/src/Tools/blocFissure/gmu/whichSideMulti.py index db013217c..74cb17218 100644 --- a/src/Tools/blocFissure/gmu/whichSideMulti.py +++ b/src/Tools/blocFissure/gmu/whichSideMulti.py @@ -1,5 +1,5 @@ # -*- coding: utf-8 -*- -# Copyright (C) 2014-2019 CEA/DEN, EDF R&D +# Copyright (C) 2014-2021 EDF R&D # # This library is free software; you can redistribute it and/or # modify it under the terms of the GNU Lesser General Public @@ -17,12 +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) +from .geomsmesh import geompy def whichSideMulti(facerefs, ifil, obj, centre, tol = 1.e-3): """ @@ -39,6 +38,7 @@ def whichSideMulti(facerefs, ifil, obj, centre, tol = 1.e-3): nbEdges = geompy.NbShapes(obj, geompy.ShapeType["EDGE"]) # --- attention ! pour une seule edge presente, renvoie 2 logging.debug(" nbEdges %s", nbEdges) vertices = geompy.ExtractShapes(obj, geompy.ShapeType["VERTEX"], False) + if nbEdges > 0 : # --- edges if nbEdges <= 2: point = geompy.MakeVertexOnCurve(obj, 0.5) @@ -50,6 +50,7 @@ def whichSideMulti(facerefs, ifil, obj, centre, tol = 1.e-3): vertices.append(point) else: # --- vertices vertices = [obj] + nbOnRef = 0 nbOnOther = 0 for vertex in vertices: @@ -70,20 +71,21 @@ def whichSideMulti(facerefs, ifil, obj, centre, tol = 1.e-3): break else: nbMiss += 1 - pass # peut-être inside, tester les autres faces + # peut-être inside, tester les autres faces else: if i == ifil: nbOnRef +=1 # le point est sur la face de référence, on continue avec les autres points break else: nbOnOther += 1 # le point est sur une autre face, mais il peut aussi être sur la face de référence... - pass # on peut tester les autres faces + # on peut tester les autres faces if nbMiss == len(facerefs): side = 1 # inside if side != 0: break + if side == 0 and nbOnRef < len(vertices): - side = 1 # inside + side = 1 # inside logging.debug(" side %s", side) - return side + return side