#NEWLINE# def test_random_tall_right( self ) : #NEWLINE# #INDENT# m = 200 #NEWLINE# n = 100 #NEWLINE# for k in range ( 2 ) : #NEWLINE# #INDENT# a = random ( [ m , n ] ) #NEWLINE# ( q , r ) = qr ( a , mode = 'economic' ) #NEWLINE# c = random ( [ m ] ) #NEWLINE# ( cq , r ) = qr_multiply ( a , c ) #NEWLINE# assert_array_almost_equal ( dot ( c , q ) , cq ) #NEWLINE# ( cq , r ) = qr_multiply ( a , identity ( m ) ) #NEWLINE# assert_array_almost_equal ( cq , q )