Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improve the spectral periodization #28

Open
savmickael opened this issue Aug 8, 2018 · 1 comment
Open

Improve the spectral periodization #28

savmickael opened this issue Aug 8, 2018 · 1 comment

Comments

@savmickael
Copy link
Contributor

No description provided.

@dumasl
Copy link
Contributor

dumasl commented Aug 9, 2018

As the doc says the periodization is to be used with a low pass filter (user kernel interpolator). Otherwise (zero padding to enlarge the spectral domain), the user kernel interpolator is actually convolve with a sinc interpolator which is not what the user wants.

We are using periodization then because this solution provides us some continuity in spectral domain when increasing the number of frequencies by a factor N. The problem is if the user kernel interpolator is actually not a good low pass filter (say it has rebounds) then some high frequencies coming from the periodization might be kept by these rebounds.
I remember having read a paper on variational method to emulate what the spectrum should actually be if we enlarge its domain. I shall try to look for it but in the mean time there are some tests that could be done on the periodization to try to lesser the magnitude of the added frequencies so that if there are rebounds they will only kept small frequencies. Or we could decide to remove the rebounds but it comes with a sinc convolution in spatial domain so we want to only cut the far away rebound (high frequency) so that the sinc convolution is as close as possible to a dirac convolution.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants