Implement Moulinec-Suquet basic scheme solver
Requires to add the \Gamma
operator as a new projection and a new solver.
We'll need this to make sure that the ringing artefacts discussed in #18 (closed) are not worse in the Fourier-Galerkin method than the original one.