code: def _sumw(y, j_l, j_u, logWmax, p, phi):
      j = np.arange(j_l, j_u + 1)
      sumw = np.sum(np.exp(_logWj(y, j, p, phi) - logWmax))
      return sumw
  
file_path: statsmodels/statsmodels/genmod/_tweedie_compound_poisson.py
imports: import numpy
  from scipy._lib._util import _lazywhere
  from scipy.special import gammaln
label_desc: Use image node_id 5 for calling a global function with example usage: _sumw(y, j_l, j_u, logWmax, p, phi) and returns: sumw