Intel® oneAPI Threading Building Blocks
Ask questions and share information about adding parallelism to your applications when using this threading library.
Announcements
Welcome to the Intel Community. If you get an answer you like, please mark it as an Accepted Solution to help others. Thank you!

Uisng MKL_lapack function with TBB

Sunny_Beast
Beginner
92 Views
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

{Details: http://software.intel.com/en-us/forums/showthread.php?t=85278}

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;

Sincerly!
Sunny



0 Kudos
0 Replies
Reply