Skip to content

getDataAsArray() convert the sparse matrix to a full one. This results in a...

Benoit Fabreges requested to merge fix_mosek_sparse_dense into master

getDataAsArray() convert the sparse matrix to a full one. This results in a bad_alloc when there are a lot of particles.

This MR remove the call to gemv of Mosek and compute the sparse Matrix/Vector product with a simple loop using the COO sparse structure.

Also, the vector u has a wrong size, remove the +1 to stats from the beginning of m_Xlvl, otherwise the size is not right and the Matrix/Vector product makes no sense.

Merge request reports