- [ ] banded SPD can be implemented as multiplication of two banded triangular matrices - [ ] symmetric orthogonal can be implemented as: QDQ' where D is a diagonal matrix of eigenvalues
