[WP34s] Bug or feature?



#23

I just realized that the 34s rounds displayed values in a way that is different from classic HPs (at least the ones I know of). If in FIX mode the value in X is sufficiently low, the display switches to SCI (or ENG) notation. That's fine. But the thresholds are different:

  FIX 2                Display 
34s HP
80 [1/x] => 0,01 0,01
2 [/] => 6,25 -3 0,01
2 [/] => 3,13 -3 3,13 -3
 
100,1 [1/x] => 9,99 -3 0,01
100,01 [1/x] => 1,00 -2 0,01
The 34s switches as soon as the value itself falls below 0,01 - while HPs own devices do so as soon as the rounded value passes this limit (i.e. below 0,005). Is this behaviour intentional? At least a subsequent ROUND seems to return the same results in both cases.

While we're at it: The way SLV works still does not match the manual. The latter claims the user interface is the same as on the 15C, but it isn't. On the 34s, the display setting affects the results, while on the 15C (and 34C and 35s and...) the function always returns the best possible result. Walter, you recently asked for an example and I posted one, but there was no reaction.

Dieter


#24

Dieter,

regarding SLV it depends on Pauli: if we want SLV to operate as in the HP-15C (as advertized) then apparently the SW has to be adapted. Else we have to specify how it operates ;-)

d:-)


#25

Honestly - if we really wanted the 34s to behave exactly as the 15C, this would require the same code of the Solve function. Many examples presented in the 34C and 15C manual cannot be duplicated on the 34s, simply because it uses a different algorithm. Finding poles or local extrema can be done with HP's original SOLVE implementation, but this will not work with SLV on the 34s.

That's why I would suggest to remove the 15C reference from the 34s manual. The only thing both implementations have in common is the stack layout on entry and on exit. And even this is not exactly the same in both cases.

Dieter


#26

Quote:
That's why I would suggest to remove the 15C reference from the 34s manual.

Done. Will be published in due course.

d:-)

#27

The FIX to SCI change over is a troublesome part of probably the hairiest piece of code in the device, I'm not surprised something like this did appear. I've long held a suspicion that HP's display code contained similar quirks that became accepted as normal behaviour.

The code in question is lines 945 - 976 of trunk/display.c The rounding step would have to be performed before the mode switching checks, but it certainly isn't a simple matter of rearranging the code. The question is how important is make a change like this? An awful lot of testing would have to follow -- there are many edge and obscure cases the display code must deal with.


As for the SLV command. Making it not depend on display settings is easy. Change the five x[approx]? tests in trunk/xrom/solve.wp34s to x=? ones instead. Is this desirable behaviour or not?


- Pauli


#28

Quote:
As for the SLV command. Making it not depend on display settings is easy. Change the five x[approx]? tests in trunk/xrom/solve.wp34s to x=? ones instead. Is this desirable behaviour or not?

No, at least not for me. I find it ok that the user is able to determine the solver accuracy with the display setting - it's not always necessary to get a solution with full precision, especially as the solving process may not be the fastest for some complicated expressions.

Franz


#29

Add a simple round command as the final step in your user routine, and you're done. Easy as that.

But there's a much more important point to this: what exactly are we talking about here? Is it the accuracy of the returned root or is it the tolerance of the function result? In other words: Would FIX 2 in a format-dependent SLV routine return a result with at least two valid decimals, or would this return a root for which the function is approximately zero (when rounded to two decimal places)? That's quite a difference.

Dieter

Edited: 30 Dec 2012, 5:27 p.m.


#30

Quote:
Add a simple round command as the final step in your user routine, and you're done. Easy as that.

It seems you didn't understand. Rounding _after_ solving won't make the solver faster at all. Why should I calculate 12 or 16 digits when I only need 3?

Franz


#31

Quote:
It seems you didn't understand. Rounding _after_ solving won't make the solver faster at all. Why should I calculate 12 or 16 digits when I only need 3?

Oh, Franz.... #-)

The rounding is done within the user routine, not after SLV has returned a full-accuracy result.

Determine the root of x^2-5 where f(x)...
...is exactly 0 ...is zero if rounded to 3 digits
 
