| #!/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 ) |