LAPACK Archives

[Lapack] Regarding clapack 3.1

Dear Piotr,

Thanks for the response.

In the meantime I have done some investigating myself in this problem, 
but did not find a conclusive solution yet.

Nonetheless I could confirm that the unnecessary static storage class 
definitions cause are one of the main problems to the usage of CLAPACK 
3.0 in OPENMP.  Two different threads running the same CLAPACK routine, 
could corrupt each others static variables.

I did not investigate further, but imagine that this could lead to the 
memory corruption I mentioned, there I think that the static variables 
are indirectly related to a bunch of pointer operations that could 
corrupt my memory if all of a sudden in the middle of a pointer 
operation, the needed variable changes.

The memory corruption I am speaking of is given by an error during 
execution of the program. (the program used routines are zpptrf, zpptri, 
zhptrf, zhptri)  Each run gives a different error, all of them listed 

FATAL: exception not rethrown

*** glibc detected *** ./MatrixElements-pbc.test.out: free(): invalid 
next size (fast): 0x000000000062c4a0 ***

*** glibc detected *** ./MatrixElements-pbc.test.out: free(): invalid 
pointer: 0x00002aaaaab06ff0 ***

terminate called without an active exception

*** glibc detected *** ./MatrixElements-pbc.test.out: malloc(): memory 
corruption: 0x0000000000608130 ***

complex division by zero

Segmentation fault

CLAPACK has been compiled with g++ in the standard way, grabbed straight 
from the site.  (somehow i have troubles with compilation with icc 
(intel c compiler))

The problems occur in a program that does roughly this.  In a loop 
running over the dimension of square hermitian complex matrices, 
matrices are filled with data and inverted with a function (i give a 
rough non working sketch of the function below).  This loop is 
parallelised with openmp.  This makes it possible that 2 function calls 
happen at the same time, and thus corrupt each others data trough the 
static variables.

This program is compiled with icpc (intel c++) with -openmp.
The program is currently solved a bit, because i have placed all clapack 
routines in critical omp pragma directives (only one can run at a time).

Before putting the critical directives, I found out that f2c contains a 
-a flag which changes the default storage class from static to 
automatic.  Would this not be a problem?

I hope this is not to messy and is of any help to you and me both.

With kind regards,

Klaas Vantournhout

(rough sketch of a function)

MatrixCH Inverse(MatrixCH A) {  // MatrixCH is complex hermitian matrix

   lapack::doublecomplex *temp1;

   temp1 = (lapack::doublecomplex *) malloc((n*(n+1))/2 * 
sizeof(lapack::doublecomplex));      // here n is the size of A

// PUT HERE DATA FROM A IN temp1  basically conversion from std::complex
// to lapack::doublecomplex

   long int N = n, info=0;
   char UPLO = 'U';

   lapack::zpptrf_(&UPLO, &N, temp1, &info);

   if (info == 0) {  // zpptrf
     lapack::zpptri_(&UPLO, &N, temp1, &info);
     if (info != 0) exit(1);
   } else if (info > 0) { // zpptrf

     long int *ipiv;
     ipiv = (long int *) malloc (n * sizeof(long int));  // n size of A

     lapack::zhptrf_(&UPLO, &N, temp1, ipiv, &info);
     if (info == 0) {
       lapack::doublecomplex *work;
       work = (lapack::doublecomplex *) malloc (n * 
       lapack::zhptri_(&UPLO, &N, temp1, ipiv, work, &info);

       if (info != 0) exit(1);

       free(work); work = NULL;   // C variable malloc <-> free

     } else if (info != 0) exit(1);

     free(ipiv); ipiv = NULL; // C variable malloc <-> free

   } else if (info < 0) exit (1);

   MatrixCH temp;

   // put data from temp1 in temp

   free(temp1); temp1 = NULL; // C variable malloc <-> free

   return temp;

Piotr Luszczek wrote:
Dear Klaas,

thank you for your input. We're working on CLAPACK 3.1 right now.
I'm guessing that some unnecessary 'static' attributes might have
caused the trouble.

Could you please let me know more details about your problems?
Which routines are you using? What system is it on? What compiler?



Klaas Vantournhout wrote, On 01/06/2007 09:45 AM:
Dear mighty developers of LAPACK and CLAPACK,

In the past year I started to use CLAPACK and was flabbergasted by its 
power.  At the current time I try to implement openmp directives in my 
code to parallelise some procedures and found out that at the usage of 
some CLAPACK functions created severe memory problems.  These got 
solved when stating the CLAPACK functions as critical, so that they 
ran one at a time.

Although the critical pragma directive is a solution, it is not 
preferred.  Simultaneous runs of CLAPACK functions would give me the 
best result.

Lately I read that LAPACK 3.1 was released and that it got threadsafe 
implementations into it.  Therefore I assume that CLAPACK 3.0 is not 
threadsafe and is the cause of the above mentioned problem.  Is this 
correct?  If so, are there plans to create CLAPACK 3.1, or are there 
other solutions to my problem.

With kind regards,

Klaas Vantournhout

"Several billion trillion tons of superhot
exploding hydrogen nuclei rose slowly above
the horizon and managed to look small, cold
and slightly damp."
        Douglas Adams - The Hitch Hickers
                Guide to the Galaxy

<Prev in Thread] Current Thread [Next in Thread>

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or