Salome HOME
Merge from V6_3_BR branch (Windows porting) 27/10/2011
[modules/smesh.git] / src / SMESH_SWIG / smeshDC.py
index 5dba9897fff30e2b6cb36b9e372f0d15814216c5..29b6c39fdd41dac1032feb1efd16e8267b54e51d 100644 (file)
@@ -509,6 +509,26 @@ def AssureGeomPublished(mesh, geom, name=''):
         mesh.geompyD.addToStudyInFather( mesh.geom, geom, name )
     return
 
+## Return the first vertex of a geomertical edge by ignoring orienation
+def FirstVertexOnCurve(edge):
+    from geompy import SubShapeAll, ShapeType, KindOfShape, PointCoordinates
+    vv = SubShapeAll( edge, ShapeType["VERTEX"])
+    if not vv:
+        raise TypeError, "Given object has no vertices"
+    if len( vv ) == 1: return vv[0]
+    info = KindOfShape(edge)
+    xyz = info[1:4] # coords of the first vertex
+    xyz1  = PointCoordinates( vv[0] )
+    xyz2  = PointCoordinates( vv[1] )
+    dist1, dist2 = 0,0
+    for i in range(3):
+        dist1 += abs( xyz[i] - xyz1[i] )
+        dist2 += abs( xyz[i] - xyz2[i] )
+    if dist1 < dist2:
+        return vv[0]
+    else:
+        return vv[1]
+
 # end of l1_auxiliary
 ## @}
 
@@ -677,6 +697,17 @@ class smeshDC(SMESH._objref_SMESH_Gen):
             aMeshes.append(aMesh)
         return aMeshes, aStatus
 
+    ## Creates a Mesh object(s) importing data from the given SAUV file
+    #  @return a list of Mesh class instances
+    #  @ingroup l2_impexp
+    def CreateMeshesFromSAUV( self,theFileName ):
+        aSmeshMeshes, aStatus = SMESH._objref_SMESH_Gen.CreateMeshesFromSAUV(self,theFileName)
+        aMeshes = []
+        for iMesh in range(len(aSmeshMeshes)) :
+            aMesh = Mesh(self, self.geompyD, aSmeshMeshes[iMesh])
+            aMeshes.append(aMesh)
+        return aMeshes, aStatus
+
     ## Creates a Mesh object importing data from the given STL file
     #  @return an instance of Mesh class
     #  @ingroup l2_impexp
@@ -796,6 +827,8 @@ class smeshDC(SMESH._objref_SMESH_Gen):
                      UnaryOp=FT_Undefined,
                      BinaryOp=FT_Undefined,
                      Tolerance=1e-07):
+        if not CritType in SMESH.FunctorType._items:
+            raise TypeError, "CritType should be of SMESH.FunctorType"
         aCriterion = self.GetEmptyCriterion()
         aCriterion.TypeOfElement = elementType
         aCriterion.Type = self.EnumToLong(CritType)
@@ -851,6 +884,7 @@ class smeshDC(SMESH._objref_SMESH_Gen):
             # Checks the treshold
             try:
                 aCriterion.Threshold = self.EnumToLong(aTreshold)
+                assert( aTreshold in SMESH.GeometryType._items )
             except:
                 if isinstance(aTreshold, int):
                     aCriterion.Threshold = aTreshold
@@ -1390,6 +1424,15 @@ class Mesh:
     def Projection1D(self, geom=0):
         return Mesh_Projection1D(self,  geom)
 
+    ## Creates a projection 1D-2D algorithm for faces.
+    #  If the optional \a geom parameter is not set, this algorithm is global.
+    #  Otherwise, this algorithm defines a submesh based on \a geom subshape.
+    #  @param geom If defined, the subshape to be meshed
+    #  @return an instance of Mesh_Projection2D algorithm
+    #  @ingroup l3_algos_proj
+    def Projection1D2D(self, geom=0):
+        return Mesh_Projection2D(self,  geom, "Projection_1D2D")
+
     ## Creates a projection 2D algorithm for faces.
     #  If the optional \a geom parameter is not set, this algorithm is global.
     #  Otherwise, this algorithm defines a submesh based on \a geom subshape.
