code: def _compute_intensity(ccube, bexpcube):
          
          bexp_data = np.sqrt(bexpcube.data[0:-1, 0:] * bexpcube.data[1:, 0:])
          intensity_data = ccube.data / bexp_data
          intensity_map = HpxMap(intensity_data, ccube.hpx)
          return intensity_map
docs: Compute the intensity map