-fp-trap-all=all generates SIGFPE in log10

[mfactor]

Hello all,

I'm testing the icpc compiler and found a curious case: if I compile the code below

    [cpp]

#include <math.h>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <fenv.h>

using namespace std;

int main() {
feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
cout << log10(0.0049999999992) <<endl;

return 0;
}

[/cpp]

If I compile it with icpc log10.cxx -lm -g, it runs ok and gives the right answer (-2.30103), but if I enable icpc log10.cxx -lm -g -fp-trap-all=all , I get a SIGFPE. I'm using icpc version icpc (ICC) 13.0.0 20120731 on Linux x86_64. Any ideas?

Cheers 

[SergeyKostrov]

>>...cout << log10(0.0049999999992)... Please try to do another test with CRT-function printf or separate calls to log10 and cout. Actually this is the 3rd time I see a problem when cout is used with some CRT-function from math.h header file. Thanks in advance for results of your additional test.

[mfactor]

Hello Sergey,

[cpp]

#include <math.h>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <fenv.h>

int main() {
feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
double x = 0;
x = log10(0.0049999999992) ;
printf("%f",x);
return 0;
}

[/cpp]

Exactly same behavior. With g++ 4.6.3 exactly the same code does not generate a SIGFPE.

Cheers

[mfactor]

Hello again,

I managed to narow it down a little bit more. An ANSI C version of the same code compiled with gcc, icc AND icpc gives no SIGFPE:

[cpp]

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <fenv.h>

int main() {
feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
double x = 0;
x = log10(0.0049999999992) ;
printf("%f",x);
return 0;
}

[/cpp] 

[cpp]

icc log10.cxx  -g -fp-trap-all=all - no SIGFPE - ok

icpc log10.cxx  -g -fp-trap-all=all  - no SIGFPE - ok

gcc log10.cxx -g -lm - no SIGFPE - ok

[/cpp] 

This code also compiles with all the options above and runs without SIGFPE:

[cpp]

#include <math.h>
#include <cstdio>
#include <cstdlib>
#include <fenv.h>

int main() {
feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
double x = 0;
x = log10(0.0049999999992) ;
printf("%f",x);
return 0;
}

[/cpp]

But if <iostream> is included, I get the SIGFPE with icpc and ftrap, even without any mention to cout or printf.

Cheers

[SergeyKostrov]

Thank you. Take into account that on a 64-bit platform with Intel CPUs a precision control can not be used and this is by design. You'll need to do a set of tests like: Debug Configuration - For 32-bit platform with "iostream" - For 32-bit platform without "iostream" - For 64-bit platform with "iostream" - For 64-bit platform without "iostream" Release Configuration - For 32-bit platform with "iostream" - For 32-bit platform without "iostream" - For 64-bit platform with "iostream" - For 64-bit platform without "iostream" and with all the rest variants of the test case (!). If some code in "iostream" ( possibly static which does some initialization ) tries to change a precision of the Floating-Point Unit with a CRT-function _control87 ( or similar ) on a 64-bit platform than an exception will be thrown and this is what you see. Also, it is not clear for me on what platform you're doing your tests. You've metioned ''...Linux x86_64...". So, did you get the exception on 32-bit platform or 64-bit platform?

[Bernard]

>>>You've metioned ''...Linux x86_64...".>>>

That's mean that 64-bit Linux version is used.

[TimP]

Original post did leave unspecified whether the 32-bit ia32 or 64-bit Intel64 compiler was in use.  As g++ 4.6.3 was mentioned, one might expect a recent linux distro where it's easier to install the Intel64 compiler than ia32, but this is all speculation.

I don't find -fp-trap-all listed as a g++ option, so I might suspect g++ ignored it.

As I've seen log10() implemented as a macro in terms of log(), the traditional comparison between icpc and g++ pre-processed files also might be relevant, along with the traditional icpc -V and g++ -v quotations.

[SergeyKostrov]

>>...As I've seen log10() implemented as a macro in terms of log()... I think it needs to be confirmed by 'mfactor' user if something like this ... #define log10( x ) ( log( x ) / log( 10.0L ) ) ... is used. However, it doesn't bring clarity on all the rest our questions.

