Sparse matrix format in Chrono::Flex

[melanz]

Suggestion from Radu:

Apart from the marginally more complicated specification of the index arrays, is there a reason why you use COO format as opposed to say CSR? If I am not mistaken, the cusp::multiply() algorithm should be more efficient for SPMV operations when using CSR vs. COO.