g ++ opciones de optimización afectan el valor de la función sin

Question

Tengo un problema con la función "sin" de libc.

#include <cmath> 
#include <stdio.h> 

int main(int argc, char **argv) 
{ 
    double tt = 6.28318530717958620000; // 2 * M_PI 
    double yy = ::sin(tt); 

    printf("%.32f\n", yy); 

    return 0; 
} 

Al compilar el código anterior utilizando "g ++" sin ninguna opción de optimización, lo haría de salida "-0,00000000000000024492127076447545". Pero si con la opción "-O3", mostraría "-0.00000000000000024492935982947064".

¿Por qué no devuelve "-0.00000000000000024492935982947064" sin "-O3"? Gracias de antemano.

Answer 1

Porque con "-O3" el compilador calcula previamente sin(2*pi) en tiempo de compilación, con un algoritmo. Sin "-O3" esto se calcula en tiempo de ejecución, con otro algoritmo.

Esto puede deberse a que el compilador en sí se creó con alguna biblioteca matemática, que difiere de su biblioteca matemática.

actualización

La única entidad, dando el resultado "-0,00000000000000024492127076447545" es la versión de 32 bits de libstdC++. La versión de 64 bits de la misma biblioteca y el mismo gcc producen "-0.00000000000000024492935982947064".

Por lo tanto, la actualización a una versión más nueva no ayudará. También probé varias opciones, propuestas aquí: ni -ffloat-store, ni -fno-builtin no hacen ninguna diferencia, así como long double y sinl.

32-bit libstdC++ usa 387 instrucciones de coma flotante, mientras que gcc aparentemente usa instrucciones de SSE. Aquí está la diferencia. Probablemente, la única forma de hacer que sean coherentes es reconstruir gcc a partir de las fuentes, orientándolas a usar solo 387 instrucciones internamente.