Instead of using two basically identical transfer matrix solver classes,
we merge them into a single class that is parameterized with the solver
and the transfer matrix rhs factory.