La división y el módulo enteros suelen ser operaciones costosas en los procesadores modernos, con latencias elevadas y un rendimiento modesto, mientras que la división en coma flotante es más rápida, tiene mayor throughput y suele disponer de más unidades de ejecución. El artículo demuestra que, para enteros sin signo que caben en 24 o 53 bits según se trabaje en precisión simple o doble, resulta seguro y directo calcular cociente y resto mediante aritmética de coma flotante: d = trunc(x / y) y m = -fma(d, y, -x). La conversión de vuelta a entero hace gratuito el truncamiento. La clave de la corrección es que, en redondeo al más cercano con empates al par, la división binaria no produce empates porque carece de puntos medios exactos. Esto garantiza que la parte fraccional nunca basta para redondear al entero siguiente, de modo que el cociente en coma flotante coincide con el entero. El autor verifica el resultado con divisores concretos (1, 2, 3, 4-7) y demuestra algebraicamente que el máximo fraccional posible siempre tiene menos bits significativos que el umbral de redondeo. Entre las consideraciones prácticas menciona el coste de cambiar el modo de redondeo, la utilidad de instrucciones con modo seleccionable como _mm_div_round_sd de AVX-512, la ausencia en la mayoría de CPUs x64 de operaciones nativas de conversión sin signo a coma flotante y la ventaja adicional cuando se trabaja en SIMD. Para divisores constantes el compilador ya optimiza, y para divisores reutilizados en tiempo de ejecución existen librerías específicas como libdivide. El artículo enlaza con un trabajo previamente verificado formalmente sobre el mismo método.
Noticias agregadas con IA