Здравствуйте! Есть такая проблема, что при вычислении определенного интеграла на отрезке [0;1] (), при разном количестве прямоугольников, результат не всегда является верным. 
Например при n = 10000, выдает 3,14159 и т.д
а при  n = 1000000, выдает 0,20384.

Я не пойму почему так происходит.
Я использую GTS 250 1Gb with CC 1.1 и Ubuntu 14.04 x64

Вот код
main.cpp


kernel.h


kernel.cu

Скорее всего, при сложении кучи маленьких чисел, накапливается ошибка округления.
Будьте аккуратнее!
(a+b) + c != a + (b + c)

заменить на 



А то у вас только максимум первые 65536 значений вычисляются (blockDim * gridDim). Потому и, кстати, правильно всё получается если меньшэ 65536 шагов.

Ну, и по поводу точности тожэ правильное замечание -- если складывать от бОльшых значений к мЕньшым -- то за 1000000 итэрацый можэт, при некотом невезении, до 1/17 от результирующего значения ошыбка накопиться. В общем, чтобы не париться, складывайте хотя бы в тип double. К тому жэ, сложэние тожэ можно распараллелить -- ну, в смысле, в рамках одного потока провести на GPU. Чтобы в итоге подсложыть только до 64k чисел. Но делайте через double всё равно.

У меня аппаратно не поддерживается данный тип. Карта старая gts250 на g92. На сколько я знаю, если даже double и написать, то все равно к float приводить будет.



Ок, попробую

Так ты ведь не на карте суммируешь.

tzirechnoy, Да, все верно, просто бегло прочитал и не сразу понял о чем речь.

Попробовал использовать Ваши изменения и вот что вышло

Просто поменял код ядра
n = 10000 --> result: 3.1415908
n = 1000000 --> result: 3.9995644

Поменял код ядра и тип данных на double
n = 10000 --> result:  0.4370185
n = 1000000 --> result: 1.1864045



Это сообщение отредактировал(а) ioManip - 13.11.2015, 13:02