“…As a result, it is not suitable to directly apply these algorithms to critical eigenvalues tracing, which is a procedure for the calculation of critical eigenvalues on multiple equilibriums. Implicitly restarted Arnoldi (IRA) – with proper matrix transformation [8, 9] or Kryrov subspace dimension switching strategy [10, 11] – is one such algorithm. Here repeating IRA at multiple equilibriums for critical eigenvalues tracing is referred to as the direct approach.…”