Gauss-Legendre-Lobatto points used for the definition of the Lagrange elements

[chrisrichardson]

Original report by Ben Crestel (Bitbucket: bcrestel, GitHub: bcrestel).

---

There is a general agreement that these interpolation points are better suited for high-order elements. It seems like this should be a direct modification of the Lagrange class (for someone who's familiar with FIAT.....).

[chrisrichardson]

Original comment by David Ham (Bitbucket: David_Ham, ).

---

This is a fair point. However you probably don't want to use FIAT for multiway tensor product (maybe except at the very lowest polynomial degree) because the scaling of assembly using fully tabulated elements (i.e. what FIAT does) is absolutely awful. The solution to this is to use FInAT, which avoids the naïve tabulation and facilitates the form compiler generating the optimal sum factored assembly algorithm.