This package allows the user to perform a (Linear) General Least Squares Inversion with EPIC Tikhonov regularization. EPIC stands for Equal Posterior Information Condition.
EPIC_LS includes 2 subpackages:
- EPIC: codes to compute the EPIC for the general linear least squares problem.
- LeastSquares: codes to solve the linear least squares problem and the general linear least squares problem (with linear regularization).
- The article describing the details of the methodology and notation of the codes: Ortega-Culaciati, F., Simons, M., Ruiz, J., Rivera, L., & Diaz-Salazar, N. (2021). An EPIC Tikhonov regularization: Application to quasi-static fault slip inversion. Journal of Geophysical Research: Solid Earth, 126, e2020JB021141. https://doi.org/10.1029/2020JB021141 (Also see CITATION.bib).
- Examples on how to use the EPIC_LS package will be added to this repository during August 2021.
-
The calculation time will decrease considerably if scipy is compiled against mkl openblas (or similar) cappable of using multicore (default in anaconda python scipy). Remember to check environment variables OMP_NUM_THREADS, MKL_NUM_THREADS, MKL_DOMAIN_NUM_THREADS to have the right number of cores of your machine.
-
To use this package add this folder to your PYTHONPATH environment variable.