def galaxy_fit_phase():
pixel_scales = 0.1
image_shape = (150, 150)
integration_util.reset_paths(test_name=test_name, output_path=output_path)
grid = al.grid.uniform(shape_2d=image_shape, pixel_scales=pixel_scales, sub_size=4)
galaxy = al.Galaxy(
redshift=0.5,
light=al.lp.SphericalExponential(
centre=(0.0, 0.0), intensity=1.0, effective_radius=0.5
),
)
image = galaxy.profile_image_from_grid(galaxies=[galaxy], grid=grid)
noise_map = al.array.manual_2d(
array=np.ones(image.shape), pixel_scales=pixel_scales
)
data = al.GalaxyData(image=image, noise_map=noise_map, pixel_scales=pixel_scales)
phase1 = al.PhaseGalaxy(
phase_name=test_name + "/",
galaxies=dict(
gal=al.GalaxyModel(redshift=0.5, light=al.lp.SphericalExponential)
),
use_image=True,
sub_size=4,
optimizer_class=af.MultiNest,
)
phase1.run(galaxy_data=[data])
def galaxy_fit_phase():
pixel_scales = 0.1
image_shape = (150, 150)
integration_util.reset_paths(test_name=test_name, output_path=output_path)
grid = al.Grid.uniform(shape_2d=image_shape, pixel_scales=pixel_scales, sub_size=4)
galaxy = al.Galaxy(
redshift=0.5,
mass=al.mp.SphericalIsothermal(centre=(0.0, 0.0), einstein_radius=1.0),
)
deflections = galaxy.deflections_from_grid(galaxies=[galaxy], grid=grid)
noise_map = al.Array.manual_2d
sub_array_1d=np.ones(deflections[:, 0].shape), pixel_scales=pixel_scales