#15535: LinBox: you are running out of primes. 1000 coprime primes found
----------------------------------+------------------------
Reporter: vbraun | Owner:
Type: defect | Status: new
Priority: major | Milestone: sage-6.4
Component: linear algebra | Resolution:
Keywords: random_fail | Merged in:
Authors: | Reviewers:
Report Upstream: N/A | Work issues:
Branch: | Commit:
Dependencies: | Stopgaps:
----------------------------------+------------------------
Comment (by jmantysalo):
Following works for n=2..6. After several hours it fails for n=7 with same
error message "1000 coprime primes". Failing function is eigenvalues(),
but is this actually same bug?
{{{
from itertools import product
S = [0,1]
for n in range(2,10):
smallest_value=1
for m in product(*[S if i>j else [0] if i<j else [1] for i in range(n)
for j in range(n)]):
M = Matrix(ZZ,n,n,m)
t = min( (M*M.transpose()).eigenvalues() )
if t < smallest_value:
smallest_value = t
smallest_matrix=M
print "N=", n, " Value", smallest_value
print smallest_matrix
print '-'*10
}}}
This was tested on Sage 6.3. Also n=8 fails.
--
Ticket URL: <http://trac.sagemath.org/ticket/15535#comment:5>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
