If you need to store a lot of diagonal 3x3 matrices then an array of diagonal static matrices looks like a good choice, something like [Diagonal(SA[i, 1.0, 2.0]) for i in 1:10]. Size information should be correctly propagated, and if it’s not for some operation, please report it
. By the way, HybridArray doesn’t currently support representing this kind of structure.
1 Like