def METHOD_NAME():
    import time
    t1 = time.time()
    imgfile = 'fits_files/blankimg.fits'
    orig_nx = 147
    orig_ny = 124
    orig_seed = 151241
    # first, make the base image
    orig_img = galsim.ImageF(orig_nx, orig_ny, scale=1.)
    gal = galsim.Gaussian(sigma=2.5, flux=100.)
    gal.drawImage(orig_img, method='no_pixel')
    for iter in range(n_iter):
        # make it into an InterpolatedImage padded with imgfile
        int_im = galsim.InterpolatedImage(orig_img)
        # do it again with a particular seed
        int_im = galsim.InterpolatedImage(orig_img, rng = galsim.GaussianDeviate(orig_seed))
        # repeat
        int_im = galsim.InterpolatedImage(orig_img, rng = galsim.GaussianDeviate(orig_seed))
    t2 = time.time()
    print('time for %s = %.2f'%(funcname(),t2-t1))