1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
=================================
Test of QR factorization routines
=================================
QR of Real matrix
A = [[0.764873 0.355365 -0.946328 -0.0494359 1.02093]
[-1.13073 -0.556704 0.403829 -0.235643 1.25461]
[0.452336 1.00977 -0.129499 -0.268905 1.00032]
[-0.451817 0.251912 0.28002 -1.53793 -1.12474]
[1.36468 -0.902459 -0.439783 1.67876 0.369076]]
norm(A - Q*R) = 0
A = [[0.611692 -1.37943]
[-0.0799272 -0.534532]
[-0.101654 1.51576]
[1.02917 -0.116042]]
norm(A - Q*R) = 0
A = [[1.18937 0.103149 -0.420697 -0.286478]
[-0.105928 -0.86184 0.297823 0.14747]]
norm(A - Q*R) = 0
QR of Real matrix with pivoting
A = [[-1.17376 0.445168 0.995892 0.185442 -0.627135]
[0.114017 0.361577 -0.0862101 1.09564 -1.64426]
[0.522963 -0.741776 -0.401752 -0.561107 1.41366]
[-0.556702 0.485322 0.945912 0.155321 -1.02616]
[-0.166006 0.252425 -0.153037 1.05107 0.0790642]]
norm(A*mat(P) - Q*R) = 0
A = [[0.303477 -1.13464]
[-0.228092 1.10821]
[0.788165 -1.23059]
[-0.943667 1.58406]]
norm(A*mat(P) - Q*R) = 0
A = [[0.100877 1.73945 0.438779 -0.205583]
[-0.185544 1.5372 0.919884 -0.0438112]]
norm(A*mat(P) - Q*R) = 0
QR of Complex matrix
A = [[-0.433019+1.0476i 0.211864-0.918654i 0.100501+1.05596i 1.01785-0.561079i -0.248716-0.630004i]
[0.653553+0.916184i -0.625162-0.943399i 0.126522-0.348387i -0.23416+0.548687i -0.183664+1.09633i]
[-0.210851-0.0526413i 0.429391-0.538272i -0.2209+0.166731i 0.319975+0.714349i -0.147774-0.213195i]
[-0.0937395+0.192222i 0.561267-0.750045i -0.540764+0.476032i 0.045842-1.01049i -0.169277-0.432241i]
[-0.209132+0.394894i 0.219214-0.141661i 0.176139-0.276432i 0.83024-1.29979i -0.291467+0.780305i]]
norm(A - Q*R) = 0
A = [[-0.675195-0.038722i 0.204412+0.926516i]
[-0.761494+1.61663i 1.29624-0.103014i]
[-0.294119-0.799369i -0.093321+0.28428i]
[0.267765-1.08489i -0.806121+0.623665i]]
norm(A - Q*R) = 0
A = [[-0.869934-0.703924i -0.588264+1.18731i 0.0682477-0.995181i -0.17659+1.26968i]
[0.704778+0.0975581i 0.849613-0.857419i -0.685392-0.413759i 0.438768-0.24119i]]
norm(A - Q*R) = 0
QR of Complex matrix with pivoting
A = [[0.000363858-1.35847i 0.738196-0.527364i 0.415023+0.38752i -0.0139722-0.183752i -0.187205+1.46323i]
[-0.00599554-0.239655i -1.12676-0.697522i 0.885893+1.49547i 0.474554-0.316045i 0.41336-0.0811581i]
[0.0145848+0.68688i 1.26072-1.13755i 0.795048-0.675515i 0.180024+0.14861i -0.151134+0.76442i]
[-0.777352-0.566061i -0.0444372-0.480627i -0.579474-0.369294i -0.117702-1.64864i -0.737179-0.994019i]
[-1.63595+0.769185i -1.55569+0.72902i 0.152014+0.190999i -0.309757+0.322777i -0.172383+0.61953i]]
norm(A*mat(P) - Q*R) = 0
A = [[0.475768-0.141198i -0.394214+0.917501i]
[0.159679+0.374208i 0.864439-1.24452i]
[0.837619-1.19847i -1.25025+0.457809i]
[0.595088+0.598712i -0.265076+0.366291i]]
norm(A*mat(P) - Q*R) = 0
A = [[-0.420762-0.553852i 0.186794-0.836022i -0.927501-0.351825i 0.379794+0.73907i]
[0.812911-0.0722411i -0.964924-0.17228i -0.615892-0.994348i -0.739438+0.079721i]]
norm(A*mat(P) - Q*R) = 0
|
