Converted to a tuple of sparse matrices for storing the values of the state variables.  Would like to test timing when compared to a conj grad solver. Profiling may also be the next step.