DirectSmoothingOperator shows unwanted behaviour on edges
By fixing the normalization of the direct/brute-force smoother, the edges now show some instabilities:
from nifty import * x = RGSpace(64, harmonic=True) p = PowerSpace(x) sm = SmoothingOperator(p, sigma = 1.) sm.domain f = Field(p, val=1.) f.val[16] = 2 plotter = plotting.PowerPlotter() plotter.title = 'test' plotter.plot(sm(f))
-> The contributions to the edge pixels don't add up anymore. This is somehow the 'adjoint problem' of what has been solved by Martin's fix.
Is there any chance to fix this?