[petsc-users] How to get the eigenvectors in Slepc
Sonya Blade
sonyablade2010 at hotmail.com
Mon Apr 15 02:22:15 CDT 2013
>Let me put it more clearly: you are not getting eigenvector entries, your printing statement is
>nonsense (you >print a pointer as a floating point number), so you cannot say the imaginary part is
>nonzero. It is indeed >zero, SLEPc gives the right solution, your program is wrong.
>Jose
Sorry and thank you for clarifying that,
One last question, I got the correct eigenvalues, now I got the
eigenvectors, but they differ from the exact solution.
For example, for the first eigenvalue(2405.247) I got the following eigenvector
set where it differ from the exact solution, what could be the possible reason of that?
Regards,
Row Exact Results SLEPC RESULTS
0 0.2255511 -0.014234
1 -5.2313502 0.330131
2 3.1352583 -0.197855
3 -4.4245184 0.279215
4 0.0898345 -0.005669
5 1.9278406 -0.121659
6 0.0033757 -0.000213
7 -0.7077308 0.044662
8 0.0687009 -0.004335
9 0.1684281 -0.010629
10 -2.81293611 0.177514
11 1.93270712 -0.121966
12 -0.00306213 0.000193
13 0.88278714 -0.055709
14 -0.70857415: 0.044715
15 0.03025516: -0.001909
16 -2.81094417: 0.177388
17 1.12005518: -0.070683
18 2.73596119: -0.172656
19 0.22734020: -0.014347
20 -4.42534221: 0.279267
21 2.22134222: -0.140181
22 -5.00448323: 0.315815
23 0.17399224: -0.01098
24 2.38934725: -0.150783
25 -3.75380226: 0.236889
26 0.09633427: -0.006079
27 0.48140228: -0.03038
28 -1.52250229: 0.09608
29 -0.00132830: 0.000084
30 1.16923331: -0.073786
31 0.09701232: -0.006122
32 0.00268833: -0.00017
33 1.59855934: -0.100879
34 -3.75642735: 0.237054
35 0.13951936: -0.008805
36 1.17316037: -0.074034
37 -1.32576838: 0.083664
38 -0.00203439: 0.000128
39 0.47943940: -0.030256
40 0.09694141: -0.006118
41 0.00071442: -0.000045
42 0.14814343: -0.009349
43 -5.23126844: 0.330126
44 2.33293345: -0.147223
45 3.01784046: -0.190445
46 -5.00416947: 0.315795
47 0.17591748: -0.011101
More information about the petsc-users
mailing list