Start with \( \hat{P}_k(0) = P(0) \).

Real and non-negative coefficients \( a_{\ell,k} \) and \( c_{\ell,k} \) are determined.

For each time instant \( n \geq 1 \), perform:

For each node \( k = 1, 2, ..., K \), perform:

\( \hat{\psi}_k(n+1) = \hat{P}_k(n) + \mu_k \sum_{\ell \in \mathcal{N}_k} c_{\ell,k} \left[ |x_\ell(n)|^2 - \hat{P}_k(n) \right] \)

\( \hat{P}_k(n+1) = \sum_{\ell \in \mathcal{N}_k} a_{\ell,k} \hat{\psi}_\ell(n+1) \)

Decide \( H_0: \hat{P}_k(n+1) < \gamma \) or \( H_1: \hat{P}_k(n+1) > \gamma \).

For threshold selection, refer to equation (32) in [42].