|
|
|
|
|
|
by enriquto
1998 days ago
|
|
|
Sure, B+ trees are cool. I was just wondering about the use case for element insertion (except maybe for the construction of the matrix). Most of my time is spent in algebraic operations with immutable matrices. Regarding the construction of the matrix, instead of insertions very often you can do it with tricks using kronecker products and the like. |
|
|
Yep, that’s precisely why I used that thing.
> Most of my time is spent in algebraic operations with immutable matrices.
CSR is probably the way to go then.
Also, if you do that math on CPU and your matrix is not randomly sparse but made of small dense areas of non-zeroes, consider less compressed SIMD-friendly variant of CSR, where instead of individual elements you keep dense SIMD vectors.
For example, if you have AVX and need fp64 precision, keep 4 elements long slices of rows there. This way you can use 256-bit loads and FMA instructions to do that math.
Or if you spending most of time multiplying matrices as opposed to matrix*vector, might be reasonable to keep 2x2 dense blocks as opposed to 4x1.