Обработка математических ошибок в ARM Опции

[cyrax0]

Здравствуйте!
Я использую IAR EWARM 6.40 и STM32F407 (с FPU). Нужно останавливать контроллер при возникновении математических ошибок, например, при log(-1). В этом случае я получаю NaN в качестве результата, но по умолчанию никакого вызова функции или exception не происходит. Нету ли для ARM функции, аналогичной стандартной _matherr (math.h)? На других архитектурах она вызывается при следующих ошибках:

Код

    DOMAIN = 1,    /* argument domain error -- log (-1)        */
    SING,          /* argument singularity  -- pow (0,-2))     */
    OVERFLOW,      /* overflow range error  -- exp (1000)      */
    UNDERFLOW,     /* underflow range error -- exp (-1000)     */
    TLOSS,         /* total loss of significance -- sin(10e70) */
    PLOSS,         /* partial loss of signif. -- not used      */
    STACKFAULT     /* floating point unit stack overflow       */

Через FPU все эти ошибки, как я понял, тоже нельзя определить, например, регистр FPSCR Cortex-M4 среди битов исключений domain error не имеет:

Bit 7 IDC: Input denormal cumulative exception bit. Cumulative exception bit for floating-point
exception.
1: Indicates that the corresponding exception occurred since 0 was last written to it.
Bit 6:5 Reserved
Bit 4 IXC: Inexact cumulative exception bit. Cumulative exception bit for floating-point exception.
1: Indicates that the corresponding exception occurred since 0 was last written to it.
Bit 3 UFC: Underflow cumulative exception bit. Cumulative exception bit for floating-point exception.
1: Indicates that the corresponding exception occurred since 0 was last written to it.
Bit 2 OFC: Overflow cumulative exception bit. Cumulative exception bit for floating-point exception.
1: Indicates that the corresponding exception occurred since 0 was last written to it.
Bit 1 DZC: Division by zero cumulative exception bit. Cumulative exception bit for floating-point
exception. 1: Indicates that the corresponding exception occurred since 0 was last written to it.
Bit 0 IOC: Invalid operation cumulative exception bit. Cumulative exception bit for floating-point
exception. 1: Indicates that the corresponding exception occurred since 0 was last written to it.

Мне даже переполнение при результате операции 0.Infinite не удалось зафиксировать в этих флагах. Может, что-то не так сделал, но вроде бы в любом случае ситуацию log(-1) ни один из этих флагов не должен определять?
Заранее спасибо.

[Xenia]

Господа! Вы что, с ума сошли? Переменную errno в отдельную секцию выносить??? Все это можно сделать много проще, а не такими драконовскими мерами.

Заметим, что ни в одной функции, и log() в том числе, значение errno не присваивается явно, а только через вызов дефайнов:

Код

#define _FERAISE_INVALID()    { errno = EDOM;   }
#define _FERAISE_DIVBYZERO()  { errno = ERANGE; }
#define _FERAISE_OVERFLOW()   { errno = ERANGE; }
#define _FERAISE_UNDERFLOW()  { errno = ERANGE; }
#define _FERAISE_INEXACT()

Так что, если вы взялись за дело настолько серьезно, что решили перекомпилировать библиотеку, то было бы куда проще переопределить эти 5 дефиниций на свои функции.

Однако я полагаю, что и в этом нет большой нужды, т.к. эти дефиниции определены в хидере xmath.h не жестко, а условно:

Код

#if defined(math_errhandling) && math_errhandling == MATH_ERRNO

  /* Optimization of the _Feraise function when the IAR default is used. */
  #define _FERAISE_INVALID()    { errno = EDOM;   }
  #define _FERAISE_DIVBYZERO()  { errno = ERANGE; }
  #define _FERAISE_OVERFLOW()   { errno = ERANGE; }
  #define _FERAISE_UNDERFLOW()  { errno = ERANGE; }
  #define _FERAISE_INEXACT()

#else
... // здесь определена функция _Feraise()

#define _FERAISE_INVALID()    _Feraise(FE_INVALID)
#define _FERAISE_DIVBYZERO()  _Feraise(FE_DIVBYZERO)
#define _FERAISE_OVERFLOW()   _Feraise(FE_OVERFLOW)
#define _FERAISE_UNDERFLOW()  _Feraise(FE_UNDERFLOW)
#define _FERAISE_INEXACT()    _Feraise(FE_INEXACT)

#endif

Т.е. достаточно определить "math_errhandling" иначе, чем MATH_ERRNO, чтобы определение этих дефиниций изменилось на альтернативное, с явным вызовом функции _Feraise(), получающей код ошибки в качестве параметра.

Саму функцию _Feraise() при желании можно легко заменить на свою (т.к. самодельный вариант линкуется раньше библиотечного), а можно оставить той, что написана авторами компилятора - она там errno присваивает, а потом исключение вызывает.

Одного только я уверенно сказать не могу, нужно ли для этого перекомпилировать библиотеку или достататочно будет определить math_errhandling перед включением хидера xmath.h.

[scifi]

Одного только я уверенно сказать не могу, нужно ли для этого перекомпилировать библиотеку или достататочно будет определить math_errhandling перед включением хидера xmath.h.

Конечно перекомпилировать.
Но перекомпиляция библиотеки сопряжена с неудобствами. К примеру, чтобы подготовить новый компьютер к сборке программ, нужно помимо установки яра колдовать с библиотекой. В общем, надо взвешивать плюсы и минусы.

[Xenia]

Конечно перекомпилировать.
Но перекомпиляция библиотеки сопряжена с неудобствами. К примеру, чтобы подготовить новый компьютер к сборке программ, нужно помимо установки яра колдовать с библиотекой. В общем, надо взвешивать плюсы и минусы.

Тогда уж гораздо проще перекомпилировать один единственный модуль, в котором функция log() находится , а потом добавить этот объектник в проект. Жалко библиотеку курочить.

[scifi]

Тогда уж гораздо проще перекомпилировать один единственный модуль, в котором функция log() находится , а потом добавить этот объектник в проект. Жалко библиотеку курочить.

Подозреваю, что топикстартера не только логарифм интересует. А если он ещё захочет ловить деление на ноль, то ему точно от Си++ и перегрузки операторов не отвертеться.

[cyrax0]

Подозреваю, что топикстартера не только логарифм интересует. А если он ещё захочет ловить деление на ноль, то ему точно от Си++ и перегрузки операторов не отвертеться.

Деление на ноль как раз ловится элементарно:
SCB->CCR |= 0x10;
Без этого (по умолчанию) деление на ноль дает ноль.
С этим - HardFault.