I\'d like to play with those traps for educational purpose.
A common problem with the default behavior in numerical calculus is that we \"miss\" the Nan (or +-inf) that
As far as I know, you have two choices for floating point exception handling in C and C++:
First, if you disable/mask floating point exceptions (which most environments do by default), you can see whether any floating point exceptions have occurred by calling fetestexcept. fetestexcept isn't available in Visual C++, but you can steal the MinGW Runtime's implementation easily enough. (It's in the public domain.) Once an exception has been flagged, it's not cleared until you call feclearexcept, so you can call fetestexcept at the end of a series of calculations to see if any of them raised an exception. This doesn't give you the traps that you asked for, but it does let you test if problems like NaN or +/-inf have occurred and react as needed.
Second, you can enable/unmask floating point exceptions by calling feenableexcept in Linux or _controlfp in Windows. How the operating system handles a processor-generated floating point exception depends on your operating system.
__try
/ __catch
block in C or try
/ catch
block in C++._MM_SET_EXCEPTION_MASK
from xmmintrin.h
, and as long as you use the default compiler options (i.e., don't disable SSE), you should be able to catch exceptions using SIGFPE. (I've written a bit more on this and other floating point issues in C and C++ in this blog posting if you're curious.)