dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

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     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help   : https://help.launchpad.net/ListHelp