News for hackers. Stuff that matters.
Re:Beware of Fast-Math by goalieca on Saturday May 31, 2025 @09:10AM (#44144065)
Floating point arithmetic is neither commutative or associative so you shouldn’t.
Re:Beware of Fast-Math by lo0dot0 on Saturday May 31, 2025 @09:51AM (#44144277)
While it technically correct to say this it also gets the wrong point across because it leaves out the fact that ordering changes create only a small difference. Other examples where arithmetic is not commutative, e.g. matrix multiplication , can create much larger differences.
Re:Beware of Fast-Math by kstrauser on Saturday May 31, 2025 @03:42PM (#44146532)
> ordering changes create only a small difference.

That can’t be assumed.

You can easily fall into a situation like:

  total = large_float_value
  for _ in range(1_000_000_000):
    total += .01
  assert total == large_float_value
Without knowing the specific situation, it’s impossible to say whether that’s a tolerably small difference.
Re:Beware of Fast-Math by layer8 on Saturday May 31, 2025 @10:42AM (#44144597)
IEEE-754 addition and multiplication is commutative. It isn't distributive, though.
Re:Beware of Fast-Math by eapriv on Saturday May 31, 2025 @10:06AM (#44144362)
Why is it not commutative?
Re:Beware of Fast-Math by layer8 on Saturday May 31, 2025 @10:40AM (#44144586)
It actually is commutative according to IEEE-754, except that in the case of a NaN result you might get a different NaN representation.
Re:Beware of Fast-Math by adgjlsfhk1 on Saturday May 31, 2025 @06:03PM (#44147229)
having multiple NaNs and no spec for how they should behave feels like such an unforced error to me
Re:Beware of Fast-Math by layer8 on Saturday May 31, 2025 @08:59PM (#44147950)
For mathematical use, NaN payloads shouldn’t matter, and behave identically (aside from quiet vs. signaling NaNs). It also doesn’t matter for equality comparison, because NaNs always compare unequal.
Re:Beware of Fast-Math by adgjlsfhk1 on Sunday June 01, 2025 @08:57AM (#44150546)
from the user perspective it's not too bad, but from the compiler perspective it is. The result of this is that LLVM has decided that trying to figure out which nan you got (e.g. by casting to an Int and comparing) is UB, which means pretty much every floating point operation becomes non-deterministic.

This also adds extra complexity to the CPU. you need special hardware for == rather than just using the perfectly good integer unit, and every fpu operation needs to devote a bunch of transistors to handling this nonsense that buys the user absolutely nothing.

there are definitely things to criticize about the design of Posits, but the thing they 100% get right is having a single NaN and sane ordering semantics

Re:Beware of Fast-Math by StefanKarpinski on Monday June 02, 2025 @09:29AM (#44158650)
Floating-point arithmetic is non-associative, but it is commutative for the operations that are algebraically commutative: x + y == y + x and x*y == y*x. And x - y = -(y - x) so subtraction is properly anti-commutative.

The only very marginal exception to this is that when both arguments are NaN, the return value will be NaN, but which NaN payload is returned can depend on argument order. But no one ever uses this because it's not specified, so it can't be used reliably for anything useful. The behavior I wish IEEE 754 had specified for this is to define a standard NaN value (or two), and when the return value of an op is NaN, and some of the arguments are non-standard NaNs, then one of those non-standard NaN values must be returned. This doesn't depend on argument order and allows NaN payloads to be reliably propagated, which would let you encode useful debugging information in NaN payloads and know that it will flow through the program.