>>>>> "Edwin" == Edwin Woollett <woollett at charter.net> writes:
Edwin> How can I speed up the use of bromberg?
Edwin> I have a simple non-singular integrand
Edwin> depending only on sqrt, over a small domain, which
Edwin> quad_qags does (using the defaults) in almost no time.
Edwin> I am only asking bromberg for 10 digit precision.
Edwin> The integrand is 1/(1 + sqrt(x), the domain is [x,0,1].
Edwin> If I don't translate the function, bromberg takes about
Edwin> 9 minutes.
Edwin> If I translate the function, bromberg takes about 8.5 min.
Edwin> If I both compile brmbrg.lisp and the function (see below)
Edwin> bromberg takes about 7.8 min.
Try the version of bromberg at
http://common-lisp.net/~rtoy/maxima/bromberg.lisp
It's much faster for me than the default bromberg, even without
compiling it. (It's basically a direct copy of romberg, but slightly
modified to use maxima's bigfloat class.)
Ray