brendanjmeade / celeri

Next generation earthquake cycle kinematics
BSD 3-Clause "New" or "Revised" License
24 stars 6 forks source link

Alternative packaging for TDE calls #72

Closed brendanjmeade closed 2 years ago

brendanjmeade commented 2 years ago

We'd love to do all TDE to OBS calculations at once, and Ben's toolbox is set up for this. This would likely be much faster than the single TDE per call version that we're doing now.

However, because we have a different nonlinear map projection for each TDE this is not possible. We could consider a different way of packaging this information so that we could do this all in one call, albeit one with a lot of redundancy.

all_projected_station_coordinates = list()

for each TDE: do optimal local projection Get local TDE coordinates Get local station coordinates all_projected_station_coordinates(local station coordinates)

At the end of this, all_projected_station_coordinates will have a size equal to len(TDEs) * len(stations). It's huge.

One giant call to Ben's TDE solver. What is the correct function name here?

This will have many more station coordinates than we need but we might be able to zero out a huge number of combinations or select only the ones we need after the fact. Not sure.