Ubuntu-gr Forum

yallou έγραψε:
Ο πίνακας Jac είναι συμμετρικός και τριδιαγώνιος δεν θα έπρεπε οι 2 υπορουτίνες να δίνουν τις ίδιες ιδιοτιμές

Είναι όντως συμμετρικός; Τσέκαρε το.

Κώδικας: Επιλογή όλων

do i=1,n
do j=i+i,n
val=jac(i,j)-jac(j,i)
print val
end do
end do


Αν αυτό σου βγάζει val περίπου 0, κάνε το εξής πείραμα : Στέιλε στην dsyev μία φορά το άνω τμήμα του πίνακα, και μία φορά το κάτω. Βγαίνουν ίδιες μεταξύ των δύο αυτών περιπτώσεων οι ιδιοτιμές;