On Mon, 21 Feb 2011 15:33:03 +0500, Andrey Bezrukov <bezrukov.andrey(a)gmail.com>
I try to use meshpy in my work to make 2D and 3D interplolations.
Unfortunately, there is not enough documentation on the web, and mailing
list is not accessible. Is it possible to hear from you ore some one
else some explanations?
this is my conceptual code:
mesh_info = MeshInfo()
print 'Building mesh'
mesh = build(mesh_info, max_volume = 10000, allow_boundary_steiner = True) #
points_out = 
for p in mesh.points:
My problem is, that I want to know exactly how initial set of points
corresponds to output points_out. In most cases there is exact connection:
points_out[0:N] == points_xy[0:N],
where N = len(points_xy).
But sometimes an extra point appears, for example K, such that
points_out[0:K] == points_xy[0:K] and points_out[K+1:N+1] == points_xy[K:N].
What is the right way to identify points?
point_markers is the right tool. They're vertex numbers that persist
through the meshing process.
In documentation it is not
described, what are point_markers should be used for. Attempt to use
point_markers to identify points fails, the program tels me that
point_markers array is not allocated.
Call mesh_obj.point_markers.setup(). (Please let me know if that
helped--if so, this needs to be added to set_points() in meshpy.common.)
Hope this helps,