def testqr(a, eps):
q,r = qr_decomposition(a,mode='e')
q,r = qr_decomposition(a,mode='r')
q,r = qr_decomposition(a,mode='full')
res = num.matrixmultiply(q,r)
b = num.ravel( num.abs(res - a) )
if num.maximum.reduce(b) > eps:
raise SelftestFailure
else:
print "OK"
def testqr(a, eps):
q, r = qr_decomposition(a, mode='e')
q, r = qr_decomposition(a, mode='r')
q, r = qr_decomposition(a, mode='full')
res = num.matrixmultiply(q, r)
b = num.ravel(num.abs(res - a))
if num.maximum.reduce(b) > eps:
raise SelftestFailure
else:
print "OK"
def cndns(m):
return Numeric.maximum.reduce(Numeric.abs(m).getflat())
def cndns(m):
return Numeric.maximum.reduce(Numeric.abs(m).getflat())
def cndns(m):
return num.maximum.reduce(num.inputarray(num.abs(m)).flat)
def cndns(m):
return num.maximum.reduce(num.inputarray(num.abs(m)).flat)
