On Complexity of Stability Analysis in Higher-order Ecological Networks through Tensor Decompositions
Anqi Dong, Can Chen
TL;DR
This work tackles the scalability challenge of stability analysis in ecological networks with higher-order interactions by focusing on the HOGLV model. It develops a tensor-based framework that recasts the HOGLV dynamics through HOSVD, CPD, and TT decompositions to enable efficient computation of the Jacobian at equilibrium and its eigenvalues. The authors derive explicit memory and computational complexity bounds for full, HOSVD-, CPD-, and TT-based representations, highlighting that CPD often minimizes complexity but can be numerically unstable, while TT-based methods offer stable, scalable performance. Numerical experiments illustrate how higher-order terms influence stability and demonstrate the practical advantages of TT representations in large systems. The approach provides a foundation for scalable stability analysis in complex ecological networks and potentially other domains with high-order interactions.
Abstract
Complex ecological networks are often characterized by intricate interactions that extend beyond pairwise relationships. Understanding the stability of higher-order ecological networks is salient for species coexistence, biodiversity, and community persistence. In this article, we present complexity analyses for determining the linear stability of higher-order ecological networks through tensor decompositions. We are interested in the higher-order generalized Lotka-Volterra model, which captures high-order interactions using tensors of varying orders. To efficiently compute Jacobian matrices and thus determine stability in large ecological networks, we exploit various tensor decompositions, including higher-order singular value decomposition, Canonical Polyadic decomposition, and tensor train decomposition, accompanied by in-depth computational and memory complexity analyses. We demonstrate the effectiveness of our framework with numerical examples.
