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

[litwr]

Программа для z80 http://www.z88dk.org/wiki/doku.php?id=libnew:examplesi 800 знаков считает за несколько минут и написана на си.
Интересно, пригодятся ли при таком расчете команды z80, отсутствующие у 8080?
Пользователи БК жалко отсутствуют. Радио-РК БК определенно должна обогнать, а вот Спектрум скорее всего нет.
В связи с таким расчетом запустил на стандартном консольном калькуляторе на Raspberry Pi c 900 МГц
time echo 'scale=1000;16*a(1/5)-4*a(1/239)'|bc -l
т.е. 1000 знаков - отсчитано с 1.5 сек.
Обратите внимание на ПРАВИЛЬНУЮ ФОРМУЛУ - в материалах выше плюс перепутан с минусом.
Рожица после ссылки вставилась сама - движок не должен делать подстановки в ссылках!

[Advertiser]

[Titus]

Рожица после ссылки вставилась сама - движок не должен делать подстановки в ссылках!

Об этом надо написать в тему о переезде форума, тут это могут не заметить.

[perestoronin]

Программа для z80 http://www.z88dk.org/wiki/doku.php?id=libnew:examples

Благодарю за найденную ошибку, за ссылку на готовую библиотеку и за результаты для ARM, ARM конечно же бесспорный лидер

Опробуйте пожалуйста на z80 при 100 знаках точных цифр. Сколько секунд получилось и при какой частоте ?
Про БК, если программу переписать на ассемблере macro11 и использовать команду целочисленного умножения (есть такая у процессора КР1801ВМ1Г и его можно впаять вместо штатного КР1801ВМ1А), то еще вопрос кто вырвется в лидеры. Также интересен результат и для процессоров ВМ2, у которых более богатые команды целочисленного умножения и деления.
Может кто на эмуляторах БК попробует ?

https://en.wikipedia.org/wiki/Machin-like_formula
Согласно вики, в файле действительно ошибка в комментариях:

Сопоставлять результаты будем когда накопится с 10 различных платформ, но уже видно, то ARM у нас 32 разрядный, а не 8 разрядный и частота в 250 раз выше чем у zx. Т.е. в 1000 раз производительнее, и результ полученный на ARM при сравнении с zx нужно делить примерно на 1000.

[hobot]

движок не должен делать подстановки в ссылках!

offtop

Просто при расширенном режиме окна редактирования в низу поставьте галочку
в пункте

Отключить смайлы в тексте
При включении, :) НЕ будет заменён на :)

[свернуть]

[ivagor]

Почитал, насколько быстро считает сишная прога на z80 и оптимизировал умножение в pirk14. Теперь 535 знаков считает на полминуты быстрее (6 мин 33 сек), кроме того прога стала короче.

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

Переделал умножение на сложение двух произведений по таблице 5*8 (16 Кб). В итоге 535 знаков стал считать еще на 45 секунд быстрее - 5 мин 47 сек. Из них примерно 1 сек - расчет таблицы. Для 535 знаков это не критично, а для 100 знаков таблицу можно и заранее отдельно посчитать.

[ivagor]

Этого Мачина на асме 8080 я бы тоже с интересом посмотрел.
Насчет сравнения скорости spigot для разного числа знаков: 535 в 28.6 раз медленнее 100 (b2m приводил подобную оценку), 800 в 2.2 раза медленнее 535 (если не учитывать необходимости увеличения разрядности части операций). Небольшие отличия вносит вывод на экран. Возможно стоит разделить расчет и вывод на экран.

[Vslav]

Метод цифра за цифрой хорош когда есть жесткие требования к размеру доступной памяти

Что-то мне кажется, что метод:

/* A spigot algorithm for the Digits of \pi, Stanley Rabinowitz
and Stan Wagon, Amer.Math.Monthly, March 1995, 195-203 */

На базе которого в этой ветке пишут тесты, не является методом "цифра за цифрой". Там именно идет добавление очередного члена ряда к длинному числу, которое находится в массиве, при этом еще идет одновременное умножение всего числа на 10, чтобы выцепить очередную цифру мантиссы в десятичном виде.

[ivagor]

Прочитал, что принятая в русской литературе транскрипция фамилии автора - Мэчин. И что это не самый быстрый способ. Но, имхо, он вполне заслуживает реализации, повторюсь, я бы с интересом посмотрел. Только не надо никого заставлять

[ivagor]

Сишная программа для вычисления с использованием арктангенсов, в т.ч. Мэчина.

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

примерно 1 сек - расчет таблицы

Это было просто неприлично. Ускорил расчет таблицы в четыре с половиной раза, теперь рк сосчитает её примерно за четверть секунды. И сократил на 28 байт. Т.е. даже для расчета 100 знаков не обязательно заранее отдельно генерировать таблицу.

[b2m]

Что-то мне кажется, что метод:
На базе которого в этой ветке пишут тесты, не является методом "цифра за цифрой". Там именно идет добавление очередного члена ряда к длинному числу, которое находится в массиве, при этом еще идет одновременное умножение всего числа на 10, чтобы выцепить очередную цифру мантиссы в десятичном виде.

Нету там никакого добавления очередного члена ряда. Там просто преобразование числа Pi из одной системы счисления в другую (в данном случае десятичную), по одной цифре (умножая на 10 и выделяя цифру). Просто не все знают, что есть система счисления в которой число Pi равно 2.222222222222...