Annotation of OpenXM_contrib/pari/src/test/32/elliptic, Revision 1.1
1.1 ! maekawa 1: realprecision = 38 significant digits
! 2: echo = 1 (on)
! 3: ? ellinit([0,0,0,-1,0])
! 4: [0, 0, 0, -1, 0, 0, -2, 0, -1, 48, 0, 64, 1728, [1.0000000000000000000000000
! 5: 000000000000, 0.E-38, -1.0000000000000000000000000000000000000]~, 2.62205755
! 6: 42921198104648395898911194136, 2.6220575542921198104648395898911194136*I, -0
! 7: .59907011736779610371996124614016193910, -1.79721035210338831115988373842048
! 8: 58173*I, 6.8751858180203728274900957798105571979]
! 9: ? ellinit([0,0,0,-17,0],1)
! 10: [0, 0, 0, -17, 0, 0, -34, 0, -289, 816, 0, 314432, 1728]
! 11: ? ellj(I)
! 12: 1728.0000000000000000000000000000000000 + 0.E-45*I
! 13: ? ellsub(ellinit([0,0,0,-17,0]),[-1,4],[-4,2])
! 14: [9, -24]
! 15: ? acurve=ellinit([0,0,1,-1,0])
! 16: [0, 0, 1, -1, 0, 0, -2, 1, -1, 48, -216, 37, 110592/37, [0.83756543528332303
! 17: 544481089907503024040, 0.26959443640544455826293795134926000404, -1.10715987
! 18: 16887675937077488504242902444]~, 2.9934586462319596298320099794525081778, 2.
! 19: 4513893819867900608542248318665252253*I, -0.47131927795681147588259389708033
! 20: 769964, -1.4354565186686843187232088566788165076*I, 7.3381327407895767390707
! 21: 210033323055881]
! 22: ? apoint=[2,2]
! 23: [2, 2]
! 24: ? elladd(acurve,apoint,apoint)
! 25: [21/25, -56/125]
! 26: ? ellak(acurve,1000000007)
! 27: 43800
! 28: ? ellan(acurve,100)
! 29: [1, -2, -3, 2, -2, 6, -1, 0, 6, 4, -5, -6, -2, 2, 6, -4, 0, -12, 0, -4, 3, 1
! 30: 0, 2, 0, -1, 4, -9, -2, 6, -12, -4, 8, 15, 0, 2, 12, -1, 0, 6, 0, -9, -6, 2,
! 31: -10, -12, -4, -9, 12, -6, 2, 0, -4, 1, 18, 10, 0, 0, -12, 8, 12, -8, 8, -6,
! 32: -8, 4, -30, 8, 0, -6, -4, 9, 0, -1, 2, 3, 0, 5, -12, 4, 8, 9, 18, -15, 6, 0
! 33: , -4, -18, 0, 4, 24, 2, 4, 12, 18, 0, -24, 4, 12, -30, -2]
! 34: ? ellap(acurve,10007)
! 35: 66
! 36: ? ellap(acurve,10007,1)
! 37: 66
! 38: ? deu=direuler(p=2,100,1/(1-ellap(acurve,p)*x+if(acurve[12]%p,p,0)*x^2))
! 39: [1, -2, -3, 2, -2, 6, -1, 0, 6, 4, -5, -6, -2, 2, 6, -4, 0, -12, 0, -4, 3, 1
! 40: 0, 2, 0, -1, 4, -9, -2, 6, -12, -4, 8, 15, 0, 2, 12, -1, 0, 6, 0, -9, -6, 2,
! 41: -10, -12, -4, -9, 12, -6, 2, 0, -4, 1, 18, 10, 0, 0, -12, 8, 12, -8, 8, -6,
! 42: -8, 4, -30, 8, 0, -6, -4, 9, 0, -1, 2, 3, 0, 5, -12, 4, 8, 9, 18, -15, 6, 0
! 43: , -4, -18, 0, 4, 24, 2, 4, 12, 18, 0, -24, 4, 12, -30, -2]
! 44: ? ellan(acurve,100)==deu
! 45: 1
! 46: ? ellisoncurve(acurve,apoint)
! 47: 1
! 48: ? acurve=ellchangecurve(acurve,[-1,1,2,3])
! 49: [-4, -1, -7, -12, -12, 12, 4, 1, -1, 48, -216, 37, 110592/37, [-0.1624345647
! 50: 1667696455518910092496975959, -0.73040556359455544173706204865073999595, -2.
! 51: 1071598716887675937077488504242902444]~, -2.99345864623195962983200997945250
! 52: 81778, -2.4513893819867900608542248318665252253*I, 0.47131927795681147588259
! 53: 389708033769964, 1.4354565186686843187232088566788165076*I, 7.33813274078957
! 54: 67390707210033323055881]
! 55: ? apoint=ellchangepoint(apoint,[-1,1,2,3])
! 56: [1, 3]
! 57: ? ellisoncurve(acurve,apoint)
! 58: 1
! 59: ? ellglobalred(acurve)
! 60: [37, [1, -1, 2, 2], 1]
! 61: ? ellheight(acurve,apoint)
! 62: 0.40889126591975072188708879805553617287
! 63: ? ellheight(acurve,apoint,1)
! 64: 0.40889126591975072188708879805553617296
! 65: ? ellordinate(acurve,1)
! 66: [8, 3]
! 67: ? ellztopoint(acurve,ellpointtoz(acurve,apoint))
! 68: [0.99999999999999999999999999999999999990 + 0.E-38*I, 2.99999999999999999999
! 69: 99999999999999998 + 0.E-38*I]
! 70: ? ellpow(acurve,apoint,10)
! 71: [-28919032218753260057646013785951999/292736325329248127651484680640160000,
! 72: 478051489392386968218136375373985436596569736643531551/158385319626308443937
! 73: 475969221994173751192384064000000]
! 74: ? ellwp(acurve)
! 75: x^-2 + 1/5*x^2 - 1/28*x^4 + 1/75*x^6 - 3/1540*x^8 + 1943/3822000*x^10 - 1/11
! 76: 550*x^12 + 193/10510500*x^14 - 1269/392392000*x^16 + 21859/34684650000*x^18
! 77: - 1087/9669660000*x^20 + 22179331/1060517858400000*x^22 - 463/124093970000*x
! 78: ^24 + 47495017/70175140035000000*x^26 - 34997918161/291117454720092000000*x^
! 79: 28 + O(x^30)
! 80: ? ellpointtoz(acurve,apoint)
! 81: 0.72491221490962306778878739838332384646 + 0.E-58*I
! 82: ? q*Ser(ellan(acurve,100),q)
! 83: q - 2*q^2 - 3*q^3 + 2*q^4 - 2*q^5 + 6*q^6 - q^7 + 6*q^9 + 4*q^10 - 5*q^11 -
! 84: 6*q^12 - 2*q^13 + 2*q^14 + 6*q^15 - 4*q^16 - 12*q^18 - 4*q^20 + 3*q^21 + 10*
! 85: q^22 + 2*q^23 - q^25 + 4*q^26 - 9*q^27 - 2*q^28 + 6*q^29 - 12*q^30 - 4*q^31
! 86: + 8*q^32 + 15*q^33 + 2*q^35 + 12*q^36 - q^37 + 6*q^39 - 9*q^41 - 6*q^42 + 2*
! 87: q^43 - 10*q^44 - 12*q^45 - 4*q^46 - 9*q^47 + 12*q^48 - 6*q^49 + 2*q^50 - 4*q
! 88: ^52 + q^53 + 18*q^54 + 10*q^55 - 12*q^58 + 8*q^59 + 12*q^60 - 8*q^61 + 8*q^6
! 89: 2 - 6*q^63 - 8*q^64 + 4*q^65 - 30*q^66 + 8*q^67 - 6*q^69 - 4*q^70 + 9*q^71 -
! 90: q^73 + 2*q^74 + 3*q^75 + 5*q^77 - 12*q^78 + 4*q^79 + 8*q^80 + 9*q^81 + 18*q
! 91: ^82 - 15*q^83 + 6*q^84 - 4*q^86 - 18*q^87 + 4*q^89 + 24*q^90 + 2*q^91 + 4*q^
! 92: 92 + 12*q^93 + 18*q^94 - 24*q^96 + 4*q^97 + 12*q^98 - 30*q^99 - 2*q^100 + O(
! 93: q^101)
! 94: ? bcurve=ellinit([0,0,0,-3,0])
! 95: [0, 0, 0, -3, 0, 0, -6, 0, -9, 144, 0, 1728, 1728, [1.7320508075688772935274
! 96: 463415058723669, 0.E-38, -1.7320508075688772935274463415058723669]~, 1.99233
! 97: 28995834907073368080310227454215, 1.9923328995834907073368080310227454215*I,
! 98: -0.78842061340415606811560792095228873037, -2.36526184021246820434682376285
! 99: 68661911*I, 3.9693903827627596663162680332564652027]
! 100: ? elllocalred(bcurve,2)
! 101: [6, 2, [1, 1, 1, 0], 1]
! 102: ? elltaniyama(bcurve)
! 103: [x^-2 - x^2 + 3*x^6 - 2*x^10 + 7*x^14 + O(x^15), -x^-3 + 3*x - 3*x^5 + 8*x^9
! 104: - 9*x^13 + O(x^14)]
! 105: ? ccurve=ellinit([0,0,-1,-1,0])
! 106: [0, 0, -1, -1, 0, 0, -2, 1, -1, 48, -216, 37, 110592/37, [0.8375654352833230
! 107: 3544481089907503024040, 0.26959443640544455826293795134926000404, -1.1071598
! 108: 716887675937077488504242902444]~, 2.9934586462319596298320099794525081778, 2
! 109: .4513893819867900608542248318665252253*I, -0.4713192779568114758825938970803
! 110: 3769964, -1.4354565186686843187232088566788165076*I, 7.338132740789576739070
! 111: 7210033323055881]
! 112: ? l=elllseries(ccurve,2)
! 113: 0.38157540826071121129371040958008663666
! 114: ? elllseries(ccurve,2,1.2)-l
! 115: 1.3224311443045735350000000000000000000 E-38
! 116: ? tcurve=ellinit([1,0,1,-19,26]);
! 117: ? ellorder(tcurve,[1,2])
! 118: 6
! 119: ? elltors(tcurve)
! 120: [12, [6, 2], [[-2, 8], [3, -2]]]
! 121: ? mcurve=ellinit([0,0,0,-17,0])
! 122: [0, 0, 0, -17, 0, 0, -34, 0, -289, 816, 0, 314432, 1728, [4.1231056256176605
! 123: 498214098559740770251, 0.E-38, -4.1231056256176605498214098559740770251]~, 1
! 124: .2913084409290072207105564235857096009, 1.2913084409290072207105564235857096
! 125: 009*I, -1.2164377440798088266474269946818791934, -3.649313232239426479942280
! 126: 9840456375802*I, 1.6674774896145033307120230298772362381]
! 127: ? mpoints=[[-1,4],[-4,2]]~
! 128: [[-1, 4], [-4, 2]]~
! 129: ? mhbi=ellbil(mcurve,mpoints,[9,24])
! 130: [-0.72448571035980184146215805860545027439, 1.307328627832055544492943428892
! 131: 1943055]~
! 132: ? ma=ellheightmatrix(mcurve,mpoints)
! 133:
! 134: [1.1721830987006970106016415566698834134 0.447697388340895169139483498064433
! 135: 13906]
! 136:
! 137: [0.44769738834089516913948349806443313906 1.75502601617295071363242692695662
! 138: 74446]
! 139:
! 140: ? matsolve(ma,mhbi)
! 141: [-1.0000000000000000000000000000000000000, 1.0000000000000000000000000000000
! 142: 000000]~
! 143: ? cmcurve=ellinit([0,-3/4,0,-2,-1])
! 144: [0, -3/4, 0, -2, -1, -3, -4, -4, -1, 105, 1323, -343, -3375, [2.000000000000
! 145: 0000000000000000000000000, -0.62500000000000000000000000000000000000 + 0.330
! 146: 71891388307382381270196920490755321*I, -0.6250000000000000000000000000000000
! 147: 0000 - 0.33071891388307382381270196920490755321*I]~, 1.933311705616811546733
! 148: 0768390298137310, 0.96665585280840577336653841951490686552 + 2.5575309899160
! 149: 994790492257969408742846*I, -0.85584863309985585256414907906657726119 + 4.98
! 150: 60612787306308740000000000000000000 E-39*I, -0.42792431654992792628207453953
! 151: 328863060 - 2.7571612171661472068745043203629203178*I, 4.9445046002825467364
! 152: 981969681843776438]
! 153: ? ellpow(cmcurve,[x,y],quadgen(-7))
! 154: [((-2 + 3*w)*x^2 + (6 - w))/((-2 - 5*w)*x + (-4 - 2*w)), ((34 - 11*w)*y*x^2
! 155: + (40 - 28*w)*y*x + (22 + 23*w)*y)/((-90 - w)*x^2 + (-136 + 44*w)*x + (-40 +
! 156: 28*w))]
! 157: ? \p96
! 158: realprecision = 96 significant digits
! 159: ? precision(cmcurve)
! 160: 38
! 161: ? getheap
! 162: [56, 4651]
! 163: ? print("Total time spent: ",gettime);
! 164: Total time spent: 813
! 165: ? \q
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to elliptic CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / pari / src / test / 32