[ublas] uBlas and Intel MKL

From: Bumbrlik (bumbrlik007_at_[hidden])
Date: 2007-07-31 19:46:00

Hello,

I am new to uBlas and I would like to know if what I am doing make sense. Basically, I want to use uBlas with Intel MKL - I realize I could get Atlas for free but I believe MKL should be faster on Intel processors than Atlas. 

This is what I do - when I can I use uBlas. When I need some functionality that is not in uBlas (like Cholesky decomp), I would do something like this:

const int dim = 3; 
// need to make sure that the matrix is stored in row order
matrix<double,column_major> varCovar(dim, dim);
....
// calling Intel MKL function 
dpotrf( &uplo, (int *) (&dim), (double *) &varCovar(0,0), (int *) (&dim), &info );

Does this seem like a reasonable approach ? Is there a better way ? What is all that talk about binding, what kind of advantage would it give me ? 

Thank you,
Bumbrlik