hydromodpy.solver.boussinesq.runtimes.stationary_picard_lscheme

Strict experimental Picard/L-scheme initializer for steady Boussinesq.

This module is intentionally isolated from the default runtime selection path. It is an investigation helper for difficult steady obstacle cases. The method keeps the original head-only problem definition: no artificial minimum saturated thickness is introduced and no drainage/surface conductance is added.

The Picard iterate solves a lagged linear transmissivity problem, adds a purely algorithmic L-scheme diagonal damping term, relaxes the update, projects the head into the physical bounds, and evaluates the strict bounded residual.

This module is the public facade for the picard sub-package, which splits the runtime in four concerns:

• picard.lscheme: the bounded relaxed Picard iterate (linear solve loop).

• picard.picard: VI cycles and strict VI obstacle assembly.

• picard.diagnostics: shared math/geometry helpers and record dataclasses.

• picard.io: JSON/CSV diagnostic writers.