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

Interpolation #23

Closed
rowanc1 opened this issue Jan 17, 2017 · 1 comment
Closed

Interpolation #23

rowanc1 opened this issue Jan 17, 2017 · 1 comment
Assignees

Comments

@rowanc1
Copy link
Member

rowanc1 commented Jan 17, 2017

From @rowanc1 on November 17, 2015 22:0

In the interp utils we should be using np.searchsorted rather than np.argmin over the vectors:

Speed up of O(log(n)) vs O(n) for each point you are interpolating, could be significant, especially on large meshes.

Copied from original issue: simpeg/simpeg#170

@jcapriot
Copy link
Member

jcapriot commented Feb 6, 2017

I think we could actually get away without having the interputil_cython extension if you used np.searchsorted as that function accepts a list of points to insert. Unless you wanted to further optimize the matrix creation routines themselves. As they stand now, there isn't any speedup from cython except very little on the _interp_point_1D function.

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

No branches or pull requests

2 participants