uafgeotools / lts_array

Least trimmed squares array processing for infrasound and seismic data
MIT License
12 stars 6 forks source link

Fix rthEllipse handling in OLS code #30

Closed jwbishop closed 1 year ago

jwbishop commented 1 year ago

Adds additional error handling for the OLS code (alpha = 1.0) specifically for the uncertainty calculation. For some array geometries and slowness vector components, the rthEllipse function may have issues due to eExtrm[0:2] = [distE.min(), distE.max()]. This update catches the resulting ValueError and leaves the uncertainty value as a NaN.