Beatrice Meini

Fabio Durastante, Beatrice Meini, Luca Saluzzi

Time-evolving networks, or temporal networks, play a crucial role in modeling dynamic interactions across various domains, including biology, social sciences, and information technology. Unlike static networks, these systems undergo continuous changes in topology and edge weights, influencing processes such as information flow, transportation efficiency, and neural activity. Understanding and controlling these networks are essential for predicting future behavior and optimizing dynamic processes. This work focuses on the problem of dynamic centrality, a measure of node importance in time-dependent networks. Specifically, we address how to steer network centrality to a desired state by making minimal modifications to the network structure. This problem is formulated as an optimal control problem for an ordinary differential equation, either matrix- or vector-based, where the control acts on network edges. The proposed framework generalizes centrality control problems studied in static networks and leverages the Pontryagin Maximum Principle for efficient solutions. For large-scale problems, the required matrix-function actions are approximated by Krylov-type techniques, avoiding the explicit formation of dense matrix functions. Numerical experiments on synthetic and real temporal networks show that the proposed framework can effectively steer receive centrality under prescribed control constraints.

D. A. Bini, S. Dendievel, G. Latouche et al.

The Erlangian approximation of Markovian fluid queues leads to the problem of computing the matrix exponential of a subgenerator having a block-triangular, block-Toeplitz structure. To this end, we propose some algorithms which exploit the Toeplitz structure and the properties of generators. Such algorithms allow to compute the exponential of very large matrices, which would otherwise be untreatable with standard methods. We also prove interesting decay properties of the exponential of a generator having a block-triangular, block-Toeplitz structure.

Dario A. Bini, Bruno Iannazzo, Beatrice Meini et al.

Component-wise accurate algorithms for computing the principal square root of an M-matrix are designed in terms of triplet representations. A triplet representation of an M-matrix $A$ is the triple $(P, {\bf u},{\bf v})$, where the matrix $P$ is such that $p_{ij}=-a_{ij}$ for $i\ne j$, $p_{ii}=0$, and ${\bf u}>0$, ${\bf v}\ge 0$ are two vectors such that $A{\bf u}={\bf v}$. It is shown that if $A$ is an M-matrix representable by a triplet, then its principal square root exists and is an M-matrix represented by a triplet as well. New versions of the Cyclic Reduction and the Incremental Newton iterations are provided in terms of triplets, to compute the principal matrix square root of $A$. It is shown that these algorithms are component-wise numerically stable independently of the singularity of $A$ and of its condition number. Numerical experiments are shown to confirm the component-wise stability.

Stefano Cipolla, Fabio Durastante, Miryam Gnazzo et al.

Reversibility is a key property of Markov chains, central to algorithms such as Metropolis-Hastings and other MCMC methods. Yet many applications yield non-reversible chains, motivating the problem of approximating them by reversible ones with minimal modification. We formulate this task as a matrix nearness problem and focus on the practically relevant case of sparse transition matrices. The resulting optimization problem is a quadratic programming problem, and numerical experiments illustrate the effectiveness of the approach. This framework provides a principled way to enforce reversibility and sparsity patterns in Markov chains with applications in MCMC, computational chemistry, and data-driven modeling.

Fabio Durastante, Miryam Gnazzo, Beatrice Meini

We address the algorithmic problem of determining the reversible Markov chain $\tilde X$ that is closest to a given Markov chain $X$, with an identical stationary distribution. More specifically, $\tilde X$ is the reversible Markov chain with the closest transition matrix, in the Frobenius norm, to the transition matrix of $X$. To compute the transition matrix of $\tilde X$, we propose a novel approach based on Riemannian optimization. Our method introduces a modified multinomial manifold endowed with a prescribed stationary vector, while also satisfying the detailed balance conditions, all within the framework of the Fisher metric. We evaluate the performance of the proposed approach in comparison with an existing quadratic programming method and demonstrate its effectiveness through a series of synthetic experiments, as well as in the construction of a reversible Markov chain from transition count data obtained via direct estimation from a stochastic differential equation.