DynamicMatrix is essentially a vector of vectors, so the memory layout is
not contiguous. We would like to be able to use algorithms from lapack and
blas, so I changed the backend to be a single std::vector plus index arithmetic