HP41: accuracy of 13digit routines

09032015, 06:22 AM
(This post was last modified: 09032015 06:23 AM by Ángel Martin.)
Post: #9




RE: HP41: accuracy of 13digit routines
(09022015 06:07 PM)Dieter Wrote: So the accuracy level depends on the function, and sometimes even on its argument. Since the four basic arithmetic functions as well as x² and √x deliver correctly rounded results they must be virtually exact internally, i.e. in all 13 internal digits, maybe the last one truncated. This would yield a max. error of 0,501 ULP. Sure, I meant that the accuracy will always be correct to the 10th decimal digit but I didn't imply that all algorithms would be equally effective in the digits beyond that. Definitely the transcendental functions are more complicated to implement, but not only those show weak points. I don't know if there were changes made to the 15C MCODE, but I can tell you that on the 41 the square root at least is far from being that accurate. If you do a chained calculation in MCODE with [SQR13] followed by a 13digit multiplication you not always get the same original argument to the 13th. digit.... and there's NOT such a thing as X^2 as individual function either, it uses the same multiplication code. (09022015 06:07 PM)Dieter Wrote: What does this mean for a Gamma implementation in 13digit MCode? The functions used in the Lanczos algorithm are mostly +, –, x and / which are (nearly) exact. On the other hand the logs may only have 11 valid digits. So this is what can be expected: about 11 exact digits, and a final error similar to that of the internal log or trig functions. the X! on the 15C is the gamma function, right? So there you also have it to test for the accuracy to the 10th decimal digit. As far as GAMMA in the SandMath I checked the accuracy for the natural arguments were ok (see the table in its manual) but didn't do any further checks for noninteger arguments. BTW there's also exponentials in the Lanczos formula, as well as sine factors in the reflection formula. Cheers, 'AM "To live or die by your own sword one must first learn to wield it aptly." 

« Next Oldest  Next Newest »

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