double frob_norm(double *A, int M, int N)
{
double norm;
double *work = NULL;
const char norm_type = 'f';
norm = dlange(&norm_type, &M, &N, A, &M, work);
return norm;
}
DLLEXPORT double d_matrix_norm(char norm, MKL_INT m, MKL_INT n, double a[], double work[])
{
return dlange(&norm, &m, &n, a, &m, work);
}