Sparse matrix support for larger data sets

[levisc8]

Only good for GEEs, cannot insert 0s into WRM since that will lead to erroneous signals (or lack thereof!)