An asymmetric weighted residual method for the solution of fluid dynamics equations is described that leads to local operators with a finite difference structure independent of the degree of approximating polynomials. It combines high accuracy on a coarse computational mesh and short computing time per space point and thus brings about sizable reductions of running time, relative to conventional methods.