void print_mat(real**,int,int);
void print_vec(real*,int);
real vec_vec(real*,real*,int);
void mat_mat(real**,int,int,real**,int,real**);
void mat_vec(real**,int,int,real*,real*);