|
|
Ublas : |
From: Jeffrey Brent McBeth (mcbeth_at_[hidden])
Date: 2007-07-31 22:23:40
On Tue, Jul 31, 2007 at 07:46:00PM -0400, Bumbrlik wrote:
> 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 ?
There are bindings available, so that all the typecasting, workspace
creation, etc are done behind the scenes. The Intel MKL is compatible with
the LAPACK bindings in SVN, and we use those on a daily basis.
Unfortunately, I didn't bring home the work laptop tonight, so I can't
answer your exact question at this time, but basically, you are correct.
Jeff
-- ---------------------------------------------------------------------------- "The man who does not read good books has no advantage over the man who cannot read them." -- Mark Twain ----------------------------------------------------------------------------