Introduced a limit on the number of ruptures per source fragment

openquake/calculators/preclassical.py

@@ -35,6 +35,7 @@
 from openquake.commonlib import readinput
 from openquake.calculators import base
 
+MAX_NUM_RUPTURES = 27_000  # so that the drouet calculation runs
 U16 = numpy.uint16
 U32 = numpy.uint32
 F32 = numpy.float32
@@ -98,7 +99,8 @@ def _filter_mag(srcs, min_mag):
     mmag = getdefault(min_mag, srcs[0].tectonic_region_type)
     out = [src for src in srcs if src.get_mags()[-1] >= mmag]
     for ss in out:
-        ss.num_ruptures
+        if ss.num_ruptures > MAX_NUM_RUPTURES:
+            raise RuntimeError('%s has too many ruptures' % ss)
     return out
 
 

openquake/hazardlib/contexts.py

@@ -1306,7 +1306,7 @@ def estimate_weight(self, src, srcfilter):
             # may happen for CollapsedPointSources
             return EPS
         src.nsites = len(sites)
-        step = 1 if src.code in b'pP' else 4
+        step = 1 if src.code in b'pP' else max(src.num_ruptures // 100, 4)
         C = sum(len(ctx) for ctx in self.get_ctxs(src, sites, step=step))
         src.dt = time.time() - t0
         if not C:

openquake/hazardlib/source/base.py

@@ -375,8 +375,8 @@ def __repr__(self):
         String representation of a source, displaying the source class name
         and the source id.
         """
-        return '<%s %s, weight=%.1f>' % (
-            self.__class__.__name__, self.source_id, self.weight)
+        return '<%s %s, rups=%d>' % (
+            self.__class__.__name__, self.source_id, self.num_ruptures)
 
 
 class ParametricSeismicSource(BaseSeismicSource, metaclass=abc.ABCMeta):