code: def image_single_point_source(self, image_model_class, kwargs_lens, kwargs_source, kwargs_lens_light,
                                    kwargs_ps):
          
          # reconstructed model with given psf
          model, error_map, cov_param, param = image_model_class.image_linear_solve(kwargs_lens, kwargs_source,
                                                                                kwargs_lens_light, kwargs_ps)
          #model = image_model_class.image(kwargs_lens, kwargs_source, kwargs_lens_light, kwargs_ps)
          data = image_model_class.Data.data
          mask = image_model_class.ImageNumerics.mask
          point_source_list = image_model_class.point_sources_list(kwargs_ps, kwargs_lens)
          n = len(point_source_list)
          model_single_source_list = []
          for i in range(n):
              model_single_source = (data - model + point_source_list[i]) * mask
              model_single_source_list.append(model_single_source)
          return model_single_source_list
docs: return model without including the point source contributions as a list (for each point source individually)
          :param image_model_class: ImageModel class instance
          :param kwargs_lens: lens model kwargs list
          :param kwargs_source: source model kwargs list
          :param kwargs_lens_light: lens light model kwargs list
          :param kwargs_ps: point source model kwargs list
          :return: list of images with point source isolated