@@ -1397,7 +1440,7 @@ class Mesh:
     #  @return an instance of Mesh_Projection2D algorithm
     #  @ingroup l3_algos_proj
     def Projection2D(self, geom=0):
-        return Mesh_Projection2D(self,  geom)
+        return Mesh_Projection2D(self,  geom, "Projection_2D")
 
     ## Creates a projection 3D algorithm for solids.
     #  If the optional \a geom parameter is not set, this algorithm is global.
@@ -1715,6 +1758,15 @@ class Mesh:
         else:
             self.mesh.ExportToMEDX(f, auto_groups, version, overwrite)
 
+    ## Exports the mesh in a file in SAUV format
+    #  @param f is the file name
+    #  @param auto_groups boolean parameter for creating/not creating
+    #  the groups Group_On_All_Nodes, Group_On_All_Faces, ... ;
+    #  the typical use is auto_groups=false.
+    #  @ingroup l2_impexp
+    def ExportSAUV(self, f, auto_groups=0):
+        self.mesh.ExportSAUV(f, auto_groups)
+
     ## Exports the mesh in a file in DAT format
     #  @param f the file name
     #  @param meshPart a part of mesh (group, sub-mesh) to export instead of the mesh
@@ -1760,6 +1812,8 @@ class Mesh:
     def ExportCGNS(self, f, overwrite=1, meshPart=None):
         if isinstance( meshPart, list ):
             meshPart = self.GetIDSource( meshPart, SMESH.ALL )
+        if isinstance( meshPart, Mesh ):
+            meshPart = meshPart.mesh
         elif not meshPart:
             meshPart = self.mesh
         self.mesh.ExportCGNS(meshPart, f, overwrite)
@@ -4566,6 +4620,42 @@ class Mesh_Algorithm:
         hyp.SetIgnoreFaces(ignoreFaces)
         return hyp
 
+    ## Transform a list of ether edges or tuples (edge 1st_vertex_of_edge)
+    #  into a list acceptable to SetReversedEdges() of some 1D hypotheses
+    #  @ingroupl3_hypos_1dhyps
+    def ReversedEdgeIndices(self, reverseList):
+        resList = []
+        geompy = self.mesh.geompyD
+        for i in reverseList:
+            if isinstance( i, int ):
+                s = geompy.SubShapes(self.mesh.geom, [i])[0]
+                if s.GetShapeType() != geompyDC.GEOM.EDGE:
+                    raise TypeError, "Not EDGE index given"
+                resList.append( i )
+            elif isinstance( i, geompyDC.GEOM._objref_GEOM_Object ):
+                if i.GetShapeType() != geompyDC.GEOM.EDGE:
+                    raise TypeError, "Not an EDGE given"
+                resList.append( geompy.GetSubShapeID(self.mesh.geom, i ))
+            elif len( i ) > 1:
+                e = i[0]
+                v = i[1]
+                if not isinstance( e, geompyDC.GEOM._objref_GEOM_Object ) or \
+                   not isinstance( v, geompyDC.GEOM._objref_GEOM_Object ):
+                    raise TypeError, "A list item must be a tuple (edge 1st_vertex_of_edge)"
+                if v.GetShapeType() == geompyDC.GEOM.EDGE and \
+                   e.GetShapeType() == geompyDC.GEOM.VERTEX:
+                    v,e = e,v
+                if e.GetShapeType() != geompyDC.GEOM.EDGE or \
+                   v.GetShapeType() != geompyDC.GEOM.VERTEX:
+                    raise TypeError, "A list item must be a tuple (edge 1st_vertex_of_edge)"
+                vFirst = FirstVertexOnCurve( e )
+                tol    = geompy.Tolerance( vFirst )[-1]
+                if geompy.MinDistance( v, vFirst ) > 1.5*tol:
+                    resList.append( geompy.GetSubShapeID(self.mesh.geom, e ))
+            else:
+                raise TypeError, "Item must be either an edge or tuple (edge 1st_vertex_of_edge)"
+        return resList
+
 # Public class: Mesh_Segment
 # --------------------------
 
