print("RNG State:") print(random.getstate()) V = lambda x, y: np.array([x, y], dtype=np.double) renderers.setDefaultLensPlaneRenderer( "threads") # threads, mpi, opencl, None or a Renderer object renderers.setDefaultMassRenderer( "threads") # threads, mpi, None, or a Renderer object inversion.setDefaultInverter("threads") # threads, mpi or an Inverter object plotutil.setDefaultAngularUnit(ANGLE_ARCSEC) z_lens = 0.45 iws = inversion.InversionWorkSpace(z_lens, 150 * ANGLE_ARCSEC, cosmology=Cosmology(0.7, 1.0, 0, 0)) iws.addImageDataToList(images.ImagesData.load("images1pointgroups.imgdata"), 2.5, "extendedimages") iws.addImageDataToList(images.ImagesData.load("null1.imgdata"), 2.5, "extendednullgrid") iws.addImageDataToList(images.ImagesData.load("images2.imgdata"), 1.5, "extendedimages") iws.addImageDataToList(images.ImagesData.load("null2.imgdata"), 1.5, "extendednullgrid") # For a quick test to see if all code executes # iws.setDefaultInversionArguments(maximumGenerations = 10) iws.setUniformGrid(15) lens1, fitness, fitdesc = iws.invert(512)
print("RNG State:") print(random.getstate()) V = lambda x, y: np.array([x, y], dtype=np.double) renderers.setDefaultLensPlaneRenderer( "threads") # threads, mpi, opencl, None or a Renderer object renderers.setDefaultMassRenderer( "threads") # threads, mpi, None, or a Renderer object inversion.setDefaultInverter("threads") # threads, mpi or an Inverter object plotutil.setDefaultAngularUnit(ANGLE_ARCSEC) z_lens = 0.7 iws = inversion.InversionWorkSpace(z_lens, 60 * ANGLE_ARCSEC, cosmology=Cosmology(0.7, 0.3, 0, 0.7)) iws.addImageDataToList(images.ImagesData.load("images_00.imgdata"), 2.74714, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_01.imgdata"), 3.07785, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_02.imgdata"), 2.85271, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_03_td.imgdata"), 3.32098, "extendedimages") #iws.addImageDataToList(images.ImagesData.load("images_03_td.imgdata"), 3.32098, "extendedimages", { "timedelay": False }) iws.addImageDataToList(images.ImagesData.load("images_04.imgdata"), 3.1543, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_12.imgdata"), 3.00706,
-43.55867768 -49.49272263 19.30124747 -16.9427861 -0.4142903709 -6.376420558 -15.50428292 -1.884038019 -60.06804983 14.08154015 37.96515147 68.83852853 -64.12416086 -82.79545034 3.2336 36.2571622 13.72493334 39.32244461 15.59162438 """ z_lens = 0.4 for useWeights in [ False, True ]: iws = inversion.InversionWorkSpace(z_lens, 250*ANGLE_ARCSEC, cosmology=Cosmology(0.7, 0.27, 0, 0.73)) imgList = images.readInputImagesFile(pointData, True) for i in imgList: iws.addImageDataToList(i["imgdata"], i["z"], "pointimages") iws.setUniformGrid(15) if showPlots: plotutil.plotSubdivisionGrid(iws.getGrid()) plt.show() lens1, fitness1, fitdesc = iws.invert(popSize, maximumGenerations=maxGen, rescaleBasisFunctions=useWeights) iws.setSubdivisionGrid(lens1, 300, 400) if showPlots: plotutil.plotSubdivisionGrid(iws.getGrid()) plt.show()
-0.4142903709 -6.376420558 -15.50428292 -1.884038019 -60.06804983 14.08154015 37.96515147 68.83852853 -64.12416086 -82.79545034 3.2336 36.2571622 13.72493334 39.32244461 15.59162438 """ z_lens = 0.4 for useWeights in [False, True]: iws = inversion.InversionWorkSpace(z_lens, 250 * ANGLE_ARCSEC, cosmology=Cosmology(0.7, 0.27, 0, 0.73)) imgList = images.readInputImagesFile(pointData, True) for i in imgList: iws.addImageDataToList(i["imgdata"], i["z"], "pointimages") iws.setUniformGrid(15) if showPlots: plotutil.plotSubdivisionGrid(iws.getGrid()) plt.show() lens1, fitness1, fitdesc = iws.invert(popSize, maximumGenerations=maxGen, rescaleBasisFunctions=useWeights, basisFunctionType=basisFunctionType) iws.setSubdivisionGrid(lens1, 300, 400)
# Write the RNG state, in case we want to reproduce the run exactly # (note that the GRALE_DEBUG_SEED environment variable will need # to be restored as well) import random print("RNG State:") print(random.getstate()) V = lambda x, y: np.array([x,y], dtype=np.double) renderers.setDefaultLensPlaneRenderer("threads") # threads, mpi, opencl, None or a Renderer object renderers.setDefaultMassRenderer("threads") # threads, mpi, None, or a Renderer object inversion.setDefaultInverter("threads") # threads, mpi or an Inverter object plotutil.setDefaultAngularUnit(ANGLE_ARCSEC) z_lens = 0.7 iws = inversion.InversionWorkSpace(z_lens, 60*ANGLE_ARCSEC, cosmology=Cosmology(0.7, 0.3, 0, 0.7)) iws.addImageDataToList(images.ImagesData.load("images_00.imgdata"), 2.74714, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_01.imgdata"), 3.07785, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_02.imgdata"), 2.85271, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_03.imgdata"), 3.32098, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_04.imgdata"), 3.1543 , "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_12.imgdata"), 3.00706, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_13.imgdata"), 3.28815, "extendedimages") iws.addImageDataToList(images.ImagesData.load("images_73.imgdata"), 3.3154 , "extendedimages") iws.setDefaultInversionArguments(sheetSearch = "genome") # For a quick test to see if all code executes # iws.setDefaultInversionArguments(maximumGenerations = 10, sheetSearch = "genome")