реклама на сайте
подробности

 
 
> Обработка математических ошибок в ARM
cyrax0
сообщение Jan 14 2015, 14:36
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 22
Регистрация: 9-12-14
Пользователь №: 84 046



Здравствуйте!
Я использую 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) ни один из этих флагов не должен определять?
Заранее спасибо.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Xenia
сообщение Jan 14 2015, 22:42
Сообщение #2


Гуру
******

Группа: Модератор FTP
Сообщений: 4 479
Регистрация: 20-02-08
Из: Москва
Пользователь №: 35 237



Цитата(cyrax0 @ Jan 14 2015, 20:02) *
А можно поподробнее, что за errno? Отмечу еще раз, что мне нужна не просто проверка на ошибки (не делать же ее после каждой операции), а что-то вроде исключения, если результат математической операции = NaN, Inf, etc.


А чем плоха проверка результата? Неприятие к этой процедуре можно было еще объяснить, если бы STM32F вычислял логарифмы аппаратно, но ведь этого-то нет! Он способен лишь на то, чтобы выполнять с плавающей точкой 4 арифметрических действия. Ну, еще два числа умеет сравнивать. Но ведь функций-то вычислять аппаратно он не умеет! А, стало быть, все эти логарифмы, экспоненты и тригонометрия - галимая эмуляция sm.gif. А раз так, то вы крайне незначительно замедлите вычисление логарифма, если решитесь на проверку результата. Опять же всегда можно создать свою функцию типа my_log(), которая вызывает нативную функцию log(), проверяя после нее результат.

Учитывая эти обстоятельства, предлагаю вам просто заглянуть в исходники и поглядеть, как там реализована функция log(). Я сама заглядывала (только у меня EWARM не 6.40, а 7.30), но нашла там лишь присвоение кода ошибки глобальной переменной errno. Хотя в хидере xmath.h есть что-то похожее на возможность установки собственного хэндлера для реакции на математические ошибки. Глубже не разбиралась.
Go to the top of the page
 
+Quote Post
cyrax0
сообщение Jan 15 2015, 07:59
Сообщение #3


Участник
*

Группа: Участник
Сообщений: 22
Регистрация: 9-12-14
Пользователь №: 84 046



Цитата(Xenia @ Jan 15 2015, 02:42) *
А чем плоха проверка результата? Неприятие к этой процедуре можно было еще объяснить, если бы STM32F вычислял логарифмы аппаратно, но ведь этого-то нет!

У нас математические алгоритмы пишут разные сотрудники, и кто-то может забыть проверить результат. Поэтому нужно проверять такие ошибки на системном уровне.

Цитата
Хотя в хидере xmath.h есть что-то похожее на возможность установки собственного хэндлера для реакции на математические ошибки. Глубже не разбиралась.

Хорошая мысль, но похоже эта _Feraise не вызывается автоматически, по крайней мере брекпоинт на нее не ставится. Либо у меня не подключена DLib, не знаю, как это делается.

Цитата
Через MPU защищаете некую область памяти от записи, линкуете туда errno и получаете прерывание MPU каждый раз при попытке записи в errno.

Прикольная идея, спасибо. Под линковкой понимается приделать к errno:
Код
@ какой-то_адрес

и как это сделать, модифицировать стандартный код?
Go to the top of the page
 
+Quote Post
jcxz
сообщение Jan 15 2015, 08:20
Сообщение #4


Гуру
******

Группа: Свой
Сообщений: 5 228
Регистрация: 3-07-08
Из: Омск
Пользователь №: 38 713



Цитата(cyrax0 @ Jan 15 2015, 13:59) *
Прикольная идея, спасибо. Под линковкой понимается приделать к errno:
Код
@ какой-то_адрес

и как это сделать, модифицировать стандартный код?

Если errno не вынесен в отдельную секцию, которую можно отдельно линковать, то придётся править исходники - вынести его в отдельную секцию.
Если вынесен (или в секции, где он находится, остальные переменные используются как read-only) - то библиотеку править не надо, только свой icf-файл.
Посмотрите свой map-файл - в какой секции находится errno.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- cyrax0   Обработка математических ошибок в ARM   Jan 14 2015, 14:36
- - scifi   Цитата(cyrax0 @ Jan 14 2015, 17:36) Я исп...   Jan 14 2015, 14:55
- - cyrax0   А можно поподробнее, что за errno? Отмечу еще раз,...   Jan 14 2015, 17:02
|- - scifi   Цитата(cyrax0 @ Jan 14 2015, 20:02) А мож...   Jan 14 2015, 19:17
- - _Pasha   Насколько я помню, Вам нужно переопределить хандле...   Jan 14 2015, 19:15
|- - cyrax0   ЦитатаНасколько я помню, Вам нужно переопределить ...   Jan 14 2015, 19:43
|- - jcxz   Цитата(cyrax0 @ Jan 15 2015, 01:43) Если ...   Jan 15 2015, 05:26
|- - scifi   Цитата(jcxz @ Jan 15 2015, 08:26) Через M...   Jan 15 2015, 06:57
- - _Pasha   Собственно можно пойти далеко в обход: написать ну...   Jan 15 2015, 04:32
- - cyrax0   Не очень понимаю, с какой стати errno может быть в...   Jan 15 2015, 08:35
|- - jcxz   Цитата(cyrax0 @ Jan 15 2015, 14:35) Не оч...   Jan 15 2015, 10:14
- - scifi   Цитата(cyrax0 @ Jan 15 2015, 10:59) У нас...   Jan 15 2015, 09:08
|- - jcxz   Цитата(scifi @ Jan 15 2015, 15:08) Не фак...   Jan 15 2015, 10:27
- - cyrax0   ЦитатаА если привлечь Си++, то можно и операторы п...   Jan 15 2015, 09:25
|- - scifi   Цитата(cyrax0 @ Jan 15 2015, 12:25) В ито...   Jan 15 2015, 10:22
|- - klen   здравствуйте. вставлю свои 5 копеек. 1. FPU вычисл...   Jan 16 2015, 07:25
- - Xenia   Господа! Вы что, с ума сошли? Переменную errn...   Jan 15 2015, 12:28
|- - scifi   Цитата(Xenia @ Jan 15 2015, 15:28) Одного...   Jan 15 2015, 13:01
||- - Xenia   Цитата(scifi @ Jan 15 2015, 16:01) Конечн...   Jan 15 2015, 13:29
||- - scifi   Цитата(Xenia @ Jan 15 2015, 16:29) Тогда ...   Jan 15 2015, 13:47
||- - cyrax0   Цитата(scifi @ Jan 15 2015, 17:47) Подозр...   Jan 15 2015, 13:53
||- - scifi   Цитата(cyrax0 @ Jan 15 2015, 16:53) Делен...   Jan 15 2015, 14:51
|- - jcxz   Цитата(Xenia @ Jan 15 2015, 18:28) Господ...   Jan 15 2015, 17:08
- - cyrax0   Цитата(scifi @ Jan 15 2015, 14:22) Там ес...   Jan 15 2015, 13:46
- - cyrax0   Цитата(jcxz @ Jan 15 2015, 21:08) Так пом...   Jan 16 2015, 06:09
|- - jcxz   Цитата(cyrax0 @ Jan 16 2015, 12:09) Я поп...   Jan 16 2015, 17:55
- - WitFed   Ещё 1 коп: обычно можно переопределить "систе...   Jan 16 2015, 09:19


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 21st July 2025 - 19:29
Рейтинг@Mail.ru


Страница сгенерированна за 0.01416 секунд с 7
ELECTRONIX ©2004-2016