search

FEniCS / fiat

FIAT is now maintained at https://github.com/firedrakeproject/
https://fenicsproject.org/
GNU Lesser General Public License v3.0
34 stars 21 forks source link

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

Open chrisrichardson opened 8 years ago

chrisrichardson commented 8 years ago

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 commented 7 years ago

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.