search

uafgeotools / rtm

A Python package for locating infrasound sources using reverse time migration
https://uaf-rtm.readthedocs.io/
MIT License
38 stars 13 forks source link

Grid search hangs if starttime and endtime clipping is requested for large stack functions #51

Closed liamtoney closed 2 years ago

liamtoney commented 3 years ago

For large stack functions S, the grid_search() function hangs on this code:

    if starttime:
        S = S.where((S.time >= np.datetime64(starttime)), drop=True)
    if endtime:
        S = S.where((S.time <= np.datetime64(endtime)), drop=True)

This appears to be a memory-related issue. It's likely we could devise a work-around for this.

liamtoney commented 2 years ago

I think this will be fixed by #65...