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
