TODO

[LeilyR]

• passing a sparse matrix by reference instead of copying it,

• changed the double to float for the output matrix,

• documentation,

• list of dependencies