001 LBL A 001 LBL A
002 x^2 002 x^2
003 5 003 5
004 - 004 -
005 RTN 005 RDP 03
006 RTN
Got the idea? This method was already mentioned in the 34C and 15C manual. And believe me, on those machines it really made sense - they were slow. Really slow.

Dieter


#32

Quote:
Oh, Franz.... #-)


Ok, ok, I've overlooked your part with the "user routine". ;-)

But this is something you've mentioned in your previous post: does SLV stop when dx~0 or f(x)~0?

Using ROUND in the user routine would influence the accuracy of f(x) (and dx only indirectly), but I guess the WP34s SLV makes its decision (when to terminate) wrt. dx-value.

Franz


#33

Quote:
...but I guess the WP34s SLV makes its decision (when to terminate) wrt. dx-value.

Take a look at the SLV code (in solve.wp34s) and no guessing is required. ;-)

Actually there are three tests that check if the function result (!) is approximately zero (applies to f(a), f(b) and f(new_estimate). And finally SLV checks whether the two last approximations agree within 5 ULP.

So, yes, the current SLV implementation seems to do its display-setting-dependent tests based on f(x). Yes, it also checks if the last two estimates are nearly identical (to all digits except the last), but this is not affected by the display setting.

In other words: the rounding done by SLV can just as well be done within the user's function routine. Simply add one of the round commands as the final step.

Dieter

Edited: 30 Dec 2012, 7:17 p.m. after one or more responses were posted


#34

The tests for f(x) being almost zero can be done via a round in the user's code. In fact, they would be better done there because they only do something useful in FIX mode. I'll make this change.

The termination code avoids rounded compares due to some problems there which resulted in bad termination conditions. I don't recollect the specifics anymore.

That just leaves two rounded comparisons. These are used at the start if a singular guess is given -- a search interval needs to be created from such. I suspect any interval creation method would be viable here. Will have a think about what might work nicely.


- Pauli


#35

Is there a chance after that SLV stops giving pole rather than root as an answer in all cases? I haven't checked it lately (are there any changes to SLV?) but I had couple of examples of doing exactly that. I know ... 'everyone is free to develop their own solver' etc., just asking. I'm not capable of doing it.
Cheers,


#36

Is there a chance asking more specific? I.e. unvealing your 'couple of examples'? That would support the quest ... TIA

d:-)


#37

Yes, there is, I'll post them here first thing in the New Year which is now less than 6 hours away. I was actually hoping PD would remember so I asked him.

All the best for the new 2013 to all!

#38

It certainly is possible to make an attempt to detect poles, it just takes time to code and test and a bit of inspiration to escape the influence of the pole.


- Pauli

#39

By the way - there also is a way to make SLV return the root (and not f(x)) with a given accuracy. The idea here is to compare the current estimate with the previous one:

   001 LBL A
002 FILL // can be omitted as SLV fills the stack anyway
003 X<> 00
004 -
005 RDP 02 // |previous - current estimate| less than 0,005?
006 x=0?
007 RTN // then return zero, make SLV believe a root was found
008 DROP
009 x^2 // usual function code follows here
010 5
011 -
012 RTN
Provide two sufficiently different guesses and store the larger one in R00.
   2 [ENTER] 3  [STO] 00
Then start SLV as usual. Even with no accuracy limitiations, e.g. in ALL mode, SLV now returns a result within the specified tolerance:
   [ALL] 03
[SLV] [A]
=> 2,24000661691
The exact solution is sqrt(5) = 2,236..., so the error in X is less than 0,005 - just as specified.

Dieter


Edited: 31 Dec 2012, 9:16 a.m.


#40

I think I'll have to add a big CAVEAT here. The method I suggested works well if the solver produces a strictly convergent series of estimates, i.e. the current estimate is as least as accurate as the difference between the last two:

  |x_new - x_exact|  <  |x_new - x_old|
This is not (always) true for the 34s solver, and it does not apply to various other algorithms eiher.

Here's an example. Let's consider the following simple function:

  ln(x) - 1  =  0
It has a root at x = e. With 2 and 5 as initial guesses The 34s solver produces the following sequence of estimates:
  estimate     abs. error
