feexceptErsatz.H

Go to the documentation of this file.

// ============================================================================
// Generic handling of floating point exceptions
//
// 2018 Alexey Matveichev
// 2021 Tatsuya Shimizu - ARM64 support
// 2023 Guanyang Xue - improve ARM64 support
// ----------------------------------------------------------------------------
//
// Original Author
//     David N. Williams
//
// License
//     Public Domain
//
// Description
//     Alternate, nondefault handling of IEEE 754 floating-point exceptions
//     in OS X, based on the GNU functions feenableexcept(), fedisableeexcept()
//     and fegetexcept() [in libm]
//
// http://www-personal.umich.edu/~williams/archive/computation/fe-handling-example.c
// ============================================================================
 
#ifndef Foam_feexceptErsatz_H
#define Foam_feexceptErsatz_H
 
#ifdef __APPLE__
#include <fenv.h>
 
// Workaround for Apple Silicon - has SIGILL (illegal instruction) not SIGFPE
#if defined __arm64__
#undef  SIGFPE
#define SIGFPE SIGILL
#endif
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
inline int feenableexcept(unsigned int excepts)
{
    static fenv_t fenv;
    unsigned int new_excepts = excepts & FE_ALL_EXCEPT;
    unsigned int old_excepts;   // previous masks
 
    if (fegetenv(&fenv))
    {
        return -1;
    }
 
#if defined __arm64__
    old_excepts = fenv.__fpsr & FE_ALL_EXCEPT;
 
    fenv.__fpsr |= new_excepts;
    fenv.__fpcr |= (new_excepts << 8);
#else
    old_excepts = fenv.__control & FE_ALL_EXCEPT;
 
    fenv.__control &= ~new_excepts;
    fenv.__mxcsr   &= ~(new_excepts << 7);
#endif
 
    return fesetenv(&fenv) ? -1 : old_excepts;
}
 
 
inline int fedisableexcept(unsigned int excepts)
{
    static fenv_t fenv;
    unsigned int new_excepts = excepts & FE_ALL_EXCEPT;
    unsigned int old_excepts;   // all previous masks
 
    if (fegetenv(&fenv))
    {
        return -1;
    }
 
#if defined __arm64__
    old_excepts = fenv.__fpsr & FE_ALL_EXCEPT;
 
    fenv.__fpsr &= ~new_excepts;
    fenv.__fpcr &= ~(new_excepts << 8);
#else
    old_excepts = fenv.__control & FE_ALL_EXCEPT;
 
    fenv.__control |= new_excepts;
    fenv.__mxcsr   |= (new_excepts << 7);
#endif
 
    return fesetenv(&fenv) ? -1 : old_excepts;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif  // __APPLE__
#endif
 
// ************************************************************************* //