|
| 1 | +# I want to use the autolens_assistant to perform Expectation Propagation |
| 2 | + |
| 3 | +Type: docs |
| 4 | +Target: PyAutoLens |
| 5 | +Repos: |
| 6 | +- PyAutoLens |
| 7 | +- autolens_assistant |
| 8 | +- autolens_workspace |
| 9 | +Difficulty: too-large |
| 10 | +Autonomy: supervised |
| 11 | +Priority: high |
| 12 | +Status: formalised |
| 13 | + |
| 14 | +I want to use the autolens_assistant to perform Expectation Propagation (EP) analysis of a Cosmology science case. |
| 15 | + |
| 16 | +The science case is inferring the Hubble constant from time delay lensed quasars, and a example package already |
| 17 | +containing EP scripts and runs, simulators and other key things is at: |
| 18 | + |
| 19 | +/mnt/c/Users/Jammy/Science/concr/scritps/cosmology and /mnt/c/Users/Jammy/Science/concr/simulator/cosmology.py |
| 20 | + |
| 21 | +This is actually quite a mature project, but its a git old so probably has some API drift (e.g. old PyAutoLens) |
| 22 | +and generally needs to be brushed up. |
| 23 | + |
| 24 | +There are also EP and other examples in autolens_workspace/scripts/guides/modeling/advanced |
| 25 | + |
| 26 | +Can you therefore make a new science project using the autolens_assistant, which uses the main scripts in thee projects |
| 27 | +project to perform the EP fit but also allows for the one by one / graphical fits. However, for now, I want to descope |
| 28 | +down from the Hubble coonstant for now, and have the only goal be to simulate N strong lens using the power-law + shear |
| 29 | +mass model, where the power-law slopes are drawn from a hierarchical distribution. The example scripts should |
| 30 | +then be built around recovering the slopes acurate but, more important the mean and scatter of the slope hierarchical |
| 31 | +parameter. |
| 32 | + |
| 33 | +The goal are: |
| 34 | + |
| 35 | +1) Show that we can do this for large samples without EP, using JAX gradient samples (SPECIFY). |
| 36 | + |
| 37 | +2) Show that we can do this with EP, ideally showing it recovers the same values and errors. |
| 38 | + |
| 39 | +3) We are going to scale up to large lens samples, so make sure all of this can run on RAL and its HPC via the HPC link. |
| 40 | + |
| 41 | +4) To test all the graphical EP diagnostics and result analysis code, especially tyhe recent EP updates we did last week. |
| 42 | + |
| 43 | +<!-- formalised by the Intake (Conception) Agent on 2026-07-16 from file:/tmp/claude-1000/-home-jammy-Code-PyAutoLabs/32708468-5918-4dbc-a763-583805364341/scratchpad/intake_ep_cosmology.md --> |
0 commit comments