def vector_mul(x, n): for i in range(n): xx = ModVector.random(x) yy = ModVector.random(x) print xx print yy print xx * yy, "\n"
def vector_mul(x,n): for i in range(n): xx = ModVector.random(x) yy = ModVector.random(x) print xx print yy print xx * yy, "\n"
def append_test(x, y, n): for i in range(n): AA = ModMatrix.random(x, y) zz = ModVector.random(x) print "AA\n", AA, "\n" print "zz\n", zz, "\n" print "append\n", AA.get_append_columns([zz, zz + zz]), "\n"
def append_test(x,y,n): for i in range(n): AA = ModMatrix.random(x,y) zz = ModVector.random(x) print "AA\n", AA, "\n" print "zz\n", zz, "\n" print "append\n", AA.get_append_columns([zz,zz + zz]), "\n"
def solve_test(x, y, n): for i in range(n): AA = ModMatrix.random(x, y) bb = ModVector.random(x) print AA.can_solve(bb) xx = AA.solve(bb) print "A\n", AA, "\n" print "rrefA\n", AA.get_rref(), "\n" print "b\n", bb, "\n" print "x\n", xx, "\n" print "Ax\n", AA * xx print AA * xx == ModMatrix([bb]).get_transpose()
def solve_test(x,y,n): for i in range(n): AA = ModMatrix.random(x,y) bb = ModVector.random(x) print AA.can_solve(bb) xx = AA.solve(bb) print "A\n", AA, "\n" print "rrefA\n", AA.get_rref(), "\n" print "b\n", bb, "\n" print "x\n", xx, "\n" print "Ax\n", AA * xx print AA * xx == ModMatrix([bb]).get_transpose()
def leading_index(x, n): for i in range(n): xx = ModVector.random(x) print xx print xx.get_leading_index() print len(xx)
def scalar_mul(x, n): for i in range(n): xx = ModVector.random(x) aa = random.randrange(opts.prime) print aa, xx print aa * xx
def leading_index(x,n): for i in range(n): xx = ModVector.random(x) print xx print xx.get_leading_index() print len(xx)
def scalar_mul(x,n): for i in range(n): xx = ModVector.random(x) aa = random.randrange(opts.prime) print aa, xx print aa * xx