AbstractThis article establishes a discrete maximum principle (DMP) for the approximate solution of convection–diffusion–reaction problems obtained from the weak Galerkin (WG) finite element method on nonuniform rectangular partitions. The DMP analysis is based on a simplified formulation of the WG involving only the approximating functions defined on the boundary of each element. The simplified weak Galerkin (SWG) method has a reduced computational complexity over the usual WG, and indeed provides a discretization scheme different from the WG when the reaction terms are present. An application of the SWG on uniform rectangular partitions yields some 5‐ and 7‐point finite difference schemes for the second order elliptic equation. Numerical experiments are presented to verify the DMP and the accuracy of the scheme, particularly the finite difference scheme.