- const LinAl::SquareMatrix<T, D+1> &get_matrix() const { return matrix; }
- operator const LinAl::SquareMatrix<T, D+1> &() const { return matrix; }
+ const LinAl::Matrix<T, D+1, D+1> &get_matrix() const { return matrix; }
+ operator const LinAl::Matrix<T, D+1, D+1> &() const { return matrix; }