Meshpy Triangle Point Attributes

[clockwork85]

Hello everyone,

I am having issues understanding how to use point_attributes with MeshPy, similar to what I'm doing with Rufat's triangle wrapper with vertex attributes. Meshpy's ability to pass a refinement function is something that I also need, and I don't think Rufat's library has this ability. I can't figure out how, or where, to set point_attributes. Here's a somewhat pseudo-code version of what I've tried to do:

from meshpy.triangle import build, MeshInfo
points3d = # some 3d numpy array of length 4310
points2d = points3d[:, :2] # some 2d numpy array of length 4310
elevation_attribs = points3d[:, 2] # some 1d numpy array of attributes - in this case z values for the points to create surface later
info = MeshInfo()
info.set_points(points2d) # can only set markers here
assert len(info.point_attributes) == len(points2d) # has correct length

info.point_attributes = elevation_attribs
# raises AttributeError: can't set attribute

info.point_attributes.resize(len(elevation_attribs))
# RuntimeError: sizes of slave arrays cannot be changed

Does anyone have any experience using point attributes with meshpy? I am a bit at a loss, and I really need both the refinement function ability with point attributes.

Thank you!

[inducer]

test_point_attributes may help:

meshpy/test/test_meshpy.py

Lines 76 to 90 in aac52c4

	def test_point_attributes():
	import meshpy.triangle as triangle
	points = [(1, 1), (-1, 1), (-1, -1), (1, -1)]
	info = triangle.MeshInfo()
	info.set_points(points)
	info.number_of_point_attributes = 2
	info.point_attributes.setup()
	for i in range(len(points)):
	info.point_attributes[i] = [0, 0]
	triangle.build(info)