Ticket #11358: 11358.sage

File 11358.sage, 631 bytes (added by was, 8 years ago)

this is a usable version of the test program in the ticket description

Line 
1for n in range(2000):
2    set_random_seed(n)
3    A = MatrixSpace(ZZ,3,2).random_element()   
4    B = MatrixSpace(ZZ,2,10).random_element()
5    try_once = A._multiply_multi_modular(B)
6    for k in range(100):
7        set_random_seed(k)
8        try_again = A._multiply_multi_modular(B)
9        if try_once != try_again:
10            print "="*60
11            print "n = %s, k = %s"%(n,k)
12            print "A = "
13            print A
14            print "B ="
15            print B
16            print "first attempt = "
17            print try_once
18            print "k-th retry = "
19            print try_again
20            raise RuntimeError