The eigensolver was only called once, leading to very bad convergences, making the inversion procedure
absolutely unusable.
The asymptotics was not properly done.

This commits also does the following:
 - Add OpenMP support
 - Refactor and rename variables for more readability