opensha / opensha-dev

Development sandbox for OpenSHA project
BSD 3-Clause "New" or "Revised" License
3 stars 5 forks source link

RJ model revision #6

Closed mbarall closed 6 years ago

mbarall commented 6 years ago

This is a revision of the Reasenberg-Jones aftershock models. Highlights:

Add extensive documentation to the R&J model classes, so it is clear what they do.

Fix bugs, precision issues, and performance issues. Improve the design of this class hierarchy; consolidate and clean up code.

Replace the buggy adaptive quadrature routine with a correct one. Auto-switch to analytic integration when possible.

Add a routine to generate simulated aftershock sequences, according to R&J statistics with Page et al. magnitude of completeness. Use it to test the sequence-specific model.

Clip the tail of the epistemic probability distribution so that time is not spent doing calculations that contribute negligably to the final result.

Fix issue which caused the automatic system to attempt forming Bayesian combinations of incompatible models.