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 [ ]:
Content source: LorenzoBi/courses
Similar notebooks: