Repository navigation

The discrete adjoint: built into the solvers, driven from the run transcript - #744

Merged

lmoresi merged 28 commits into

developmentfrom

feature/discrete-adjoint

Oct 13, 2026

Merged

The discrete adjoint: built into the solvers, driven from the run transcript#744
lmoresi merged 28 commits into
developmentfrom
feature/discrete-adjoint

Commits on Sep 15, 2026
- feat: the discrete adjoint — in the solvers, and over the run transcript
  
  lmoresi
  and
  claude
  committed
  6c3a976
- feat: a misfit read through a gradient, and the fault-segments example
  
  lmoresi
  and
  claude
  committed
  e305485
- The fault-segments example is listric: a flat decollement, a circular ramp, four strengths
  
  lmoresi
  and
  claude
  committed
  459856e
Commits on Sep 16, 2026
Commits on Sep 17, 2026
- fault_friction: PyVista renders of the true state on the mesh's own triangulation
  
  lmoresi
  and
  claude
  committed
  f7d9ccf
- The fault example records itself; adjoint solves are operators in every transcript view
  
  lmoresi
  and
  claude
  committed
  077af00
- adjoint_solve transposes the matrices and solves forwards with the solver's own KSP
  
  lmoresi
  and
  claude
  committed
  b1d0682
- sensitivity: leave numeric constants as named expressions when peeling, so the kernel is not recompiled per value
  
  lmoresi
  and
  claude
  committed
  57fcdf3
- plot_convergence: read the per-observation data files
  
  lmoresi
  and
  claude
  committed
  4568899
- The adjoint operator is assembled from the transposed pointwise kernels
  
  lmoresi
  and
  claude
  committed
  52deffe
- solver.gradient(misfit, parameters, fields): the steady adjoint as one call, and the adjoint problem in writing
  
  lmoresi
  and
  claude
  committed
  557417b
- A misfit on a boundary, and a parameter that enters through one
  
  lmoresi
  and
  claude
  committed
  30b0938
- A misfit with terms on several domains; the fault example's surface terms are boundary integrals
  
  lmoresi
  and
  claude
  committed
  3ccc2b6
- fault_friction: the optimiser sees the misfit relative to its starting value
  
  lmoresi
  and
  claude
  committed
  c3adaf8
- fault_friction: orientation is an interior observable; the surface gives the strain rate
  
  lmoresi
  and
  claude
  committed
  503989c
Commits on Sep 18, 2026
- uw.adjoint.minimise: PETSc TAO as the inversion driver
  
  lmoresi
  and
  claude
  committed
  ab4b7b0
- fault_friction: scipy stays the default driver, since the note's figures come from it; tao is the switch
  
  lmoresi
  and
  claude
  committed
  53417bb
- fault_friction: noise on the observations, a Tikhonov term, and bounds for TAO's blmvm
  
  lmoresi
  and
  claude
  committed
  183ea41
- fault_friction: the objective is a negative log posterior, with the prior given by its width
  
  lmoresi
  and
  claude
  committed
  3be3f18
Commits on Sep 19, 2026
- fault_friction: the example as a notebook, in kilometres, with the run header in the declared units
  
  lmoresi
  and
  claude
  committed
  b019c75
Commits on Sep 20, 2026
- adjoint: a parameter in a Dirichlet datum gets its reaction term
  
  lmoresi
  and
  claude
  committed
  6fb6055
- adjoint: the datum term through the transcript driver, and a refusal for data on other fields
  
  lmoresi
  and
  claude
  committed
  e6e4e61
Commits on Sep 24, 2026
- Merge development into the adjoint branch: model.step, the cell radius, and the
  lmoresi
  committed
  82e697b
Commits on Oct 12, 2026
- Merge development into the discrete-adjoint branch
  
  lmoresi
  and
  claude
  committed
  4649653
- Adjoint derivatives go through diff_wrt_field, not sympy.diff
  
  lmoresi
  and
  claude
  committed
  33a4b19