Ticket #11339: remove_polys.patch

File remove_polys.patch, 14.2 KB (added by vbraun, 10 years ago)

Remove some initialization form ell_curve_isogeny.py for debugging.

  • sage/schemes/elliptic_curves/ell_curve_isogeny.py

    # HG changeset patch
    # User Volker Braun <vbraun@stp.dias.ie>
    # Date 1309010285 -7200
    # Node ID 9e4521759128a07327e2b5161d225788f7431dc5
    # Parent  6ca08864b80c1e04599ba23ba295bfcf6a61f19b
    imported patch remove_polys.patch
    
    diff --git a/sage/schemes/elliptic_curves/ell_curve_isogeny.py b/sage/schemes/elliptic_curves/ell_curve_isogeny.py
    a b  
    39123912#Psi7 = sum([X**i * f.coefficient({x:i})(X,t) for i in range(4)]).monic()
    39133913
    39143914
    3915 Psi2 = X + (t + 16)/(t - 8)
    3916 Psi3 = X + (t**2 + 30*t + 81)/(t**2 + 18*t - 27)
    3917 Psi5 = X**2 + ((2*t**2 + 20*t + 10)/(t**2 + 4*t - 1))*X + (t**6 + 42*t**5 + 639*t**4 + 4300*t**3 + 12015*t**2 + 9450*t + 2225)/(t**6 + 30*t**5 + 315*t**4 + 1300*t**3 + 1575*t**2 - 978*t + 125)
    3918 Psi7 = X**3 + (3*t**4 + 54*t**3 + 345*t**2 + 774*t + 147)/(t**4 + 14*t**3 + 63*t**2 + 70*t -    7)*X**2 + (3*t**8 + 108*t**7 + 1614*t**6 + 12864*t**5 + 57945*t**4 + 142464*t**3 +    163710*t**2 + 51708*t + 4851)/(t**8 + 28*t**7 + 322*t**6 + 1904*t**5 + 5915*t**4 +    8624*t**3 + 4018*t**2 - 980*t + 49)*X + (t**12 + 54*t**11 + 1269*t**10 + 17186*t**9 +    148602*t**8 + 854574*t**7 + 3288621*t**6 + 8273166*t**5 + 12820026*t**4 + 10888834*t**3    + 4132533*t**2 + 697110*t + 43169)/(t**12 + 42*t**11 + 777*t**10 + 8246*t**9 +    54810*t**8 + 233730*t**7 + 628425*t**6 + 999306*t**5 + 801738*t**4 + 159838*t**3 -    93639*t**2 + 10290*t - 343)
    3919 Psi13=  X**6 + (6*t**6 + 72*t**5 + 408*t**4 + 1260*t**3 + 2136*t**2 + 1512*t + 78)/(t**6 + 10*t**5     + 46*t**4 + 108*t**3 + 122*t**2 + 38*t - 1)*X**5 + (15*t**14 + 450*t**13 + 6495*t**12     + 59124*t**11 + 376152*t**10 + 1755732*t**9 + 6148842*t**8 + 16247820*t**7 +              32110602*t**6 + 46286868*t**5 + 46339800*t**4 + 29344116*t**3 + 9588735*t**2 +            826626*t + 20943)/(t**14 + 26*t**13 + 325*t**12 + 2548*t**11 + 13832*t**10 + 54340*t**9    + 157118*t**8 + 333580*t**7 + 509366*t**6 + 534820*t**5 + 354536*t**4 + 124852*t**3 +      15145*t**2 - 982*t + 13)*X**4 + (20*t**20 + 840*t**19 + 17060*t**18 + 222056*t**17 +     2070768*t**16 + 14658344*t**15 + 81442892*t**14 + 362320920*t**13 + 1305160220*t**12 +    3824743728*t**11 + 9110443680*t**10 + 17528242992*t**9 + 26902547308*t**8 +              32277253368*t**7 + 29335252636*t**6 + 19226924776*t**5 + 8380801392*t**4 +               2097809704*t**3 + 224201908*t**2 + 10406568*t + 175396)/(t**20 + 36*t**19 + 631*t**18     + 7102*t**17 + 57192*t**16 + 348178*t**15 + 1652611*t**14 + 6221436*t**13 +               18725317*t**12 + 45094452*t**11 + 86379696*t**10 + 129896364*t**9 + 149994455*t**8 +      128453916*t**7 + 77238065*t**6 + 29760110*t**5 + 6132072*t**4 + 332258*t**3 -             50875*t**2 + 1476*t - 13)*X**3 + (15*t**28 + 900*t**27 + 26250*t**26 + 495564*t**25 +    6798111*t**24 + 72094704*t**23 + 613899300*t**22 + 4303114008*t**21 +                    25254801624*t**20 + 125568545352*t**19 + 533156573604*t**18 + 1942964145168*t**17 +      6093526623066*t**16 + 16455167665008*t**15 + 38205334804860*t**14 +                     75996871379184*t**13 + 128765242080378*t**12 + 184277137115280*t**11 +                  220161982061220*t**10 + 216123510416520*t**9 + 170576695582296*t**8 +                   105031037279064*t**7 + 48334104078756*t**6 + 15601810268784*t**5 + 3207779251023*t**4    + 367285489644*t**3 + 22695568842*t**2 + 711871524*t + 8913567)/(t**28 + 52*t**27 +      1326*t**26 + 21996*t**25 + 265785*t**24 + 2484144*t**23 + 18623020*t**22 +                114646168*t**21 + 588732456*t**20 + 2548318760*t**19 + 9358879244*t**18 +                29265789072*t**17 + 77993055062*t**16 + 176864657072*t**15 + 339914911284*t**14 +        549949737584*t**13 + 741655098470*t**12 + 822222337872*t**11 + 735158601164*t**10 +      516213585992*t**9 + 274020098664*t**8 + 103737033400*t**7 + 25289774380*t**6 +           3099363696*t**5 - 6620367*t**4 - 26498628*t**3 + 1358094*t**2 - 25532*t + 169)*X**2     + (6*t**34 + 432*t**33 + 15132*t**32 + 344052*t**31 + 5709918*t**30 + 73699992*t**29 +    769487790*t**28 + 6672103956*t**27 + 48942245532*t**26 + 307804367136*t**25 +    1676073313398*t**24 + 7958878083744*t**23 + 33126710022252*t**22 +    121277328806544*t**21 + 391336247726520*t**20 + 1113873747905400*t**19 +    2795560626375132*t**18 + 6176928558284304*t**17 + 11982419368604844*t**16 +    20324125159073592*t**15 + 29974811649927576*t**14 + 38159370086015376*t**13 +    41537814088619580*t**12 + 38195725056468384*t**11 + 29209678804898958*t**10 +    18203416490394912*t**9 + 9000165833061804*t**8 + 3405735803167284*t**7 +    939484884915270*t**6 + 177123317266200*t**5 + 21220677934758*t**4 +    1562246408532*t**3 + 68267386188*t**2 + 1624644528*t + 16229070)/(t**34 + 62*t**33 +    1892*t**32 + 37756*t**31 + 552479*t**30 + 6303400*t**29 + 58229885*t**28 +    446585620*t**27 + 2893400900*t**26 + 16033796610*t**25 + 76671686091*t**24 +    318352372832*t**23 + 1152541108082*t**22 + 3646932539236*t**21 + 10094977882824*t**20    + 24431638381728*t**19 + 51596515180206*t**18 + 94759230136176*t**17 +    150580571258898*t**16 + 205622547731472*t**15 + 239093561232552*t**14 +    233930961126044*t**13 + 189610795016734*t**12 + 124713933296992*t**11 +    64709325948813*t**10 + 25426241898798*t**9 + 7087502758996*t**8 + 1233476382836*t**7    + 88878726643*t**6 - 6438268888*t**5 - 837389711*t**4 + 75009548*t**3 - 2307692*t**2 +    31954*t - 169)*X + (t**42 + 90*t**41 + 3951*t**40 + 112976*t**39 + 2368311*t**38 +    38817366*t**37 + 517880995*t**36 + 5779780752*t**35 + 55024462326*t**34 +    453373071376*t**33 + 3268823312982*t**32 + 20799376011552*t**31 +    117570250802657*t**30 + 593413518802194*t**29 + 2684969136846585*t**28 +    10922547545780848*t**27 + 40033688945917857*t**26 + 132383684067116478*t**25 +    395225474841339419*t**24 + 1065301721692202256*t**23 + 2590965176032226004*t**22 +    5679139393727958880*t**21 + 11197539488052244692*t**20 + 19809141347016503568*t**19    + 31337344703196264683*t**18 + 44145543697033615230*t**17 +    55090950408765652881*t**16 + 60516047590917322480*t**15 + 58058828547614202153*t**14    + 48186407491493228562*t**13 + 34193160536571257393*t**12 +    20445602003818379424*t**11 + 10116440549699970582*t**10 + 4048631186832794128*t**9 +    1273256824201491894*t**8 + 303519895786350096*t**7 + 52538628166252339*t**6 +    6321448897556694*t**5 + 513273187133223*t**4 + 27378446805968*t**3 +    916904893215*t**2 + 17467984794*t + 144422161)/(t**42 + 78*t**41 + 3003*t**40 +    75920*t**39 + 1414959*t**38 + 20695506*t**37 + 246960727*t**36 + 2468218896*t**35 +    21049407990*t**34 + 155295575344*t**33 + 1001454266670*t**32 + 5689547290752*t**31 +    28649170403321*t**30 + 128443250574150*t**29 + 514438663036389*t**28 +    1845008125106224*t**27 + 5933895362321121*t**26 + 17125495350262650*t**25 +    44346853027068215*t**24 + 102951022273700016*t**23 + 213919716149439924*t**22 +    396901875045903136*t**21 + 655411304131422852*t**20 + 959204141060458800*t**19 +    1237542066720323195*t**18 + 1398213766917750042*t**17 + 1372033638327823605*t**16 +    1157344876243332400*t**15 + 828397238649616545*t**14 + 494847675800288070*t**13 +    241330300300762205*t**12 + 93195177370728576*t**11 + 27216183343932726*t**10 +    5548446154604080*t**9 + 655832005477038*t**8 + 13945989607056*t**7 -    5655555970805*t**6 - 193929669102*t**5 + 43376116875*t**4 - 2043712624*t**3 +    45287151*t**2 - 497874*t + 2197)
    3920 
    3921 
    3922 isogeny_kernel_polys = dict([(5,Psi5),(7,Psi7),(13,Psi13)])
     3915#Psi2 = X + (t + 16)/(t - 8)
     3916#Psi3 = X + (t**2 + 30*t + 81)/(t**2 + 18*t - 27)
     3917#Psi5 = X**2 + ((2*t**2 + 20*t + 10)/(t**2 + 4*t - 1))*X + (t**6 + 42*t**5 + 639*t**4 + 4300*t**3 + 12015*t**2 + 9450*t + 2225)/(t**6 + 30*t**5 + 315*t**4 + 1300*t**3 + 1575*t**2 - 978*t + 125)
     3918#Psi7 = X**3 + (3*t**4 + 54*t**3 + 345*t**2 + 774*t + 147)/(t**4 + 14*t**3 + 63*t**2 + 70*t -    7)*X**2 + (3*t**8 + 108*t**7 + 1614*t**6 + 12864*t**5 + 57945*t**4 + 142464*t**3 +    163710*t**2 + 51708*t + 4851)/(t**8 + 28*t**7 + 322*t**6 + 1904*t**5 + 5915*t**4 +    8624*t**3 + 4018*t**2 - 980*t + 49)*X + (t**12 + 54*t**11 + 1269*t**10 + 17186*t**9 +    148602*t**8 + 854574*t**7 + 3288621*t**6 + 8273166*t**5 + 12820026*t**4 + 10888834*t**3    + 4132533*t**2 + 697110*t + 43169)/(t**12 + 42*t**11 + 777*t**10 + 8246*t**9 +    54810*t**8 + 233730*t**7 + 628425*t**6 + 999306*t**5 + 801738*t**4 + 159838*t**3 -    93639*t**2 + 10290*t - 343)
     3919#Psi13=  X**6 + (6*t**6 + 72*t**5 + 408*t**4 + 1260*t**3 + 2136*t**2 + 1512*t + 78)/(t**6 + 10*t**5     + 46*t**4 + 108*t**3 + 122*t**2 + 38*t - 1)*X**5 + (15*t**14 + 450*t**13 + 6495*t**12     + 59124*t**11 + 376152*t**10 + 1755732*t**9 + 6148842*t**8 + 16247820*t**7 +              32110602*t**6 + 46286868*t**5 + 46339800*t**4 + 29344116*t**3 + 9588735*t**2 +            826626*t + 20943)/(t**14 + 26*t**13 + 325*t**12 + 2548*t**11 + 13832*t**10 + 54340*t**9    + 157118*t**8 + 333580*t**7 + 509366*t**6 + 534820*t**5 + 354536*t**4 + 124852*t**3 +      15145*t**2 - 982*t + 13)*X**4 + (20*t**20 + 840*t**19 + 17060*t**18 + 222056*t**17 +     2070768*t**16 + 14658344*t**15 + 81442892*t**14 + 362320920*t**13 + 1305160220*t**12 +    3824743728*t**11 + 9110443680*t**10 + 17528242992*t**9 + 26902547308*t**8 +              32277253368*t**7 + 29335252636*t**6 + 19226924776*t**5 + 8380801392*t**4 +               2097809704*t**3 + 224201908*t**2 + 10406568*t + 175396)/(t**20 + 36*t**19 + 631*t**18     + 7102*t**17 + 57192*t**16 + 348178*t**15 + 1652611*t**14 + 6221436*t**13 +               18725317*t**12 + 45094452*t**11 + 86379696*t**10 + 129896364*t**9 + 149994455*t**8 +      128453916*t**7 + 77238065*t**6 + 29760110*t**5 + 6132072*t**4 + 332258*t**3 -             50875*t**2 + 1476*t - 13)*X**3 + (15*t**28 + 900*t**27 + 26250*t**26 + 495564*t**25 +    6798111*t**24 + 72094704*t**23 + 613899300*t**22 + 4303114008*t**21 +                    25254801624*t**20 + 125568545352*t**19 + 533156573604*t**18 + 1942964145168*t**17 +      6093526623066*t**16 + 16455167665008*t**15 + 38205334804860*t**14 +                     75996871379184*t**13 + 128765242080378*t**12 + 184277137115280*t**11 +                  220161982061220*t**10 + 216123510416520*t**9 + 170576695582296*t**8 +                   105031037279064*t**7 + 48334104078756*t**6 + 15601810268784*t**5 + 3207779251023*t**4    + 367285489644*t**3 + 22695568842*t**2 + 711871524*t + 8913567)/(t**28 + 52*t**27 +      1326*t**26 + 21996*t**25 + 265785*t**24 + 2484144*t**23 + 18623020*t**22 +                114646168*t**21 + 588732456*t**20 + 2548318760*t**19 + 9358879244*t**18 +                29265789072*t**17 + 77993055062*t**16 + 176864657072*t**15 + 339914911284*t**14 +        549949737584*t**13 + 741655098470*t**12 + 822222337872*t**11 + 735158601164*t**10 +      516213585992*t**9 + 274020098664*t**8 + 103737033400*t**7 + 25289774380*t**6 +           3099363696*t**5 - 6620367*t**4 - 26498628*t**3 + 1358094*t**2 - 25532*t + 169)*X**2     + (6*t**34 + 432*t**33 + 15132*t**32 + 344052*t**31 + 5709918*t**30 + 73699992*t**29 +    769487790*t**28 + 6672103956*t**27 + 48942245532*t**26 + 307804367136*t**25 +    1676073313398*t**24 + 7958878083744*t**23 + 33126710022252*t**22 +    121277328806544*t**21 + 391336247726520*t**20 + 1113873747905400*t**19 +    2795560626375132*t**18 + 6176928558284304*t**17 + 11982419368604844*t**16 +    20324125159073592*t**15 + 29974811649927576*t**14 + 38159370086015376*t**13 +    41537814088619580*t**12 + 38195725056468384*t**11 + 29209678804898958*t**10 +    18203416490394912*t**9 + 9000165833061804*t**8 + 3405735803167284*t**7 +    939484884915270*t**6 + 177123317266200*t**5 + 21220677934758*t**4 +    1562246408532*t**3 + 68267386188*t**2 + 1624644528*t + 16229070)/(t**34 + 62*t**33 +    1892*t**32 + 37756*t**31 + 552479*t**30 + 6303400*t**29 + 58229885*t**28 +    446585620*t**27 + 2893400900*t**26 + 16033796610*t**25 + 76671686091*t**24 +    318352372832*t**23 + 1152541108082*t**22 + 3646932539236*t**21 + 10094977882824*t**20    + 24431638381728*t**19 + 51596515180206*t**18 + 94759230136176*t**17 +    150580571258898*t**16 + 205622547731472*t**15 + 239093561232552*t**14 +    233930961126044*t**13 + 189610795016734*t**12 + 124713933296992*t**11 +    64709325948813*t**10 + 25426241898798*t**9 + 7087502758996*t**8 + 1233476382836*t**7    + 88878726643*t**6 - 6438268888*t**5 - 837389711*t**4 + 75009548*t**3 - 2307692*t**2 +    31954*t - 169)*X + (t**42 + 90*t**41 + 3951*t**40 + 112976*t**39 + 2368311*t**38 +    38817366*t**37 + 517880995*t**36 + 5779780752*t**35 + 55024462326*t**34 +    453373071376*t**33 + 3268823312982*t**32 + 20799376011552*t**31 +    117570250802657*t**30 + 593413518802194*t**29 + 2684969136846585*t**28 +    10922547545780848*t**27 + 40033688945917857*t**26 + 132383684067116478*t**25 +    395225474841339419*t**24 + 1065301721692202256*t**23 + 2590965176032226004*t**22 +    5679139393727958880*t**21 + 11197539488052244692*t**20 + 19809141347016503568*t**19    + 31337344703196264683*t**18 + 44145543697033615230*t**17 +    55090950408765652881*t**16 + 60516047590917322480*t**15 + 58058828547614202153*t**14    + 48186407491493228562*t**13 + 34193160536571257393*t**12 +    20445602003818379424*t**11 + 10116440549699970582*t**10 + 4048631186832794128*t**9 +    1273256824201491894*t**8 + 303519895786350096*t**7 + 52538628166252339*t**6 +    6321448897556694*t**5 + 513273187133223*t**4 + 27378446805968*t**3 +    916904893215*t**2 + 17467984794*t + 144422161)/(t**42 + 78*t**41 + 3003*t**40 +    75920*t**39 + 1414959*t**38 + 20695506*t**37 + 246960727*t**36 + 2468218896*t**35 +    21049407990*t**34 + 155295575344*t**33 + 1001454266670*t**32 + 5689547290752*t**31 +    28649170403321*t**30 + 128443250574150*t**29 + 514438663036389*t**28 +    1845008125106224*t**27 + 5933895362321121*t**26 + 17125495350262650*t**25 +    44346853027068215*t**24 + 102951022273700016*t**23 + 213919716149439924*t**22 +    396901875045903136*t**21 + 655411304131422852*t**20 + 959204141060458800*t**19 +    1237542066720323195*t**18 + 1398213766917750042*t**17 + 1372033638327823605*t**16 +    1157344876243332400*t**15 + 828397238649616545*t**14 + 494847675800288070*t**13 +    241330300300762205*t**12 + 93195177370728576*t**11 + 27216183343932726*t**10 +    5548446154604080*t**9 + 655832005477038*t**8 + 13945989607056*t**7 -    5655555970805*t**6 - 193929669102*t**5 + 43376116875*t**4 - 2043712624*t**3 +    45287151*t**2 - 497874*t + 2197)
     3920
     3921
     3922#isogeny_kernel_polys = dict([(5,Psi5),(7,Psi7),(13,Psi13)])
    39233923
    39243924def isogenies_prime_degree_genus_0(E, l=None):
    39253925    """