Який теоретичний коефіцієнт конвергенції для вирішувача FFT Poison?
Я розв'язую рівняння Пуассона: з n ( x , y , z ) = 3
Ось програма з використанням NumPy, яка робить розрахунок.
from numpy import empty, pi, meshgrid, linspace, sum
from numpy.fft import fftn, fftfreq
E_exact = 128/(35*pi)
print "Hartree Energy (exact): %.15f" % E_exact
f = open("conv.txt", "w")
for N in range(3, 384, 10):
print "N =", N
L = 2.
x1d = linspace(0, L, N)
x, y, z = meshgrid(x1d, x1d, x1d)
nr = 3 * ((x-1)**2 + (y-1)**2 + (z-1)**2 - 1) / pi
ng = fftn(nr) / N**3
G1d = N * fftfreq(N) * 2*pi/L
kx, ky, kz = meshgrid(G1d, G1d, G1d)
G2 = kx**2+ky**2+kz**2
G2[0, 0, 0] = 1 # omit the G=0 term
tmp = 2*pi*abs(ng)**2 / G2
tmp[0, 0, 0] = 0 # omit the G=0 term
E = sum(tmp) * L**3
print "Hartree Energy (calculated): %.15f" % E
f.write("%d %.15f\n" % (N, E))
f.close()
А ось графік конвергенції (просто побудуйте графік conv.txt
з вищевказаного сценарію, ось ноутбук, який робить це, якщо ви хочете грати з цим самостійно):
Як бачите, конвергенція лінійна, що було для мене несподіванкою, я думав, що FFT конвергується набагато швидше, ніж це.
Оновлення :
Рішення має зусилля на кордоні (я цього раніше не усвідомлював). Для того щоб FFT швидко конвергувався, рішення повинно мати всі похідні гладкими. Тому я також спробував наступну праву частину:
nr = 3*pi*sin(pi*x)*sin(pi*y)*sin(pi*z)/4
Хтось знає якийсь тест у 3D, щоб я міг бачити швидше конвергенцію, ніж лінійний?