Investigate sparse RHS

The numerical RHS (inhomogeneity) is sparse, but stored as dense vector/matrix, which wastes memory and potentially runtime. Replace by sparse vectors and investigate runtime benefit.