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)]], )