Non-monotone methods with Shishkin meshes are employed in obtaining finite difference schemes for solving a linear two-dimensional steady state convection-diffusion problem. Preconditioners are used that significantly reduce the number of iterations of the linear solver. Computational results for a Galerkin method are presented which indicate parameter robust, super-linear orders of convergence. (C) 2003 Elsevier B.V. All rights reserved.