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

 
 
 
Reply to this topicStart new topic
> Деление с остатком
d7d1cd
сообщение Jun 25 2013, 18:44
Сообщение #1


Местный
***

Группа: Участник
Сообщений: 442
Регистрация: 26-11-10
Пользователь №: 61 199



Привет всем. У меня снова вопрос по написанию программы для MSP430 на ассемблере.

Для выполнения определенных задач мне необходимо 32-х разрядное число умножить на 8-ми разрядное, результат разделить на 100, получить результат, а так же остаток от деления. Для умножения я использую аппаратный умножитель. А вот с делением проблема.

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

Однако данный алгоритм деления долгий. Подскажите, может можно быстрее выполнить деление на 100 с остатком?

Go to the top of the page
 
+Quote Post
kovigor
сообщение Jun 25 2013, 20:00
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 5 273
Регистрация: 30-03-10
Пользователь №: 56 295



Цитата(d7d1cd @ Jun 25 2013, 21:44) *
Подскажите, может можно быстрее выполнить деление на 100 с остатком?

Так ли уж необходимо делить на 100 ? Если это для индикации, то просто сместите запятую влево на два разряда. А если без деления никак, то у Atmel есть апп. ноут:

http://www.atmel.com/Images/doc0936.pdf
http://zalil.ru/34605572

Кстати, а почему бы такие вещи на Си не написать ? Зачем вы используете ассемблер ?
Go to the top of the page
 
+Quote Post
rezident
сообщение Jun 25 2013, 21:00
Сообщение #3


Гуру
******

Группа: Свой
Сообщений: 10 920
Регистрация: 5-04-05
Пользователь №: 3 882



Цитата(d7d1cd @ Jun 25 2013, 23:44) *
Для выполнения определенных задач мне необходимо 32-х разрядное число умножить на 8-ми разрядное, результат разделить на 100, получить результат, а так же остаток от деления.
Определитесь более конкретно с диапазонами чисел и требуемой точностью результата деления. Если умножить любое большое 32-х разрядное число на 8-ми разрядное, то получится переполнение 32-х разрядной сетки. Значит нужно считать в 64-х разрядных числах?
Есть способ быстрого деления на константы с использованием умножения на "магические" числа и последующего деления с помощью операций сдвига. Но для подбора "магического" числа нужны критерии оглашенные выше. Пока могу только намекнуть.
Математическое преобразование X/100 = X * 1/100 является тождественным
А X * 1/100 ≈ X * 655/65536 = (X * 655) >> 16 - приближенным с заданной точностью (около -0,054%). В нем присутствует одно умножение и сдвиг на 16 бит вправо, либо в качестве результата деления можно просто взять старшее 16-и разрядное слово из 32-х разрядного числа.
Go to the top of the page
 
+Quote Post
rx3apf
сообщение Jun 25 2013, 21:36
Сообщение #4


Гуру
******

Группа: Участник
Сообщений: 3 834
Регистрация: 14-06-06
Из: Moscow, Russia
Пользователь №: 18 047



Цитата(d7d1cd @ Jun 25 2013, 22:44) *
Однако данный алгоритм деления долгий. Подскажите, может можно быстрее выполнить деление на 100 с остатком?

А когда-то даже в школе учили делить "в столбик" wink.gif

Есть делимое, делитель, выделим регистры под остаток (в данном случае нужно 8 битов, так что один регистр, и тот наполовину), обнулим их. Сдвигаем влево делимое (младший бит обнуляем), старшие биты уползают в остаток. И пробуем из этого остатка вычесть делитель. Получилось ? Запишем в младший бит делимого (там же набирается частное) "1". Не получилось - оставляем 0. И так 32 раза. Частное будет там, где было делимое, остаток там, куда сдвигали.

Сообщение отредактировал rx3apf - Jun 25 2013, 21:38
Go to the top of the page
 
+Quote Post
d7d1cd
сообщение Jun 26 2013, 06:13
Сообщение #5


Местный
***

Группа: Участник
Сообщений: 442
Регистрация: 26-11-10
Пользователь №: 61 199



