Cool. In fact my deleted post was a suggestion to allow arrays to do exactly that (a[:UpperTriangular]), but then I thought I hadn’t thought it all the way through. That is a nice and elegant solution IMO.

Four years late, but here’s an answer anyway. Yes, it must be possible to manipulate triangular matrices using BLAS, because that is what LINPACK does in Fortran.