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

 
 
> Как правильно в плис делить
sergey sva
сообщение Jul 15 2015, 03:26
Сообщение #1


Гуру
******

Группа: Свой
Сообщений: 2 546
Регистрация: 23-05-07
Из: Самарская область Сызрань
Пользователь №: 27 923



Как правильно в плис делить, очень медленная операция, на высокой частоте появляются слаки. Если сдвигом то можно поделить на 2 4 8 16.... а если нужно поделить на 9 или 5?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
sergey sva
сообщение Jul 15 2015, 06:32
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 2 546
Регистрация: 23-05-07
Из: Самарская область Сызрань
Пользователь №: 27 923



Да очень благодарен Maverick по ссылкам есть то что искал раньше.
Как оптимизировать что бы получить min latency max FREQ ?
Битовый сдвиг выполняется за одну операцию, или на каждый бит по такту?
Go to the top of the page
 
+Quote Post
Maverick
сообщение Jul 15 2015, 06:49
Сообщение #3


я только учусь...
******

Группа: Модераторы
Сообщений: 3 447
Регистрация: 29-01-07
Из: Украина
Пользователь №: 24 839



Цитата(sergey sva @ Jul 15 2015, 09:32) *


быстрее алгритм Radix 4 division

CODE

#include <stdio.h>
#include <stdlib.h>

/* This code assumes that unsigned short is 16 bits, and unsigned int is 32 bits */

unsigned short radix4_division_u16(unsigned short divisor, unsigned short dividend)
{
unsigned short quotient = 0;
unsigned char i = 16; /* Number of bits to process */

if(dividend == 0)
{
printf("Divide by zero");
exit(0);
}

while(i > 0)
{
unsigned char j;
unsigned int d1,d2,d3;

i -= 2;

d1 = dividend << i;
d2 = dividend << (i+1);
d3 = d1+d2;

if(divisor < d1)
j = 0;
else if(divisor < d2) {
j = 1;
divisor -= d1;
} else if(divisor < d3) {
j = 2;
divisor -= d2;
} else {
j = 3;
divisor -= d3;
}
quotient = (quotient << 2)+j;
}
return quotient;
}

int main(int c, char *v[])
{
int i,j;
for(i = 0; i < 256*256; i++)
{
for(j = 1; j < 256*256; j++)
{
unsigned k = radix4_division_u16(i,j);
if(i/j != k)
{
printf("Error with %i/%i != %i\n",i,j,k);
exit(0);
}
}
if(i%650 == 0)
printf("%2.2f%% tested\n",i*100.0/256/256);
}
return 0;
}


алгоритм по шагам:

CODE
Each step

Inputs:

Existing quotient
Current dividend
Divisor
Clock signal
Bit offset of bits to generate

Constant (generic) information also needed:

Width of Divisor
Width of Quotient
Bit offset of bits to generate

Outputs:

updated quotient
updated dividend
Divisor

It is possible to reduce the length of comparisons - rather than comparing n2-2 bits only 4 bits need subtraction - the lowwer order bits can just be ORed together.

Equations for first step of a/b:

Comparing against divisor
x1 <= (0 & a(n-1 ... n-2)))
y1 <= (OR(b(n-1 ... 2)) & b(1 ... 0))
out1 <= (x-y)(1 ... 0) & a(n-3 ... 0)

Comparing against 2*divisor
x2 <= (0 & a(n-1))
y2 <= (OR(b(n-1 ... 1)) & b(0))
out2 <= (x-y)(0) & a(n-2 ... 0)

Comparing against 3*divisor
sum <= (0 & b(1 ... 0)) + (0 & b(0) & 0)
x3 <= (0 & a(n-1) & a(n-2))
y3 <= (OR(d(n-1 ... 1),sum(2)) & sum(1 ... 0))
out3 <= (x-y)(1 ... 0) & a(n-3 ... 0)

which 'out' is passed to the next stage is selected based on which of out1, out2 or out3 were positive (if any).

If you have n bits, you need 'n/2' stages, and if you number stages from (n/2-1) (dealing to the high bits) to 0 (dealing to bits 1 and 0) the equations become more generic:

