dolfin team mailing list archive

[marc spiegelman <question101686@xxxxxxxxxxxxxxxxxxxxx>]

New question #101686 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/101686

Hi All,  
   I've implemented a rather nifty semi-lagrangian advection scheme in dolfin that requires rapid checking to see if a point intersects a cell (without necessarily firing up the full intersection operator).  Pre 0.9.7 this was implemented as cell.intersects(point).  

Everything works fine in 0.9.6   with the new CGAL intersection operator replacing the GTS intersection operator (with about a 20% speedup).  However, in 0.9.7 the code breaks because  cell.intersects() has been removed as a GTS dependency (actually, looking at TriangleCell.intersects(point), it looks like it simply reproduces GTS code but wasn't actually dependent on GTS) 

The question is... is there is any equivalent functionality in CGAL and is anyone interested in implementing it (either as cell.intersects(??) or some other interface).  The functionality should be straightforward (and I'd be surprised if it didn't exist), given a *known cell* and a point, does the cell intersect the point.  This would be useful for particle based methods as well.  

On a related note...are there any plans for a parallel implementation of Function::eval(Array<double>& values, const Array<double>& x)
by release 1.0?  

All assistance greatly appreciated
marc



-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.