A diagonal splitting method for solving semidiscretized parabolic partial differential equations