BUG, SIMD: Fix rounding large numbers >= 2^52 on SSE2

  Before SSE41, there were no native instructions for rounding
  operations on double precision. We usually emulate it by assuming
  that the `MXCR` register is set to rounding, adding a
  large number `2^52` to `X` and then subtracting it back to
  eliminate any excess precision as long as `|X|` is less than `2^52`
  otherwise returns `X.`

  The current emulated intrinics `npyv_[rint,floor, ceil, trunc]_f64`
  was not checking whether `|x|` equal or large `2^52` which leads
  to losing accuracy on large numbers.

0 files changed, 0 insertions, 0 deletions