Дискретизация оператора Лапласа и быстрое решение уравнения Пуассона для внешности тела вращения

При расчете обтекания тела вращения под углом атаки пучком заряженных частиц (плазмы), движущихся в самосогласованном электрическом поле, требуется на каждом шаге по времени пересчитывать потенциал электрического поля, т. е. решать уравнение Пуассона. Для того чтобы эта задача решалась за приемлемое время, необходим быстрый алгоритм решения уравнения Пуассона. Построен алгоритм без насыщения, который требует $O(m^2n^2L\log_3L)$ операций, где $m$, $n$ и $L$ – число узлов сетки по переменным $r$, $\zeta$ и $\varphi$ некоторой криволинейной системы координат.