Showing results for 
Search instead for 
Did you mean: 

Uisng MKL_lapack function with TBB

Dear Support Provider

I working on an application which uses MKL_LAPACK routine from inside of a C++ program on MAC-Ox 10.6. using icpc

My data is complex type 3D data of A[4x4xn] and B[4x1xn] dimensions;

{n ranges upto 2^48};

I have to find out least square estimation for equation B[4x4xi] = A[4x4xi].X;

I am using zgelss function for performing this routine;

now I have directed by MKL support provider as the task size of matrises are small so internal optimization of routines can not be used

So I have to paralleize the data over the routine and gave example openMP for doing so


I have not used openMP or TBB before. I spent some time reading openMP and find out it is not possible to safely parrallelize a routine which is used during link time.
{Please comment: right/wrong/other}

I belive My sitution is best for using Data Parallel model.

the sequential statement is

for (int loop = 0 ; loop < 1000000; loop++){
zgelss_( &m, &n, &nrhs, a, &lda, b, &ldb, s, &rcond, &rank, work, &lwork, rwork, &info );

My query --> is it possible to parallelize the statement using TBB.
If possible How it can be done;


0 Kudos
0 Replies