dolfin team mailing list archive

[Garth Wells <gnw20@xxxxxxxxx>]

Has this been fixed?

-- 
BarycenterQuadrature gives wrong results for tetrahedron triangle intersection
https://bugs.launchpad.net/bugs/575493
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.

Status in DOLFIN: New

Bug description:
If a tetrahedron and triangle (converted from MeshEntities to CGAL Nef_polyhedron_3) are intersected and resulting polygon has non-vanishing area, BarycenterQuadrature computes quadrature wrong. Actually it gives zero area, since some loops are not entered in the BarycenterQuadrature::compute_quadrature function. (Such a intersection object has surprisingly *NO* half_facets). Check against CGAL is needed.