Numerical analysis of a first-order computational algorithm for reaction-diffusion equations via the primal-dual hybrid gradient method