For stage i (need to verify!):

Comparing against divisor
x1 <= (0 & a(n-1 ... 2i) )
y1 <= (OR(b(n-1 ... n-2i)) & b(n-2i ... 0))
out1 <= (x-y)(n/2-i-1 ... 0) & a(2i-1 ... 0)

Comparing against 2*divisor
x2 <= (0 & a(n-1 ... 2i+1) )
y2 <= (OR(b(n-1 ... n-2i+1)) & b(0))
out2 <= (x-y)(n/2-2i-2..0) & a(2i-1 ... 0)

Comparing against 3*divisor
sum <= (0 & b(n-2i-1 .. 0)) + (0 & b(n-2i-2 ... 0) & 0)
x3 <= (0 & a(n-1 ... 2i))
y3 <= (OR(d(n-1 ... n-2i),sum(2)) & sum(1 ... 0))
out3 <= (x-y)(n/2-i-1 ... 0) & a(2i-1 ... 0)


PS насчет скорости и кол-во ресурсов не знаю не реализовывал, но вроде считает раза в 2 быстрее (если не изменяет память) предложенного ранее реализации

PS PS есть Radix более высокого прядка , например Radix8

PS PS PS алгоритмы деления на константу - рекомендую к просмотру


--------------------
If it doesn't work in simulation, it won't work on the board.

"Ты живешь в своих поступках, а не в теле. Ты — это твои действия, и нет другого тебя" Антуан де Сент-Экзюпери повесть "Маленький принц"
Go to the top of the page
 
+Quote Post
Fat Robot
сообщение Jul 15 2015, 07:10
Сообщение #4


ʕʘ̅͜ʘ̅ʔ
*****

Группа: Свой
Сообщений: 1 008
Регистрация: 3-05-05
Пользователь №: 4 691



Метод Ньютона-Рафсона и его реализация: Goldschmidt division

The Goldschmidt method is used in AMD Athlon CPUs and later models.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- sergey sva   Как правильно в плис делить   Jul 15 2015, 03:26
- - Barktail   Мне вот если делители известны заранее нравится ум...   Jul 15 2015, 03:36
- - Golikov A.   Делить правильно делителем. Видов делителей велик...   Jul 15 2015, 03:59
- - sergey sva   Вариантов много, но нужно еще и быстрый. То что си...   Jul 15 2015, 04:15
|- - des00   Цитата(sergey sva @ Jul 15 2015, 12:15) В...   Jul 15 2015, 05:05
- - dima32rus   А что если задействовать аппаратный делитель? У Ал...   Jul 15 2015, 04:43
- - Maverick   Цитата(sergey sva @ Jul 15 2015, 06:26) К...   Jul 15 2015, 05:04
- - sergey sva   Сейчас попробую, у циклона 3 есть dsp, если их зад...   Jul 15 2015, 05:18
|- - Maverick   Цитата(sergey sva @ Jul 15 2015, 08:18) С...   Jul 15 2015, 05:22
- - sergey sva   Если на нем сделать сдвиг, что бы сдвиг выполнился...   Jul 15 2015, 05:42
|- - Maverick   Цитата(sergey sva @ Jul 15 2015, 08:42) Е...   Jul 15 2015, 05:45
- - sergey sva   Благодарю за информацию, ссылок много сейчас смотр...   Jul 15 2015, 06:14
- - Kuzmi4   2 Maverick я так понял, что человеку надо min late...   Jul 15 2015, 06:16
- - serjj   И снова CORDIC стр 6 документа - умножение с нако...   Jul 15 2015, 06:39
- - likeasm   Цитата(Maverick @ Jul 15 2015, 09:22) зач...   Jul 15 2015, 12:36
|- - Kuzmi4   Цитата(likeasm @ Jul 15 2015, 14:36) Можн...   Jul 15 2015, 14:04
|- - D Mike   если делимое не более 12 разрядов, делитель=конста...   Aug 18 2015, 09:43
- - Golikov A.   да правда быстро, особенно если знаешь сколько буд...   Aug 18 2015, 11:14


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

 


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


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