Main result
For primitive matrices (i.e., connected networks, no dead-ends/sinks)
\[\mathbf{h} \propto AA^T \mathbf{h} = \lambda_h AA^T \mathbf{h}\]
\[\mathbf{au} \propto A^TA \mathbf{au} = \lambda_{au} A^TA \mathbf{au}\]
convergence is assured;
dominant \(\lambda\) is unique and
values for \(\mathbf{h}\) and \(\mathbf{au}\) will be all positive, as desired.
(negative values have no interpretation here)