AtomDiag()
Iterate over flat indices (block_name, n1, n2) for every matrix-valued block.
(block_name, n1, n2)
(block_name, n1, n2) triple.