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
| #!/usr/bin/env python3
|
| import sys
| import numpy as np
| import scipy.linalg as la
|
| n = 3
|
| if sys.argv.__len__()!=n*n+1:
| print( "bad # args" )
| sys.exit()
|
|
| A = np.empty([n,n], dtype='f8')
|
| for j in range(n):
| for i in range(n):
| A[j][i] = sys.argv[1+j*n+i]
|
| print( "A:" )
| print( A )
|
| Q,R,P = la.qr(A, pivoting=True)
|
| print( "Q:" )
| print( Q )
|
| print( "R:" )
| print( R )
|
| print( "P:" )
| print( P )
|
|
