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")