Цитата(kovigor @ Jun 26 2013, 00:00) *
Так ли уж необходимо делить на 100 ? Если это для индикации, то просто сместите запятую влево на два разряда. А если без деления никак, то у Atmel есть апп. ноут:

http://www.atmel.com/Images/doc0936.pdf
http://zalil.ru/34605572

Кстати, а почему бы такие вещи на Си не написать ? Зачем вы используете ассемблер ?

Да, именно необходим результат (целая часть и остаток) от деления на 100. Это не для индикации. Причину выбора ассемблера, если Вам интересно, могу написать в ЛС.

Цитата(rezident @ Jun 26 2013, 01:00) *
Определитесь более конкретно с диапазонами чисел и требуемой точностью результата деления. Если умножить любое большое 32-х разрядное число на 8-ми разрядное, то получится переполнение 32-х разрядной сетки. Значит нужно считать в 64-х разрядных числах?
Есть способ быстрого деления на константы с использованием умножения на "магические" числа и последующего деления с помощью операций сдвига. Но для подбора "магического" числа нужны критерии оглашенные выше. Пока могу только намекнуть.
Математическое преобразование X/100 = X * 1/100 является тождественным
А X * 1/100 ≈ X * 655/65536 = (X * 655) >> 16 - приближенным с заданной точностью (около -0,054%). В нем присутствует одно умножение и сдвиг на 16 бит вправо, либо в качестве результата деления можно просто взять старшее 16-и разрядное слово из 32-х разрядного числа.

Достоверно известно, что результат умножения не будет превышать 32-х разрядов, так как исходное 32-х разрядное число на самом деле будет занимать не более 24-х разрядов. Запас разрядов нужен для хранения результата умножения.
Деление способом умножения на "магическое" число я рассматривал, но в этом случае будет потеря какой-то части результата, а мне необходимо учитывать все. Остаток от деления - идеальный вариант для моей задачи.

Цитата(rx3apf @ Jun 26 2013, 01:36) *
А когда-то даже в школе учили делить "в столбик" wink.gif

Есть делимое, делитель, выделим регистры под остаток (в данном случае нужно 8 битов, так что один регистр, и тот наполовину), обнулим их. Сдвигаем влево делимое (младший бит обнуляем), старшие биты уползают в остаток. И пробуем из этого остатка вычесть делитель. Получилось ? Запишем в младший бит делимого (там же набирается частное) "1". Не получилось - оставляем 0. И так 32 раза. Частное будет там, где было делимое, остаток там, куда сдвигали.

Я проходил в школе деление столбиком wink.gif Однако применить это знание в моей задаче не додумался.
Когда Вы пишете "пробуем из этого остатка вычесть делитель", то вы имеете ввиду, что вычитание "получается" если результат не отрицательный? И еще: если вычитание "получилось", то результат вычитания должен быть записан в регистр для остатка?
Go to the top of the page
 
+Quote Post
rx3apf
сообщение Jun 26 2013, 07:56
Сообщение #6


Гуру
******

Группа: Участник
Сообщений: 3 834
Регистрация: 14-06-06
Из: Moscow, Russia
Пользователь №: 18 047



Цитата(d7d1cd @ Jun 26 2013, 10:13) *
Когда Вы пишете "пробуем из этого остатка вычесть делитель", то вы имеете ввиду, что вычитание "получается" если результат не отрицательный?

Совершенно верно.
Цитата
И еще: если вычитание "получилось", то результат вычитания должен быть записан в регистр для остатка?

Да.
Go to the top of the page
 
+Quote Post
d7d1cd
сообщение Jun 26 2013, 09:37
Сообщение #7


Местный
***

Группа: Участник
Сообщений: 442
Регистрация: 26-11-10
Пользователь №: 61 199



rx3apf, спасибо за "школьный" алгоритм. Все работает как надо!
Go to the top of the page
 
+Quote Post
d7d1cd
сообщение Oct 29 2013, 17:35
Сообщение #8


Местный
***

Группа: Участник
Сообщений: 442
Регистрация: 26-11-10
Пользователь №: 61 199



