Pattern for triangular matrices

[devreal]

Some (numerical) algorithms do not require the full matrix but work on the upper or lower triangular matrix. Thus, it would be sufficient to store half of the matrix. It's important that it's not strictly a triangular matrix in that the blocks on the diagonal should remain intact.