Вопрос по двумерному преобразованию Фурье через fft2 на Python

Решаю двумерное стационарное уравнение дифракции, используя двумерное преобразование Фурье. Делаю на языке Python средствами fft2. Для получения решения сначала выполняется прямое преобразование Фурье только от граничного условия (а оно в виде exp(-(x^2+y^2))). Затем умножаю на некоторую функцию, зависящую от i и k, где i соответствует x, k - y (то есть уже на сетке). А затем делаю обратное преобразование Фурье. Вопрос в следующем. Выводится амплитуда поля (как корень из суммы квадратов действительной и мнимой частей) при конкретном z = const. В итоге все равно должна получиться поверхность в гауссовой форме, которая при увеличении z должна расплываться. Но вот это не получается. Появляются дополнительные "бугры", а исходный "бугор" смещается. Не понимаю, почему это происходит. Предполагаю, что что-то на стадии умножения на функцию (перед взятием обратного Фурье). Может быть фаза какая-то дополнительная. Буду благодарен за любые помощь и подсказки.