Under some conditions, the new domegadk feature, which extrapolates the new initial guess based on the previous solution in of the dispersion relation, diverges. The code then loses a given mode and presents a diverging solution for which the frequency goes to +/- infinity until the code crashes.

We should include a safeguard to check that the derivative domegadk does not lead to such a diverging solution.

We could also add a user command that allows to turn the feature on or off as needed. In some examples, the code works much better with this feature, especially when crossing through gamma=0.

The feature has been disabled for now.