Re^2: How to get better exponentiation?

Thanks, haukex, POSIX::cbrt is just what I need! But I’m disappointed to find that POSIX::pow has the same limitations as **. So the more general question remains open.

Thanks again,

Athanasius <°(((>< contra mundum Iustus alius egestas vitae, eros Piratica,

Comment on Re^2: How to get better exponentiation? Select or Download Code

Replies are listed 'Best First'.
Re^3: How to get better exponentiation? (updated) by haukex (Archbishop) on Jan 22, 2022 at 11:58 UTC
But I’m disappointed to find that `POSIX::pow` has the same limitations as ``. So the more general question remains open. I like the module Math::Prime::Util::GMP a lot because of its many functions that also go beyond prime numbers. `$ perl -MMath::Prime::Util::GMP=rootreal -le 'print rootreal(-8, 3)' -2.000000000000000000000000000000000000000` [download] Use `0+rootreal(...)` to get just -2 in this case. `powreal(-8, 1/3)` suffers a bit from the floating-point inaccuracies in `1/3`. Note both functions support an optional third argument to specify the number of significant digits, e.g. `powreal(-8, 1/3, 5)` is `"-2.0000"`. Update: However, note the issues raised by vr and syphilis here!	[reply] [d/l] [select]

Replies are listed 'Best First'.

Re^3: How to get better exponentiation? (updated)
by haukex (Archbishop) on Jan 22, 2022 at 11:58 UTC

But I’m disappointed to find that POSIX::pow has the same limitations as **. So the more general question remains open.

I like the module Math::Prime::Util::GMP a lot because of its many functions that also go beyond prime numbers.

$ perl -MMath::Prime::Util::GMP=rootreal -le 'print rootreal(-8, 3)'
-2.000000000000000000000000000000000000000
[download]

Use 0+rootreal(...) to get just -2 in this case. powreal(-8, 1/3) suffers a bit from the floating-point inaccuracies in 1/3. Note both functions support an optional third argument to specify the number of significant digits, e.g. powreal(-8, 1/3, 5) is "-2.0000".

Update: However, note the issues raised by vr and syphilis here!

[reply]
[d/l]
[select]