docs: def _resample_residuals(self, stars, epsf):
          """
          Compute normalized residual images for all the input stars.
  
          Parameters
          ----------
          stars : `EPSFStars` object
              The stars used to build the ePSF.
  
          epsf : `EPSFModel` object
              The ePSF model.
  
          Returns
          -------
          star_imgs : 3D `~numpy.ndarray`
              A 3D cube containing the resampled residual images.
          """
  
          shape = (stars.n_good_stars, epsf.shape[0], epsf.shape[1])
          star_imgs = np.zeros(shape)
          for i, star in enumerate(stars.all_good_stars):
              star_imgs[i, :, :] = self._resample_residual(star, epsf)
  
          return star_imgs