@Jonas I guess that you can try to use a Julia wrapper to scipy.sparse.linalg.lobpcg — SciPy v1.9.1 Manual
which includes many extra tricks to improve stability compared to the current native Julia code.
1 Like