[mfactor]

Hi all,

Some informations: running Ubuntu 12.10 with kernel

[cpp]

Linux voyager 3.5.0-25-generic #38-Ubuntu SMP Mon Feb 18 23:27:42 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux

[/cpp]

So I confirm this is a 64 bits platform. g++ -v prints out

[cpp]

COLLECT_GCC=g++
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/4.7/lto-wrapper
Destino: x86_64-linux-gnu
Configurado com: ../src/configure -v --with-pkgversion='Ubuntu/Linaro 4.7.2-2ubuntu1' --with-bugurl=file:///usr/share/doc/gcc-4.7/README.Bugs --enable-languages=c,c++,go,fortran,objc,obj-c++ --prefix=/usr --program-suffix=-4.7 --enable-shared --enable-linker-build-id --with-system-zlib --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --with-gxx-include-dir=/usr/include/c++/4.7 --libdir=/usr/lib --enable-nls --with-sysroot=/ --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --enable-gnu-unique-object --enable-plugin --enable-objc-gc --disable-werror --with-arch-32=i686 --with-tune=generic --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Modelo de thread: posix
versão do gcc 4.7.2 (Ubuntu/Linaro 4.7.2-2ubuntu1)

[/cpp]

icpc -V prints

[cpp]

Intel(R) C++ Compiler XE for applications running on IA-32, Version 13.0.0.079 Build 20120731
Copyright (C) 1985-2012 Intel Corporation. All rights reserved.

[/cpp]

and

[cpp]

Intel(R) C++ Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 13.0.0.079 Build 20120731
Copyright (C) 1985-2012 Intel Corporation. All rights reserved.

[/cpp]

Depending the "source compilevars.sh" I use (ia32 or intel64).

And now, for the tests:

1) g++ -m32 -g log10.cxx -lm with iostream: OK

2) g++ -m64 -g log10.cxx -lm with iostream: OK

3) g++ -m32 -g log10.cxx -lm without iostream: OK

4) g++ -m64 -g log10.cxx -lm without iostream: OK

Tests with Intel compiler:

1) icpc -m64 -fp-trap-all=all -g log10.cxx with iostream: *** SIGFPE ***

2) icpc -m64 -fp-trap-all=all -g log10.cxx without iostream: OK

3) icpc -m64  -g log10.cxx with iostream: OK

4) icpc -m64  -g log10.cxx without iostream: OK

5) icpc -m32 -g log10.cxx with iostream: OK

6) icpc -m32 -g log10.cxx without iostream: OK

7) icpc -m32 -fp-trap-all=all -g log10.cxx with iostream: *** SIGFPE ***

8) icpc -m32 -fp-trap-all=all -g log10.cxx without iostream: OK

I'm not using any macros whatsoever. The test code follows:

[cpp]

#include <math.h>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <fenv.h>

using namespace std;

int main() {
feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
double x = 0.0;
x = log10(0.0049999999992) ;
printf("%f",x);
return 0;
}

[/cpp]

Hope this helps.

[SergeyKostrov]

>>Tests with Intel compiler: >> >>1) icpc -m64 -fp-trap-all=all -g log10.cxx with iostream: *** SIGFPE *** >>... >>7) icpc -m32 -fp-trap-all=all -g log10.cxx with iostream: *** SIGFPE *** Thank you for these test results.

[SergeyKostrov]

There is some inconsistency with your tests for g++ because you didn't use option fp-trap-all. If that option is not supported than another two could be used to verify that the problem with SIGFPE exception doesn't exist for g++: -ftrapping-math - Assume floating-point operations can trap -ftrapv - Trap for signed overflow in addition, subtraction and multiplication Thanks in advance.

[mfactor]

Hello Sergey,

The function feenableexcept should catch this kind of errors in C++, with these options:

1) g++ -m32 -ftrapping-math -ftrapv -fsignaling-nans -g log10.cxx with with iostream: OK

2) g++ -m64 -ftrapping-math -ftrapv -fsignaling-nans -g log10.cxx with with iostream: OK

