Вид для печати
Вот из вышеупомянутой книги, как ни странно - работает :)
Код:DIV32: MOV A,B ; DE = HLDE/BC, HL = HLDE%BC
CMA
MOV B,A
MOV A,C
CMA
MOV C,A
INX B
XRA A
DIV321: DAD H
RAR
XCHG
DAD H
XCHG
JNC $+4
INX H
PUSH H
DAD B
JNC DIV322
RAL
DIV323: INX D
INX SP
INX SP
ADI 10h
JNC DIV321
RET
DIV322: RAL
JC DIV323
POP H
ADI 10h
JNC DIV321
RET
Подобные вещи и в pirk катят, но убрал, чтобы работало и при разрешенных прерываниях (там где они есть).Цитата:
Сообщение от b2m
- - - Добавлено - - -
Это я к тому, что желательно этот фрагмент потом заменить
- - - Добавлено - - -
Пардон, здесь проблем не будет (у меня несколько другая ситуация была).Цитата:
Сообщение от ivagor
Я тоже не согласен с 45000 :biggrin:Цитата:
Сообщение от ivagor
Взял Пи вот тут
и сравнил с вычисленным по 4-цифровому шпигату (по ссылке что давал litwr и в моем листинге), оно навернулось на 16278 цифре - начиная с нее не совпадает с эталонными цифрами. То есть, найденная мной ошибка - там не единственная.
Update:
там банальное целочисленное переполнение происходит, то есть - не совсем ошибка, просто неприменимость в данных условиях. Поменял типы на 64-битные - доехало до 50К, потом опять ERROR
Update2:
Правильная последовательность с 54920-ой цифры:
61527146690058147000026330
Вычисленная по 4-цифровому шпигату (с int64):
61527146690058146000026330
То есть - перенос 10000 не отработал, ошибка в наличии.
Скопировал оттуда, сравнил с ранее полученным результатомЦитата:
Сообщение от Vslav
- эти самые 54932 знака совпалиЦитата:
Сообщение от ivagor
Да, а потом встречается перенос за четыре десятичных разряда и он игнорится, на 57936-ой позиции имеем липовую цифру 6 вместо верной 7, я заслужил с полки пирожок :)Цитата:
Сообщение от ivagor
Нехватает однако 16*16=32 и 32/16=16. Наткнулся на итерацию, в которой после деления 17 бит получается, дальше можно не проверять. Нужны нормальные 32-битные процедуры, но тут регистров уже не хватит. Да и на первый взгляд, медленновато пока получается (пусть даже и неверный результат). Точно не засекал, но даже на глаз - медленнее. Так что первый шпигат - для восьмибиток пока лучше.
Я так понимаю, что нужна процедура переваривающая 32 битные: делимое, делитель и частное. Остаток м.б. 16 битным
Cделал код для 6502. На Коммодоре+4 800 знаков отсчитало за 6 мин 32 сек. Там по cсылке есть указания на алгоритм, по которому миллиарды знаков отсчитал некто Чудновский, с 32-разрядными базовыми константами. Но кода не искал, может где и есть в сети.
Для затравки процедура 32 битного деления для 8080 - Вложение 54857
Все 32битное - делитель, делимое, частное и остаток. Очень медленно. Единственный плюс - легко масштабировать хоть до 256 бит.
С делением я разобрался, доделал до 32-битного частного и 16-битного остатка. Сделал так: если надо, деление будет выполнено дважды (для старших двух байт и для младших двух байт). Умножение оставил как есть, но сделал цикл сложением для старшего байта 24-битных чисел. Вобщем, пока считает 100 знаков за 12,5 сек.