Nonlinear solver: automatic warm start, model-advertised yield homotopy, and an inf-safe power-mean yield law by lmoresi · Pull Request #441 · underworldcode/underworld3

feat(constitutive): power-mean yield soft-min smoother + rampable delta constant

lmoresi

committed

fix(constitutive): DP consistent-Newton recursion on Max yield + smooth floors

lmoresi

committed

feat(systems): yield_continuation driver — multi-solve delta-continuation (reference impl)

lmoresi

committed

docs(design): nonlinear-solver automatic warm-start + single-parameter yield homotopy

lmoresi

committed

feat(solver): has_solution warm-start status + cold consistent-Newton Picard warm-up (Layer 1a)

lmoresi

committed

docs(solver): nonlinear-solver skill (recipe + config trap list) + mark Layer 1a landed

lmoresi

committed

feat(solver): default the FMG velocity smoother to gmres+sor, fixed-cost V-cycle (Layer 3)

lmoresi

committed

docs(solver): record the measured Layer 3 result and the multigrid-depth method

lmoresi

committed

feat(solver): model-advertised yield homotopy and stokes.solve(homotopy=True) (Layer 2)

lmoresi

committed

docs(solver): Layer 2 landed — one-call homotopy, and the cold-start NaN correction

lmoresi

committed

fix(constitutive): power-mean yield NaN at zero strain rate (rigid/unyielded points)

lmoresi

committed

docs(solver): correct the cold-start diagnosis — inf-safe algebra, not a strain-rate floor

lmoresi

committed

feat(solver): zero_init_guess becomes tri-state and auto-detects cold vs warm (Layer 1b)

lmoresi

committed

docs(solver): all four layers landed; drop the last stale cold-start claims

lmoresi

committed

fix(constitutive): yield_stress_min defaults to 0, floored via symbolic smooth_max

lmoresi

committed

fix(solver): close the critical adversarial-review findings in the homotopy march

lmoresi

committed

fix(constitutive/solver): close review findings M6, M7, M9

lmoresi

committed

fix(solver): close the remaining review minors and add parallel coverage

lmoresi

committed

feat(constitutive): apply the rate regularisation on the DP model; pin the homotopy capability

lmoresi

committed

chore: drop the constitutive_models.py.bak copy committed by accident

lmoresi

committed

fix(test): pin the JIT-constant ramp test to a cold initial guess

lmoresi

committed

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Uh oh!

Nonlinear solver: automatic warm start, model-advertised yield homotopy, and an inf-safe power-mean yield law - #441

Nonlinear solver: automatic warm start, model-advertised yield homotopy, and an inf-safe power-mean yield law#441
lmoresi merged 21 commits into
developmentfrom
feature/nonlinear-warmstart-homotopy

Commits on Jul 28, 2026

feat(constitutive): power-mean yield soft-min smoother + rampable delta constant

fix(constitutive): DP consistent-Newton recursion on Max yield + smooth floors

feat(systems): yield_continuation driver — multi-solve delta-continuation (reference impl)

docs(design): nonlinear-solver automatic warm-start + single-parameter yield homotopy

feat(solver): has_solution warm-start status + cold consistent-Newton Picard warm-up (Layer 1a)

docs(solver): nonlinear-solver skill (recipe + config trap list) + mark Layer 1a landed

feat(solver): default the FMG velocity smoother to gmres+sor, fixed-cost V-cycle (Layer 3)

docs(solver): record the measured Layer 3 result and the multigrid-depth method

feat(solver): model-advertised yield homotopy and stokes.solve(homotopy=True) (Layer 2)

docs(solver): Layer 2 landed — one-call homotopy, and the cold-start NaN correction

fix(constitutive): power-mean yield NaN at zero strain rate (rigid/unyielded points)

docs(solver): correct the cold-start diagnosis — inf-safe algebra, not a strain-rate floor

feat(solver): zero_init_guess becomes tri-state and auto-detects cold vs warm (Layer 1b)

docs(solver): all four layers landed; drop the last stale cold-start claims

fix(constitutive): yield_stress_min defaults to 0, floored via symbolic smooth_max

fix(solver): close the critical adversarial-review findings in the homotopy march

fix(constitutive/solver): close review findings M6, M7, M9

fix(solver): close the remaining review minors and add parallel coverage

feat(constitutive): apply the rate regularisation on the DP model; pin the homotopy capability

chore: drop the constitutive_models.py.bak copy committed by accident

fix(test): pin the JIT-constant ramp test to a cold initial guess