HP Forums

Full Version: Extended precision library
You're currently viewing a stripped down version of our content. View the full version with proper formatting.
This is a library of functions to allow unlimited floating point precision for all arithmetic and most transcendental functions on the Prime.

Explanations can be found (in French...) in the 11th issue of the french journal "La Gazette des Pocketicaires" available for free download at :
http://silicium.org/site/index.php/telec...keticaires.

To start type xinit(7) for example where 7 is the number of slices of 6 digits that you want to have available.
Here 7 means you calculate on 42 digits.
Values are stored as lists.
To convert from/to the readable representation which is a string, you use two functions.
Example :
xtof("1256.23589") returns a list
xtostr(that list) gives you back the string "1256.23589".

Then you can use the lists created this way to calculate, see the examples in the test file attached.

Feedback welcome.
Just an idea on the user interface
I am wondering whether users might find it more natural to specify the number of digits they want to use, rather than the number of slices.
As xinit(42) seems more intuitive than xinit(7) when the user wants 42 digits.
Hello,
Sure good idea, that's fairly simple.

List of changes to do :
- use Newton methode for division
- change xinit so that it takes a number of digits and not slices
- normalize angles as input to sin/cos/tan
- (hard) fix precision issues in cordic functions

G.E.
Maybe I'm doing something wrong?

xinit(5) -> 5
xtostr(xlog(xtof("0.01"))) -> "-999997697417.20959.....

while

LN(0.01) -> -4.6051....
Reference URL's
• HP Forums: https://www.hpmuseum.org/forum/index.php
• :