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