-
Notifications
You must be signed in to change notification settings - Fork 15
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
Add routine to Fourier transform collection of (x, y, flux) sampled randomly #195
base: master
Are you sure you want to change the base?
Conversation
Slow because it uses brute force to find which triangle each point is in.
Also further fixes like passing a pointer to the Delauney triangulation and adding some timing tests.
Also added further timing tests and make sure to clean up.
On my end, the code seems to work, but needs to be restructured to better fit the organization/layout of the galario code. The changes I think I need to make (though input/feedback welcome):
|
and one helper: - triangulate_h - bin_triangles_h - interpolate_or_bin_to_image_h - interpolate_on_triangle_h
Also cleaned up some commented out code and some small rearranging.
Ok, I spent some time today renaming most of the new functions to have an _h suffix to indicate that they are run on the host. I also split interpolate_to_image into four smaller functions ( Also added a few test scripts that I've been using to test that the code works properly. A few questions I still have/places I still need help that might require some discussion:
|
Per our discussion over e-mail, this pull request is meant to add functionality to galario so that it can take the Fourier transform of an "unstructured image", i.e. a list of (x, y, flux) where (x,y) do not need to fall on a regularly sampled grid. It does so by triangulating the collection of points and using the triangulation to interpolate onto a regular grid. Once that is done, the regular routines in galario can be used to do the Fourier transform and then to sample onto the relevant baselines.