code: def _differential_to_integral(hpx_map):
          
          ebins = hpx_map.hpx.ebins
          ratio = ebins[1:] / ebins[0:-1]
          half_log_ratio = np.log(ratio) / 2.
          int_map = ((hpx_map.data[0:-1].T * ebins[0:-1]) +
                     (hpx_map.data[1:].T * ebins[1:])) * half_log_ratio
          return HpxMap(int_map.T, hpx_map.hpx)
docs: Convert a differential map to an integral map
  
          Here we are using log-log-quadrature to compute the integral quantities.