arXiv research

Method estimates M-matrices in graphical models with improved accuracy.

DAGMA learns DAGs faster and more accurately using log-determinant acyclicity.

Consider a random vector with finite second moments. If its precision matrix is an M-matrix, then all partial correlations are non-negative. If that random vector is additionally Gaussian, the corresponding Markov random field (GMRF) is called attractive. We study estimation of M-matrices taking the role of inverse sec…

In 1981 Masur proved the existence of a dense geodesic in the moduli space for a Teichmüller space. We prove an analogue theorem for reduced Outer Space endowed with the Lipschitz metric. We also prove two results possibly of independent interest: we show Brun's unordered algorithm weakly converges and from this prove …

We study the following generalized matrix rank estimation problem: given an $n \times n$ matrix and a constant $c \geq 0$, estimate the number of eigenvalues that are greater than $c$. In the distributed setting, the matrix of interest is the sum of $m$ matrices held by separate machines. We show that any deterministic…

Nonnegative matrix factorization (NMF) has become a very popular technique in machine learning because it automatically extracts meaningful features through a sparse and part-based representation. However, NMF has the drawback of being highly ill-posed, that is, there typically exist many different but equivalent facto…

We study rank-1 {L1-norm-based TUCKER2} (L1-TUCKER2) decomposition of 3-way tensors, treated as a collection of $N$ $D \times M$ matrices that are to be jointly decomposed. Our contributions are as follows. i) We prove that the problem is equivalent to combinatorial optimization over $N$ antipodal-binary variables. ii)…

This work is motivated by numerical solutions to Hamilton-Jacobi-Bellman quasi-variational inequalities (HJBQVIs) associated with combined stochastic and impulse control problems. In particular, we consider (i) direct control, (ii) penalized, and (iii) semi-Lagrangian discretization schemes applied to the HJBQVI proble…

We examine the action of the fundamental group $Γ$ of a Riemann surface with $m$ punctures on the middle dimensional homology of a regular fiber in a Lefschetz fibration, and describe to what extent this action can be recovered from the intersection numbers of vanishing cycles. Basis changes for the vanishing cycles re…

We define in the space of n by m matrices of rank n, n less or equal than m, the condition Riemannian structure as follows: For a given matrix A the tangent space of A is equipped with the Hermitian inner product obtained by multiplying the usual Frobenius inner product by the inverse of the square of the smallest sing…

Diagonal Frog: High-order positivity-preserving FD schemes for anisotropic Fokker-Planck equations