search

The-ReSolver / Fields.jl

0 stars 0 forks source link

Modification of spectral fields in IFFT #9

Open tb6g16 opened 2 years ago

tb6g16 commented 2 years ago

Turns out that one (of possibly more) errors with the local residual was the modification of the spectral field during the IFFT.

To remedy I have added a couple methods with an extra auxiliary spectral field that can be modified during the transform without problem. This obviously isn't particularly desirable since it necessitates an extra spectral array. The ultimate solution will probably come from digging through the docs in extreme detail (I've had a cursory look already and found nothing).

https://github.com/The-ReSolver/Fields.jl/blob/0741e73f1eb6772909d98ef0aa20d6056658874d/src/fft.jl#L92-L102

gasagna commented 2 years ago

As referenced in #11 (which can be closed), it's fine to have this method. Only use it where necessary.