Вычисление числа Пи на ассемблере

[b2m]

Вот из вышеупомянутой книги, как ни странно - работает

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

[Advertiser]

[ivagor]

INX SP INX SP

Подобные вещи и в pirk катят, но убрал, чтобы работало и при разрешенных прерываниях (там где они есть).

- - - Добавлено - - -

Это я к тому, что желательно этот фрагмент потом заменить

- - - Добавлено - - -

желательно этот фрагмент потом заменить

Пардон, здесь проблем не будет (у меня несколько другая ситуация была).

[Vslav]

У меня пока нет оснований не доверять тому архиву, т.ч. с 45000 я не согласен.

Я тоже не согласен с 45000
Взял Пи вот тут
и сравнил с вычисленным по 4-цифровому шпигату (по ссылке что давал litwr и в моем листинге), оно навернулось на 16278 цифре - начиная с нее не совпадает с эталонными цифрами. То есть, найденная мной ошибка - там не единственная.

Update:
там банальное целочисленное переполнение происходит, то есть - не совсем ошибка, просто неприменимость в данных условиях. Поменял типы на 64-битные - доехало до 50К, потом опять ERROR

Update2:
Правильная последовательность с 54920-ой цифры:
61527146690058147000026330
Вычисленная по 4-цифровому шпигату (с int64):
61527146690058146000026330
То есть - перенос 10000 не отработал, ошибка в наличии.

[ivagor]

Взял Пи вот тут

Скопировал оттуда, сравнил с ранее полученным результатом

правильно считает до 54932 знака включительно

- эти самые 54932 знака совпали

[Vslav]

- эти самые 54932 знака совпали

Да, а потом встречается перенос за четыре десятичных разряда и он игнорится, на 57936-ой позиции имеем липовую цифру 6 вместо верной 7, я заслужил с полки пирожок

[b2m]

Нехватает однако 16*16=32 и 32/16=16. Наткнулся на итерацию, в которой после деления 17 бит получается, дальше можно не проверять. Нужны нормальные 32-битные процедуры, но тут регистров уже не хватит. Да и на первый взгляд, медленновато пока получается (пусть даже и неверный результат). Точно не засекал, но даже на глаз - медленнее. Так что первый шпигат - для восьмибиток пока лучше.

[ivagor]

Я так понимаю, что нужна процедура переваривающая 32 битные: делимое, делитель и частное. Остаток м.б. 16 битным

[litwr]

Cделал код для 6502. На Коммодоре+4 800 знаков отсчитало за 6 мин 32 сек. Там по cсылке есть указания на алгоритм, по которому миллиарды знаков отсчитал некто Чудновский, с 32-разрядными базовыми константами. Но кода не искал, может где и есть в сети.

[ivagor]

Для затравки процедура 32 битного деления для 8080 - div32vslow.zip
Все 32битное - делитель, делимое, частное и остаток. Очень медленно. Единственный плюс - легко масштабировать хоть до 256 бит.

[b2m]

С делением я разобрался, доделал до 32-битного частного и 16-битного остатка. Сделал так: если надо, деление будет выполнено дважды (для старших двух байт и для младших двух байт). Умножение оставил как есть, но сделал цикл сложением для старшего байта 24-битных чисел. Вобщем, пока считает 100 знаков за 12,5 сек.