hydromodpy.solver.boussinesq.runtimes.picard.lscheme#

Picard/L-scheme iterate: bounded relaxed linear system + projection.

Functions

accept_relaxed_candidate(*, mesh, head_old, ...)

Halve omega until the candidate residual fits the growth factor.

bounded_picard_lscheme(inputs, *[, ...])

Run the experimental strict bounded Picard/L-scheme initializer.

iteration_record(*, mesh, result, iteration, ...)

Build one PicardIterationRecord from a Picard iterate result.

lagged_internal_tau(mesh, *, thickness, ...)

Return one lagged internal-edge transmissibility coefficient.

picard_linear_system(mesh, *, head_old, ...)

Assemble A h_new = rhs with strict lagged physical transmissivity.

picard_summary(*, mesh, result, options, ...)

Return a JSON-friendly summary of one Picard/L-scheme run.

resolve_lscheme_L(mesh, *, base_diag, lscheme_L)

Resolve the L-scheme diagonal stabilization.

solve_sparse_system(matrix, rhs)

Solve a sparse linear system; raise on warnings to detect singular solves.