First commit. Edited function using l2-projection method. In 2D: create a linear system of the projections on edges. In 3D: create a linear system first on edges, then on faces (using the residuals of the edge values)