Initial implementation of the sparse matrix.

@bsoetaer 
Tests pass when the contracts are disabled.
again contract coupling