Sorry for the delay - I didn't see your original reply
in the Mailing List until today - been rather
over-committed...
I am getting a Lisp error message when I use your version:
--------------------------------------------
Maxima 5.28.0-2 http://maxima.sourceforge.net
using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
2012-11-19
(%i1) load("c:/work2/bromberg.lisp");
(%o1) "c:/work2/bromberg.lisp"
(%i2) [brombergtol,brombergabs,brombergit,brombergmin];
(%o2) [1.0E-4,0.0,11,0]
(%i3) (brombergtol:0.0b0,brombergit:100)$
(%i4) brombergabs:1.0b-10;
(%o4) 1.0b-10
(%i5) g(x):= 1/(1+sqrt(x))$
(%i6) bromberg(g,0,1);
Maxima encountered a Lisp error:
Error in BIGFLOAT::TWO-ARG-< [or a callee]: No matching method for the
generic-function #<compiled-closure BIGFLOAT::TWO-ARG-<>,
when called with arguments (((BIGFLOAT SIMP 56) 61897001964269014 -33)
+3.407640327429691b-2).
Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
(%i7)
------------------------------------
Ted
======================
=====================
----- Original Message -----
From: Raymond Toy
To: Edwin Woollett
Sent: Friday, November 16, 2012 7:13 PM
Subject: Re: how can I speed up bromberg?
>>>>> "Raymond" == Raymond Toy <toy.raymond at gmail.com> writes:
>>>>> "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.
Raymond> Try the version of bromberg at
Raymond> http://common-lisp.net/~rtoy/maxima/bromberg.lisp
Did you try this out? Does it work better?
Ray