In [220]:
import numpy as np
import numpy.linalg as LA
import matplotlib.pyplot as plt
from pandas import *
%matplotlib inline

In [296]:
def Arnodi_algo(A, b, k):
    N = A.shape[0]
    Q = np.zeros((N, k), dtype='complex')
    H = np.zeros((k, k), dtype='complex')
    Q[:, [0]] = b / LA.norm(b)
    
    for j in range(k - 1):
        r = A.dot(Q[:,[j]])
        for i in range(j + 1):
            H[i][j] = Q[:, [i]].T.conj().dot(r)[0][0]
            r -= Q[:, [i]] * H[i, j]
        H[j + 1, j] = LA.norm(r)
        Q[:, [j + 1]] = r / H[j + 1, j]
    for i in range(k):
        H[i, k - 1] = Q[:, [i]].T.conj().dot(r)[0][0]
        r -= Q[:, [i]] * H[i, k - 1]
    return Q, H

def Lanczos_algo(A, b, k):
    N = A.shape[0]
    Q = np.zeros((N, k), dtype='complex')
    alpha = np.zeros(k, dtype='complex')
    beta = np.zeros(k - 1, dtype='complex')
    Q[:, [0]] = b / LA.norm(b)
    r = np.dot(A, Q[:,[0]])
    for i in range(k - 1):
        alpha[i] = np.dot(Q[:, [i]].T.conj(), r)[0][0]
        r -= Q[:, [i]] * alpha[i]
        beta[i] = LA.norm(r)
        Q[:, [i + 1]] = r / beta[i]
        r = np.dot(A, Q[:,[i + 1]]) - beta[i] * Q[:,[i]]
    alpha[k - 1] = np.dot(Q[:, [k - 1]].T.conj(), r)[0][0]
    return Q, alpha, beta

In [302]:
np.random.seed(13)
N = 200
k = 200
A = np.random.randn(N, N) + 1j * np.random.randn(N, N)
A = (A + A.T.conj()) / np.sqrt(2)
b = (np.random.randn(N, 1) + 1j * np.random.randn(N, 1)) / np.sqrt(2)

In [303]:
Q_ar, H_ar = Arnodi_algo(A, b, k)
my_id_ar = np.dot(Q_ar.T.conj(), Q_ar)
my_id_ar[np.absolute(my_id_ar) < 0.0001] = 0
DataFrame(np.absolute(my_id_ar));

In [304]:
Q_lan, alpha, beta = Lanczos_algo(A, b, k)
H_lan = np.diag(beta, -1) + np.diag(alpha) + np.diag(beta, 1)
my_id_lan = np.dot(Q_lan.T.conj(), Q_lan)
my_id_lan[np.absolute(my_id_lan) < 0.0001] = 0
DataFrame(np.absolute(my_id_lan));

In [305]:
lam_true = LA.eigvalsh(A)
lam_ar = LA.eigvals(H_ar)
lam_lan = LA.eigvals(H_lan)

In [306]:
plt.plot(lam_true, 'x')
plt.plot(sorted(lam_ar.real))
plt.plot(sorted(lam_lan.real))


Out[306]:
[<matplotlib.lines.Line2D at 0x7f38f03a9b00>]

In [307]:
A[np.random.random(A.shape) + np.eye(N) < .7] = 0

In [308]:
Q_ar, H_ar = Arnodi_algo(A, b, k)
my_id_ar = np.dot(Q_ar.T.conj(), Q_ar)
my_id_ar[np.absolute(my_id_ar) < 0.0001] = 0
DataFrame(np.absolute(my_id_ar));

In [309]:
Q_lan, alpha, beta = Lanczos_algo(A, b, k)
H_lan = np.diag(beta, -1) + np.diag(alpha) + np.diag(beta, 1)
my_id_lan = np.dot(Q_lan.T.conj(), Q_lan)
my_id_lan[np.absolute(my_id_lan) < 0.0001] = 0
DataFrame(np.absolute(my_id_lan))


