Related papers: Computation of the nearest structured matrix tripl…
Training deep neural networks (DNNs) can be difficult due to the occurrence of vanishing/exploding gradients during weight optimization. To avoid this problem, we propose a class of DNNs stemming from the time discretization of Hamiltonian…
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…
We consider problems of estimation of structured covariance matrices, and in particular of matrices with a Toeplitz structure. We follow a geometric viewpoint that is based on some suitable notion of distance. To this end, we overview and…
Explicit and semi-explicit geometric integration schemes for dissipative perturbations of Hamiltonian systems are analyzed. The dissipation is characterized by a small parameter $\epsilon$, and the schemes under study preserve the…
We consider the problem of computing the closest stable/unstable non-negative matrix to a given real matrix. This problem is important in the study of linear dynamical systems, numerical methods, etc. The distance between matrices is…
Linear discriminant analysis (LDA) is a classical method for dimensionality reduction, where discriminant vectors are sought to project data to a lower dimensional space for optimal separability of classes. Several recent papers have…
Optimization tasks are crucial in statistical machine learning. Recently, there has been great interest in leveraging tools from dynamical systems to derive accelerated and robust optimization methods via suitable discretizations of…
The structural invariant subspaces of the discrete-time singular Hamiltonian system are used in 1] to give an analytic nonrecursive expression of all the admissible trajectories. A deeper insight into the features of these subspaces,…
Several important algorithms for machine learning and data analysis use pairwise distances as input. On Riemannian manifolds these distances may be prohibitively costly to compute, in particular for large datasets. To tackle this problem,…
Riemannian optimization uses local methods to solve optimization problems whose constraint set is a smooth manifold. A linear step along some descent direction usually leaves the constraints, and hence retraction maps are used to…
Given a structured matrix $A$ we study the problem of finding the closest normal matrix with the same structure. The structures of our interest are: Hamiltonian, skew-Hamiltonian, per-Hermitian, and perskew-Hermitian. We develop a…
For a general formulation of linearised hybrid inverse problems in impedance tomography, the qualitative properties of the solutions are analysed. Using an appropriate scalar pseudo-differential formulation, the problems are shown to permit…
Let $S(A)$ denote the orbit of a complex or real matrix $A$ under a certain equivalence relation such as unitary similarity, unitary equivalence, unitary congruences etc. Efficient gradient-flow algorithms are constructed to determine the…
In this note we discuss several questions concerning the implementation of overdetermined least-squares collocation methods for higher-index differential algebraic equations (DAEs). Since higher-index DAEs lead to ill-posed problems in…
In this paper implicit and explicit exact difference schemes (EDS) for system $\textbf{x}' = A\textbf{x}$ of three linear differential equations with constant coefficients are constructed. Numerical simulations for stiff problem and for…
Many discrete optimization problems amount to selecting a feasible set of edges of least weight. We consider in this paper the context of spatial graphs where the positions of the vertices are uncertain and belong to known uncertainty sets.…
Discrete gradient methods are a powerful tool for the time discretization of dynamical systems, since they are structure-preserving regardless of the form of the total energy. In this work, we discuss the application of discrete gradient…
Consider $n \times n$ matrix $A$ and a set $\Lambda$ consisting of $k \le n$ prescribed complex numbers. Lippert (2010) in a challenging article, studied geometrically the spectral norm distance from $A$ to the set $\Lambda$ and constructed…
We analyze infinite-dimensional Hamiltonian systems corresponding to partial differential equations on one-dimensional spatial domains formulated with formally skew-adjoint Hamiltonian operators and nonlinear Hamiltonian density. In various…
For a prototype quadratic Hamiltonian describing a driven, dissipative system, exact matrix elements of the reduced density matrix are obtained from a generating function in terms of the normal characteristic functions. The approach is…