-------------------------
2,00000000000 -7,18E-1
5,00000000000 2,28E 0
3,00465761171 2,86E-1 // here the estimate is 0,286 off
2,99797706383 2,80E-1 // but the solver adjusts it by only 0,0068
2,49898853191 -2,19E-1 // which may lead to the conclusion that
2,71681463363 -1,47E-3 // the estimate is already that exact
2,85739584873 1,39E-1
2,71828330193 1,47E-6
2,71828215475 3,26E-7
2,71754839419 -7,33E-4
2,71828182846 -9,00E-15 // here the estimate is almost exact
2,71828199161 1,63E-7 //
2,71828182846 0,00E+0 // and here the estimate evnb is dead on
2,71828191003 8,16E-8 // but due to roundoff the function result
2,71828186925 4,08E-8 // is not exactly zero and the solver again
2,71828184885 2,04E-8 // moves away from the exact solution.
2,71828183866 1,02E-8 // Starting with an estimate that's 8,16E-8 off
2,71828183356 5,10E-9 // it now switches to bisection
2,71828183101 2,55E-9 // so that it takes more than 20 iterations
2,71828182973 1,27E-9 // to find the final digits
2,71828182910 6,37E-10
2,71828182878 3,19E-10
2,71828182862 1,59E-10
2,71828182854 7,97E-11
2,71828182850 3,98E-11
2,71828182848 1,99E-11
2,71828182847 9,96E-12
2,71828182846 4,98E-12
2,71828182846 2,49E-12
2,71828182846 1,25E-12
2,71828182846 6,23E-13
2,71828182846 3,11E-13
2,71828182846 1,55E-13
2,71828182846 7,70E-14
2,71828182846 3,90E-14
2,71828182846 1,90E-14
2,71828182846 9,00E-15
2,71828182846 5,00E-15
2,71828182846 3,00E-15
So the current solver eventually comes to a solution, but it may take a while (which is fine if on the other hand a solution is reliably found in more complex cases). Line 3...5 show why the proposed method of comparing the current estimate with the previous one does not work here.

Here is the sequence returned with 1 and 5 as initial guesses:

  estimate     abs. error
-------------------------
1,00000000000 -1,72E+0
5,00000000000 2,28E+0
3,48533973824 7,67E-1
2,24266986912 -4,76E-1
2,68993254763 -2,83E-2
3,08763614293 3,69E-1 // here the solver moves away from the solution
2,71846362445 1,82E-4 // ...comes closer
2,71839999019 1,18E-4
2,70416626891 -1,41E-2 // ...moves away again
2,71828182711 -1,35E-9 // ...comes very close
2,71834090865 5,91E-5 // ...but then moves away again
2,71828182846 0,00E+0 // Here x is exact but f(x) is not zero
2,71828182846 2,00E-15 // so the solver prefers this solution
Finally I tried the same with the 35s and 34C solvers.
         HP 35s                         HP-34C
estimate abs. error estimate abs. error
-------------------------------------------------------
1,00000000000 -1,72E+0 1,000000000 -1,72E+0
5,00000000000 2,28E+0 5,000000000 2,28E+0
3,48533973824 7,67E-1 3,485339739 7,67E-1
2,75051645695 3,22E-2 2,442055911 -2,76E-1
2,71892340384 6,42E-4 2,756337820 3,81E-2
1,85946170192 -8,59E-1 2,720245459 1,96E-3
2,71814296203 -1,39E-4 2,718268115 -1,37E-5
2,71828182743 1,03E-9 2,718281833 5,00E-9
2,71828182846 0,00E+0 2,718281828 0,00E+0
2,71828182846 0,00E+0 2,718281829 1,00E-9
Which leads to the question if there are some deeper insights on how HP Solve (resp. a particular implementation) works in detail. I found an interesting discussion here in the forum back in 2003, cf. this thread.

Dieter


#41

Hallo Dieter,

Quote:
I found an interesting discussion here in the forum back in 2003, cf. this thread.

Thanks for making us aware of that old thread. Interesting read.

d:-)


#42

Yes, the HP solver (here: the version used in the 34C and 15C) is known to produce reliable results, so maybe we can find out more about the way it works and use this for the 34s.

