getDataAsArray() convert the sparse matrix to a full one. This results in a...
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.