One branch's equation, and the derivatives the Jacobian needs.
Two forms, chosen per branch per iteration:
``dp`` ``(p_up - p_dn) - dp(mdot) = 0``. The usual one.
``choked`` ``K (mdot_ceiling - mdot) = 0``. Written when a component
reports that the available drop has taken it past choking, at
which point mass flow no longer depends on downstream pressure
and the dp form is not invertible. ``K`` carries the row into
pressure units so it is scaled like its neighbours.
Definition at line 602 of file steady.py.