Refactoring Solution::solve(), phases 2 and 3: broke out the population of stiffness matrix and load vector, as well as the solve itself.  Tests continue to pass.