On Wednesday April 20 2011 15:07:07 Neilen Marais wrote: > Hi, > > I'm running a dorsal compiled development version of dolfin, and get a > "not implemented" error for mesh.closeset_point(). Is this the current > status of tet mesh support, or did I build it wrong?
It is no implemented. The comment above the method gives the cause: // Partial special for 3D since the nearest_point_3 which is internally used // in CGAL can not yet handles tetrahedrons. // Have to supply myself :) myself here is Andre Massing, I guess ;). Johan > Just want to > check before I dig out an old Bio.KDTree based python implementation I > used in my old code. > > from dolfin import * > > # Define mesh > mesh = UnitCube(1,1,1) > mesh.coordinates()[:] *= [29,23,19] > > p = Point(1,1,1) > cp = mesh.closest_point(p) > > /tmp/python-1810z9J.py in <module>() > 4 mesh = UnitCube(1,1,1) > 5 mesh.coordinates()[:] *= [29,23,19] > 6 > 7 p = Point(1,1,1) > ----> 8 cp = mesh.closest_point(p) > > RuntimeError: *** Error: The function 'compute' has not been > implemented (in > /home/nmarais/src/fenics/dolfin/dolfin/mesh/IntersectionOperatorImplementat > ion.h line 218). > > Thanks > Neilen > > _______________________________________________ > Mailing list: https://launchpad.net/~dolfin > Post to : [email protected] > Unsubscribe : https://launchpad.net/~dolfin > More help : https://help.launchpad.net/ListHelp _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
