The origins of the finite element method can be traced back to the 1950s when engineers started to solve numerically structural mechanics problems in aeronautics. Since then, the field of applications has widened steadily and nowadays encompasses nonlinear solid mechanics, fluid/structure interactions, flows in industrial or geophysical settings, multicomponent reactive turbulent flows, mass transfer in porous media, viscoelastic flows in medical sciences, electromagnetism, wave scattering problems, and option pricing (to cite a few examples). Numerous commercial and academic codes based on the finite element method have been developed over the years. The method has been so successful to solve Partial Differential Equations (PDEs) that the term 'Finite Element Method' nowadays refers not only to the mere interpolation technique it is, but also to a fuzzy set of PDEs and approximation techniques. The efficiency of the finite element method relies on two distinct ingredi ents: the
