#9353: Sage mvrank.pyx doctest failure
-----------------------+----------------------------------------------------
Reporter: retry | Owner: mvngu
Type: defect | Status: new
Priority: major | Milestone: sage-4.5
Component: doctest | Keywords:
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
-----------------------+----------------------------------------------------
[re...@retry-laptop ~]$ sage -t
"devel/sage/sage/libs/mwrank/mwrank.pyx"sage -t
"devel/sage/sage/libs/mwrank/mwrank.pyx"
**********************************************************************File
"/opt/sage/sage/devel/sage/sage/libs/mwrank/mwrank.pyx", line 340:
sage: E.silverman_bound()Expected: 6.5222617951910102Got:
6.5222617951910111**********************************************************************File
"/opt/sage/sage/devel/sage/sage/libs/mwrank/mwrank.pyx", line 372:
sage: E.silverman_bound()Expected: 6.5222617951910102Got:
6.5222617951910111P1 = ![0:1:0] is torsion point, order 1P1 =
![-3:0:1] is generator number 1saturating up to 20...Checking
2-saturation Points have successfully been 2-saturated (max q used =
7)Checking 3-saturation Points have successfully been 3-saturated (max q
used = 7)Checking 5-saturation Points have successfully been 5-saturated
(max q used = 23)Checking 7-saturation Points have successfully been
7-saturated (max q used = 41)Checking 11-saturation Points have
successfully been 11-saturated (max q used = 17)Checking
13-saturation Points have successfully been 13-saturated (max q used =
43)Checking 17-saturation Points have successfully been 17-saturated (max
q used = 31)Checking 19-saturation Points have successfully been
19-saturated (max q used = 37)doneP2 = ![-2:3:1] is generator number
2saturating up to 20...Checking 2-saturation possible kernel vector =
[1,1]This point may be in 2E(Q): ![14:-52:1]...and it is! Replacing old
generator !#1 with new generator ![1:-1:1]
Points have successfully been 2-saturated (max q used = 7)
Index gain = 2!^1
Checking 3-saturation
Points have successfully been 3-saturated (max
q used = 13)
Checking
5-saturation
Points have successfully been 5-saturated (max q used =
67)
Checking 7-saturation
Points have successfully been 7-saturated (max q used = 53)
Checking 11-saturation
Points have
successfully been 11-saturated (max q used = 73)
Checking 13-saturation
Points have successfully been
13-saturated (max q used = 103)
Checking 17-saturation
Points have successfully been 17-saturated
(max q used = 113)
Checking
19-saturation
Points have successfully been 19-saturated (max q used =
47)
done (index = 2).Gained index
2, new generators = [ ![1:-1:1] ![-2:3:1] ]P3 = ![-14:25:8] is generator
number 3saturating up to 20...Checking 2-saturation Points have
successfully been 2-saturated (max q used = 11)Checking
3-saturation Points have successfully been 3-saturated (max q used =
13)Checking 5-saturation Points have successfully been 5-saturated (max q
used = 71)Checking 7-saturation Points have successfully been 7-saturated
(max q used = 101)Checking 11-saturation Points have successfully been
11-saturated (max q used = 127)Checking 13-saturation Points have
successfully been 13-saturated (max q used = 151)Checking
17-saturation Points have successfully been 17-saturated (max q used =
139)Checking 19-saturation Points have successfully been 19-saturated (max
q used = 179)done (index = 1).P4 = ![-1:3:1] = -1*P1 + -1*P2 + -1*P3
(mod torsion)P4 = ![0:2:1] = 2*P1 + 0*P2 + 1*P3 (mod torsion)P4 =
![2:13:8] = -3*P1 + 1*P2 + -1*P3 (mod torsion)P4 = ![1:0:1] = -1*P1
+ 0*P2 + 0*P3 (mod torsion)P4 = ![2:0:1] = -1*P1 + 1*P2 + 0*P3 (mod
torsion)P4 = ![18:7:8] = -2*P1 + -1*P2 + -1*P3 (mod torsion)P4 =
![3:3:1] = 1*P1 + 0*P2 + 1*P3 (mod torsion)P4 = ![4:6:1] = 0*P1 +
-1*P2 + -1*P3 (mod torsion)P4 = ![36:69:64] = 1*P1 + -2*P2 + 0*P3 (mod
torsion)P4 = ![68:-25:64] = -2*P1 + -1*P2 + -2*P3 (mod torsion)P4
= ![12:35:27] = 1*P1 + -1*P2 + -1*P3 (mod torsion)Computed 78519 primes,
largest is 1000253reading primes from file
/home/retry/.sage//temp/retry_laptop/2378//PRIMESread extra prime
10000019finished reading primes from file
/home/retry/.sage//temp/retry_laptop/2378//PRIMESExtra primes in list:
10000019 Saturation index bound = 10WARNING: saturation at primes p > 2
will not be done; points may be unsaturated at primes between 2 and index
boundFailed to saturate MW basis at primes [ ]Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Basic pair: I=336,
J=-10800disc=350922242-adic index bound = 2By Lemma 5.1(a), 2-adic index =
12-adic index = 1One (I,J) pairLooking for quartics with I = 336, J =
-10800Looking for Type 2 quartics:Trying positive a from 1 up to 5 (square
a first...)(1,0,-24,52,-20) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Searching for points (bound =
8)...done: found points of rank 3 and regulator
0.41714355875838396981711954461809339674981010609846Processing points
found during 2-descent...done: now regulator =
0.41714355875838396981711954461809339674981010609846No saturation being
doneBasic pair: I=336, J=-10800disc=350922242-adic index bound = 2By Lemma
5.1(a), 2-adic index = 12-adic index = 1One (I,J) pairLooking for quartics
with I = 336, J = -10800Looking for Type 2 quartics:Trying positive a from
1 up to 5 (square a first...)(1,0,-24,52,-20)
--nontrivial...(!x:y:z) = (1 : 1 : 0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Searching for points (bound =
8)...done: found points of rank 3 and regulator
0.41714355875838396981711954461809339674981010609846Processing points
found during 2-descent...done: now regulator =
0.41714355875838396981711954461809339674981010609846No saturation being
doneBasic pair: I=336, J=-10800disc=350922242-adic index bound = 2By Lemma
5.1(a), 2-adic index = 12-adic index = 1One (I,J) pairLooking for quartics
with I = 336, J = -10800Looking for Type 2 quartics:Trying positive a from
1 up to 5 (square a first...)(1,0,-24,52,-20)
--nontrivial...(!x:y:z) = (1 : 1 : 0)Point = ![4:6:1] height =
1.4667784839930708031879744782308494153040088690549Rank of B=im(eps)
increases to 1(1,0,-12,4,16) --nontrivial...(!x:y:z) = (1 : 1 : 0)Point =
![2:0:1] height =
0.76704335533154620579545064655221715456242461918851Rank of B=im(eps)
increases to 2Trying positive a from 1 up to 5 (...then non-square
a)Trying negative a from -1 down to -1Finished looking for Type 2
quartics.Looking for Type 1 quartics:Trying positive a from 1 up to 6
(square a first...)(1,0,6,28,25) --nontrivial...(!x:y:z) = (1 : 1 :
0)Point = ![-1:3:1] height =
1.2050811041858521515551130942606110675347928058166Rank of B=im(eps)
increases to 3 (The previous point is on the egg)Exiting search for Type 1
quartics after finding one which is globally soluble.Mordell rank
contribution from B=im(eps) = 3Selmer rank contribution from B=im(eps) =
3Sha rank contribution from B=im(eps) = 0Mordell rank contribution
from A=ker(eps) = 0Selmer rank contribution from A=ker(eps) = 0Sha
rank contribution from A=ker(eps) = 0Searching for points (bound =
8)...done: found points of rank 3 and regulator
0.41714355875838396981711954461809339674981010609846Processing points
found during 2-descent...done: now regulator =
0.41714355875838396981711954461809339674981010609846No saturation being
done**********************************************************************2
items had failures: 1 of 6 in !__main!__.example_10 1 of 6 in
!__main!__.example_11***Test Failed*** 2 failures.For whitespace errors,
see the file /home/retry/.sage//tmp/.doctest_mwrank.py [25.3
s] ----------------------------------------------------------------------The
following tests failed: sage -t
"devel/sage/sage/libs/mwrank/mwrank.pyx"Total time for all tests: 25.4
seconds
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9353>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
