Hi all,
I am trying to find vertices that are near an edge within some kind of
tolerance.
Imagine lots of short edges intersecting a surface. The edges are
parallel to the surface normal and evenly distributed over the surface.
The surface has a mesh and I have a dictionary of said edges and a
dictionary of mesh nodes (vertices).
So what I'm trying to do is assign all nodes that are within ~3mm of an
edge to that edge.
I used BRepExtrema_DistShapeShape which seems to do the job quite well
for small numbers of points and edges (600 points, 72 edges).
For larger numbers however, the process consumes all available memory
(16GB) and then hangs the machine.
So my questions are:
* Is there a better way to achieve what i want?
* How can I reduce memory consumption?
Unrelated question: why does the conversion from TopoDS_Vertex to gp_Pnt
fail when i write TopoDS_Vertex instead of TopoDS_vertex? What exactly
is the difference? If it's in the docs I must have missed it.
Thanks,
Uwe
#---------------------------------------
pntDic = {} # some dictionary of gp_Pnt
curveDic = {} # some dictionary of TopoDS_Edge
tol = 3
curveToNode = {}
for k in pntDic:
for l in curveDic:
v = BRepBuilderAPI_MakeVertex(pntDic[k])
dist = BRepExtrema_DistShapeShape(v.Vertex(), curveDic[l])
if dist.Value() <= tol:
# assign curve to node
# get normalized edge direction
# get start and end vertex of edge
ex = TopExp_Explorer(curveDic[l], TopAbs_VERTEX)
# no need for a loop since we know for a fact that
# the edge has only one start and one end
c = ex.Current()
cv = TopoDS_vertex(c)
v0 = BRep_Tool_Pnt(cv)
ex.Next()
c = ex.Current()
cv = TopoDS_vertex(c)
v1 = BRep_Tool_Pnt(cv)
dx = v1.X() - v0.X()
dy = v1.Y() - v0.Y()
dz = v1.Z() - v0.Z()
edgeScaling = sqrt(dx ** 2 + dy ** 2 + dz ** 2)
edgeDirection = [dx / edgeScaling,
dy / edgeScaling,
dz / edgeScaling]
curveToNode[k] = (l, edgeDirection)
del ex, v1, v0, dx, dy, dz, c, cv
del edgeScaling, edgeDirection
else:
pass
del v, dist
#---------------------------------------
--
INTALES GmbH Engineering Solutions
Uwe Schlifkowitz -- Software Engineer
Innsbrucker Str. 1, A-6161 Natters
Phone: +43 512 54611120 -- schlifkow...@intales.com
_______________________________________________
Pythonocc-users mailing list
Pythonocc-users@gna.org
https://mail.gna.org/listinfo/pythonocc-users