3) g++ -m32 -ftrapping-math -ftrapv -fsignaling-nans -g log10.cxx with without iostream: OK

4) g++ -m64 -ftrapping-math -ftrapv -fsignaling-nans -g log10.cxx with without iostream: OK

I currently don't have another platform to do the tests. Can you compile the code that I submitted? Does it run? Does it give SIGFPE?

Cheers

[SergeyKostrov]

>>...Can you compile the code that I submitted? Yes. >>...Does it run? Does it give SIGFPE? I'll provide a report for two Intel C++ compilers, versions 12 and 13 ( 32-bit and 64-bit ) on two Windows platforms ( XP and 7 ).

[SergeyKostrov]

Here are results for 32-bit test application on Windows XP: [ Windows XP 32-bit - Debug - with 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 12.1.3.300 Build 20120130 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 8.00.50727.762 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows XP 32-bit - Debug - without 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 12.1.3.300 Build 20120130 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 8.00.50727.762 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows XP 32-bit - Release - with 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 12.1.3.300 Build 20120130 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 8.00.50727.762 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows XP 32-bit - Release - without 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 12.1.3.300 Build 20120130 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 8.00.50727.762 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions

[SergeyKostrov]

I'll continue investigation with Intel C++ compiler ( v13 / 32-bit and 64-bit ) on Windows 7 Professional ( 64-bit ) some time later this week.

[SergeyKostrov]

Could you upload cstdio, cstdlib and iostream header files ( of course GCC versions ) for review?

[mfactor]

I found them in the GNU software compilation:

http://gcc.gnu.org/onlinedocs/libstdc++/libstdc++-html-USERS-3.4/iostream-source.html

http://gcc.gnu.org/onlinedocs/libstdc++/libstdc++-html-USERS-3.4/classstd_1_1basic__iostream.html

http://gcc.gnu.org/onlinedocs/libstdc++/libstdc++-html-USERS-3.3/cstdlib-source.html

http://gcc.gnu.org/onlinedocs/gcc-4.7.2/libstdc++/api/a01063_source.html

http://gcc.gnu.org/onlinedocs/gcc-4.7.2/libstdc++/api/

Cheers

[mfactor]

I posted the links, put it appeared that my comment was flagged for review. 

[SergeyKostrov]

Here are results for a 32-bit test application on Windows 7 Professional ( 64-bit ): [ Windows 7 Professional - 32-bit application - Debug - with 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 32-bit application - Debug - without 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 32-bit application - Release - with 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 32-bit application - Release - without 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Compiler XE for applications running on IA-32, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions

[SergeyKostrov]

Here are results for a 64-bit test application on Windows 7 Professional ( 64-bit ): [ Windows 7 Professional - 64-bit application - Debug - with 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 64-bit application - Debug - without 'iostream' ] ..\FpTrapAll>icl /MDd /Qfp-trap-all=all /D_DEBUG Main.cpp Intel(R) C++ Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 64-bit application - Release - with 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions [ Windows 7 Professional - 64-bit application - Release - without 'iostream' ] ..\FpTrapAll>icl /MD /Qfp-trap-all=all /DNDEBUG Main.cpp Intel(R) C++ Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 13.0.0.089 Build 20120731 Copyright (C) 1985-2012 Intel Corporation. All rights reserved. Main.cpp Microsoft (R) Incremental Linker Version 9.00.30729.01 Copyright (C) Microsoft Corporation. All rights reserved. -out:Main.exe Main.obj ..\FpTrapAll>Main.exe -2.301030 Note: No Exceptions

[SergeyKostrov]

My conclusion: I couldn't reproduce the problem with SIGFPE exception on Windows 32-bit and 64-bit platforms with Intel C++ compilers ( v12 and v13 ). It would be nice to have additional verification by Intel software engineers. [ Source code of All tests ] #include "math.h" #include "cstdio" #include "cstdlib" #include "iostream" #include "fenv.h" using namespace std; int main( void ) { feenableexcept( FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW ); double x = 0.0L; x = log10( 0.0049999999992 ); printf( "%f", x ); return ( int )0; } Note: I'll take a look at sources for the STL header files you've provided. Thank you.