Improve `find_k` algorithm efficiency

[doddgray]

Previously find_k for several eigenmodes used the n-1 mode solutions as inital guesses for the solve_k_single routine for the n-th mode. This strategy was removed in 13a4e53 to improve stability near mode crossings.

Let's fix this