Efficient computation of coherent multimode instabilities in