pub type SMatrix<S> = Matrix<f64, Dynamic, Dynamic, S>;
matrix with arbitrary storage S: Storage<f64, Dynamic, Dynamic>