For subsonic flows and upwind-discretized, linearized 1-D Euler equations, the smoothing behavior of multigrid-accelerated point Gauss-Seidel relaxation is analyzed. Error decay by convection across domain boundaries is also discussed. A fix to poor convergence rates at low Mach numbers is sought in replacing the point relaxation applied to unconditioned Euler equations, by locally implicit "time"-stepping applied to preconditioned Euler equations. The locally implicit iteration step is optimized for good damping of high-frequency errors. Numerical inaccuracy at low Mach numbers is also addressed.