Solving highly-oscillatory NLS with SAM: numerical efficiency and geometric properties