docs: def _get_faulting_style_term(self, C, rup):
          """
          Compute and return faulting style term, that is the sum of the second
          and third terms in equation 1, page 74.
          """
          # this implements equations 5 and 6 at page 1032. f7 is the
          # coefficient for reverse mechanisms while f8 is the correction
          # factor for normal ruptures
          if rup.mag > 5.0:
              f7 = C['a11']
              f8 = C['a12']
          elif rup.mag >= 4:
              f7 = C['a11'] * (rup.mag - 4.)
              f8 = C['a12'] * (rup.mag - 4.)
          else:
              f7 = 0.0
              f8 = 0.0
          # ranges of rake values for each faulting mechanism are specified in
          # table 2, page 1031
          return (f7 * float(rup.rake > 30 and rup.rake < 150) +
                  f8 * float(rup.rake > -150 and rup.rake < -30))