mostly done with Matrix.h
SymMatrix