Find eigenvalues and eigenvectors

[anacondgame]

I Have a 8x8 symmetric matrix A and I wont compute the eigenvalues and eigenvectors with dsyevd function.

The parameter are following:

...
// Var dichiaration

job = 'V';

uplo = 'L';  

n = 8;

lda = 8;

lwork = 3*pow(n,2)+(5+2*3)*n+1;

liwork = 5*n+2;

// alloc the vectors

work = (double*)calloc(lwork,sizeof(double));

w = (double*)calloc(n,sizeof(double));

iwork = (int*)calloc(liwork,sizeof(int));

// call dsyevd for compute the eigenvalues and eigenvectors

dsyevd(&job,&uplo,&n,A,&lda,w,work,&lwork,iwork,&liwork,&info);
...

but the output program is incorrect !!? why?
help me, please!!!

[Zhanghong_T_]

(float *)--->(double *)

or

dsyevd--->ssyevd

[anacondgame]

sorry, I have been wrong to write...
but the problem remains!

[anacondgame]

the "A" matrix is a two-dimensional or one-dimensional array?

A = (double**)malloc(n*double*);

for(i=0; i

  A=(double*)calloc(n,sizeof(double));

 }

or

A = (double*)calloc(n,sizeof(double));

where n = 8 ?????

[TimP]

It's a 2-dimensional Fortran array, not a **.

[Intel_C_Intel]

You can declare an array in C of size n*n of type double for !, and be sure to load in in Fortran style, i.e., column order, as I had stated in an earlier response.

lda will be set to n in this case.

Bruce