High‐order methods for the numerical solution of the BiGlobal linear stability eigenvalue problem in complex geometries