OpenQuake Engine 3.17.0

[Michele Simionato (@micheles)]

• Raise an early error if the user disable both realizations and statistics
• Internal: added dbserver.user in openquake.cfg
• Solved the hanging of classical calculations due to large zmq packets
  and reduced memory occupation

[Marco Pagani (@mmpagani), Michele Simionato (@micheles)]

• Added support for the "direct" method in MRD calculations

[Michele Simionato (@micheles)]

• Reimplemented the conditional_spectrum calculator as a post-processor

[Paolo Tormene (@ptormene), Michele Simionato (@micheles)]

• Raise an error when the same parameter is set in different sections
  of the job.ini file

[Michele Simionato (@micheles)]

• Fixed tau-phi inversion in LanzanoLuzi2019
• Fixed another bug with conditioned GMFs appearing as the error
  array[m, n, bad] = exp(mean_covs[0, g, m, n], im) TypeError: list indices must be integers or slices, not tuple

[Claudio Schill (@claudio525)]

• Fixed sorting bug in the sigma_mu adjustment factor in
  the Kuehn et al. (2020) GMM

[Michele Simionato (@micheles)]

• Fixed a subtle filtering bug in conditioned GMFs
• Fixed shakemap calculations to take the risk IMTs and not the hazard IMTs
• Extended shakemap2gmfs to amplify the GMFs based on user-provided vs30
  values
• Removed sorting of the sections in the multiFaultSource reader to keep
  the rup_idxs consistent with the order in the files
• GMFs too small to produce losses now give a warning rather than an error
• Fixed bug in get_ry0_distance breaking conditioned GMFs
• Made sure lon and lat are rounded to 5 digits in the site collection

[Graeme Wheaterill]

• Fixes to bchydro_2016_epistemic, eshm20_craton and lanzano_2019

[Michele Simionato (@micheles)]

• Extended sensitity analysis to file parameters
• Introduced new loss type residents

[Marco Pagani (@mmpagani)]

• Fixed bug in get_profiles_from_simple_fault_data affecting KiteSurfaces

[Michele Simionato (@micheles)]

• Fixed the rupture exporter in presence of GMPETables
• Optimized the postclassical phase
• Fixed a serious issue in the classical calculator with tiles
• Vectorized the association assets <-> hazard sites (up to 32x speedup)
• Fixed bug in disaggregation calculations associated to an error
  in the line [mag] = np.unique(np.round(ctx.mag, 6))
• Optimized the calculation of geohashes by using numba

[Paolo Tormene (@ptormene)]

• Extended total_losses to include business_interruption

[Michele Simionato (@micheles)]

• Recognized occupants_avg field in the exposure
• Optimized reading the exposure (both speed and memory)
• Optimized the rupture sampling for MultiFaultSources and improved
  the parallelization for MultiFaultSources and MultiPointSources
• Made the parameter truncation_level mandatory
• Fixed the damage state header in the aggrisk outputs
• Changed the order of the IMTs to be by period and not lexicographic
• Fixed the realizations extractor for scenario calculations reading
  ruptures with multiple different tectonic region types

[Kendra Johnson (@kejohnso)]

• Fixed a bug in disaggregation by Lon-Lat-TRT when there are multifault
  ruptures and multiple sites

[Michele Simionato (@micheles)]

• Deprecated the conversion in the exposure
• Added new loss types "area" and "number" for usage in risk calculations
• Made the contexts immutable and fixed mutability bugs potentially
  affecting the modules akkar_bommer_2010, bindi_2011_ipe, can_shm6_inslab,
  lanzano_2019, youd_etal_2002, zhang_zhao_2005, cauzzi_faccioli_2008_swiss,
  chiou_youngs_2008_swiss campbell_bozorgnia_2014
• Changed the disaggregation algorithm (see manual)
• Added a view risk_by_rup
• Internal: added command oq compare med_gmv <imt>
• Extended extract/ruptures to extract a single rupture given the rup_id
• Fixed a bug in classical_damage causing a ValueError: could not broadcast
  input array from shape (X,Y) into shape (X,)
• Added check for missing site parameters

[Paolo Tormene (@ptormene)]

• Internal: reduced the size of the gsim column in the view
  required_params_per_trt in the .rst report

[Michele Simionato (@micheles)]

• Optimized the generation of events
• Fixed the GMFs CSV importer so that it can import exported GMFs even in
  the presence of a custom_site_id
• Implemented sampling for mutex sources with grp_probability
• Introduced extreme_gmv parameter in the job.ini
• Introduced a new output "Annual Frequency of Events" and refined the command
  oq plot event_based_mfd?
