Post Reply 
WP-34S Integration speed
02-21-2014, 03:08 AM
Post: #1
WP-34S Integration speed
Is there a way to balance accuracy with speed of integration on the 34S as with HP's that can integrate? The HP's give the user the ability to require, say, 4 significant digits of accuracy vs. the full 12 digits in order to arrive at a result in a significantly shorter amount of time. Can this be accomplished on the wonderful 34S?
Find all posts by this user
Quote this message in a reply
02-22-2014, 11:44 AM
Post: #2
RE: WP-34S Integration speed
The integration routine is modeled after the IG routine from the HP 41 PPC ROM. It uses the display accuracy to determine the number of iterations by means of the x≈ Y comparison command. You can improve the integration speed by setting a different display mode such as FIX 4 or SCI 4.

Marcus von Cube
Wehrheim, Germany
http://www.mvcsys.de
http://wp34s.sf.net
http://mvcsys.de/doc/basic-compare.html
Find all posts by this user
Quote this message in a reply
02-22-2014, 12:35 PM
Post: #3
RE: WP-34S Integration speed
However, an integral that is approximately zero will hit a snag and will continue iterating long after it should.

Nobody like the Gauss-Kronrod quadrature I initially implemented which was fast and good for smooth functions but in no way adaptive.


- Pauli
Find all posts by this user
Quote this message in a reply
02-23-2014, 02:02 PM
Post: #4
RE: WP-34S Integration speed
(02-22-2014 12:35 PM)Paul Dale Wrote:  However, an integral that is approximately zero will hit a snag and will continue iterating long after it should.

Pauli, I remember we had a discussion on this very topic some time ago, and I suggested a possible fix: compare the current value of \(\int f(x)\) with an approximation of \(\int |f(x)|\). It worked well with the examples we discussed then. Details can be found here, cf. message 6 ff.

This change was supposed to get implemented in v. 3460 resp. 3462. So I thought this issue had been resolved.

Dieter
Find all posts by this user
Quote this message in a reply
02-23-2014, 09:42 PM
Post: #5
RE: WP-34S Integration speed
Sorry, Dieter is correct. We implemented a fix for the integral approaching zero. I'd just forgotten about it.


- Pauli
Find all posts by this user
Quote this message in a reply
02-27-2014, 12:14 AM
Post: #6
RE: WP-34S Integration speed
This thread was moved from the General forum by request.
Visit this user's website Find all posts by this user
Quote this message in a reply
02-27-2014, 06:43 AM (This post was last modified: 02-27-2014 06:44 AM by walter b.)
Post: #7
RE: WP-34S Integration speed
χαιρε, Δρακων. Now it's where ο δημος will look for it.

d:-)
Find all posts by this user
Quote this message in a reply
02-27-2014, 02:28 PM (This post was last modified: 02-27-2014 02:29 PM by Dieter.)
Post: #8
RE: WP-34S Integration speed
(02-23-2014 09:42 PM)Paul Dale Wrote:  We implemented a fix for the integral approaching zero.

I see the fix got implemented, but one final improvement I posted back then still is missing (cf. last message in the mentioned thread). It's one single line that should be added:

Code:
        ...
        FC?C f_final        /* final iteration done? */
        x[approx]? Y        /* or convergence? */
             JMP int_done
        RCL/ r_SkAbs
        RCL/ r_deltau
        RCL/ r_ba4        <==  add this line
        ABS
        ...

Maybe it can be added with the next update.

Dieter
Find all posts by this user
Quote this message in a reply
02-27-2014, 09:28 PM
Post: #9
RE: WP-34S Integration speed
Dieter, I must have missed that final improvement. Added now but not built.

- Pauli
Find all posts by this user
Quote this message in a reply
03-02-2014, 01:15 PM
Post: #10
RE: WP-34S Integration speed
(02-27-2014 09:28 PM)Paul Dale Wrote:  Dieter, I must have missed that final improvement. Added now but not built.
My Job! Build 3530 is Committed to SVN... Smile

Marcus von Cube
Wehrheim, Germany
http://www.mvcsys.de
http://wp34s.sf.net
http://mvcsys.de/doc/basic-compare.html
Find all posts by this user
Quote this message in a reply
Post Reply 




User(s) browsing this thread: 1 Guest(s)