A stable fast solver for quasi-Helmholtz decomposition methods