On Fri, Apr 17, 2009 at 03:40:37PM -0500, Robert Kirby wrote: > Are there any 2d quadrature rules built into dolfin? All I can seem to find > are 1d rules. > Even the dumb collapsed-coordinate ones are fine.
We had this a *long* time ago, see for example http://www.fenics.org/hg/dolfin?f=e26c8159ce5e;file=src/kernel/quadrature/dolfin/TriangleMidpointQuadrature.h But we removed this when FIAT was introduced. :-) The 1D rules are used by the ODE solvers. If you want quadrature rules on a triangle or tetrahedron, you can generate code for a single finite element. Put something like this in a .ufl file: element = FiniteElement("Quadrature", "triangle", 5) Then compile and include the generated code. You can call this function from the generated ufc::dof_map: void tabulate_coordinates(double** coordinates, const cell& c) const = 0; You'll get out the quadrature points you've put into FIAT. :-) -- Anders
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
