That sort of makes sense :-)
Is there a module that covers over these sort of cracks? It would seem to be a routine problem. | [reply] |
Is there a module that covers over these sort of cracks?
The crack that has perl's print() function often lie to us about a floating point value is rather annoying.
You just need to remember that the floating point value that print() provides will be the same as that provided by printf "%.14e", whereas, if you want a practically useful value, you need to printf "%.16e".
Outputting more than 17 decimal digits of a double is generally of no use. If you want to see the exact value, use printf "%a" (which will display that exact value in hex).
I gather that there was once a $# special variable (see perldoc perlvar) which allowed one to specify the decimal precision that print() would use.
But it has been abandoned and disabled owing to difficulties in its implementation.
You can, of course, write your own subroutine for printing floats in the format of your choice:
sub p17 {
printf "%.16e\n", $_;
}
sub ph {
printf "%a\n", $_;
}
The breaking of the Associative Law of Addition that you experienced is something that can happen when floating point operations overflow the fixed precision.
To be guaranteed of avoiding that, I think you'd need to work in an environment that avoids overflow by allowing the precision to increase as needed.
Math::BigFloat does appear to provide that.
Cheers,
Rob | [reply]
[d/l]
[select] |
| [reply] |
