Source code for openquake.hazardlib.gsim.multi

#  -*- coding: utf-8 -*-
#  vim: tabstop=4 shiftwidth=4 softtabstop=4

#  Copyright (c) 2018, GEM Foundation

#  OpenQuake is free software: you can redistribute it and/or modify it
#  under the terms of the GNU Affero General Public License as published
#  by the Free Software Foundation, either version 3 of the License, or
#  (at your option) any later version.

#  OpenQuake is distributed in the hope that it will be useful,
#  but WITHOUT ANY WARRANTY; without even the implied warranty of
#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#  GNU Affero General Public License for more details.

#  You should have received a copy of the GNU Affero General Public License
#  along with OpenQuake.  If not, see <http://www.gnu.org/licenses/>.

"""
Module exports :class:`MultiGMPE`, which can create a composite of
multiple GMPEs for different IMTs when passed a dictionary of ground motion
models organised by IMT type or by a string describing the association
"""
import collections
from openquake.hazardlib import const
from openquake.hazardlib.gsim.base import GMPE
from openquake.hazardlib.imt import from_string

uppernames = '''
DEFINED_FOR_INTENSITY_MEASURE_TYPES
DEFINED_FOR_STANDARD_DEVIATION_TYPES
REQUIRES_SITES_PARAMETERS
REQUIRES_RUPTURE_PARAMETERS
REQUIRES_DISTANCES
'''.split()


[docs]class MultiGMPE(GMPE, collections.Mapping): """ The MultiGMPE can call ground motions for various IMTs when instantiated with a dictionary of ground motion models organised by IMT or a string describing the association. In the case of spectral accelerations the period of the IMT must be defined explicitly and only SA for that period will be computed. """ #: Supported tectonic region type is undefined DEFINED_FOR_TECTONIC_REGION_TYPE = "" #: Supported intensity measure types are not set DEFINED_FOR_INTENSITY_MEASURE_TYPES = set() #: Supported intensity measure component is horizontal #: :attr:`~openquake.hazardlib.const.IMC.HORIZONTAL`, DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = const.IMC.HORIZONTAL #: Supported standard deviation type DEFINED_FOR_STANDARD_DEVIATION_TYPES = set([const.StdDev.TOTAL]) #: Required site parameters will be set be selected GMPES REQUIRES_SITES_PARAMETERS = set() #: Required rupture parameter is magnitude, others will be set later REQUIRES_RUPTURE_PARAMETERS = set(('mag', )) #: Required distance metrics will be set by the GMPEs REQUIRES_DISTANCES = set() def __init__(self, gsim_by_imt): """ Instantiate with a dictionary of GMPEs organised by IMT """ super().__init__() self.gsim_by_imt = gsim_by_imt for name in uppernames: setattr(self, name, getattr(self, name).copy()) for imt, gsim in gsim_by_imt.items(): imt_class = from_string(imt).__class__ if imt_class not in gsim.DEFINED_FOR_INTENSITY_MEASURE_TYPES: raise ValueError("IMT %s not supported by %s" % (imt, gsim)) for name in uppernames: getattr(self, name).update(getattr(gsim, name)) def __iter__(self): yield from self.gsim_by_imt def __getitem__(self, imt): return self.gsim_by_imt[imt] def __len__(self): return len(self.gsim_by_imt) def __hash__(self): items = tuple((imt, str(gsim)) for imt, gsim in sorted(self.gsim_by_imt.items())) return hash(items)
[docs] def get_mean_and_stddevs(self, sctx, rctx, dctx, imt, stddev_types): """ Call the get mean and stddevs of the GMPE for the respective IMT """ return self.gsim_by_imt[str(imt)].get_mean_and_stddevs( sctx, rctx, dctx, imt, stddev_types)