def METHOD_NAME(self):
    subpixel_x, subpixel_y = self._adaptive_grid._high_res_coordinates
    x, y = self._adaptive_grid._x_low_res, self._adaptive_grid._x_low_res
    model = LightModel(light_model_list=["GAUSSIAN"])
    kwargs_light = [{"center_x": 0, "center_y": 0, "sigma": 1, "amp": 1}]
    subgrid_values = model.surface_brightness(subpixel_x, subpixel_y, kwargs_light)
    image1d = model.surface_brightness(x, y, kwargs_light)
    image_added = self._adaptive_grid._merge_low_high_res(image1d, subgrid_values)
    added_array = util.image2array(image_added)
    supersampled_values = self._adaptive_grid._average_subgrid(subgrid_values)
    assert (
        added_array[util.image2array(self._supersampling_indexes)]
        == supersampled_values
    )
    image_high_res = self._adaptive_grid._high_res_image(subgrid_values)
    assert len(image_high_res) == self.nx * self._supersampling_factor