Related papers: Index of a matrix, complex logarithms, and multidi…
Consider the $2n$-by-$2n$ matrix $M=(m_{i,j})_{i,j=1}^{2n}$ with $m_{i,j} = 1$ for $i,j$ satisfying $|2i-2n-1|+|2j-2n-1| \leq 2n$ and $m_{i,j} = 0$ for all other $i,j$, consisting of a central diamond of 1's surrounded by 0's. When $n \geq…
An efficient method is described to handle mesh indexes in multidimensional problems like numerical integration of partial differential equations, lattice model simulations, and determination of atomic neighbor lists. By creating an…
We analyze an algorithm for computing a skew-Hermitian logarithm of a unitary matrix. This algorithm is very easy to implement using standard software and it works well even for unitary matrices with no spectral conditions assumed. Certain…
In this paper we study the problem -\Delta u =\left(\frac{2+\alpha}{2}\right)^2\abs{x}^{\alpha}f(\lambda,u), & \hbox{in}B_1 \\ u > 0, & \hbox{in}B_1 u = 0, & \hbox{on} \partial B_1 where $B_1$ is the unit ball of $\R^2$, $f$ is a smooth…
We revisit the Callias index formula for Dirac-type operators $L$ in odd space dimension $n$, and prove that \begin{align} \text{ind} \, (L) =\bigg(\frac{i}{8\pi}\bigg)^{\frac{n-1}{2}}\frac{1}{2(\frac{n-1}{2})!} \lim_{\Lambda…
Given a symmetric matrix $M$ and a vector $\lambda$, we present new bounds on the Frobenius-distance utility of the Gaussian mechanism for approximating $M$ by a matrix whose spectrum is $\lambda$, under $(\varepsilon,\delta)$-differential…
We present an algorithm computing the determinant of an integer matrix A. The algorithm is introspective in the sense that it uses several distinct algorithms that run in a concurrent manner. During the course of the algorithm partial…
The classic graphical Cheeger inequalities state that if $M$ is an $n\times n$ symmetric doubly stochastic matrix, then \[ \frac{1-\lambda_{2}(M)}{2}\leq\phi(M)\leq\sqrt{2\cdot(1-\lambda_{2}(M))} \] where…
A multiset $\Lambda=\{\lambda_1,\ldots,\lambda_n\}$ of complex numbers is said to be realizable whenever there exists a nonnegative matrix of order $n$ with spectrum $\Lambda$. One of the broadest criterion that guarantees realizability is…
Let M be a complete n-dimensional Riemannian spin manifold, partitioned by q two-sided hypersurfaces which have a compact transverse intersection N and which in addition satisfy a certain coarse transversality condition. Let E be a…
Let $A(x)=A\_0+x\_1A\_1+...+x\_nA\_n$ be a linear matrix, or pencil, generated by given symmetric matrices $A\_0,A\_1,...,A\_n$ of size $m$ with rational entries. The set of real vectors x such that the pencil is positive semidefinite is a…
For $\boldsymbol{\large {\lambda}} = (\lambda_1,\ldots,\lambda_d)$ with $\lambda_i > 0$, it is proved that \begin{equation*} \prod_{i=1}^d \frac{ 1}{(1- r x_i)^{\lambda_i}} = \frac{\Gamma(|\boldsymbol{\large {\lambda}}|)}{\prod_{i=1}^{d}…
We study the problem of approximating the eigenspectrum of a symmetric matrix $\mathbf A \in \mathbb{R}^{n \times n}$ with bounded entries (i.e., $\|\mathbf A\|_{\infty} \leq 1$). We present a simple sublinear time algorithm that…
Given a matrix $D$ describing the pairwise dissimilarities of a data set, a common task is to embed the data points into Euclidean space. The classical multidimensional scaling (cMDS) algorithm is a widespread method to do this. However,…
In this paper, we propose the tensor Noda iteration (NI) and its inexact version for solving the eigenvalue problem of a particular class of tensor pairs called generalized $\mathcal{M}$-tensor pairs. A generalized $\mathcal{M}$-tensor pair…
In this paper, we study matrix scaling and balancing, which are fundamental problems in scientific computing, with a long line of work on them that dates back to the 1960s. We provide algorithms for both these problems that, ignoring…
The limiting distribution of eigenvalues of N x N random matrices has many applications. One of the most studied ensembles are real symmetric matrices with independent entries iidrv; the limiting rescaled spectral measure (LRSM)…
The notion of the higher rank numerical range $\Lambda_{k}(L(\lambda))$ for matrix polynomials $L(\lambda)=A_{m}\lambda^{m}+...+A_{1}\lambda+A_{0}$ is introduced here and some fundamental geometrical properties are investigated. Further,…
Let $m\ge 2$ be an integer, $K$ an algebraic number field and $\alpha\in K\setminus \{0,-1\}$ with sufficiently small absolute value. In this article, we provide a new lower bound for linear form in…
In this paper, we study the multiple integral $ \displaystyle I= \int_0^1 \int_0^1 \dots \int_0^1 f(x_1+x_2 + \dots +x_n) \, dx_1 \, dx_2 \, \dots \, dx_n$. A general formula of $I$ is presented. As an application, the integral $I$ with…