Вычисление Pi с использованием массивов Float или Integer

Question

Вычисление Pi с использованием массивов Float или Integer

Я использую Aparapi, чтобы сделать число хруст в моей программе Java на графическом процессоре. Из того, что я понимаю, Апарапи отлично играет с плавающими массивами.

Я хочу вычислить Pi до N-го десятичного числа, используя Aparapi. Я думаю об использовании метода Лейбница, но я не уверен, как бы я работал с представлением и хранением длинных десятичных знаков в форме float или integer.

Будет ли массив целых чисел работать, причем размер массива равен N числу десятичных знаков?

int[] digits = new int[N];

Если бы я использовал это с методом Лейбница, мне нужно было бы вычислить массив из N целых чисел для M терминов, которые я нахожу (Liebniz говорит, что pi/4 = 1 - 1/3 + 1/5 - 1/7 + 1/9 ....), а затем добавьте их вместе и умножьте полученное число на 4. Но это означало бы, что мне нужно было бы выделить M много целых чисел для КАЖДОГО термина, который я вычислил, что бы скомпоновать и действительно отразиться на памяти.

tl; dr: Как я могу вычислить Pi, используя циклы операций с float, чтобы я мог сделать это с Aparapi?

Огромное спасибо!

1

добавлено 28 Декабрь 2012 в 09:04 автор K. Barresi

Просмотры: 2

источник

nl ja de

Просто создайте случайные точки на карте апарапи и в конце нарисуйте круг и просто разделите общее количество точек на количество точек внутри круга, вы найдете 4/pi

добавлено 24 Январь 2013 в 01:23, автор huseyin tugrul buyukisik, источник

Я пропустил метод Лейбница, извините за то, что нажал «Монтекарло».

добавлено 25 Май 2013 в 12:18, автор huseyin tugrul buyukisik, источник

Это медленный способ действительно для вычисления pi.

добавлено 28 Декабрь 2012 в 09:10, автор Blender, источник

Вы все равно будете хранить номера, независимо от вашего алгоритма. Вот гораздо более быстрый алгоритм: stackoverflow.com/questions/14065929/…

добавлено 28 Декабрь 2012 в 09:15, автор Blender, источник

Извините, вот что я пытался вставить: mathworld.wolfram.com/Brent-SalaminFormula.html </а>

добавлено 28 Декабрь 2012 в 09:21, автор Blender, источник

Вот почему я надеюсь, что у кого-то есть лучший способ сделать это, все еще используя меньшие типы данных

добавлено 28 Декабрь 2012 в 09:14, автор K. Barresi, источник

Я не думаю, что вы вставили правильный URL :)

добавлено 28 Декабрь 2012 в 09:20, автор K. Barresi, источник

Благодаря! И теперь проблема заключается в том, как хранить эти данные ...

добавлено 28 Декабрь 2012 в 09:22, автор K. Barresi, источник

1 ответы

автор huseyin tugrul buyukisik · Answer 1 · 2013-05-25T12:11Z

Вам нужно уменьшить, чтобы частично добавить m элементов в аннотированный массив @Local (используется для локальной памяти в aparapi). Допустим, что существует N общих элементов, поэтому вам нужны N/m частичные суммы, которым требуется более высокая пропускная способность (поэтому в локальной памяти).

Также вы должны использовать localBarrier (); для синхронизации ядер в вычислительных единицах или рабочих единицах в местных рабочих группах. Затем отправьте данные в основную память.

Очень короткий пример:

Вам нужно 1 + 2 + 3 + 4 + 5 + 6 + ..... и у вас есть 3 ядра на вычислительную единицу gpu.

1 + 2 выполняется в ядре № 1

3 + 4 выполняется в ядре номер 2

5 + 6 выполняется в ядре № 3

то вы добавляете выходы ядра в локальной памяти, которые разделяются всеми соседними ядрами этого вычислительного блока.

Наконец, вы получите 3,7,11 в массиве, который суммируется до 21 в локальной памяти.

Все эти данные всех вычислительных единиц в основной памяти, такие как 21,57, ... тогда вы можете добавить их все просто в CPU.

Конечно, будет (+/-) (1/(2 * n + 1)) вместо 1,2,3,4,5