feat: combined strong+weak lensing features example (weak series step 8)#251
Merged
Conversation
scripts/weak/features/strong_lensing/{simulator,fit,modeling}.py: one tracer
simulated into an Imaging dataset (central arcs, no lens light) and a
WeakDataset shear field to 10" (space-like noise 0.1, 400 galaxies); fit.py
shows the joint log likelihood is the per-dataset sum; modeling.py samples it
with Nautilus via the AnalysisFactor/FactorGraphModel API with fully shared
priors (NumPy graph loop — the weak factor is NumPy). Science framing:
weak shear as the large-radius complement to strong lensing around clusters
(Niemiec 2020 hybrid-Lenstool) and groups (Oguri 2012 SGAS). Navigator
catalogue regenerated in-branch.
Step 8 of the weak-lensing series (autolens_workspace#247).
Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
1 participant
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Summary
Adds
scripts/weak/features/strong_lensing/{simulator,fit,modeling}.py— one tracer simulated into anImagingdataset (central arcs) and a surroundingWeakDatasetshear field (to 10", σ=0.1, 400 galaxies); the joint likelihood is sampled with one Nautilus search via theAnalysisFactor/FactorGraphModelAPI with fully shared priors. Science framing: weak shear as the large-radius complement to strong lensing around clusters (Niemiec 2020 hybrid-Lenstool) and groups (Oguri 2012 SGAS). Companion to PyAutoLabs/PyAutoLens#587 (mixed-graph visualization fix) — merge the library PR first. Tracking issue: #247.Scripts Changed
scripts/weak/features/strong_lensing/simulator.py(new) — both datasets from one tracer.scripts/weak/features/strong_lensing/fit.py(new) —FitImaging+FitWeakwith a shared tracer; joint LL = exact sum (6653.7 + 676.5 = 7330.2); shear profile demo.scripts/weak/features/strong_lensing/modeling.py(new) — factor-graph joint fit; full real Nautilus run recoveredeinstein_radius = 1.5999 (1.5993–1.6005)vs truth 1.6.llms-full.txt/workspace_index.json— catalogue regenerated in-branch.Validation checklist (--auto run)
🤖 Generated with Claude Code