X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FrestreintFaceFissure.py;h=0f90e4a04b439a21408928e2e3a32b6d4ee09f92;hb=d439ade3e9fbb1e08c4f3d2277479898a1d59a3f;hp=b512bbae2804dc8716e70c0cef65cddfadcf822a;hpb=09eeca0267acf6f14581efdd342e8b7a87d09cf0;p=modules%2Fsmesh.git diff --git a/src/Tools/blocFissure/gmu/restreintFaceFissure.py b/src/Tools/blocFissure/gmu/restreintFaceFissure.py index b512bbae2..0f90e4a04 100644 --- a/src/Tools/blocFissure/gmu/restreintFaceFissure.py +++ b/src/Tools/blocFissure/gmu/restreintFaceFissure.py @@ -1,5 +1,5 @@ # -*- coding: utf-8 -*- -# Copyright (C) 2014-2020 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,14 +17,18 @@ # # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # +"""Restriction de la face de fissure au domaine solide""" import logging + +import traceback + from .geomsmesh import geompy from .geomsmesh import geomPublish -from .geomsmesh import geomPublishInFather + from . import initLog + from .sortFaces import sortFaces -import traceback from .fissError import fissError def restreintFaceFissure(shapeDefaut, facesDefaut, pointInterne): @@ -33,13 +37,16 @@ def restreintFaceFissure(shapeDefaut, facesDefaut, pointInterne): partition face fissure étendue par fillings """ logging.info('start') + partShapeDefaut = geompy.MakePartition([shapeDefaut], facesDefaut, [], [], geompy.ShapeType["FACE"], 0, [], 0) geomPublish(initLog.debug, partShapeDefaut, 'partShapeDefaut') facesPartShapeDefaut = geompy.ExtractShapes(partShapeDefaut, geompy.ShapeType["FACE"], False) + if pointInterne is not None: distfaces = [(geompy.MinDistance(face,pointInterne), i, face) for i, face in enumerate(facesPartShapeDefaut)] distfaces.sort() - logging.debug("selection de la face la plus proche du point interne, distance={}".format(distfaces[0][0])) + texte = "selection de la face la plus proche du point interne, distance={}".format(distfaces[0][0]) + logging.debug(texte) facesPortFissure = distfaces[0][2] else: try: @@ -52,8 +59,10 @@ def restreintFaceFissure(shapeDefaut, facesDefaut, pointInterne): texte += "
  • le prémaillage de la face de fissure est trop grossier, les mailles à enlever dans le maillage sain " texte += "n'ont pas toutes été détectées.
  • " raise fissError(traceback.extract_stack(),texte) - logging.debug("surfaces faces fissure étendue, min {}, max {}".format(minSurf, maxSurf)) + texte = "surfaces faces fissure étendue, min {}, max {}".format(minSurf, maxSurf) + logging.debug(texte) facesPortFissure = facesPartShapeDefautSorted[-1] geomPublish(initLog.debug, facesPortFissure, "facesPortFissure") + return facesPortFissure