The standard code to compute elliptic integrals in sage appears to be _very_ slow numerically. To compute the complete elliptic integral on 100 points in the unit interval takes about 16 seconds for me.
sage: %time pts=[ elliptic_ec(m) for m in srange(0,1,0.01,include_endpoint=True)] CPU times: user 0.59 s, sys: 0.06 s, total: 0.65 s Wall time: 16.44 s By contrast, Octave computes these in less than 0.2 seconds on the same machine octave:37> tic; [k,e]=ellipke([0:0.01:1]'); toc Elapsed time is 0.01437 seconds. Am I doing something wrong here? I realise sage isn't a numerical system, but this code appears to be quite slow. -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
