Stable nite di erence approximations of convection-di usion equations lead to large sparse linear systems of equations whose coe cient matrix is an M{matrix, which is highly non symmetric when the convection dominates. For an e cient iterative solution of such systems, it is proposed to consider in the non symmetric case an algebraic multilevel preconditioning method formerly proposed for pure ...