To Do List

operator for sparse triangular solves with implicit unit diagonal (e.g., for LLDL)
block operators, e.g., [A B ; C D] (in e942fc0)
restriction operator, e.g., A[1:n,1:m].
normest()
condest()