def _get_stddevs(self, C, stddev_types, num_sites, mag, c1_rrup,
                     log_phi_ss, mean_phi_ss):
        """
        Return standard deviations
        """
        phi_ss = _compute_phi_ss(C, mag, c1_rrup, log_phi_ss, mean_phi_ss)

        stddevs = []
        for stddev_type in stddev_types:
            assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES

            if stddev_type == const.StdDev.TOTAL:
                stddevs.append(np.sqrt(
                               C['tau'] * C['tau'] +
                               phi_ss * phi_ss) +
                               np.zeros(num_sites))

            elif stddev_type == const.StdDev.INTRA_EVENT:
                stddevs.append(phi_ss + np.zeros(num_sites))

            elif stddev_type == const.StdDev.INTER_EVENT:
                stddevs.append(C['tau'] + np.zeros(num_sites))
        return stddevs
def _get_stddevs ( self , C , stddev_types , num_sites , mag , c1_rrup , log_phi_ss , mean_phi_ss ) : phi_ss = _compute_phi_ss ( C , mag , c1_rrup , log_phi_ss , mean_phi_ss ) stddevs = [ ] for stddev_type in stddev_types : assert stddev_type in self . DEFINED_FOR_STANDARD_DEVIATION_TYPES if stddev_type == const . StdDev . TOTAL : stddevs . append ( np . sqrt ( C [ 'tau' ] * C [ 'tau' ] + phi_ss * phi_ss ) + np . zeros ( num_sites ) ) elif stddev_type == const . StdDev . INTRA_EVENT : stddevs . append ( phi_ss + np . zeros ( num_sites ) ) elif stddev_type == const . StdDev . INTER_EVENT : stddevs . append ( C [ 'tau' ] + np . zeros ( num_sites ) ) return stddevs
Return standard deviations