def METHOD_NAME(self):
    "Test the bleed_channel function."
    channel = (np.ones(self.nypix, dtype=int) *
               np.random.poisson(lam=self.sky_level))
    star_center = 1000
    star_size = 20
    channel[star_center-star_size:star_center+star_size] = 2*self.full_well
    total_count = sum(channel)
    # Run the bleed channel function.
    bled_channel = imsim.bleed_channel(channel, self.full_well)
    # Check that charge is conserved in this case.
    self.assertEqual(total_count, sum(bled_channel))
    # Check that maximum pixel value is full well.
    self.assertEqual(self.full_well, max(bled_channel))
    # Check that an un-interrupted bleed trail is produced
    index = np.where(bled_channel == self.full_well)
    nmin = max(0, min(index[0])-1)
    nmax = max(index[0])
    expected_sum = (nmax*(nmax+1))//2 - (nmin*(nmin+1))//2
    self.assertEqual(sum(index[0]), expected_sum)
    # Check that bleed trail is centered on original star.
    self.assertLessEqual(np.abs(star_center - (nmin + nmax)/2), 1)