A direct solver for the gradient equation

A new nite element discretization of the equation grad p = g is introduced. This discretization gives rise to an invertible system that can be directly solved, taking a number of operations that is proportional to the number of unknowns. Assuming that g is such that the continuous system has a solution, we obtain an optimal error estimate. We discuss a number of applications related to the Stokes equations