Related papers: Shifted Lanczos method for quadratic forms with He…
Inhomogeneous dynamical mean-field theory has been employed to solve many interesting strongly interacting problems from transport in multilayered devices to the properties of ultracold atoms in a trap. The main computational step,…
We study approximations to the Moreau envelope -- and infimal convolutions more broadly -- based on Laplace's method, a classical tool in analysis which ties certain integrals to suprema of their integrands. We believe the connection…
The matrix logarithm, when applied to Hermitian positive definite matrices, is concave with respect to the positive semidefinite order. This operator concavity property leads to numerous concavity and convexity results for other matrix…
Variational procedure is developed that yields lowest frequencies of small-amplitude oscillations of classical Hamiltonian systems. Genuine Lanczos recursion is generalized to treat related non-Hermitian eigenvalue problems.
The $k$-step Lanczos bidiagonalization reduces a matrix $A\in\mathbb{R}^{m\times n}$ into a bidiagonal form $B_k\in\mathbb{R}^{(k+1)\times k}$ while generates two orthonormal matrices $U_{k+1}\in\mathbb{R}^{m\times (k+1)}$ and…
The concept of Gauss quadrature can be generalized to approximate linear functionals with complex moments. Following the existing literature, this survey will revisit such generalization. It is well known that the (classical) Gauss…
A fundamental class of matrix optimization problems that arise in many areas of science and engineering is that of quadratic optimization with orthogonality constraints. Such problems can be solved using line-search methods on the Stiefel…
We extend all cohomological invariants of similarity classes of quadratic forms to anti-hermitian forms over a quaternion algebra. This uses the fact that such invariants can be lifted to Witt invariants, which can be described as…
We consider the approximation of $B^T (A+sI)^{-1} B$ for large s.p.d. $A\in\mathbb{R}^{n\times n}$ with dense spectrum and $B\in\mathbb{R}^{n\times p}$, $p\ll n$. We target the computations of Multiple-Input Multiple-Output (MIMO) transfer…
We explore the Jacobi Last Multiplier as a means for deriving the Lagrangian of a fourth-order differential equation. In particular we consider the classical problem of the Pais-Uhlenbeck oscillator and write down the accompanying…
We review the numerical analysis' understanding of Krylov subspace methods for solving (non-hermitian) systems of equations and discuss its implications for lattice gauge theory computations using the example of the Wilson fermion matrix.…
Quadratic descent of hermitian and skew hermitian forms over division algebras with involution of the first kind in arbitrary characteristic is investigated and a criterion, in terms of systems of quadratic forms, is obtained. A refined…
Reduction theory has played a major role in the study of Hamiltonian systems. On the other hand, the Hamilton-Jacobi theory is one of the main tools to integrate the dynamics of certain Hamiltonian problems and a topic of research on its…
We obtain an expression for the error in the approximation of $f(A) \boldsymbol{b}$ and $\boldsymbol{b}^T f(A) \boldsymbol{b}$ with rational Krylov methods, where $A$ is a symmetric matrix, $\boldsymbol{b}$ is a vector and the function $f$…
Shifted Laplacian multigrid preconditioner has become a tool du jour for solving highly indefinite Helmholtz equations. The idea is to add a complex damping to the original Helmholtz operator and then apply a multigrid processing to the…
We present an efficient method to solve the impurity Hamiltonians involved in Dynamical Mean-Field Theory at low but finite temperature, based on the extension of the Lanczos algorithm from ground state properties alone to excited states.…
In this paper, we discuss a general approach to find periodic solutions bifurcating from equilibrium points of classical Vlasov systems. The main access to the problem is chosen through the Hamiltonian representation of any Vlasov system,…
A new iterative method is developed to numerically calculate the periodic, matched beam envelope solution of the coupled Kapchinskij-Vladimirskij (KV) equations describing the transverse evolution of a beam in a periodic, linear focusing…
Gutzwiller-projected fermionic states can be efficiently implemented within quantum Monte Carlo calculations to define extremely accurate variational wave functions for Heisenberg models on frustrated two-dimensional lattices, not only for…
For several classes of mathematical models that yield linear systems, the splitting of the matrix into its Hermitian and skew Hermitian parts is naturally related to properties of the underlying model. This is particularly so for…