Классический алгоритм деления без восстановления остатка реализован на 8-разрядном контроллере PIC16F690 (код написан на ассемблере). Делится 32-разрядное число на 16-разрядное. При некоторых значениях аргументов программа выдаёт неправильные результаты. Проверил на бумажке - получаю то же самое.
Например, деление 0x8000 на 0x1AAE даёт результат 5, а должно получаться 4.
Что я делал (и, соответственно, делает проверяемая программа):
1. Нормализуем делимое и делитель, т.е. сдвигаем влево до получения 1 в старшем разряде. При этом счётчик сдвигов инкрементируется при сдвиге делимого и декрементируется при сдвиге делителя.
В итоге делимое - 0x80000000, делитель - 0xD570 счётчик сдвигов - ZZ=0x0D.
2. Вычитаем делитель из старших разрядов делимого и инкрементируем счётчик сдвигов:
0x8000-0xВ570=0xAA90, ZZ=0x0E, C=0 => следующая операция - сложение делимого и удвоенного остатка (как я понимаю, старший бит при сдвиге делимого влево не учитывается? или именно тут я и ошибаюсь?)
2*0xAA90=0x5520, 0x5520+0xD570=0x2A90, ZZ=0x0F, C=1 => следующая операция - вычитание делимого из удвоенного остатка.
2*0x2A90=0x5520, 0x5520-D570=0x7FB0, ZZ=0x10, C=0 => следующая операция - сложение делимого и удвоенного остатка.
2*0x7FB0=0xFF60, 0xFF60+0xD570=0xD4D0, ZZ=0x11 (признак окончания деления), C=1.
Итак, получается результат 0x0101. А виндовый калькулятор и простая логика даёт результат 0x0100.
Где тут грабли?