A Fourier method solution for the time dependent Schrödinger equation: A study of the reaction $\bfH^+ + H_2,  D^+ + HD$, and $\bfD^+ + H_2$