Подскажите, а возможно ли быстро выполнить деление 16 разрядного числа на 3? Не используя принцип деления столбиком. Умножение на 3 сделать легко: Х*3 = Х*2+Х. А деление - это операция обратная умножению. Но я не могу додуматься. И возможно ли это?
Go to the top of the page
 
+Quote Post
_pv
сообщение Oct 29 2013, 17:51
Сообщение #9


Гуру
******

Группа: Свой
Сообщений: 2 563
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



при наличии аппаратного умножения делать его сдвигами будет только медленнее.
деления на константы можно заменять на умножения и сдвиги x/3 = x * 21845 / 65536 = (x * 21845) >> 16

Go to the top of the page
 
+Quote Post
◠◡◠
сообщение Oct 29 2013, 19:58
Сообщение #10





Группа: Новичок
Сообщений: 3
Регистрация: 28-10-13
Пользователь №: 78 944



вот быстрое деление на 10:
Код
u8  divmod10_u32_rem;
u32 divmod10_u32(u32 n) {
  u32 quot = n >> 1;
  quot  += quot >> 1;
  quot  += quot >> 4;
  quot  += quot >> 8;
  quot  += quot >> 16;
  u32 qq = quot & ~7ul;
  quot >>= 3;
  divmod10_u32_rem = n - ((quot << 1) + qq);
  if (divmod10_u32_rem > 9) {
    divmod10_u32_rem -= 10;
    quot++;
  }
  return quot;
}


на 100 - поделить два раза по 10.
Go to the top of the page
 
+Quote Post
_pv
сообщение Oct 30 2013, 09:16
Сообщение #11


Гуру
******

Группа: Свой
Сообщений: 2 563
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



унможить два 32х разрядных числа (на 2^32/10 = 429496730) и взять старшие 32 разряда (поделить на 2^32) будет всё равно быстрее, тем более что в новых 5 и 6 серии умножитель 32х разрядный.
Go to the top of the page
 
+Quote Post
Harvester
сообщение Oct 30 2013, 15:07
Сообщение #12


Местный
***

Группа: Участник
Сообщений: 338
Регистрация: 1-02-06
Из: Королев, М.О.
Пользователь №: 13 846



Цитата(d7d1cd @ Oct 29 2013, 21:35) *
Подскажите, а возможно ли быстро выполнить деление 16 разрядного числа на 3? Не используя принцип деления столбиком. Умножение на 3 сделать легко: Х*3 = Х*2+Х. А деление - это операция обратная умножению. Но я не могу додуматься. И возможно ли это?

x/3 = x/2 - x/4 + x/8 - x/16 + x/32 - x/64 + ...


--------------------
-Да как так-то?/-Да как-то так/-Ну так-то да
Go to the top of the page
 
+Quote Post
◠◡◠
сообщение Oct 31 2013, 07:50
Сообщение #13





Группа: Новичок
Сообщений: 3
Регистрация: 28-10-13
Пользователь №: 78 944



Цитата
Подскажите, а возможно ли быстро выполнить деление 16 разрядного числа на 3?
Код
U16 div_by3_U16_soft(U16 data_in) {
    U32 U32_01, result;
    result = U32_01 = data_in;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    result += 0x5555; //correction
    return(result >> 16);
}

U16 div_by3_U16_hmul(U16 data_in) {
    return((((U32)data_in * 0x5555) + 0x5555) >> 16);
}
Go to the top of the page
 
+Quote Post
d7d1cd
сообщение Oct 31 2013, 14:51
Сообщение #14


Местный
***

Группа: Участник
Сообщений: 442
Регистрация: 26-11-10
Пользователь №: 61 199



Цитата(◠◡◠ @ Oct 31 2013, 11:50) *
Код
U16 div_by3_U16_soft(U16 data_in) {
    U32 U32_01, result;
    result = U32_01 = data_in;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    U32_01 <<= 2; result += U32_01;
    result += 0x5555; //correction
    return(result >> 16);
}

U16 div_by3_U16_hmul(U16 data_in) {
    return((((U32)data_in * 0x5555) + 0x5555) >> 16);
}

Ну вообще то на С можно и просто написать result = data_in / 3;
Я же пишу на асме. Краткость и быстрота кода - вот цель!
Go to the top of the page
 
+Quote Post

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

 


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


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