Consider implementing a 'strict/"most determinstic" floating point' mode

[axic]

This would mean to return the canonical NaN (which is subset of the arithmetic NaNs) in every case where the spec allows for non-canonical NaNs.

[chfast]

Precise: return positive canonical NaN.