Anyway, I always thought the 34s solver used some kind of implementation of Brent's algorithm. But in this case it would not require 30 to 40 iterations (!) for a simple function like ln(x) = 1. As shown in my previous post, the current method changes to bisection very early, so it converges very slowly.

I tried a simple Brent-style implementation in VBA for Excel, and here's the sequence of approximations I got:

 Estimate            abs. error
-------------------------------
2,0000 00000 00000 -7,18E-1
5,0000 00000 00000 2,28E+0
3,0046 57611 71378 2,86E-1
2,7574 19085 23714 3,91E-2
2,7180 91677 30546 -1,90E-4
2,7182 83194 09203 1,37E-6
2,7182 81828 50681 4,78E-11
2,7182 81828 45905 0,00E+0
2,7182 81828 45905
That's 15+ digits after just 9 function calls (including those for the two initial guesses) in a nice convergent series of approximations.

So what's going on in the 34s? Why does the current solver converge so slowly?

Dieter


#43

The 34S solver uses a variety of different methods. It certainly did use Brent at one stage, I'm not how much of that I removed as it evolved. Brent style quadratic interpolation and secant steps are used as is bisection and sometimes Ridder's method (but only after a bisection step). Bisection is used when the root is bracketed and which ever one of the other methods being used appears to produce a divergent x value -- this happens quite a lot unfortunately.

I'm not overly attached to the current solver and would certainly consider replacing it with something better.


- Pauli


#44

As far as I understand there are various methods that make sure that the interval narrows with each iteration step, i.e. the root (if there is one) is found within the interval initially given - provided the two initial guesses bracket the solution.

I just tried a modified regula falsi algorithm (with Anderson/Björck adjustment) - it finds a root within the given interval quite fast and since the calculation is simple and straightforward (as opposed to the Brent method) one or two more iterations do not matter much.

So I see two major points that should be resolved: first, if the initial guesses do not bracket a root, define a method to find such an interval. Second, find a foolproof way to decide whether a) a solution is found and, more important, b) no solution can be found. I mean, something better than "max_iterations exceeded". ;-)

Dieter


#45

That max iteration message is a reside of the C implementation -- I was extremely memory constrained and couldn't afford space to do much else.

- Pauli

#46

Have any of the 41 machine code gurus figured out how the advantage pack solver works?

If the algorithm used there is understood, it shouldn't be too hard to implement in the 34S.


- Pauli


#47

While the mcode specialists analyze the 41's advantage ROM code, here's an interesting document I found. It mentions two HP journal articles by W. Kahan and P. McClellan on the HP solvers built into the 34C and 18C. The proposed solver sounds very promising to me. The PDF can be found here.

What do the experts say?

Dieter


#48

Interesting, I hadn't seen that one.


- Pauli


Possibly Related Threads…
Thread Author Replies Views Last Post
  Prime feature request Stefan Dröge (Germany) 0 1,028 11-06-2013, 11:06 AM
Last Post: Stefan Dröge (Germany)
  HP Prime "Notes" feature request Charles Bennett 0 921 09-27-2013, 12:14 PM
Last Post: Charles Bennett
  wp34s binomial bug Andrew Nikitin 4 1,835 09-22-2013, 05:20 PM
Last Post: Paul Dale
  Expon bug in wp34s Andrew Nikitin 7 2,284 07-14-2013, 03:23 AM
Last Post: Marcus von Cube, Germany
  another wp34s bug Andrew Nikitin 8 2,559 06-26-2013, 01:01 AM
Last Post: Paul Dale
  weird statistics bug in wp34s Andrew Nikitin 5 2,159 06-20-2013, 01:54 PM
Last Post: Namir
  WP-34s feature suggestion: "Follow jump" shortcut Andrew Nikitin 3 1,528 06-12-2013, 01:42 AM
Last Post: Walter B
  [WP34S] A funny bug in Pi (prod) Eduardo Duenez 3 1,442 01-28-2013, 03:41 AM
Last Post: Walter B
  [WP34S] WP34S firmware on the AT91SAM7L-STK dev kit? jerome ibanes 1 1,226 10-04-2012, 04:59 PM
Last Post: Paul Dale
  Bug HP39GII....or feature? Bunuel66 14 3,735 08-11-2012, 02:59 PM
Last Post: Gilles Carpentier

Forum Jump: