Hello again, I noticed two problems with my reply.
• I was speaking about currents while you interpolate densities. (This does not matter much.) • I was assuming that your output field is 5000^2 while it’s rather (20*9)^2 because you keep the n parameter to interpolate_density at its default value of 9. Can you verify this? The second point modifies my analysis quite a lot in that the kernel size is only 9*9 and not 250*250. Now it’s something like 32000 CPU cycles per kernel evaluation, which seems horrible. Could you send me a minimal script that reproduces your problem, i.e. that interpolates some (whatever) density with the same parameters that you use? Then I can check things directly instead of guessing. Thanks Christoph
smime.p7s
Description: S/MIME cryptographic signature