[an error occurred while processing this directive]
А вот еще последний гвоздь (+)
(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)
|
Отправлено
homekvn 31 октября 2006 г. 19:27
В ответ на: Ответ (+) отправлено
<font color=gray>homekvn</font> 31 октября 2006 г. 19:14
|
 |
|
|
Если Вы начнете суммировать квадраты "в лоб" в фиксед пойнте, предварительно их не пронормировав, то вы рискуете получить переполнение. Например, пусть у Вас есть вектор:
x = [2^15 - 1, -2^15, 1, -2^15 , 2^15 - 1, -2^15];
(Я здесь умалчиваю, что данные Ваши могут и превосходить 16 битную сетку). Попробуйте отыскать средний квадрат их "в лоб". Ваше счастье, если длина вектора кратна степени двойки N=2^n. В этом случае, прежде чем возводить в квадрат и суммировать, Вы сдвините каждое число на n бит вправо (это правильный алгоритм для фиксед-пойнт процессора; вот тут-то малые величины и улетучатся). А в случае произвольной длины, Вы должны будете каждое значение разделить на N, что, согласитесь не добавляет привлекательности Блекфина для решения этой задачи.
Составить ответ
|||
Конференция
|||
Архив
Ответы
- Вообще-то я ничего не сдвигаю - просто складываю раздельно сумму и число отсчетов... (+) — ASergej_R19 (31.10.2006 19:40 80.250.160.170, 165 байт)
- Уж не знаю как там у блэкфина, но делить каждое число перед возведением в квадрат - это какое-то особенный шик, больше напоминающий извращение. Вот Вам навскидку метОд. Отдельно суммируете старшие 16-ти разрядные чапсти квадратов, отдельно младшие. Потом суммы сшиваете и результат делите. Или сначала делите, потом сшиваете, пофигу. И гвоздь Ваш совсем ржавый в результате получился. — -=ВН=- (31.10.2006 19:36 193.125.71.140, пустое)
- Можно эту операцию расписать для того примера, что я привел ранее (+) — homekvn (31.10.2006 19:53 212.185.161.237, 552 байт)
- А я Вам ответил, что точный результат получается за время ~2*N. И если нужно деление (считается средний квадрат), то делится означенный точный результат, а вовсе не каждый операнд. И годится это дело для любого вектора, хорошего ли, плохого ли. 2*N гораздо хуже, чем N, гораздо, но много лучше деления каждого отсчета, много лучше. — -=ВН=- (31.10.2006 20:02 193.125.71.140, пустое)
- А никто и не говорил, что у Блекфина производительность больше чем у Шарка... (+) — ASergej_R19 (31.10.2006 20:02 80.250.160.170, 222 байт)
- А если использовать оба аккумулятора A0 и A1, то посчитает за время вдвое меньшее.. — quark (31.10.2006 20:06 62.140.241.223, пустое)
- С последним утверждением я не спорю и не возражал против него. Просто пример Вы привели, на мой взгляд, не очень удачный для демонстрации этого утверждения. Вот там-то я и завозражал. — homekvn (31.10.2006 20:05 212.185.161.237, пустое)
- А я начал возражать Вам, потому что Вы возражали неправильным применением fixed point... :-) (+) — ASergej_R19 (31.10.2006 20:21 80.250.160.170, 398 байт)
- Ну, про точность я готов продолжить. Приведите, если Вам не трудно, пример наихудшего для floating-point вычислений вектора, где бы начались проблемы с точностью, но при этом 32-битный fixed-point подобных проблем бы не имел. — homekvn (31.10.2006 20:29 212.185.161.237, пустое)
- Не в состоянии прочитать всю дискуссию, но сейчас особенно моден формат ieee754. Каковой для флоата обычного оговаривает на все прол все 32-х разряда. С учетом этого сложите во float 2 числа и доложите результат. Первое число=1073741823, второе =1073741822. — -=ВН=- (31.10.2006 20:38 193.125.71.140, пустое)
- Ну-у! Так и я могу. А теперь Вы сложите серию чисел вроде 8.21e+48 или (e-48), которые также соответствуют приведенному Вами стандарту. И потом Вы мне исходные числа приведите, которые надо в квадрат возвести, а потом просуммировать. А потом сами на блекфине сделайте то же самое. — homekvn (31.10.2006 21:08 212.185.161.237, пустое)
- Вы подменили вопрос. 8.21e+48 не входят в ДИАПАЗОН 32-х разрядных чисел целых. Не входят. А приведенные мной числа в ДИАПАЗОН флоата одинарной точности входят. Я ведь не просил Вас складывать в одинарном флоате числа порядка 8.21e+16799. Нехорошо, нехорошо, Вы пытаетесь меня обмануть :-((( — -=ВН=- (31.10.2006 21:18 193.125.71.140, пустое)
- Да признаюсь, перебрал. Докладываю результат 3,9999997 И что? Кстати (+) — homekvn (31.10.2006 21:41 212.185.161.237, 635 байт)
- Да где Вы видели, что кто-то Ваш Шарк принижал? :-) Это ж разные процессоры, заточены под разные вещи... (+) — ASergej_R19 (31.10.2006 22:05 80.250.160.170, 339 байт)
- Вы знаете, я как-то догадывался, что Шарик и с целыми числами умеет, я даже с шариком работал, правда это было давно, шарик старый был, 21065, подох наверное уже. И мне в общем по барабану, кто кого сделает, шарик блэкфина или наоборот. Вы задаете вопрос или приводите пример, я Вам в ответ на это и только на это и демонстрирую всякие ноу и даже хау. — -=ВН=- (31.10.2006 21:48 193.125.71.140, пустое)
- Отлично! :-) (+) — ASergej_R19 (31.10.2006 20:42 80.250.160.170, 289 байт)
- Вот кстати и так тоже можно... (-) — ASergej_R19 (31.10.2006 19:41 80.250.160.170, пустое)
- В силу последнего замечания для нахождения RMS Блекфин даже и точнее-то ну никак не будет (ввиду обязательной нормировки входного вектора). — homekvn (31.10.2006 19:33 212.185.161.237, пустое)
Перейти к списку ответов
|||
Конференция
|||
Архив
|||
Главная страница
|||
Содержание
