Use the discrete inner product , with endpoints . The centered first-difference matrix is skew-Hermitian, while summation by parts gives
Since is real, the centered advection contribution has zero real part. Therefore
This proves the energy contraction for centered drift-diffusion, and hence stability for every fixed real and every mesh width, in the discrete L2 norm. The extra inequality would make the off-diagonal coefficients nonnegative and support a discrete maximum principle, but is not needed for this energy proof. No time-stepping method is involved here.