def METHOD_NAME(skymodel, nxy_min=10, nxy_max=300, npts=8, pixel_scale=0.2):
    timing = OrderedDict()
    nphot = OrderedDict()
    for nxy in [int(x) for x in np.logspace(np.log10(nxy_min),
                                            np.log10(nxy_max), npts)]:
        image = galsim.Image(nxy, nxy, scale=pixel_scale)
        t0 = time.clock()
        skymodel.addNoiseAndBackground(image, photParams=skymodel.photParams)
        timing[nxy] = time.clock() - t0
        try:
            nphot[nxy] = skymodel.photon_array.size()
        except AttributeError:
            nphot[nxy] = 0
    return timing, nphot