Fortran Solving the Fortran LAPACK DSYEV Error with a 3x3 Matrix

  • Thread starter Thread starter tarekph
  • Start date Start date
  • Tags Tags
    Fortran Lapack
AI Thread Summary
The discussion revolves around a beginner's attempt to use the LAPACK DSYEV function in Fortran to compute the eigenvalues of a 3x3 matrix. The user encounters a compilation error indicating that one of the parameters passed to DSYEV has an illegal value. Key issues identified include a syntax error in the matrix declaration, where a comma should replace a period, and the incorrect passing of the matrix variable (Q) instead of the defined matrix (A) in the DSYEV call. Additionally, a resource link to DSYEV documentation is provided to clarify parameter requirements. The conversation emphasizes the importance of correct syntax and parameter usage in Fortran programming.
tarekph
Messages
1
Reaction score
0
Hello
I am a beginner in programming with Fortran and I want to use the LAPACK DSYEV in my program.
I am an example to determine the eigenvalues ​​of a matrix (3x3) after I use it in my program but when I compile I get this error
"" "" "** One entry to DSYEV parameter number 5 Had an illegal value" "" ""
Here is my example:
Program eigenvalue

implicit none

! statements
double precision :: A (3.3), w (3), Work (3)
integer :: i INFO, LWORK, LDA = 3
character :: N, U

! A matrix
A (1 : ) = (/ 3, 1, 3 /)
A (2, : ) = (/ 1, 5, 6 /)
A (3 : ) = (/ 3, 6, 5 /)




CALL DSYEV ('N', 'U', 3, Q, W, LDA, WORK, LWORK, INFO)

! print the solution
do i = 1, 3
write (*, 9) i, w (i)
end do

9 format ('x [' i1 '] =', f5.2)
end programs eigenvalue
I know not how I met the LDA or WORK .. LWORK...
thank you
 
Technology news on Phys.org


tarekph said:
Hello
I am a beginner in programming with Fortran and I want to use the LAPACK DSYEV in my program.
I am an example to determine the eigenvalues ​​of a matrix (3x3) after I use it in my program but when I compile I get this error
"" "" "** One entry to DSYEV parameter number 5 Had an illegal value" "" ""
Here is my example:
Program eigenvalue

implicit none

! statements
double precision :: A (3.3), w (3), Work (3)
integer :: i INFO, LWORK, LDA = 3
character :: N, U

! A matrix
A (1 : ) = (/ 3, 1, 3 /)
A (2, : ) = (/ 1, 5, 6 /)
A (3 : ) = (/ 3, 6, 5 /)




CALL DSYEV ('N', 'U', 3, Q, W, LDA, WORK, LWORK, INFO)

! print the solution
do i = 1, 3
write (*, 9) i, w (i)
end do

9 format ('x [' i1 '] =', f5.2)
end programs eigenvalue
I know not how I met the LDA or WORK .. LWORK...
thank you

I see a couple of mistakes in your code.
1) In this line -- double precision :: A (3.3) -- there should be a comma between the two 3s, not a period.
2) You passed Q in the call to DSYEV. That's the parameter that should be the array A.

I found a link to some documentation on DSYEV - it explains what the parameters need to be - http://www.netlib.org/lapack/double/dsyev.f
 
  • Like
Likes GadhaR
Thread 'Is this public key encryption?'
I've tried to intuit public key encryption but never quite managed. But this seems to wrap it up in a bow. This seems to be a very elegant way of transmitting a message publicly that only the sender and receiver can decipher. Is this how PKE works? No, it cant be. In the above case, the requester knows the target's "secret" key - because they have his ID, and therefore knows his birthdate.

Similar threads

Replies
4
Views
7K
Replies
4
Views
2K
Replies
5
Views
1K
Replies
18
Views
6K
Replies
6
Views
2K
Replies
22
Views
5K
Replies
2
Views
2K
Replies
7
Views
4K
Replies
8
Views
4K
Back
Top