Δημοσιεύτηκε: 27 Μάιος 2010, 20:22
Μια ετεροχρονισμένη ερώτηση:Λύνοντας το ίδιο πρόβλημα στην Octave και καλώντας την εντολή j_eig=eig(jac) οι ιδιοτιμές που υπολογίζονται είναι διαφορετικές απο αυτές που έχω υπολογίσει με την gfortran/Lapack. 
Γιατί συμβαίνει αυτό.Αν έχω καταλάβει καλά η Octave χρησιμοποιεί τις ρουτινές της Lapack για αυτούς τους υπολογισμούς.Ενδεικτικά οι τιμές των ιδιοτιμών με την Octave είναι:
και με την fortran:
Γιατί συμβαίνει αυτό.Αν έχω καταλάβει καλά η Octave χρησιμοποιεί τις ρουτινές της Lapack για αυτούς τους υπολογισμούς.Ενδεικτικά οι τιμές των ιδιοτιμών με την Octave είναι:- Κώδικας: Επιλογή όλων
16373.817
16343.292
16292.502
16221.573
16130.681
16020.052
15889.962
15740.733
15572.737
15386.391
15182.159
14960.548
14722.109
14467.436
14197.161
13911.956
13612.53
13299.628
12974.028
12636.539
12288
11929.278
11561.264
11184.874
10801.043
10410.725
10014.891
9614.5259
9210.6236
8804.1889
8396.2322
7987.7678
7579.8111
7173.3764
6769.4741
6369.1085
5973.2749
5582.9574
5199.1263
4822.7361
4454.7223
10.183301
40.707888
91.497872
162.42698
253.31887
363.94757
494.03805
643.26687
811.26303
997.60887
1201.8411
1423.452
1661.8905
1916.5639
2186.8391
2472.044
2771.4696
3084.3715
3409.9718
3747.461
4096
-1.0
-1.0
και με την fortran:
- Κώδικας: Επιλογή όλων
16408.44313358766158
16374.36435669779348
16325.90236972711682
16257.95101013651401
16170.61500033698030
16064.18017674126349
15938.93573956556429
15795.19616984482855
15633.31140972829053
15453.67050850839041
15256.70265539303728
15042.87708480910260
14812.70243109937837
14566.72577411609927
14305.53148520734612
14029.73992740050562
13740.00603858639442
13437.01781467786532
13121.49470387844849
12794.18592022473058
12455.86868303007759
12107.34638806607654
11749.44671592963641
11383.01968287524869
11008.93563935048769
10628.08322150335880
10241.36726100321175
9849.70665861724046
9454.03222710077353
9055.28450908944251
8654.41157582096275
8252.36681266608139
7850.10669761243025
7448.58857902604723
7048.76845921635777
6651.59879055659985
6258.02629116909702
5868.98978747757974
5485.41809126300632
5108.22791923642671
4738.32186356064994
4376.58642220119236
4023.89009844151133
3681.08157931005462
3348.98800295331148
3028.41332501155193
2720.13679359562502
2424.91154119934117
2143.46329935270342
1876.48923742124362
1624.65691998182024
1388.60336704601173
1168.93418806815271
966.22274606003941
781.00929895507397
613.80008014964028
465.06636705391105
335.24387317190332
224.73364008856848
133.90790149036312
63.12622798152661
12.58800685344813
-1.00000000000000
-22.31511396855747
-2049.10568108890038