Le opzioni di ottimizzazione g ++ hanno effetto sul valore della funzione sin

Question

Ho un problema con la funzione "sin" di 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; 
} 

Quando compilare il codice di cui sopra con "g ++", senza alcuna opzione di ottimizzazione, sarebbe uscita "-0,00000000000000024492127076447545". Ma se con l'opzione "-O3", emetterebbe "-0,00000000000000024492935982947064".

Perché non restituisce "-0,00000000000000024492935982947064" senza "-O3"? Grazie in anticipo.

Answer 1

Perché con "-O3" il compilatore precompila sin(2*pi) in fase di compilazione, con un algoritmo. Senza "-O3" questo è calcolato in runtime, con un altro algoritmo.

Ciò può essere dovuto al fatto che il compilatore stesso è stato creato con alcune librerie matematiche, che differiscono dalla libreria matematica.

Aggiornamento

L'unica entità, dando il risultato "-0,00000000000000024492127076447545" è la versione a 32 bit di libstdC++. La versione a 64 bit della stessa libreria e di gcc produce "-0,00000000000000024492935982947064".

Quindi l'aggiornamento alla versione più recente non sarà di aiuto. Ho anche provato varie opzioni, proposte qui: né -ffloat-store, né -fno-builtin non fanno alcuna differenza, così come long double e sinl.

libstdC++ a 32 bit utilizza 387 istruzioni in virgola mobile, mentre gcc utilizza apparentemente le istruzioni SSE. Ecco la differenza. Probabilmente, l'unico modo per renderli coerenti è ricostruire gcc dai sorgenti, indirizzandoli a utilizzare internamente solo 387 istruzioni.