Thank you very much sincerely Julie,

I'm about to complete the compilation, if there doesn't occur any further hassle, of Primme. I though that all of these codes can be used in any platform without any code changing simply by compiling the code with compiler suit terminal targeted for specific platforms. Now I learn that it's not such an easy path .

When I try to solve the eigenvalue problem I encounter, the following errors in output is printed. I can not find any real part in none of the eigenvalues.

My problem dimension is 96 and I'm sure that both of A and B matrices are positive definite and they don't have the zero value in their diagonal

Is there something misconfigured or misinterpreted in the given code?

Eigenvalue( 1) & is numerically infinite or undetermined

ALPHAR( 1) = 2.1902E+06 , ALPHAI( 1) = 0.0000E+00 , BETA( 1) = 1.1978E+01

Eigenvalue( 2) & is numerically infinite or undetermined

ALPHAR( 2) = 2.0576E+06 , ALPHAI( 2) = 0.0000E+00 , BETA( 2) = 1.2556E+01

Eigenvalue( 3) & is numerically infinite or undetermined

ALPHAR( 3) = 1.8441E+06 , ALPHAI( 3) = 0.0000E+00 , BETA( 3) = 1.4328E+01

Eigenvalue( 4) & is numerically infinite or undetermined

ALPHAR( 4) = 1.7180E+06 , ALPHAI( 4) = 0.0000E+00 , BETA( 4) = 1.4965E+01

- Code: Select all
`NB = 64; LWORK = test*(7+NB)`

LDVR =test; N = test; Lda=test; LDb =test;

allocate(A(test,test), AlphaI(test), AlphaR(test), B(test,test), BETA(test), dummy(1,1), VR(LDVR,test), Work(LWork) )

A = assem; B = mass_m;

CALL DGGEV('No left vectors','Vectors (right)',N,A,LDA,B,LDB,ALPHAR,ALPHAI,BETA,DUMMY,1,VR,LDVR,WORK,LWORK,INFO)

Your help will be appreciated,