I built the matrix in LIL then converted to CSR like was done in SHO