4
#include "lapack_wrapper.h"
7
template<unsigned int M, unsigned int N, class T>
8
void invert(Mat<M,N,T>& matrix) {
15
LAPACK::invert_general_matrix(M, matrix.a, M, ipiv, info);
19
template<unsigned int M, unsigned int N, class T>
20
void least_squares(Mat<M,N,T>&matrix, Vec<M,T>&rhs) {
22
LAPACK::solve_least_squares('N', M, N, 1, matrix.a, M, rhs.v, max(M,N), info);