Hi, I'm writing to a sparse 2500x180000 matrix in a column wise random access pattern and I'm facing very strong performance issues which I'm not facing with the dense implementation ( where I'm facing main memory issues ) Is there another way to solve this? Best regards Romeo [[alternative HTML version deleted]]