hydromodpy.solver.boussinesq.runtimes.petsc_partition#

PETSc runtime for the regularized partition surface-interaction model.

This backend keeps the current head-only Boussinesq residual and solves it with PETSc SNES. Surface interaction follows the regularized partition law

q_ex = G_r(theta) R(balance)

instead of the mixed complementarity (h, q_ex, q_dry) formulation, and its Jacobian is assembled analytically almost everywhere.

Functions

solve_steady_problem(inputs)

Solve one steady nonlinear balance with PETSc on the head-only system.

solve_transient_step(inputs)

Solve one transient implicit step with PETSc on the head-only system.