Out[309]:
0 1 2 3 4 5 6 7 8 9 ... 190 191 192 193 194 195 196 197 198 199
0 1.000000 0.000000 0.571994 0.013005 0.311583 0.016762 0.190299 0.019272 0.123731 0.021760 ... 0.084552 0.063836 0.084558 0.063911 0.084565 0.063985 0.084571 0.064059 0.084578 0.064131
1 0.000000 1.000000 0.000000 0.872911 0.008553 0.710043 0.019407 0.582462 0.028572 0.485571 ... 0.083648 0.088499 0.083831 0.088644 0.084013 0.088783 0.084193 0.088916 0.084370 0.089043
2 0.571994 0.000000 1.000000 0.000000 0.936829 0.007024 0.849464 0.015045 0.769944 0.022126 ... 0.151481 0.098382 0.152118 0.098558 0.152748 0.098731 0.153370 0.098902 0.153985 0.099069
3 0.013005 0.872911 0.000000 1.000000 0.000000 0.961473 0.005283 0.897872 0.011328 0.835696 ... 0.072522 0.081537 0.072888 0.081251 0.073250 0.080995 0.073609 0.080767 0.073964 0.080567
4 0.311583 0.008553 0.936829 0.000000 1.000000 0.000000 0.978385 0.003668 0.937600 0.008031 ... 0.192181 0.083597 0.191823 0.083857 0.191487 0.084113 0.191173 0.084367 0.190879 0.084617
5 0.016762 0.710043 0.007024 0.961473 0.000000 1.000000 0.000000 0.983857 0.002582 0.952615 ... 0.061048 0.176936 0.061494 0.175561 0.061934 0.174223 0.062370 0.172921 0.062802 0.171654
6 0.190299 0.019407 0.849464 0.005283 0.978385 0.000000 1.000000 0.000000 0.988941 0.001852 ... 0.262297 0.068625 0.261005 0.068928 0.259745 0.069227 0.258517 0.069524 0.257319 0.069817
7 0.019272 0.582462 0.015045 0.897872 0.003668 0.983857 0.000000 1.000000 0.000000 0.991525 ... 0.053442 0.264228 0.053908 0.262366 0.054369 0.260542 0.054826 0.258755 0.055277 0.257004
8 0.123731 0.028572 0.769944 0.011328 0.937600 0.002582 0.988941 0.000000 1.000000 0.000000 ... 0.331155 0.056583 0.329304 0.056904 0.327488 0.057222 0.325706 0.057538 0.323958 0.057851
9 0.021760 0.485571 0.022126 0.835696 0.008031 0.952615 0.001852 0.991525 0.000000 1.000000 ... 0.048220 0.337372 0.048681 0.335233 0.049138 0.333131 0.049590 0.331066 0.050038 0.329038
10 0.084356 0.035628 0.700564 0.016693 0.892673 0.005686 0.965416 0.001321 0.993330 0.000000 ... 0.392095 0.047115 0.389904 0.047442 0.387747 0.047766 0.385625 0.048087 0.383537 0.048407
11 0.024230 0.410535 0.028030 0.779693 0.012161 0.917738 0.004207 0.973362 0.000994 0.994841 ... 0.044244 0.399134 0.044691 0.396824 0.045135 0.394551 0.045575 0.392315 0.046010 0.390115
12 0.061514 0.040960 0.640348 0.021105 0.848477 0.008593 0.937142 0.002982 0.978203 0.000719 ... 0.445037 0.039648 0.442630 0.039973 0.440258 0.040296 0.437920 0.040617 0.435614 0.040935
13 0.026594 0.350958 0.032925 0.730035 0.015844 0.882991 0.006585 0.951281 0.002337 0.983014 ... 0.040945 0.452039 0.041378 0.449623 0.041806 0.447243 0.042230 0.444898 0.042651 0.442588
14 0.050230 0.045000 0.587922 0.024635 0.806726 0.011090 0.907522 0.004593 0.959035 0.001653 ... 0.491088 0.033715 0.488544 0.034034 0.486033 0.034352 0.483556 0.034668 0.481111 0.034982
15 0.028806 0.302586 0.037029 0.685963 0.019086 0.849666 0.008834 0.927793 0.003769 0.967787 ... 0.038062 0.497947 0.038478 0.495466 0.038891 0.493020 0.039300 0.490608 0.039707 0.488230
16 0.046878 0.048097 0.542024 0.027431 0.767944 0.013153 0.878144 0.006011 0.938100 0.002578 ... 0.531407 0.028966 0.528777 0.029278 0.526179 0.029588 0.523614 0.029898 0.521080 0.030206
17 0.030853 0.262568 0.040525 0.646630 0.021953 0.818174 0.010919 0.904081 0.005191 0.950847 ... 0.035468 0.538214 0.035870 0.535697 0.036268 0.533214 0.036663 0.530763 0.037055 0.528345
18 0.047980 0.050506 0.501597 0.029639 0.732170 0.014819 0.849755 0.007195 0.916617 0.003400 ... 0.566978 0.025141 0.564298 0.025445 0.561649 0.025748 0.559030 0.026050 0.556442 0.026350
19 0.032739 0.228947 0.043554 0.611295 0.024511 0.788581 0.012848 0.880728 0.006567 0.933122 ... 0.033103 0.573850 0.033490 0.571317 0.033874 0.568817 0.034256 0.566348 0.034635 0.563910
20 0.050971 0.052413 0.465778 0.031385 0.699251 0.016145 0.822688 0.008151 0.895255 0.004082 ... 0.598597 0.022047 0.595891 0.022342 0.593215 0.022636 0.590569 0.022930 0.587951 0.023222
21 0.034474 0.200338 0.046223 0.579349 0.026821 0.760824 0.014639 0.858027 0.007889 0.915150 ... 0.030930 0.605624 0.031304 0.603090 0.031676 0.600587 0.032045 0.598114 0.032411 0.595671
22 0.054484 0.053949 0.433865 0.032771 0.668959 0.017191 0.797060 0.008901 0.874381 0.004621 ... 0.626896 0.019534 0.624182 0.019820 0.621498 0.020105 0.618841 0.020389 0.616213 0.020672
23 0.036071 0.175739 0.048608 0.550299 0.028930 0.734787 0.016314 0.836119 0.009157 0.897251 ... 0.028928 0.634133 0.029290 0.631610 0.029649 0.629115 0.030006 0.626650 0.030361 0.624213
24 0.057929 0.055210 0.405289 0.033879 0.641051 0.018009 0.772879 0.009472 0.854187 0.005024 ... 0.652380 0.017484 0.649673 0.017761 0.646993 0.018037 0.644340 0.018312 0.641714 0.018586
25 0.037542 0.154406 0.050765 0.523743 0.030875 0.710341 0.017887 0.815063 0.010373 0.879619 ... 0.027080 0.659854 0.027430 0.657349 0.027777 0.654872 0.028123 0.652423 0.028466 0.650000
26 0.061085 0.056266 0.379585 0.034772 0.615293 0.018645 0.750100 0.009894 0.834768 0.005305 ... 0.675457 0.015807 0.672766 0.016074 0.670101 0.016341 0.667462 0.016606 0.664848 0.016872
27 0.038900 0.135779 0.052737 0.499356 0.032683 0.687357 0.019375 0.794873 0.011542 0.862370 ... 0.025371 0.683170 0.025710 0.680690 0.026047 0.678236 0.026382 0.675809 0.026714 0.673408
28 0.063893 0.057169 0.356370 0.035500 0.591469 0.019139 0.728647 0.010195 0.816158 0.005481 ... 0.696456 0.014428 0.693789 0.014685 0.691147 0.014942 0.688529 0.015199 0.685935 0.015455
29 0.040156 0.119426 0.054556 0.476867 0.034376 0.665717 0.020787 0.775537 0.012668 0.845573 ... 0.023791 0.704394 0.024119 0.701943 0.024445 0.699518 0.024770 0.697118 0.025092 0.694743
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
170 0.084503 0.081689 0.144734 0.068651 0.197213 0.056329 0.277206 0.048506 0.351759 0.043341 ... 0.998883 0.000104 0.998664 0.000124 0.998427 0.000145 0.998175 0.000167 0.997906 0.000191
171 0.063008 0.086631 0.096431 0.086323 0.080792 0.192918 0.065407 0.285139 0.053207 0.361048 ... 0.000000 0.999013 0.000105 0.998820 0.000125 0.998612 0.000147 0.998389 0.000170 0.998153
172 0.084507 0.081895 0.145438 0.069056 0.196574 0.056824 0.275537 0.049023 0.349513 0.043852 ... 0.999106 0.000000 0.998909 0.000102 0.998694 0.000121 0.998463 0.000142 0.998216 0.000164
173 0.063097 0.086856 0.096642 0.085666 0.081090 0.191121 0.065745 0.282845 0.053559 0.358479 ... 0.000000 0.999211 0.000000 0.999037 0.000103 0.998848 0.000123 0.998644 0.000144 0.998427
174 0.084512 0.082099 0.146135 0.069457 0.195969 0.057313 0.273911 0.049535 0.347310 0.044357 ... 0.999302 0.000000 0.999126 0.000000 0.998933 0.000100 0.998723 0.000119 0.998497 0.000139
175 0.063185 0.087072 0.096850 0.085053 0.081384 0.189371 0.066079 0.280599 0.053907 0.355957 ... 0.000000 0.999384 0.000000 0.999229 0.000000 0.999059 0.000101 0.998874 0.000121 0.998675
176 0.084516 0.082301 0.146827 0.069854 0.195396 0.057797 0.272326 0.050042 0.345151 0.044857 ... 0.999472 0.000000 0.999318 0.000000 0.999146 0.000000 0.998957 0.000000 0.998752 0.000117
177 0.063271 0.087279 0.097054 0.084482 0.081674 0.187668 0.066410 0.278399 0.054252 0.353482 ... 0.000000 0.999534 0.000000 0.999398 0.000000 0.999247 0.000000 0.999081 0.000000 0.998900
178 0.084521 0.082500 0.147512 0.070247 0.194853 0.058276 0.270781 0.050543 0.343033 0.045352 ... 0.999616 0.000000 0.999484 0.000000 0.999333 0.000000 0.999165 0.000000 0.998980 0.000000
179 0.063356 0.087477 0.097254 0.083951 0.081960 0.186009 0.066737 0.276246 0.054594 0.351052 ... 0.000000 0.999662 0.000000 0.999545 0.000000 0.999412 0.000000 0.999264 0.000000 0.999102
180 0.084525 0.082697 0.148191 0.070636 0.194340 0.058751 0.269275 0.051039 0.340957 0.045842 ... 0.999737 0.000000 0.999625 0.000000 0.999495 0.000000 0.999348 0.000000 0.999183 0.000000
181 0.063439 0.087666 0.097451 0.083459 0.082243 0.184393 0.067060 0.274137 0.054933 0.348667 ... 0.000000 0.999768 0.000000 0.999669 0.000000 0.999555 0.000000 0.999426 0.000000 0.999281
182 0.084530 0.082891 0.148863 0.071021 0.193856 0.059220 0.267807 0.051530 0.338920 0.046327 ... 0.999833 0.000000 0.999743 0.000000 0.999633 0.000000 0.999506 0.000000 0.999362 0.000000
183 0.063521 0.087848 0.097644 0.083005 0.082521 0.182821 0.067380 0.272072 0.055269 0.346325 ... 0.000000 0.999853 0.000000 0.999773 0.000000 0.999677 0.000000 0.999565 0.000000 0.999439
184 0.084536 0.083084 0.149528 0.071402 0.193399 0.059684 0.266377 0.052015 0.336922 0.046808 ... 0.999907 0.000000 0.999837 0.000000 0.999748 0.000000 0.999641 0.000000 0.999517 0.000000
185 0.063602 0.088021 0.097833 0.082587 0.082795 0.181290 0.067696 0.270050 0.055602 0.344025 ... 0.000000 0.999918 0.000000 0.999857 0.000000 0.999778 0.000000 0.999684 0.000000 0.999575
186 0.084541 0.083274 0.150186 0.071779 0.192968 0.060144 0.264982 0.052496 0.334963 0.047283 ... 0.999959 0.000000 0.999909 0.000000 0.999841 0.000000 0.999754 0.000000 0.999649 0.000000
187 0.063681 0.088187 0.098019 0.082204 0.083066 0.179799 0.068009 0.268069 0.055932 0.341767 ... 0.000000 0.999964 0.000000 0.999920 0.000000 0.999860 0.000000 0.999783 0.000000 0.999691
188 0.084547 0.083462 0.150837 0.072153 0.192562 0.060598 0.263623 0.052971 0.333041 0.047754 ... 0.999990 0.000000 0.999960 0.000000 0.999911 0.000000 0.999844 0.000000 0.999759 0.000000
189 0.063759 0.088346 0.098202 0.081855 0.083333 0.178348 0.068319 0.266129 0.056259 0.339550 ... 0.000000 0.999991 0.000000 0.999965 0.000000 0.999922 0.000000 0.999863 0.000000 0.999788
190 0.084552 0.083648 0.151481 0.072522 0.192181 0.061048 0.262297 0.053442 0.331155 0.048220 ... 1.000000 0.000000 0.999990 0.000000 0.999961 0.000000 0.999913 0.000000 0.999847 0.000000
191 0.063836 0.088499 0.098382 0.081537 0.083597 0.176936 0.068625 0.264228 0.056583 0.337372 ... 0.000000 1.000000 0.000000 0.999991 0.000000 0.999966 0.000000 0.999924 0.000000 0.999866
192 0.084558 0.083831 0.152118 0.072888 0.191823 0.061494 0.261005 0.053908 0.329304 0.048681 ... 0.999990 0.000000 1.000000 0.000000 0.999990 0.000000 0.999962 0.000000 0.999915 0.000000
193 0.063911 0.088644 0.098558 0.081251 0.083857 0.175561 0.068928 0.262366 0.056904 0.335233 ... 0.000000 0.999991 0.000000 1.000000 0.000000 0.999992 0.000000 0.999966 0.000000 0.999925
194 0.084565 0.084013 0.152748 0.073250 0.191487 0.061934 0.259745 0.054369 0.327488 0.049138 ... 0.999961 0.000000 0.999990 0.000000 1.000000 0.000000 0.999991 0.000000 0.999963 0.000000
195 0.063985 0.088783 0.098731 0.080995 0.084113 0.174223 0.069227 0.260542 0.057222 0.333131 ... 0.000000 0.999966 0.000000 0.999992 0.000000 1.000000 0.000000 0.999992 0.000000 0.999967
196 0.084571 0.084193 0.153370 0.073609 0.191173 0.062370 0.258517 0.054826 0.325706 0.049590 ... 0.999913 0.000000 0.999962 0.000000 0.999991 0.000000 1.000000 0.000000 0.999991 0.000000
197 0.064059 0.088916 0.098902 0.080767 0.084367 0.172921 0.069524 0.258755 0.057538 0.331066 ... 0.000000 0.999924 0.000000 0.999966 0.000000 0.999992 0.000000 1.000000 0.000000 0.999992
198 0.084578 0.084370 0.153985 0.073964 0.190879 0.062802 0.257319 0.055277 0.323958 0.050038 ... 0.999847 0.000000 0.999915 0.000000 0.999963 0.000000 0.999991 0.000000 1.000000 0.000000
199 0.064131 0.089043 0.099069 0.080567 0.084617 0.171654 0.069817 0.257004 0.057851 0.329038 ... 0.000000 0.999866 0.000000 0.999925 0.000000 0.999967 0.000000 0.999992 0.000000 1.000000

