Add Sym/BiSym HierarchicalOperator

[MikaelSlevinsky]

This would be useful for fractal screens of constant coefficient self-adjoint PDOs.

(Only half the inner products are needed and a simpler factorization can occur as Schur complements are symmetric.)

[MikaelSlevinsky]

Hmm, maybe bisymmetry isn't the right term.