Agregar dos números de coma flotante
Me gustaría calcular la suma, redondeada, de dos números binarios IEEE 75464. Con ese fin, escribí el programa C99 a continuación:
#include <stdio.h>
#include <fenv.h>
#pragma STDC FENV_ACCESS ON
int main(int c, char *v[]){
fesetround(FE_UPWARD);
printf("%a\n", 0x1.0p0 + 0x1.0p-80);
}
Sin embargo, si compilo y ejecuto mi programa con varios compiladores:
$ gcc -v … gcc version 4.2.1 (Apple Inc. build 5664) $ gcc -Wall -std=c99 add.c && ./a.out add.c:3: warning: ignoring #pragma STDC FENV_ACCESS 0x1p+0 $ clang -v Apple clang version 1.5 (tags/Apple/clang-60) Target: x86_64-apple-darwin10 Thread model: posix $ clang -Wall -std=c99 add.c && ./a.out add.c:3:14: warning: pragma STDC FENV_ACCESS ON is not supported, ignoring pragma [-Wunknown-pragmas] #pragma STDC FENV_ACCESS ON ^ 1 warning generated. 0x1p+0
¡No funciona! (Esperaba el resultado0x1.0000000000001p0
)
De hecho, el cálculo se realizó en tiempo de compilación en el modo predeterminado de redondeo a más cercano:
$ clang -Wall -std=c99 -S add.c && cat add.s add.c:3:14: warning: pragma STDC FENV_ACCESS ON is not supported, ignoring pragma [-Wunknown-pragmas] #pragma STDC FENV_ACCESS ON ^ 1 warning generated. … LCPI1_0: .quad 4607182418800017408 … callq _fesetround movb $1, %cl movsd LCPI1_0(%rip), %xmm0 leaq L_.str(%rip), %rdx movq %rdx, %rdi movb %cl, %al callq _printf … L_.str: .asciz "%a\n"
Sí, vi la advertencia emitida por cada compilador. Entiendo que activar o desactivar las optimizaciones aplicables a escala de línea puede ser complicado. Todavía me gustaría, si fuera posible, desactivarlos a la escala del archivo, lo que sería suficiente para resolver mi pregunta.
Mi pregunta es: ¿qué opciones de línea de comandos debo usar con GCC o Clang para compilar una unidad de compilación C99 que contiene código destinado a ejecutarse con un modo de redondeo de FPU que no sea el predeterminado?
DigresiónMientras investigaba esta pregunta, encontré estoPágina de cumplimiento de GCC C99, que contiene la entrada a continuación, que dejaré aquí en caso de que alguien más lo encuentre divertido. Grrrr.
floating-point | | environment access | N/A | Library feature, no compiler support required. in <fenv.h> | |