#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: |
-----------------------+----------------------------------------------------
Description changed by retry:
Old description:
> [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
New description:
Sage mvrank.pyx doctest fails.
Sage compiled manually with GCC 4.5.0
Distribution: Arch Linux
--
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9353#comment:1>
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.