200 rows × 200 columns


In [310]:
lam_true = LA.eigvalsh(A)
lam_ar = LA.eigvals(H_ar)
lam_lan = LA.eigvals(H_lan)

In [312]:
plt.plot(lam_true, 'x')
plt.plot(sorted(lam_ar.real))
plt.plot(sorted(lam_lan.real))


Out[312]:
[<matplotlib.lines.Line2D at 0x7f38f026ce10>]

In [275]:
#DataFrame(A.real)


Out[275]:
0 1 2 3 4 5 6 7 8 9 ... 390 391 392 393 394 395 396 397 398 399
0 -1.007473 0.000000 0.000000 0.000000 1.558179 0.000000 0.000000 0.000000 0.812338 -0.476996 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.121949 0.318505 0.596888 -0.544198
1 0.397489 0.684237 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.000000 -1.545368 0.070965 0.000000 0.000000 0.325414 -0.119133
2 0.000000 0.000000 0.371328 1.026638 0.000000 0.000000 0.000000 0.000000 0.000000 -0.908022 ... 0.000000 0.000000 0.000000 0.000000 0.524016 0.000000 0.000000 0.000000 0.000000 0.000000
3 0.000000 0.196665 1.026638 0.106577 0.000000 0.000000 0.658652 0.000000 0.000000 0.079267 ... 0.000000 0.000000 0.000000 -0.359880 0.000000 0.000000 0.703978 0.000000 0.000000 0.453294
4 0.000000 -1.023011 0.000000 0.000000 1.901062 0.000000 0.000000 0.000000 0.000000 0.099173 ... 0.000000 0.000000 0.000000 0.000000 0.000000 -0.172301 0.000000 0.000000 0.000000 0.000000
5 0.000000 1.191557 0.000000 0.000000 -0.307071 -1.541704 0.000000 0.000000 0.000000 0.000000 ... 0.000000 -1.238743 -0.780350 0.000000 0.000000 0.000000 0.000000 0.000000 1.601553 0.000000
6 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.706914 0.000000 -0.335366 0.000000 ... 0.000000 0.000000 0.000000 -0.100597 0.000000 0.000000 0.000000 0.000000 0.000000 -0.862522
7 0.381292 0.000000 0.000000 -0.056006 -0.596711 -1.898772 0.000000 -0.374591 0.000000 0.000000 ... 0.000000 0.000000 -0.405805 0.000000 0.000000 0.000000 0.000000 -1.148865 -0.428393 0.773866
8 0.000000 0.000000 0.000000 0.000000 0.000000 -0.700757 0.000000 0.000000 -0.678923 -0.637436 ... 0.000000 0.000000 0.116828 0.000000 0.036805 0.000000 0.000000 0.000000 0.000000 0.000000
9 0.000000 0.591590 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.637436 0.436354 ... -0.365377 0.000000 0.000000 0.000000 0.000000 0.243423 0.000000 0.000000 1.068770 0.000000
10 -0.902057 -0.120655 0.000000 0.000000 0.000000 0.616563 0.000000 0.000000 -1.529611 -0.491985 ... -0.378176 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.647088 0.952636 0.000000
11 0.000000 -2.267805 -2.086632 0.000000 0.000000 2.294077 0.000000 0.000000 0.000000 0.467775 ... 0.812258 0.319959 0.000000 0.000000 0.000000 0.000000 1.234920 0.999286 0.000000 0.000000
12 0.000000 0.000000 0.000000 0.000000 0.000000 -0.192434 -0.719029 0.510087 0.000000 -0.035702 ... -0.177383 0.000000 0.000000 0.770974 1.590150 0.000000 0.000000 0.000000 0.000000 0.000000
13 0.000000 0.000000 -1.825278 0.674983 0.000000 0.000000 0.000000 -0.949681 0.000000 0.000000 ... 1.934794 0.000000 1.922980 0.000000 0.000000 0.000000 0.000000 0.000000 -0.619418 0.000000
14 0.743784 0.000000 0.000000 0.000000 -0.026338 0.000000 -2.064548 2.678204 0.000000 0.000000 ... -1.007783 0.151818 0.000000 0.000000 0.000000 -2.276282 0.000000 0.530824 0.000000 0.000000
15 0.946454 0.000000 0.000000 0.000000 0.000000 0.000000 0.404244 0.000000 0.000000 -0.330448 ... 1.220831 -0.143522 1.406572 0.389183 1.570672 0.873995 0.750307 0.736997 0.000000 -0.144838
16 0.000000 0.000000 1.343569 0.000000 0.000000 -0.642843 0.000000 0.109701 0.000000 0.000000 ... 1.268191 0.000000 0.160080 0.000000 0.000000 -0.333274 0.000000 0.000000 0.279655 0.263185
17 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.164605 0.000000 -0.175728 0.000000 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.416976 0.000000
18 0.000000 1.536856 0.000000 0.000000 -0.656226 0.000000 0.000000 -1.911210 0.000000 0.000000 ... 0.000000 0.000000 0.000000 -0.041144 0.000000 -0.579381 0.000000 0.000000 0.000000 0.000000
19 0.357996 0.000000 0.000000 0.000000 0.000000 0.000000 -1.941407 0.000000 0.000000 -0.129671 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.981091 0.000000
20 0.000000 0.000000 2.029514 -0.617855 0.000000 -1.118708 0.000000 -0.543274 0.000000 -0.119184 ... 0.000000 -0.316984 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
21 1.290392 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.065644 ... 0.016202 1.009938 0.115320 -0.414784 2.690379 0.779504 0.000000 0.000000 0.000000 -1.370439
22 0.000000 -0.865589 0.000000 0.000000 0.000000 0.000000 0.201773 0.000000 0.401044 0.336854 ... 0.000000 -1.839831 0.300046 0.000000 0.000000 0.399042 0.000000 0.000000 0.000000 0.000000
23 0.000000 0.112785 -1.924978 0.000000 0.136069 0.000000 0.000000 -0.612775 0.000000 0.000000 ... 0.161727 -0.579980 0.000000 1.619896 0.000000 0.000000 0.000000 0.000000 0.000000 1.605283
24 0.000000 0.000000 0.000000 0.000000 0.000000 -0.757277 0.000000 -1.089632 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
25 0.000000 0.000000 0.000000 0.000000 0.197128 0.000000 0.390439 -2.350113 0.000000 2.304953 ... 0.000000 0.456982 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.591682 0.000000
26 -1.366564 0.000000 1.133727 0.000000 0.515547 0.000000 0.000000 1.525873 -1.092390 0.000000 ... 0.000000 0.000000 1.683429 0.000000 -0.756256 0.880976 0.000000 0.000000 0.000000 0.000000
27 0.000000 -0.289998 0.000000 0.000000 0.000000 0.000000 -3.105419 0.000000 0.000000 0.000000 ... 0.000000 -1.513790 0.000000 0.000000 0.000000 0.000000 0.000000 0.440952 0.000000 -0.254144
28 0.000000 0.466033 0.000000 -0.985593 0.000000 0.000000 0.134139 0.000000 0.253170 0.000000 ... 0.000000 0.000000 0.000000 0.252844 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
29 -1.509591 0.000000 0.000000 0.000000 0.000000 0.350887 0.000000 0.000000 0.000000 1.098560 ... 1.888286 0.681064 0.000000 0.000000 -0.709111 0.000000 0.000000 0.000000 0.000000 0.699920
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
370 0.000000 0.000000 -0.174791 0.000000 -0.919253 0.000000 -0.981417 0.000000 0.000000 -0.674539 ... -1.212194 0.000000 0.856122 0.000000 0.000000 0.000000 -0.251064 0.000000 0.000000 1.694221
371 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.000000 1.113272 -0.558146 0.000000 0.000000 -0.809528 0.782712
372 0.000000 1.444730 0.400549 0.000000 0.000000 0.000000 0.000000 0.000000 -0.246775 -0.666793 ... 0.000000 0.000000 0.000000 0.000000 0.406487 2.049441 0.000000 0.000000 -0.548821 0.000000
373 0.000000 0.017115 0.000000 0.729891 -2.086957 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.244962 -0.235946 -0.789702 -0.159370 -0.550853 -0.373280 0.000000 0.000000 -0.743714
374 -0.286686 0.000000 0.000000 0.000000 1.422147 0.000000 -0.851903 0.000000 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
375 0.000000 0.000000 0.000000 -0.350327 -2.210006 0.000000 0.000000 0.000000 0.999765 -0.658310 ... 0.000000 0.000000 0.000000 0.000000 -1.027861 0.000000 0.000000 0.000000 1.367930 -0.403378
376 0.000000 1.363124 0.000000 -0.916621 -0.700371 0.000000 -0.089477 -0.446579 0.899793 0.000000 ... 1.359667 0.000000 0.000000 0.000000 -0.808402 0.000000 0.000000 0.000000 -0.430210 0.000000
377 0.000000 0.000000 0.094079 0.805008 0.000000 0.000000 -0.018525 0.000000 0.000000 0.000000 ... 0.000000 0.000000 0.063996 0.000000 0.000000 0.000000 0.000000 0.000000 0.405337 0.000000
378 0.000000 0.000000 0.000000 -0.533984 0.000000 0.000000 0.000000 -1.238723 1.944544 0.237655 ... 0.000000 0.000000 0.000000 0.000000 0.000000 -0.710796 0.000000 0.000000 -0.642080 0.000000
379 -0.302441 0.158566 0.076208 0.000000 0.000000 0.000000 -1.133025 0.000000 0.000000 -0.836873 ... 0.000000 0.000000 -0.364643 -0.118761 0.384125 -0.224762 -0.936028 0.000000 0.000000 0.000000
380 -0.233551 -0.371163 0.000000 0.000000 0.000000 0.000000 0.000000 -0.728337 0.000000 -1.456232 ... 0.000000 0.000000 0.000000 0.000000 0.058762 0.410779 0.000000 0.000000 0.000000 -0.809283
381 -0.014124 1.399663 -1.586426 1.623769 0.210021 0.000000 0.000000 0.000000 -0.691970 0.000000 ... 0.000000 0.000000 -0.663298 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
382 0.000000 0.000000 -0.176483 0.000000 0.000000 0.000000 0.000000 0.000000 0.331362 0.000000 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.792786 0.000000 0.000000 0.000000 0.000000
383 -1.788750 0.663222 0.000000 0.000000 0.000000 -0.509345 0.000000 0.000000 1.036188 0.000000 ... 0.000000 -0.666414 0.000000 0.000000 1.976116 -0.859110 0.000000 0.000000 -1.162042 0.000000
384 -0.011933 0.000000 0.000000 0.679848 0.000000 1.408628 0.253139 -0.157454 -0.380721 0.000000 ... 0.000000 0.000000 0.340170 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
385 0.587745 1.114174 -0.429772 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.000000 0.000000 0.670511 0.000000 -0.162448 0.000000 0.000000
386 0.572156 0.000000 0.000000 0.000000 0.000000 0.000000 -2.200651 2.461827 0.000000 0.000000 ... 0.000000 0.000000 0.000000 0.348188 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
387 0.000000 0.000000 0.000000 0.797878 -1.611046 -0.412712 1.312795 0.005786 0.000000 0.000000 ... 0.000000 -0.636194 0.849977 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
388 0.384433 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.166455 ... -0.786908 -0.411594 0.838133 -1.149155 0.000000 0.000000 0.103973 0.000000 0.000000 0.000000
389 0.000000 0.000000 -0.330495 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.413120 0.000000 0.000000 0.000000 1.237159 0.000000 0.000000 -0.567580 0.000000
390 0.570070 0.295093 0.000000 -0.954378 0.000000 0.000000 0.000000 0.000000 -1.658608 0.000000 ... 0.785321 0.000000 -0.086981 0.000000 -0.040645 0.000000 -0.976915 0.000000 0.000000 0.265144
391 1.701736 -0.275047 1.000713 -1.049357 0.077675 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 -2.461728 0.813706 0.000000 1.609673 0.000000 0.000000 1.171391 0.000000 0.000000
392 0.000000 0.000000 0.123966 0.156651 1.030053 0.000000 0.000000 0.000000 0.116828 0.000000 ... 0.000000 0.000000 2.166448 0.000000 -0.717920 0.000000 0.000000 0.000000 0.000000 0.000000
393 0.000000 0.000000 -2.057547 0.000000 -0.084543 0.000000 0.000000 -0.185175 0.000000 0.000000 ... -0.687468 -1.041476 0.000000 -1.903664 -1.027905 0.000000 0.000000 0.000000 0.065811 1.312126
394 0.000000 0.000000 0.000000 0.000000 0.000000 -0.715674 -0.179254 1.414403 0.000000 0.000000 ... 0.000000 1.609673 0.000000 0.000000 -0.943360 0.867056 0.000000 0.000000 0.000000 0.000000
395 0.000000 0.000000 0.000000 -1.122359 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.612005 0.000000 0.000000 0.000000 0.046425 -0.543320 1.552151 0.000000 0.561114
396 0.000000 1.004994 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -1.684963 0.000000 ... 0.000000 0.000000 -1.013793 0.000000 0.000000 -0.543320 0.009184 1.103873 0.000000 0.000000
397 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.773370 0.000000 0.803754 0.000000 ... -2.529054 0.000000 0.000000 0.000000 -1.367054 0.000000 0.000000 2.439419 0.000000 0.000000
398 0.000000 0.325414 1.850208 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.068770 ... 0.000000 -1.296475 0.000000 0.065811 0.000000 0.000000 0.000000 0.000000 -0.928730 0.000000
399 -0.544198 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.773866 0.000000 2.321825 ... 0.265144 0.000000 0.000000 0.000000 -0.719529 0.000000 0.115816 0.208657 0.000000 1.704518

400 rows × 400 columns


In [ ]: