Matrix-oriented discretization methods for reaction–diffusion PDEs: Comparisons and applications