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 :

Very shortly, all function names start with the letter x.
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.
Sure good idea, that's fairly simple.

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

Maybe I'm doing something wrong?

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


LN(0.01) -> -4.6051....
Reference URL's