Numerical considerations in spectral multidomain methods for biglobal instability analysis of open cavity configurations