Fast FEM/BEM computation of SAW harmonic admittance and slowness curves