Fix compilation error. Define lexical order for complex numbers.
Also try to clarify the mapping between alpha, eta and Gamma_DP/ST
Reusing same preconditioner for the continuation is a little bit faster.
Add some docstrings. Remove results.g, which is unused (and was wrong too).