Related papers: An Adaptive Shifted Power Method for Computing Gen…
This paper is concerned with the design and analysis of a fully adaptive eigenvalue solver for linear symmetric operators. After transforming the original problem into an equivalent one formulated on $\ell_2$, the space of square summable…
In this paper we develop two new Tensor Alternating Steepest Descent algorithms for tensor completion in the low-rank $\star_{M}$-product format, whereby we aim to reconstruct an entire low-rank tensor from a small number of measurements…
The Jastrow-modified antisymmetric geminal power (JAGP) ansatz in Hilbert space successfully overcomes two key failings of other pairing theories, namely a lack of inter-pair correlations and a lack of multiple resonance structures, while…
Tensor decomposition is a fundamental unsupervised machine learning method in data science, with applications including network analysis and sensor data processing. This work develops a generalized canonical polyadic (GCP) low-rank tensor…
We apply the asymptotic iteration method (AIM) [J. Phys. A: Math. Gen. 36, 11807 (2003)] to solve new classes of second-order homogeneous linear differential equation. In particular, solutions are found for a general class of eigenvalue…
Nonlinear eigenvalue problems (NEPs) present significant challenges due to their inherent complexity and the limitations of traditional linear eigenvalue theory. This paper addresses these challenges by introducing a nonlinear…
The Generalized Eigenvalue Problem (GEVP) has been used extensively in the past in order to reliably extract energy levels from time-dependent Euclidean correlators calculated in Lattice QCD. We propose a formulation of the GEVP in…
Some numerical algorithms for elliptic eigenvalue problems are proposed, analyzed, and numerically tested. The methods combine advantages of the two-grid algorithm, two-space method, the shifted inverse power method, and the polynomial…
We study extensions of compressive sensing and low rank matrix recovery (matrix completion) to the recovery of low rank tensors of higher order from a small number of linear measurements. While the theoretical understanding of low rank…
We compute the arithmetic L-invariants (of Greenberg-Benois) of twists of symmetric powers of p-adic Galois representations attached to Iwahori level Hilbert modular forms (under some technical conditions). Our method uses the automorphy of…
We present an efficient method for finding the independent invariant tensors of a gauge theory. Our method uses a theorem relating invariant tensors and D-flat directions in field space. We apply our method to several examples-- SO(3) with…
We consider the problem of learning high dimensional polynomial transformations of Gaussians. Given samples of the form $p(x)$, where $x\sim N(0, \mathrm{Id}_r)$ is hidden and $p: \mathbb{R}^r \to \mathbb{R}^d$ is a function where every…
An exact charged solution with axial symmetry is obtained in the teleparallel equivalent of general relativity (TEGR). The associated metric has the structure function $G(\xi)=1-{\xi}^2-2mA{\xi}^3-q^2A^2{\xi}^4$. The fourth order nature of…
We find the set of generalized symmetries associated with the free graviton theory in four dimensions. These are generated by gauge invariant topological operators that violate Haag duality in ring-like regions. As expected from general QFT…
In a recent work [Manucci, Unger, ArXiv e-print 2404.10511, 2024], the authors propose using two generalized Lyapunov equations (GLEs) to derive a balancing-based model order reduction~(MOR) method for a general class of switched…
We present a greedy algorithm for computing selected eigenpairs of a large sparse matrix $H$ that can exploit localization features of the eigenvector. When the eigenvector to be computed is localized, meaning only a small number of its…
Generalized inverses of tensors play increasingly important roles in computational mathematics and numerical analysis. It is appropriate to develop the theory of generalized inverses of tensors within the algebraic structure of a ring. In…
Tensor permutation is a fundamental operation widely applied in AI, tensor networks, and related fields. However, it is extremely complex, and different shapes and permutation maps can make a huge difference. SIMD permutation began to be…
Generalized approximate message passing (GAMP) is a promising technique for unknown signal reconstruction of generalized linear models (GLM). However, it requires that the transformation matrix has independent and identically distributed…
We present a new power method to obtain solutions of eigenvalue problems. The method can determine not only the dominant or lowest eigenvalues but also all eigenvalues without the need for a deflation procedure. The method uses a functional…