LAPACK/ScaLAPACK Development

Hi i am using CLAPACK to find an inverse of the matrix.
My matrix is symmetric and real-definite positive.

Now, when i use the method sgetrf_ and then sgetri_ then i am able to get the inverse of the matrix. As this method is for general matrices.

But since i already have the symmetric and real-definite positive matrix so i wanted to use the special function so that i can have performance boost ( since my matrix is so size 3500 * 3500).
So for this i am trying to use function spotrf_ and spotri_. But i am getting error in the first function i.e. spotrf_

My test code goes like this:-

Code: Select all

    char UPLO='L';
   integer   n=4;
   integer LDA=4;
   integer info;
   int result=-1;
   int v=0;
   real *A=(real*)malloc(sizeof(real)*16);
   A[0]=4.16; A[1]=-3.12; A[2]=0.56; A[3]=-0.10;
   A[4]=-3.12; A[5]=5.03; A[6]=-0.83; A[7]=1.18;
   A[8]=0.56; A[9]=-0.83; A[10]=0.76; A[11]=0.34;
   A[12]=-0.10; A[13]=1.18; A[14]=0.34; A[15]=1.18;

[b]result=spotrf_(UPLO,&n,A,&LDA,&info);[/b]

It compiles, but when i execute it, it gives me following error:-

Unhandled exception at 0x0108464a in CLAPACK-EXAMPLE.exe: 0xC0000005: Access violation reading location 0x0000004c.

Help me find the error.

Thank You
miztaken

Try &UPLO instead of UPLO in the spotrf call.

thanks that did the trick.. :)