Add different find_grand_plus methods (computation_variant).

This introduces "computation_variant" (module parametre)
which chooses between different method of find_grand_plus
(module matrice_s) for estimating numerical stability.