Fast identification of oscillatory stability margin and damping margin using continuation of invariant subspaces with sensitivity