use compile-time sized matrices where possible to improve understability of algorithm
Introduce a check for whether a particular decomposition may be safely used.  Seems like good practice.
add another QoR metric for debug: transfer function eigenvalues