Skip to content

Instantly share code, notes, and snippets.

@asmeurer

asmeurer/gist:993982

Created May 26, 2011 20:18
Show Gist options
  • Save asmeurer/993982 to your computer and use it in GitHub Desktop.
Save asmeurer/993982 to your computer and use it in GitHub Desktop.
Download ZIP
poly vs. Poly
Raw
gistfile1.txt
In [25]: %timeit Poly(p)
1 loops, best of 3: 735 ms per loop
In [26]: %timeit poly(p)
10 loops, best of 3: 97.2 ms per loop
In [27]: print p
_z0*(-6*_z0 + 6)**2*(-3*_z0 + 3)*(6*_z0 - 6)*(-8*_z0**2 + 16*_z0 - 8)*(8*_z0**2 - 16*_z0 + 8)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280)**2*(-12800*2**(1/3)*_z0**2*(-634*_z0 - 104)*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640) + 7680*2**(1/3)*_z0**2*(-23*_z0 - 4)*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**3*(144*_z0**2 - 288*_z0 + 144)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 19200*2**(1/3)*_z0**2*(6*_z0 - 6)**3*(8*_z0 - 8)*(-4*_z0**2 + 8*_z0 - 4)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 10*2**(1/3)*_z0**2*(6*_z0 - 6)**3*(-317*_z0**2 - 104*_z0 + 16)*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(-4*_z0**3 + 12*_z0**2 - 12*_z0 + 4)*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 883200*2**(1/3)*_z0**2*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) + 76800*2**(1/3)*_z0*(-23*_z0 - 4)*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 12800*2**(1/3)*_z0*(6*_z0 - 6)**3*(-317*_z0**2 - 104*_z0 + 16)*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 38400*2**(1/3)*_z0*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**2*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 230400*2**(1/3)*_z0*(6*_z0 - 6)*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) + 38400*2**(1/3)*(6*_z0 - 6)**2*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280)) - 12800*_z0*(-3*_z0 + 3)*(6*_z0 - 6)**4*(-8*_z0**2 + 16*_z0 - 8)*(-4*_z0**2 + 8*_z0 - 4)**3*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**3*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280)**2*(-2**(1/3)*_z0*(-6*_z0 + 6)*(-_z0 + 4)*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 2**(1/3)*_z0*(-6*_z0 + 6)*(-5*_z0**2 + 24*_z0 + 8)*(8*_z0**2 - 16*_z0 + 8)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 2**(1/3)*_z0*(-6*_z0 + 6)*(8*_z0**2 - 16*_z0 + 8)*(-55*_z0**3 + 276*_z0**2 + 216*_z0 - 32)*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80) + 3*2**(1/3)*_z0*(8*_z0**2 - 16*_z0 + 8)*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) + 2**(1/3)*(-6*_z0 + 6)*(8*_z0**2 - 16*_z0 + 8)*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280)) + 12800*(-6*_z0 + 6)**2*(2*_z0 - 4)*(6*_z0 - 6)**3*(-4*_z0**2 + 8*_z0 - 4)**3*(8*_z0**2 - 16*_z0 + 8)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)**2*(-48*_z0**3 + 144*_z0**2 - 144*_z0 + 48)**2*(80*_z0**3 - 240*_z0**2 + 240*_z0 - 80)*(3*2**(1/3)*_z0**2*(-23*_z0 - 4)*(6*_z0 - 6)*(-8*_z0**2 + 16*_z0 - 8)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) - 2**(1/3)*_z0**2*(6*_z0 - 6)*(-317*_z0**2 - 104*_z0 + 16)*(-8*_z0**2 + 16*_z0 - 8)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240) - 3*2**(1/3)*_z0**2*(6*_z0 - 6)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280) + 3*2**(1/3)*_z0*(-8*_z0**2 + 16*_z0 - 8)*(-240*_z0**3 + 720*_z0**2 - 720*_z0 + 240)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280))*(_z0**4 - 4*_z0**3 + 6*_z0**2 - 4*_z0 + 1)**2*(640*_z0**4 - 2560*_z0**3 + 3840*_z0**2 - 2560*_z0 + 640)*(1280*_z0**4 - 5120*_z0**3 + 7680*_z0**2 - 5120*_z0 + 1280)**2
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment