On Oct 5, 2009, at 1:26 PM, mark mcclure wrote: > > On Oct 5, 1:19 pm, Robert Bradshaw <rober...@math.washington.edu> > wrote: >> >> you might want to call these functions in SciPy directly as C >> functions. (You'd have to look up the SciPy headers to see what to >> call them.) > > That was my first thought, in fact. I'm not certain, but it appears > that > SciPy is ultimately calling Fortran code, which I'm not eager to delve > into. > > Your suggestions did provide a further modest speedup but not another > 50%. I really think that the airy function is the likely culprit.
Try replacing the airy function with both math.sqrt, and math.h's sqrt, and see what the timings are, which for benchmarking purposes should give you a good idea if it's really the airy function. If it is, I doubt we're going to be writing a faster one ourselves... - Robert --~--~---------~--~----~------------~-------~--~----~ 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 -~----------~----~----~----~------~----~------~--~---