@@ -4638,7 +4728,8 @@ class Mesh_Segment(Mesh_Algorithm):
     ## Defines "NumberOfSegments" hypothesis to cut an edge in a fixed number of segments
     #  @param n for the number of segments that cut an edge
     #  @param s for the scale factor (optional)
-    #  @param reversedEdges is a list of edges to mesh using reversed orientation
+    #  @param reversedEdges is a list of edges to mesh using reversed orientation.
+    #                       A list item can also be a tuple (edge 1st_vertex_of_edge)
     #  @param UseExisting if ==true - searches for an existing hypothesis created with
     #                     the same parameters, else (default) - create a new one
     #  @return an instance of StdMeshers_NumberOfSegments hypothesis
@@ -4647,20 +4738,19 @@ class Mesh_Segment(Mesh_Algorithm):
         if not isinstance(reversedEdges,list): #old version script, before adding reversedEdges
             reversedEdges, UseExisting = [], reversedEdges
         entry = self.MainShapeEntry()
-        if reversedEdges and isinstance(reversedEdges[0],geompyDC.GEOM._objref_GEOM_Object):
-            reversedEdges = [ self.mesh.geompyD.GetSubShapeID(self.mesh.geom, e) for e in reversedEdges ]
+        reversedEdgeInd = self.ReversedEdgeIndices(reversedEdges)
         if s == []:
