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

 
 
 
Reply to this topicStart new topic
> Перемножить два двоичных вектора С=A^B, Как это эффективно сделать для ARM
sergvks
сообщение Jul 26 2008, 05:30
Сообщение #1


Местный
***

Группа: Свой
Сообщений: 251
Регистрация: 26-07-05
Пользователь №: 7 117



Стоит задача перемножить в GF(2) два двоичных вектора, которые хранятся в двух массивах 32-битных слов. Вопрос заключается в том, как оптимизировать код под ARM7,ARM9 и CORTEX M3 с учётом их особенностей, возможно, с использованием ассемблера.

Пока в голову пришло вот это:

unsigned int multAB(unsigned int *A,unsigned int *B, unsigned int size)
{
unsigned int y=0;
for (i=0; i<size; i++)y^=A[i]^B[i];
y=y^(y>>1);
y=y^(y>>2);
y=y^(y>>4);
y=y^(y>>8);
y=y^(y>>16);
return y;
}
Go to the top of the page
 
+Quote Post
sergvks
сообщение Jul 27 2008, 08:57
Сообщение #2


Местный
***

Группа: Свой
Сообщений: 251
Регистрация: 26-07-05
Пользователь №: 7 117



Увидел пару неточностей smile.gif
Вот так правильно:
unsigned int multAB(unsigned int *A,unsigned int *B, unsigned int size)
{
unsigned int y=0,n;
for (i=0; i<size; i++)y^=A[i]&B[i];
y=y^(y>>1);
y=y^(y>>2);
y=y^(y>>4);
y=y^(y>>8);
y=y^(y>>16);
return y&1;
}
Go to the top of the page
 
+Quote Post
zltigo
сообщение Jul 27 2008, 09:44
Сообщение #3


Гуру
******

Группа: Свой
Сообщений: 13 372
Регистрация: 27-11-04
Из: Riga, Latvia
Пользователь №: 1 244



Цитата(sergvks @ Jul 27 2008, 10:57) *
Увидел пару неточностей smile.gif

И еще одна осталась sad.gif
unsigned int y=0,n;
надо:
unsigned int y=0;
unsigned int i;

А оптимизировать в перемножении векторов произвольного размера особо нечего sad.gif. Едиственно, написать цикл по человечески - практически наверняка будет покороче и в цикле, и при инициализации цикла:
Код
    do{
        y ^= *A++ & *B++;
    }
    while( --size );


--------------------
Feci, quod potui, faciant meliora potentes
Go to the top of the page
 
+Quote Post
Alex03
сообщение Jul 28 2008, 09:26
Сообщение #4


Местный
***

Группа: Свой
Сообщений: 359
Регистрация: 9-12-05
Пользователь №: 12 034



Цитата(zltigo @ Jul 27 2008, 15:44) *
Едиственно, написать цикл по человечески - практически наверняка будет покороче и в цикле, и при инициализации цикла: ....

А не факт... Зависит от возможностей адресации проца... Проверять надо, плюс у sergvks предварительная проверка есть на size==0.

sergvks Я людь тёмный, GF(2) двоичных вектороф не знаю, что в результате то надо, признак нечётного количества попарно единичных битов в исходных массивах?
Go to the top of the page
 
+Quote Post
zltigo
сообщение Jul 28 2008, 10:00
Сообщение #5


Гуру
******

Группа: Свой
Сообщений: 13 372
Регистрация: 27-11-04
Из: Riga, Latvia
Пользователь №: 1 244



Цитата(Alex03 @ Jul 28 2008, 11:26) *
А не факт...

Слова "практически наверняка" видели?
Цитата
Зависит от возможностей адресации проца...

Речь идет про АРМ, а зависит в данном случае в первую очередь от способностей компилятора к оптимизации. Но в общем случае написание "по-сишному", без явного указания лишних переменных и лишних (зачем перемножать два отсутствующих вектора) проверок в if() благотворно скажется на результате.
Цитата
Проверять надо...

Проверьте, если интересно, на все про все минут 15 smile.gif работы. Думаю, что одну команду в цикле и 3-4 вне цикла сэкономить реально.


--------------------
Feci, quod potui, faciant meliora potentes
Go to the top of the page
 
+Quote Post
sergvks
сообщение Jul 29 2008, 05:34
Сообщение #6


Местный
***

Группа: Свой
Сообщений: 251
Регистрация: 26-07-05
Пользователь №: 7 117



Цитата(Alex03 @ Jul 28 2008, 13:26) *
sergvks Я людь тёмный, GF(2) двоичных вектороф не знаю, что в результате то надо, признак нечётного количества попарно единичных битов в исходных массивах?

Ага, возможно и так сказать.

Я поэкспериментировал с асмом, но лучше чем у кейла у меня не получилось, так что наверно оставлю всё в таком виде.
Go to the top of the page
 
+Quote Post

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

 


RSS Текстовая версия Сейчас: 20th July 2025 - 09:51
Рейтинг@Mail.ru


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