Skip to content

Stable inversion

Maximilian Schanner requested to merge stable_inversion into master

When inverting records with locations and times that are close, the kernel matrix is ill-conditioned. To address this the np.linalg.inv call is replaced by an eigen-decomposition and (numerically) negative eigenvalues are regularized.

Merge request reports