self.galaxies.lens.mass.centre_0 = af.GaussianPrior(mean=4.0, sigma=0.1)
self.galaxies.lens.mass.centre_1 = af.GaussianPrior(mean=4.0, sigma=0.1)
self.galaxies.source.light.centre_0 = af.GaussianPrior(mean=4.0, sigma=0.1)
self.galaxies.source.light.centre_1 = af.GaussianPrior(mean=4.0, sigma=0.1)
phase1 = LensPhase(
phase_name="phase_1",
phase_folders=phase_folders,
galaxies=dict(
lens=al.GalaxyModel(redshift=0.5, mass=al.mp.SphericalIsothermal),
source=al.GalaxyModel(redshift=1.0, light=al.lp.EllipticalSersic),
),
positions_threshold=0.5,
non_linear_class=non_linear_class,
)
phase1.optimizer.const_efficiency_mode = True
phase1.optimizer.n_live_points = 30
phase1.optimizer.sampling_efficiency = 0.8
return al.PipelineDataset(name, phase1)
if __name__ == "__main__":
import sys
runner.run(
sys.modules[__name__],
positions=[[(5.6, 4.0), (4.0, 5.6), (2.4, 4.0), (4.0, 2.4)]],
)
def make_pipeline(name, phase_folders, non_linear_class=af.MultiNest):
class MMPhase(al.PhaseImaging):
pass
phase1 = MMPhase(
phase_name="phase_1",
phase_folders=phase_folders,
galaxies=dict(lens=al.GalaxyModel(
redshift=0.5,
light_0=al.lp.EllipticalSersic,
light_1=al.lp.EllipticalSersic,
mass=al.mp.EllipticalIsothermal,
align_centres=True,
)),
non_linear_class=non_linear_class,
)
phase1.optimizer.const_efficiency_mode = True
phase1.optimizer.n_live_points = 20
phase1.optimizer.sampling_efficiency = 0.8
return al.PipelineDataset(name, phase1)
if __name__ == "__main__":
import sys
runner.run(sys.modules[__name__])
data_type = "lens_sis__source_smooth"
data_resolution = "lsst"
def make_pipeline(name, phase_folders, optimizer_class=af.MultiNest):
phase1 = al.PhaseImaging(
phase_name="phase_1",
phase_folders=phase_folders,
galaxies=dict(
lens=al.GalaxyModel(redshift=0.5, mass=al.mp.SphericalIsothermal),
source=al.GalaxyModel(redshift=1.0, light=al.lp.EllipticalSersic),
),
positions_threshold=0.0,
optimizer_class=optimizer_class,
)
phase1.optimizer.const_efficiency_mode = True
phase1.optimizer.n_live_points = 30
phase1.optimizer.sampling_efficiency = 0.8
return al.PipelineDataset(name, phase1)
if __name__ == "__main__":
import sys
runner.run(
sys.modules[__name__],
positions=[[(1.6, 0.0), (0.0, 1.6), (-1.6, 0.0), (0.0, -1.6)]],
)
