An efficient, memory-saving approach for the Loewner framework

TL;DR

This paper introduces a memory-efficient, scalable method for the Loewner framework that exploits matrix structure to reduce computational costs from quadratic to near-linear, enabling model reduction on large data sets.

Contribution

The paper develops a novel approach using hierarchically semiseparable matrices to efficiently compute rational models without explicitly forming large matrices.

Findings

Reduces computational complexity to O(N log N)

Significantly lowers memory requirements

Enables handling of large, highly-sampled data sets

Abstract

The Loewner framework is one of the most successful data-driven model order reduction techniques. If $N$ is the cardinality of a given data set, the so-called Loewner and shifted Loewner matrices $\mathbb{L}\in\mathbb{C}^{N\times N}$ and $\mathbb{S}\in\mathbb{C}^{N\times N}$ can be defined by solely relying on information encoded in the considered data set and they play a crucial role in the computation of the sought rational model approximation. In particular, the singular value decomposition of a linear combination of $\mathbb{S}$ and $\mathbb{L}$ provides the tools needed to construct accurate models which fulfill important approximation properties with respect to the original data set. However, for highly-sampled data sets, the dense nature of $\mathbb{L}$ and $\mathbb{S}$ leads to numerical difficulties, namely the failure to allocate these matrices in certain memory-limited environments or excessive computational costs. Even though they do not possess any sparsity pattern, the Loewner and shifted Loewner matrices are extremely structured and, in this paper, we show how to fully exploit their Cauchy-like structure to reduce the cost of computing accurate rational models while avoiding the explicit allocation of $\mathbb{L}$ and $\mathbb{S}$. In particular, the use of the \emph{hierarchically semiseparable} format allows us to remarkably lower both the computational cost and the memory requirements of the Loewner framework obtaining a novel scheme whose costs scale with $N \log N$.