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

[litwr]

Первое, что приходит в голову - деление. Тут и команды вычитания регистровых пар пригодятся (для ускорения восстановления остатка) и дополнительные регистры.

Слышал от знатоков z80, что быстрые программы IX, IY не используют. 16-битная разность почти бесполезна из-за отсутствия беспереносного вычитания, а ставить перенос на z80 морочно - проще складывать дополнительные коды. Система команд z80 определенный лидер конкурса некрасоты. Разве только второй набор регистров... Но их лучше возможно к прерываниям прикрепить.
Вот, например, деление на z80 - http://cpcwiki.eu/index.php/Programm...32bit_division - другого там раньше не было. Оно раз в 5-6 (!) медленнее DIV32 в программе для Радио-РК (pirk)! Кстати последний вариант pirk - это просто шедевр какой-то.
Пора посмотреть на сверхскоростной pi.rar... Но он, как понял, в отличие от pirk, 1000 знаков не считает? И где же он этот сверхскоростной на 7 сек?
Кстати программка на AMSTRAD CPC6128 с DIV32 и с умножением с таблицами на 16 KB сто знаков отсчитывает за 5.3 сек, 1000 соответственно 5.3*10^2 = 8м 50с - явно можно сделать получше, все же процессор примерно в два раза помощнее...
Возник вопрос. Зачем генерируете таблицы? Проще же готовыми их вставить или кто-то с кассеты грузит? ;-)

Для каждой точности есть свой вариант уточненной формулы Мэчина.

Эта формула вроде уже считается устаревшей, даже у Гаусса была лучше. А кроме того, в 90-е пооткрывали много поцифирных формул, сходящися по 10^n и лучше...
И кроме того, если соревнования между архитектурами, то зачем разные алгоритмы? Ясно, что система с худшим железом легко обгонит лучшую с негодным алгоритмом. Например, Радио-РК легко обгонит лучший суперкомпьютер, который будет считать число Фибоначчи по определению. А если объявить соревнования алгоритмов, то зачем тут ретроплатформы? Есть уже, например, http://benchmarksgame.alioth.debian.org/

[Advertiser]

[b2m]

И где же он этот сверхскоростной на 7 сек?

Вот вариант без таблиц, но он 8.89 сек. pi.rar

[ivagor]

Насчет переносов могу привести такой пример. Несколько лет назад с помощью рекомпилятора Tim0xи адаптнул распаковщик megalz с z80 на 8080. Потом оптимизировал его и все такое. А потом b2m по сишному распаковщику написал распаковщик для 8080, который (если не изменяет память) в разы быстрее. Да, там тоже можно пару байтов и пару тактов скинуть, но в разы его уже не улучшить. Это я к тому, что с си вполне можно написать нормальную программу для 8080, что уже не раз демонстрировал b2m, и, с другой стороны, с асма на асм можно довольно неудачно перенести.

[ivagor]

pi.rar - зачетно. Я собирался модифицировать первоначальный spigot подобным образом, а ты уже сделал.

[ivagor]

b2m ленится оптимизировать "по-мелкому", а мне было интересно. Заменил в последнем pi.rar умножение и деление:
100 цифр - 12868015 тактов - 7.24 сек
536 цифр (из них 535 верно) - 351958021 тактов - 3 мин 18 сек
Просто охренеть. И там еще есть резервы, я их не трогал, тупо заменил процедуры. А так явно можно 100 знаков меньше чем за 7 секунд посчитать. Как бороться с b2mом?

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

Уточню - это без учета времени расчета таблицы умножения

[balu_dark]

b2m ленится оптимизировать "по-мелкому", а мне было интересно. Как бороться с b2mом?

Самый простой вариант - добавить ему пурген в чай.

А вообше - очень интересная тема - с удовольствием слежу. пока нет возможности проверить или поправить на 8ми битке.

[b2m]

Как бороться с b2mом?

Исключительно личным примером

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

добавить ему пурген в чай.

Э... проверено на личном опыте?

[balu_dark]

Э... проверено на личном опыте?

Не Прочитано в книжках. Как ненасильственный метод борьбы Ибо согласно нашему уголовному кодексу - "публичные призывы к ...." малость наказуемы
А ivagor к тому моменту не уточнил - что он хочет побороть.

[ivagor]

Во избежание недоразумений заменю формулировку своего риторического вопроса на более корректную - Как соревноваться с b2mом? Для себя у меня есть отмазка - я (практически) не программист.
Интересно было бы увидеть реализацию для 8080 более эффективного алгоритма.

[ivagor]

Насчет выкладывания семисекундного+ варианта - я тут обещал больше не заниматься (по крайней мере в этой теме ) паразитической оптимизацией программ b2mа, так что выкладывать очередную версию было бы непоследовательно и нехорошо. Правильнее было бы даже и не писать про нее, но эмоции переполняли. При желании быструю версию можно собрать самостоятельно, нужно только заменить умножение и деление. Умножение можно взять из pirk20 (там еще можно таблицу расположить с 0000h и избавиться от лишнего inr). Деление проще взять из предыдущего pirk, в 20 уже немного иначе сделано. Если вдруг сам реализую какой-нибудь другой вариант расчета, то выложу, но это очень маловероятно.

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

Кстати, насчет реализации умножения и деления на z80 - не новая, но очень хорошая статья. Думаю про нее большинство знает, но может кто не читал. И, конечно, то, что там написано, вполне можно адаптировать для других процессоров.