The finite element method is a popular method for approximating numerically on a computer the solution of partial differential equations. It is based on a variational (weak) formulation of the PDE, followed by discretization on a finite-dimensional ansatz space which reduces the problem to a sparse linear algebra one.