def setUp() :
rowSize = 0
colSize = 0
matrixLocation = input("What's the filename of the test file? (please have the file in the same directory as the code)")
matrix = common.matrix_from_file(matrixLocation)
return matrix
if j >= i:
# if not at position pertaining to current x (i.e. a diagonal entry)
if i != j:
d[i] = d[i] - (r[i, j] * d[j])
elif i == j:
x[i] = d[i] / r[i, j]
return (x, common.error(common.mult_mat(a, x), b))
if __name__ == "__main__":
invalid = True
while invalid:
fname = raw_input("Please enter a filename: ")
try:
matrix = common.matrix_from_file(fname)
invalid = False
except:
print "Invalid filename, please try again"
print "Householder Xsol:"
result = solve_qr_b_househ(matrix[:, :-1], matrix[:, -1])
print result[0]
print "\nHouseholder Error:"
print result[1]
print "\nGivens Xsol:"
result = solve_qr_b_househ(matrix[:, :-1], matrix[:, -1])
print result[0]
print "\nGivens Error:"
print result[1]