eigen_wrapper.h

void matrix_multiply(double**, int, int, double**, int, int, double**);

double check_inversion(double**,double**,int);

double invert_lapack(double**,double**,int,int);

double get_determinant(double*,int*,int*,int);