Προγραμματισμός σε fortran

[idomeneas]

Ο πίνακας 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 μία φορά το άνω τμήμα του πίνακα, και μία φορά το κάτω. Βγαίνουν ίδιες μεταξύ των δύο αυτών περιπτώσεων οι ιδιοτιμές;

[yallou]

To val είναι ίσο με μηδέν εκτός απο δυο τιμές η μια είναι:-4064.67153393674 και η άλλη:61.9955927398473 .Προφανώς ο πίνακας δεν είναι συμμετρικός.Καλώντας τέλος την dsyev και για τις δυο επιλογές του uplo('u','l') δίνει διαφορετικές ιδιοτιμές.Idomeneas σε ευχαριστώ πάρα πολύ για μια ακόμη μια φορά.

[pc_magas]

@gman
Σου προτεινω να χρισιμοποιειϛ κατευθειαν τον gfortran και όχι κάποιο IDE διότι και εγώ που χριοσιμοποιο τον gcc έχει ποιο πολλέϛ επιλογεϛ στο compiling από ότι ένα IDE

[logari81]

!

Μήνυμα συντονισμού πρός τον/την pc_magas Σε παρακαλούμε να διαβάσεις τις οδηγίες προς νεοεισερχόμενους προσεκτικά. Όλες οι δημοσιεύσεις πρέπει να είναι σε λιτά πεζά ελληνικά, δηλ. όχι greeklish (εκτός αν δεν λειτουργούν τα ελληνικά), όχι κεφαλαία (φαίνεται ότι φωνάζεις), όχι ξένη γλώσσα (εκτός αν υπάρχει ανάγκη), όχι χωρίς λόγο πολλά χρώματα, bold, κλπ. Για αυτούς τους λόγους, κάνε τις απαραίτητες αλλαγές στην δημοσίευση σου.

[pc_magas]

@logari81
Off topic:
Ok ένα λάθος κάναμε!

[yallou]

Έχει επιχειρήσει κανείς να εγκαταστήσει και να δουλέψει με το interface της GSL για την fortran;Ρωτάω γιατί διαβάζοντας το readme σε αυτήν την διεύθυνση http://www.lrz.de/services/software/mathematik/gsl/fortran/ αναφέρει ότι η gfortran 4.4.1 τα "καταφέρνει" με το build αλλά αποτυγχάνει σε κάποιο τεστ.

[kalakouentin]

Εγώ δουλεύω εκτενώς GSL σε C/C++. Από αυτά που διαβάζω φαίνεται να δουλεύει κανονικά με το gfortran και απλά να αποτυγχάνει στην χρήση κάποιων άλλων compilers.