Help improving the performance of my implementation of a lower diagonal storage

See also the discussion in this thread: Triangular Arrays in julia? - #27 by mkborregaard

My feeling is that the factor of 2 storage benefit is not going to be worth it for the performance degradation of giving up BLAS, if you want to do linear algebra (e.g. triangular solves) with the matrix. YMMV depending on what you want to do with the array.

1 Like