• Moved the minimum magnitude filtering after the rupture sampling
• Changed the internals of rupture sampling and fixed the generated MFD
  in the case of a same source appearing in multiple logic tree branches
• Fixed a bug in scenario calculations causing a mysterious IndexError in the
  line rctx.clon = c.lons[rctx.sids]
• Changed the disaggregation exporter to add a column with the IML
  associated to the PoE for the mean disaggregation outputs
• Internal: added a parameter override_vs30 in the job.ini

[Marco Pagani (@mmpagani), Michele Simionato (@micheles)]

• Fixed a bug in the source writer when writing NonParametricSources

[Michele Simionato (@micheles)]

• Changed the conditional spectrum storage and fixed the exporter to export
  all PoEs

[Marco Pagani (@mmpagani)]

• Fixed conditional spectrum with investigation time different from 1 year

[Michele Simionato (@micheles)]

• Fixed documentation of the Extractor for hazard curves and maps
• Fixed bug in oq zip: when starting from a directory different from the
  job.ini directory, .hdf5 files for nonparametric ruptures were not included
• Honored minimum_magnitude for MultiFaultSources
• Overhauled the experimental postprocessing framework and fixed many
  bugs in the compute_mrd postprocessor

[Anirudh Rao (@raoanirudh)]

• Added a verification test suite for the Conditioned GMFs calculator

[Michele Simionato (@micheles)]

• Added output rates_by_src
• Fixed method '/v1/ini_defaults' returning NaNs for the site parameters
• Added output mean_disagg_bysrc for single site calculations with use_rates
• Added command oq reduce_smlt
• Fixed reinsurance calculations when there is a single loss type
• Changed the source seed generation to not depend on the source ID suffix
• Tracebacks were printed 3 times in the console: this is finally fixed
• Forbidden duplicate sources in atomic groups
• Forbidden duplicate mutually exclusive sources
• Fixed conditional spectrum calculator with multiple sites and implemented
  parallelization
• Fixed the source ID generation for CollapsedPointSources
• Fixed the error message for wrong weights in the gsim logic tree file
• Fixed KeyError in calculations using the EasternCan15Mid with PGA
• Replaced hard-coded limit to 25 levels in the MRD postprocessor
  with a simple warning

[Kendra Johnson (@kejohnso)]

• Added a job parameter 'reqv_ignore_sources', which indicates which
  sources - specified as a list of thier source ids - should not be collapsed
  using the equivalent distance

[Michele Simionato (@micheles)]

• Extended (re)insurance calculations to manage insurance deductibles at asset
  level
• Internal: added a parameter use_rates to compute statistical hazard curves
  in terms of rates, not probabilities
• Internal: changed disagg_by_src to work in terms of rates
• Added command oq runsite <lon>,<lat> for mosaic PSHA calculations
• Enhanced --reuse-input to reuse preclassical calculations if possible
• Normalized IMTs in the minimum_intensity dictionary (i.e. SA(1.)=>SA(1.0))
• Fixed event based calculations running out of memory due to hazard curve
  arrays being instantiated without need

[Julián Santiago Montejo Espitia (@SantiMon1417)]

• Contributed the GMPE Arteta et al. (2023) for crustal events on northern
  South America

[Michele Simionato (@micheles)]

• Setting both sites and site_model_file now works also for models
  with nonstandard site parameters, like EUR

[Paolo Tormene (@ptormene)]

• Added AELO mode for the engine webui, providing a web form to insert input
  values and launch a calculation. Input values are validated and the
  progress of the calculation can be monitored through the webui. On
  completion, the server sends an email to the job owner, providing feedback
  about failure or completion, with a link to retrieve a traceback or to
  browse the calculation outputs. Corresponding tests were added and the
  existing tests for the webui were improved, cleaning up calculations after
  they have been checked.

[Michele Simionato (@micheles)]

• Changed the default: now all realizations are consider in disaggregation
  calculations, not only the one closest to the mean hazard curve
• Added feature infer_occur_rates for multiFault sources
• Sampling the ruptures was raising a BlockingIOError sometimes

[Kendra Johnson (@kejohnso)]

• Made it possible to use extendModel for multiple branchSets in the
  source model logic tree

[Michele Simionato (@micheles)]

• Added an exporter for disagg_by_src
• Implemented mean disaggregation outputs, with the means computed in terms
  of rates, not probabilities

[Kendra Johnson (@kejohnso)]

• Corrected a bug when using the reqv feature: all sources were
  collapsed instead of the ones with the specified tectonic region types

[Michele Simionato (@micheles)]

• Changed the lon-lat binning algorithm to produce 4 times less bins
• Fixed install.py: it must not set the shared_dir
• Optimized disaggregation in performance, memory consumption and data
  transfer; also changed the task distribution algorithm