class GridPhase(af.as_grid_search(al.PhaseImaging, parallel=True)):
@property
def grid_priors(self):
return [
self.model.galaxies.subhalo.mass.centre_0,
self.model.galaxies.subhalo.mass.centre_1,
]
class GridPhase(
af.as_grid_search(phase_class=toy.PhaseImaging,
parallel=parallel)):
@property
def grid_priors(self):
return [
self.model.gaussians.sub_gaussian.centre_0,
self.model.gaussians.sub_gaussian.centre_1,
]
class GridPhase(
af.as_grid_search(phase_class=al.PhaseImaging,
parallel=slam.setup_subhalo.parallel)):
@property
def grid_priors(self):
return [
self.model.galaxies.subhalo.mass.centre_0,
self.model.galaxies.subhalo.mass.centre_1,
]
class GridPhase(
af.as_grid_search(
phase_class=al.PhaseInterferometer,
number_of_cores=slam.setup_subhalo.number_of_cores,
number_of_steps=slam.setup_subhalo.number_of_steps,
)):
@property
def grid_priors(self):
return [
self.model.galaxies.subhalo.mass.centre_0,
self.model.galaxies.subhalo.mass.centre_1,
]
class GridPhase(af.as_grid_search(al.PhaseImaging)):
@property
def grid_priors(self):
return [
self.variable.galaxies.subhalo.mass.centre_0,
self.variable.galaxies.subhalo.mass.centre_1,
]
def customize_priors(self, results):
### Lens Mass, PL -> PL, Shear -> Shear ###
# self.galaxies.lens = results.from_phase('phase_1').\
# constant.galaxies.lens
### Lens Subhalo, Adjust priors to physical masses (10^6 - 10^10) and concentrations (6-24)
self.galaxies.subhalo.mass.mass_at_200 = af.LogUniformPrior(
lower_limit=10.0e6, upper_limit=10.0e9)
self.galaxies.subhalo.mass.centre_0 = af.UniformPrior(
lower_limit=-2.0, upper_limit=2.0)
self.galaxies.subhalo.mass.centre_1 = af.UniformPrior(
lower_limit=-2.0, upper_limit=2.0)
### Source Light, Sersic -> Sersic ###
self.galaxies.source.light.centre = (
results.from_phase("phase_1").variable_absolute(
a=0.05).galaxies.source.light.centre)
self.galaxies.source.light.intensity = (
results.from_phase("phase_1").variable_relative(
r=0.5).galaxies.source.light.intensity)
self.galaxies.source.light.effective_radius = (
results.from_phase("phase_1").variable_relative(
r=0.5).galaxies.source.light.effective_radius)
self.galaxies.source.light.sersic_index = (
results.from_phase("phase_1").variable_relative(
r=0.5).galaxies.source.light.sersic_index)
self.galaxies.source.light.axis_ratio = (
results.from_phase("phase_1").variable_absolute(
a=0.1).galaxies.source.light.axis_ratio)
self.galaxies.source.light.phi = (
results.from_phase("phase_1").variable_absolute(
a=20.0).galaxies.source.light.phi)
class GridPhase(af.as_grid_search(al.PhaseImaging, parallel=True)):
@property
def grid_priors(self):
return [
self.variable.galaxies.subhalo.mass.centre_0,
self.variable.galaxies.subhalo.mass.centre_1,
]
def customize_priors(self, results):
### Lens Subhalo, Adjust priors to physical masses (10^6 - 10^10) and concentrations (6-24)
self.galaxies.subhalo.mass.kappa_s = af.UniformPrior(
lower_limit=0.0005, upper_limit=0.2)
self.galaxies.subhalo.mass.centre_0 = af.UniformPrior(
lower_limit=-2.0, upper_limit=2.0)
self.galaxies.subhalo.mass.centre_1 = af.UniformPrior(
lower_limit=-2.0, upper_limit=2.0)
class GridPhase(af.as_grid_search(al.PhaseImaging)):
@property
def grid_priors(self):
return [self.model.lens.bulge.sersic_index]
def customize_priors(self, results):
self.galaxies.lens.disk = results.from_phase(
"phase_1").instance.lens.disk
self.galaxies.lens.bulge.centre_0 = af.UniformPrior(
lower_limit=-0.01, upper_limit=0.01)
self.galaxies.lens.bulge.centre_1 = af.UniformPrior(
lower_limit=-0.01, upper_limit=0.01)
self.galaxies.lens.bulge.axis_ratio = af.UniformPrior(
lower_limit=0.79, upper_limit=0.81)
self.galaxies.lens.bulge.phi = af.UniformPrior(lower_limit=-1.0,
upper_limit=1.0)
self.galaxies.lens.bulge.intensity = af.UniformPrior(
lower_limit=0.99, upper_limit=1.01)
self.galaxies.lens.bulge.effective_radius = af.UniformPrior(
lower_limit=1.25, upper_limit=1.35)
