Decided that integrating the IR is a bad idea.
So now explicitly computes the SR using matrix exponential -- much
more accurate.