Almost finished sparsity pattern solution
Finished matrix storage allocation