Andre Massing wrote:
> Hi,
> 
> I am just wondering why GTS has been only optional in the past, although 
> the eval function
> 
> void Function::eval(double* values, const double* x) const
> http://www.fenics.org/hg/dolfin/file/32bcc3b91dd7/dolfin/function/Function.cpp#l243
> 
> in the Function class depends on the IntersectionDetector, which in turn 
> depends on GTS.
> 

GTS was added before eval(double* values, const double* x) for finite 
element functions was added and we haven't bothered reconsidering 
whether or not GTS should be an optional dependency sicne. Most demos 
should run without GTS.

> Furthermore a 1D function can not evaluated at arbitrary points for 
> functions  on 1D grids since the the IntervalCell class
> 
> http://www.fenics.org/hg/dolfin/file/32bcc3b91dd7/dolfin/mesh/IntervalCell.cpp#l183
> 
> does not provide a implementation.
> 
> So I wondering how these things work instead? At least I could not find 
> it in the docs and my benchmarks throw an exception when I tried to run 
> the old IntersectionDetector for 1D meshes and points intersection.
> 

Probably never been tested before.

Garth

> (Of course the forthcoming CGAL interface will remedy these issues ;))
> 
> Greetings,
> Andre
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@fenics.org
> http://www.fenics.org/mailman/listinfo/dolfin-dev


_______________________________________________
DOLFIN-dev mailing list
DOLFIN-dev@fenics.org
http://www.fenics.org/mailman/listinfo/dolfin-dev

Reply via email to