To avoid underflow in solvers, both the source and initial solution vectors are rescaled such that the source vector has unit magnitude.  Closes #116.