-            hyp = self.Hypothesis("NumberOfSegments", [n, reversedEdges, entry],
+            hyp = self.Hypothesis("NumberOfSegments", [n, reversedEdgeInd, entry],
                                   UseExisting=UseExisting,
                                   CompareMethod=self.CompareNumberOfSegments)
         else:
-            hyp = self.Hypothesis("NumberOfSegments", [n,s, reversedEdges, entry],
+            hyp = self.Hypothesis("NumberOfSegments", [n,s, reversedEdgeInd, entry],
                                   UseExisting=UseExisting,
                                   CompareMethod=self.CompareNumberOfSegments)
             hyp.SetDistrType( 1 )
             hyp.SetScaleFactor(s)
         hyp.SetNumberOfSegments(n)
-        hyp.SetReversedEdges( reversedEdges )
+        hyp.SetReversedEdges( reversedEdgeInd )
         hyp.SetObjectEntry( entry )
         return hyp
 
@@ -4683,7 +4773,8 @@ class Mesh_Segment(Mesh_Algorithm):
     ## Defines "Arithmetic1D" hypothesis to cut an edge in several segments with increasing arithmetic length
     #  @param start defines the length of the first segment
     #  @param end   defines the length of the last  segment
-    #  @param reversedEdges is a list of edges to mesh using reversed orientation
+    #  @param reversedEdges is a list of edges to mesh using reversed orientation.
+    #                       A list item can also be a tuple (edge 1st_vertex_of_edge)
     #  @param UseExisting if ==true - searches for an existing hypothesis created with
     #                     the same parameters, else (default) - creates a new one
     #  @return an instance of StdMeshers_Arithmetic1D hypothesis
@@ -4691,15 +4782,14 @@ class Mesh_Segment(Mesh_Algorithm):
     def Arithmetic1D(self, start, end, reversedEdges=[], UseExisting=0):
         if not isinstance(reversedEdges,list): #old version script, before adding reversedEdges
             reversedEdges, UseExisting = [], reversedEdges
-        if reversedEdges and isinstance(reversedEdges[0],geompyDC.GEOM._objref_GEOM_Object):
-            reversedEdges = [ self.mesh.geompyD.GetSubShapeID(self.mesh.geom, e) for e in reversedEdges ]
+        reversedEdgeInd = self.ReversedEdgeIndices(reversedEdges)
         entry = self.MainShapeEntry()
-        hyp = self.Hypothesis("Arithmetic1D", [start, end, reversedEdges, entry],
+        hyp = self.Hypothesis("Arithmetic1D", [start, end, reversedEdgeInd, entry],
                               UseExisting=UseExisting,
                               CompareMethod=self.CompareArithmetic1D)
         hyp.SetStartLength(start)
         hyp.SetEndLength(end)
-        hyp.SetReversedEdges( reversedEdges )
+        hyp.SetReversedEdges( reversedEdgeInd )
         hyp.SetObjectEntry( entry )
         return hyp
 
@@ -4721,7 +4811,8 @@ class Mesh_Segment(Mesh_Algorithm):
     # values are equals 1
     #  @param points defines the list of parameters on curve
     #  @param nbSegs defines the list of numbers of segments
-    #  @param reversedEdges is a list of edges to mesh using reversed orientation
+    #  @param reversedEdges is a list of edges to mesh using reversed orientation.
+    #                       A list item can also be a tuple (edge 1st_vertex_of_edge)
     #  @param UseExisting if ==true - searches for an existing hypothesis created with
     #                     the same parameters, else (default) - creates a new one
     #  @return an instance of StdMeshers_Arithmetic1D hypothesis
@@ -4729,15 +4820,14 @@ class Mesh_Segment(Mesh_Algorithm):
     def FixedPoints1D(self, points, nbSegs=[1], reversedEdges=[], UseExisting=0):
         if not isinstance(reversedEdges,list): #old version script, before adding reversedEdges
             reversedEdges, UseExisting = [], reversedEdges
-        if reversedEdges and isinstance(reversedEdges[0],geompyDC.GEOM._objref_GEOM_Object):
-            reversedEdges = [ self.mesh.geompyD.GetSubShapeID(self.mesh.geom, e) for e in reversedEdges ]
+        reversedEdgeInd = self.ReversedEdgeIndices(reversedEdges)
         entry = self.MainShapeEntry()
-        hyp = self.Hypothesis("FixedPoints1D", [points, nbSegs, reversedEdges, entry],
+        hyp = self.Hypothesis("FixedPoints1D", [points, nbSegs, reversedEdgeInd, entry],
                               UseExisting=UseExisting,
                               CompareMethod=self.CompareFixedPoints1D)
         hyp.SetPoints(points)
         hyp.SetNbSegments(nbSegs)
-        hyp.SetReversedEdges(reversedEdges)
+        hyp.SetReversedEdges(reversedEdgeInd)
         hyp.SetObjectEntry(entry)
         return hyp
 
@@ -4757,7 +4847,8 @@ class Mesh_Segment(Mesh_Algorithm):
     ## Defines "StartEndLength" hypothesis to cut an edge in several segments with increasing geometric length
     #  @param start defines the length of the first segment
     #  @param end   defines the length of the last  segment
-    #  @param reversedEdges is a list of edges to mesh using reversed orientation
+    #  @param reversedEdges is a list of edges to mesh using reversed orientation.
+    #                       A list item can also be a tuple (edge 1st_vertex_of_edge)
     #  @param UseExisting if ==true - searches for an existing hypothesis created with
     #                     the same parameters, else (default) - creates a new one
     #  @return an instance of StdMeshers_StartEndLength hypothesis
@@ -4765,15 +4856,14 @@ class Mesh_Segment(Mesh_Algorithm):
     def StartEndLength(self, start, end, reversedEdges=[], UseExisting=0):
         if not isinstance(reversedEdges,list): #old version script, before adding reversedEdges
             reversedEdges, UseExisting = [], reversedEdges
-        if reversedEdges and isinstance(reversedEdges[0],geompyDC.GEOM._objref_GEOM_Object):
-            reversedEdges = [ self.mesh.geompyD.GetSubShapeID(self.mesh.geom, e) for e in reversedEdges ]
+        reversedEdgeInd = self.ReversedEdgeIndices(reversedEdges)
         entry = self.MainShapeEntry()
-        hyp = self.Hypothesis("StartEndLength", [start, end, reversedEdges, entry],
+        hyp = self.Hypothesis("StartEndLength", [start, end, reversedEdgeInd, entry],
                               UseExisting=UseExisting,
                               CompareMethod=self.CompareStartEndLength)
         hyp.SetStartLength(start)
         hyp.SetEndLength(end)
-        hyp.SetReversedEdges( reversedEdges )
+        hyp.SetReversedEdges( reversedEdgeInd )
         hyp.SetObjectEntry( entry )
         return hyp
 
@@ -5096,35 +5186,34 @@ class Mesh_Triangle(Mesh_Algorithm):
             #  Parameter of BLSURF algo
             self.params.SetVerbosity(level)
 
-    ## To optimize the CAD (merges edges and removes nano edges).
+    ## To optimize merges edges.
     #  @ingroup l3_hypos_blsurf
-    def SetPreCADOptimCAD(self, toOptimizeCAD=False):
+    def SetPreCADMergeEdges(self, toMergeEdges=False):
         if self.Parameters():
             #  Parameter of BLSURF algo
-            self.params.SetPreCADOptimCAD(toOptimizeCAD)
+            self.params.SetPreCADMergeEdges(toMergeEdges)
 
-    ## To compute topology from scratch
+    ## To remove nano edges.
     #  @ingroup l3_hypos_blsurf
-    def SetPreCADDiscardInput(self, toDiscardInput=False):
+    def SetPreCADRemoveNanoEdges(self, toRemoveNanoEdges=False):
         if self.Parameters():
             #  Parameter of BLSURF algo
-            self.params.SetPreCADDiscardInput(toDiscardInput)
+            self.params.SetPreCADRemoveNanoEdges(toRemoveNanoEdges)
 
-    ## To help PreCAD treat some very dirty cases.
-    #  If the treated object is manifold.
+    ## To compute topology from scratch
     #  @ingroup l3_hypos_blsurf
-    def SetPreCADManifoldGeom(self, manifold=False):
+    def SetPreCADDiscardInput(self, toDiscardInput=False):
         if self.Parameters():
             #  Parameter of BLSURF algo
-            self.params.SetPreCADManifoldGeom(manifold)
+            self.params.SetPreCADDiscardInput(toDiscardInput)
 
-    ## To help PreCAD treat some very dirty cases.
-    #  If the object is also closed (imagine a shell).
+    ## Sets the length below which an edge is considered as nano 
+    #  for the topology processing.
     #  @ingroup l3_hypos_blsurf
-    def SetPreCADClosedGeom(self, closed=False):
+    def SetPreCADEpsNano(self, epsNano):
         if self.Parameters():
             #  Parameter of BLSURF algo
-            self.params.SetPreCADClosedGeom(closed)
+            self.params.SetPreCADEpsNano(epsNano)
 
     ## Sets advanced option value.
     #  @ingroup l3_hypos_blsurf
@@ -5133,6 +5222,16 @@ class Mesh_Triangle(Mesh_Algorithm):
             #  Parameter of BLSURF algo
             self.params.SetOptionValue(optionName,level)
 
+    ## Sets advanced PreCAD option value.
+    #  Keyword arguments:
+    #  optionName: name of the option
+    #  optionValue: value of the option
+    #  @ingroup l3_hypos_blsurf
+    def SetPreCADOptionValue(self, optionName, optionValue):
+        if self.Parameters():
+            #  Parameter of BLSURF algo
+            self.params.SetPreCADOptionValue(optionName,optionValue)
+
     ## Sets GMF file for export at computation
     #  @ingroup l3_hypos_blsurf
     def SetGMFFile(self, fileName):
@@ -5191,75 +5290,38 @@ class Mesh_Triangle(Mesh_Algorithm):
     #  @param x            : x coordinate
     #  @param y            : y coordinate
     #  @param z            : z coordinate
-    #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertex(self, theFace, x, y, z):
-        if self.Parameters():
-            #  Parameter of BLSURF algo
-            AssureGeomPublished( self.mesh, theFace )
-            return self.params.SetEnforcedVertex(theFace, x, y, z)
-
-    ## To set an enforced vertex as SetEnforcedVertex. The created enforced vertex is identified by a name.
-    #  @param theFace      : GEOM face (or group, compound) on which to define an enforced vertex
-    #  @param x            : x coordinate
-    #  @param y            : y coordinate
-    #  @param z            : z coordinate
-    #  @param vertexName   : name of the enforced vertex
-    #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertexNamed(self, theFace, x, y, z, vertexName):
-        if self.Parameters():
-            #  Parameter of BLSURF algo
-            AssureGeomPublished( self.mesh, theFace )
-            return self.params.SetEnforcedVertexNamed(theFace, x, y, z, vertexName)
-
-    ## To set an enforced vertex on a face (or group, compound) given a GEOM vertex, group or compound.
-    #  @param theFace      : GEOM face (or group, compound) on which to define an enforced vertex
-    #  @param theVertex    : GEOM vertex (or group, compound) to be projected on theFace.
-    #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertexGeom(self, theFace, theVertex):
-        if self.Parameters():
-            #  Parameter of BLSURF algo
-            AssureGeomPublished( self.mesh, theFace )
-            AssureGeomPublished( self.mesh, theVertex )
-            return self.params.SetEnforcedVertexGeom(theFace, theVertex)
-
-    ## To set an enforced vertex as SetEnforcedVertex and add it in the group "groupName".
-    #  @param theFace      : GEOM face (or group, compound) on which to define an enforced vertex
-    #  @param x            : x coordinate
-    #  @param y            : y coordinate
-    #  @param z            : z coordinate
-    #  @param groupName    : name of the group
-    #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertexWithGroup(self, theFace, x, y, z, groupName):
-        if self.Parameters():
-            #  Parameter of BLSURF algo
-            AssureGeomPublished( self.mesh, theFace )
-            return self.params.SetEnforcedVertexWithGroup(theFace, x, y, z, groupName)
-
-    ## To set an enforced vertex as SetEnforcedVertexNamed and add it in the group "groupName".
-    #  @param theFace      : GEOM face (or group, compound) on which to define an enforced vertex
-    #  @param x            : x coordinate
-    #  @param y            : y coordinate
-    #  @param z            : z coordinate
     #  @param vertexName   : name of the enforced vertex
     #  @param groupName    : name of the group
     #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertexNamedWithGroup(self, theFace, x, y, z, vertexName, groupName):
+    def SetEnforcedVertex(self, theFace, x, y, z, vertexName = "", groupName = ""):
         if self.Parameters():
             #  Parameter of BLSURF algo
             AssureGeomPublished( self.mesh, theFace )
-            return self.params.SetEnforcedVertexNamedWithGroup(theFace, x, y, z, vertexName, groupName)
+            if vertexName == "":
+              if groupName == "":
+                return self.params.SetEnforcedVertex(theFace, x, y, z)
+              else:
+                return self.params.SetEnforcedVertexWithGroup(theFace, x, y, z, groupName)
+            else:
+              if groupName == "":
+                return self.params.SetEnforcedVertexNamed(theFace, x, y, z, vertexName)
+              else:
+                return self.params.SetEnforcedVertexNamedWithGroup(theFace, x, y, z, vertexName, groupName)
 
-    ## To set an enforced vertex as SetEnforcedVertexGeom and add it in the group "groupName".
+    ## To set an enforced vertex on a face (or group, compound) given a GEOM vertex, group or compound.
     #  @param theFace      : GEOM face (or group, compound) on which to define an enforced vertex
     #  @param theVertex    : GEOM vertex (or group, compound) to be projected on theFace.
     #  @param groupName    : name of the group
     #  @ingroup l3_hypos_blsurf
-    def SetEnforcedVertexGeomWithGroup(self, theFace, theVertex, groupName):
+    def SetEnforcedVertexGeom(self, theFace, theVertex, groupName = ""):
         if self.Parameters():
             #  Parameter of BLSURF algo
             AssureGeomPublished( self.mesh, theFace )
             AssureGeomPublished( self.mesh, theVertex )
-            return self.params.SetEnforcedVertexGeomWithGroup(theFace, theVertex,groupName)
+            if groupName == "":
+              return self.params.SetEnforcedVertexGeom(theFace, theVertex)
+            else:
+              return self.params.SetEnforcedVertexGeomWithGroup(theFace, theVertex,groupName)
 
     ## To remove an enforced vertex on a given GEOM face (or group, compound) given the coordinates.
     #  @param theFace      : GEOM face (or group, compound) on which to remove the enforced vertex
@@ -5836,80 +5898,79 @@ class Mesh_Tetrahedron(Mesh_Algorithm):
             self.params.SetToRemoveCentralPoint(toRemove)
 
     ## To set an enforced vertex.
+    #  @param x            : x coordinate
+    #  @param y            : y coordinate
+    #  @param z            : z coordinate
+    #  @param size         : size of 1D element around enforced vertex
+    #  @param vertexName   : name of the enforced vertex
+    #  @param groupName    : name of the group
     #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedVertex(self, x, y, z, size):
-        #  Advanced parameter of GHS3D
-        if self.Parameters():
-            return self.params.SetEnforcedVertex(x, y, z, size)
-
-    ## To set an enforced vertex and add it in the group "groupName".
-    #  Only on meshes w/o geometry
-    #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedVertexWithGroup(self, x, y, z, size, groupName):
-        #  Advanced parameter of GHS3D
-        if self.Parameters():
-            return self.params.SetEnforcedVertexWithGroup(x, y, z, size,groupName)
-
-    ## To remove an enforced vertex.
-    #  @ingroup l3_hypos_ghs3dh
-    def RemoveEnforcedVertex(self, x, y, z):
+    def SetEnforcedVertex(self, x, y, z, size, vertexName = "", groupName = ""):
         #  Advanced parameter of GHS3D
         if self.Parameters():
-            return self.params.RemoveEnforcedVertex(x, y, z)
+          if vertexName == "":
+            if groupName == "":
+              return self.params.SetEnforcedVertex(x, y, z, size)
+            else:
+              return self.params.SetEnforcedVertexWithGroup(x, y, z, size, groupName)
+          else:
+            if groupName == "":
+              return self.params.SetEnforcedVertexNamed(x, y, z, size, vertexName)
+            else:
+              return self.params.SetEnforcedVertexNamedWithGroup(x, y, z, size, vertexName, groupName)
 
     ## To set an enforced vertex given a GEOM vertex, group or compound.
+    #  @param theVertex    : GEOM vertex (or group, compound) to be projected on theFace.
+    #  @param size         : size of 1D element around enforced vertex
+    #  @param groupName    : name of the group
     #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedVertexGeom(self, theVertex, size):
+    def SetEnforcedVertexGeom(self, theVertex, size, groupName = ""):
         AssureGeomPublished( self.mesh, theVertex )
         #  Advanced parameter of GHS3D
         if self.Parameters():
+          if groupName == "":
             return self.params.SetEnforcedVertexGeom(theVertex, size)
+          else:
+            return self.params.SetEnforcedVertexGeomWithGroup(theVertex, size, groupName)
 
-    ## To set an enforced vertex given a GEOM vertex, group or compound
-    #  and add it in the group "groupName".
-    #  Only on meshes w/o geometry
+    ## To remove an enforced vertex.
+    #  @param x            : x coordinate
+    #  @param y            : y coordinate
+    #  @param z            : z coordinate
     #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedVertexGeomWithGroup(self, theVertex, size, groupName):
-        AssureGeomPublished( self.mesh, theVertex )
+    def RemoveEnforcedVertex(self, x, y, z):
         #  Advanced parameter of GHS3D
         if self.Parameters():
-            return self.params.SetEnforcedVertexGeomWithGroup(theVertex, size,groupName)
+          return self.params.RemoveEnforcedVertex(x, y, z)
 
     ## To remove an enforced vertex given a GEOM vertex, group or compound.
+    #  @param theVertex    : GEOM vertex (or group, compound) to be projected on theFace.
     #  @ingroup l3_hypos_ghs3dh
     def RemoveEnforcedVertexGeom(self, theVertex):
         AssureGeomPublished( self.mesh, theVertex )
         #  Advanced parameter of GHS3D
         if self.Parameters():
-            return self.params.RemoveEnforcedVertexGeom(theVertex)
-
-    ## To set an enforced mesh.
-    #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedMesh(self, theSource, elementType):
-        #  Advanced parameter of GHS3D
-        if self.Parameters():
-            return self.params.SetEnforcedMesh(theSource, elementType)
-
-    ## To set an enforced mesh and add the enforced elements in the group "groupName".
-    #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedMeshWithGroup(self, theSource, elementType, groupName):
-        #  Advanced parameter of GHS3D
-        if self.Parameters():
-            return self.params.SetEnforcedMeshWithGroup(theSource, elementType, groupName)
-
-    ## To set an enforced mesh with given size.
-    #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedMeshSize(self, theSource, elementType, size):
-        #  Advanced parameter of GHS3D
-        if self.Parameters():
-            return self.params.SetEnforcedMeshSize(theSource, elementType, size)
+          return self.params.RemoveEnforcedVertexGeom(theVertex)
 
     ## To set an enforced mesh with given size and add the enforced elements in the group "groupName".
+    #  @param theSource    : source mesh which provides constraint elements/nodes
+    #  @param elementType  : SMESH.ElementType (NODE, EDGE or FACE)
+    #  @param size         : size of elements around enforced elements. Unused if -1.
+    #  @param groupName    : group in which enforced elements will be added. Unused if "".
     #  @ingroup l3_hypos_ghs3dh
-    def SetEnforcedMeshSizeWithGroup(self, theSource, elementType, size, groupName):
+    def SetEnforcedMesh(self, theSource, elementType, size = -1, groupName = ""):
         #  Advanced parameter of GHS3D
         if self.Parameters():
-            return self.params.SetEnforcedMeshSizeWithGroup(theSource, elementType, size, groupName)
+          if size >= 0:
+            if groupName != "":
+              return self.params.SetEnforcedMesh(theSource, elementType)
+            else:
+              return self.params.SetEnforcedMeshWithGroup(theSource, elementType, groupName)
+          else:
+            if groupName != "":
+              return self.params.SetEnforcedMeshSize(theSource, elementType, size)
+            else:
+              return self.params.SetEnforcedMeshSizeWithGroup(theSource, elementType, size, groupName)
 
     ## Sets command line option as text.
     #  @ingroup l3_hypos_ghs3dh
@@ -6061,9 +6122,9 @@ class Mesh_Projection1D(Mesh_Algorithm):
 class Mesh_Projection2D(Mesh_Algorithm):
 
     ## Private constructor.
-    def __init__(self, mesh, geom=0):
+    def __init__(self, mesh, geom=0, algoName="Projection_2D"):
         Mesh_Algorithm.__init__(self)
-        self.Create(mesh, geom, "Projection_2D")
+        self.Create(mesh, geom, algoName)
 
     ## Defines "Source Face" hypothesis, specifying a meshed face, from where
     #  a mesh pattern is taken, and, optionally, the association of vertices
@@ -6389,7 +6450,7 @@ class Mesh_UseExistingElements(Mesh_Algorithm):
     #  @param UseExisting if ==true - searches for the existing hypothesis created with
     #                     the same parameters, else (default) - creates a new one
     def SourceEdges(self, groups, toCopyMesh=False, toCopyGroups=False, UseExisting=False):
-        if self.algo.GetName() == "Import_2D":
+        if self.algo.GetName() != "Import_1D":
             raise ValueError, "algoritm dimension mismatch"
         for group in groups:
             AssureGeomPublished( self.mesh, group )