building up to sparse matrix stuff
some more bits n bobs
Added threshold to pre-run estimation of gamma