Implemented new sparse representation of soc scaling matrices.

Works fine, but quite a number of iterative refinement steps are needed.

Therefore we compute the error on the linear system and do refinement only if necessary.