Ubuntu-gr Forum

Idomeneas για μια ακόμη φορά είχες δίκιο . Χρησιμοποιώντας την dgeev έχω τις ίδιες ιδιοτιμές με την Octave.
Στην 1η στήλη οι ιδιοτιμες απο την Octave και στην 2η απο την Fortran

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


16375.529 16375.529
16345.97 16345.97
16296.752 16296.752
16228.009 16228.009
16139.909 16139.909
16032.663 16032.663
15906.53 15906.53
15761.815 15761.815
15598.866 15598.866
15418.075 15418.075
15219.879 15219.879
15004.754 15004.754
14773.219 14773.219
14525.831 14525.831
14263.188 14263.188
13985.92 13985.92
13694.697 13694.697
13390.22 13390.22
13073.223 13073.223
12744.468 12744.468
12404.749 12404.749
12054.883 12054.883
11695.714 11695.714
11328.106 11328.106
10952.945 10952.945
10571.136 10571.136
10183.597 10183.597
9791.2632 9791.2632
9395.0787 9395.0787
8995.9982 8995.9982
8594.9833 8594.9833
8193.0000 8193.0000
7791.0167 7791.0167
7390.0018 7390.0018
6990.9213 6990.9213
6594.7368 6594.7368
6202.4027 6202.4027
5814.864 5814.864
5433.0545 5433.0545
5057.894 5057.894
4690.2862 4690.2862
4331.1167 4331.1167
3981.2509 3981.2509
3641.5315 3641.5315
3312.777 3312.777
2995.7794 2995.7794
2691.3024 2691.3024
2400.0794 2400.0794
2122.8121 2122.8121
1860.1683 1860.1683
1612.781 1612.781
1381.2459 1381.2459
1166.1209 1166.1209
967.92434 967.92434
787.13355 787.13355
624.18411 624.18411
479.46856 479.46856
353.3355 353.3355
246.0887 246.0887
157.98637 157.98637
89.24023 89.24023
40.013854 40.013854
10.411118 10.411118
-1.00000 -1.00000
-1.00000 -1.00000


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

Off topic:
Η αντίστοιχη υπορουτίνα με την devcrg της IMSL που καλούσα στο εργαστήριο της σχολής τελικά είναι η dgeev και όχι η dseyv...