Inefficient quadrature rules

Issue #20 new
Martin Sandve Alnæs created an issue

An alternative family of quadrature points is the economical Gauss rules found here (reference at top of this file):

http://bazaar.launchpad.net/~syfi-core/fenics-syfi/main/view/611/site-packages/sfc/quadrature/quad_tables.py#L2941

For some choices of degrees it seems much more efficient.

Comments (3)

  1. Lawrence Mitchell

    I note that many of these rules appear to have both negative weights, and points outside the reference cell. This may or may not be desirable. If you want symmetric, all positive weights, points inside the domain, for high order rules, I recommend the rules in Witherden and Vincent (https://arxiv.org/abs/1409.1865).

  2. Log in to comment