We compute the independence number, zero-error capacity, and the values of the Lov\'asz function and the quantum Lov\'asz function for the quantum graph associated to the partial trace quantum channel $\operatorname{Tr}_n\otimes\mathrm{id}_k\colon\operatorname{B}(\mathbb{C}^n\otimes\mathbb{C}^k)\to\operatorname{B}(\mathbb{C}^k)$.
相關內容
Approximating invariant subspaces of generalized eigenvalue problems (GEPs) is a fundamental computational problem at the core of machine learning and scientific computing. It is, for example, the root of Principal Component Analysis (PCA) for dimensionality reduction, data visualization, and noise filtering, and of Density Functional Theory (DFT), arguably the most popular method to calculate the electronic structure of materials. Given Hermitian $H,S\in\mathbb{C}^{n\times n}$, where $S$ is positive-definite, let $\Pi_k$ be the true spectral projector on the invariant subspace that is associated with the $k$ smallest (or largest) eigenvalues of the GEP $HC=SC\Lambda$, for some $k\in[n]$. We show that we can compute a matrix $\widetilde\Pi_k$ such that $\lVert\Pi_k-\widetilde\Pi_k\rVert_2\leq \epsilon$, in $O\left( n^{\omega+\eta}\mathrm{polylog}(n,\epsilon^{-1},\kappa(S),\mathrm{gap}_k^{-1}) \right)$ bit operations in the floating point model, for some $\epsilon\in(0,1)$, with probability $1-1/n$. Here, $\eta>0$ is arbitrarily small, $\omega\lesssim 2.372$ is the matrix multiplication exponent, $\kappa(S)=\lVert S\rVert_2\lVert S^{-1}\rVert_2$, and $\mathrm{gap}_k$ is the gap between eigenvalues $k$ and $k+1$. To achieve such provable "forward-error" guarantees, our methods rely on a new $O(n^{\omega+\eta})$ stability analysis for the Cholesky factorization, and a smoothed analysis for computing spectral gaps, which can be of independent interest. Ultimately, we obtain new matrix multiplication-type bit complexity upper bounds for PCA problems, including classical PCA and (randomized) low-rank approximation.
In contrast to regular ordinary differential equations, the problem of accurately setting initial conditions just emerges in the context of differential-algebraic equations where the dynamic degree of freedom of the system is smaller than the absolute dimension of the described process, and the actual lower-dimensional configuration space of the system is deeply implicit. For linear higher-index differential-algebraic equations, we develop an appropriate numerical method based on properties of canonical subspaces and on the so-called geometric reduction. Taking into account the fact that higher-index differential-algebraic equations lead to ill-posed problems in naturally given norms, we modify this approach to serve as transfer conditions from one time-window to the next in a time stepping procedure and combine it with window-wise overdetermined least-squares collocation to construct the first fully numerical solvers for higher-index initial-value problems.
We show that the parameters of a $k$-mixture of inverse Gaussian or gamma distributions are algebraically identifiable from the first $3k-1$ moments, and rationally identifiable from the first $3k+2$ moments. Our proofs are based on Terracini's classification of defective surfaces, careful analysis of the intersection theory of moment varieties, and a recent result on sufficient conditions for rational identifiability of secant varieties by Massarenti--Mella.
Previous papers give accounts of quests for satisfactory formalizations of the classical informal notion of an algorithm and the contemporary informal notion of an interactive algoritm. In this paper, an attempt is made to generalize the results of the former quest to the contemporary informal notion of a concurrent algorithm. The notion of a concurrent proto-algorithm is introduced. The thought is that concurrent algorithms are equivalence classes of concurrent proto-algorithms under an appropriate equivalence relation. Three equivalence relations are defined. Two of them are deemed to be bounds for an appropriate equivalence relation and the third is likely an appropriate one. The connection between concurrency and non-determinism in the presented setting is also addressed.
We propose a fast method for computing the eigenvalue decomposition of a dense real normal matrix $A$. The method leverages algorithms that are known to be efficiently implemented, such as the bidiagonal singular value decomposition and the symmetric eigenvalue decomposition. For symmetric and skew-symmetric matrices, the method reduces to calling the latter, so that its advantages are for orthogonal matrices mostly and, potentially, any other normal matrix. The method relies on the real Schur decomposition of the skew-symmetric part of $A$. To obtain the eigenvalue decomposition of the normal matrix $A$, additional steps depending on the distribution of the eigenvalues are required. We provide a complexity analysis of the method and compare its numerical performance with existing algorithms. In most cases, the method is as fast as obtaining the Hessenberg factorization of a dense matrix. Finally, we evaluate the method's accuracy and provide experiments for the application of a Karcher mean on the special orthogonal group.
We determine the exact error and strong converse exponent for entanglement-assisted classical-quantum channel simulation in worst case input purified distance. The error exponent is expressed as a single-letter formula optimized over sandwiched R\'enyi divergences of order $\alpha \in [1, \infty)$, notably without the need for a critical rate--a sharp contrast to the error exponent for classical-quantum channel coding. The strong converse exponent is expressed as a single-letter formula optimized over sandwiched R\'enyi divergences of order $\alpha\in [\frac{1}{2},1]$. As in the classical work [Oufkir et al., arXiv:2410.07051], we start with the goal of asymptotically expanding the meta-converse for channel simulation in the relevant regimes. However, to deal with non-commutativity issues arising from classical-quantum channels and entanglement-assistance, we critically use various properties of the quantum fidelity, additional auxiliary channel techniques, approximations via Chebyshev inequalities, and entropic continuity bounds.
We focus on a family of nonlinear continuity equations for the evolution of a non-negative density $\rho$ with a continuous and compactly supported nonlinear mobility $\mathrm{m}(\rho)$ not necessarily concave. The velocity field is the negative gradient of the variation of a free energy including internal and confinement energy terms. Problems with compactly supported mobility are often called saturation problems since the values of the density are constrained below a maximal value. Taking advantage of a family of approximating problems, we show the existence of $C_0$-semigroups of $L^1$ contractions. We study the $\omega$-limit of the problem, its most relevant properties, and the appearance of free boundaries in the long-time behaviour. This problem has a formal gradient-flow structure, and we discuss the local/global minimisers of the corresponding free energy in the natural topology related to the set of initial data for the $L^\infty$-constrained gradient flow of probability densities. Furthermore, we analyse a structure preserving implicit finite-volume scheme and discuss its convergence and long-time behaviour.
The incompressible Euler equations are an important model system in computational fluid dynamics. Fast high-order methods for the solution of this time-dependent system of partial differential equations are of particular interest: due to their exponential convergence in the polynomial degree they can make efficient use of computational resources. To address this challenge we describe a novel timestepping method which combines a hybridised Discontinuous Galerkin method for the spatial discretisation with IMEX timestepping schemes, thus achieving high-order accuracy in both space and time. The computational bottleneck is the solution of a (block-) sparse linear system to compute updates to pressure and velocity at each stage of the IMEX integrator. Following Chorin's projection approach, this update of the velocity and pressure fields is split into two stages. As a result, the hybridised equation for the implicit pressure-velocity problem is reduced to the well-known system which arises in hybridised mixed formulations of the Poisson- or diffusion problem and for which efficient multigrid preconditioners have been developed. Splitting errors can be reduced systematically by embedding this update into a preconditioned Richardson iteration. The accuracy and efficiency of the new method is demonstrated numerically for two time-dependent testcases that have been previously studied in the literature.
We prove that multilevel Picard approximations are capable of approximating solutions of semilinear heat equations in $L^{p}$-sense, ${p}\in [2,\infty)$, in the case of gradient-dependent, Lipschitz-continuous nonlinearities, in the sense that the computational effort of the multilevel Picard approximations grow at most polynomially in both the dimension $d$ and the reciprocal $1/\epsilon$ of the prescribed accuracy $\epsilon$.
Unique continuation principles are fundamental properties of elliptic partial differential equations, giving conditions that guarantee that the solution to an elliptic equation must be uniformly zero. Since finite-element discretizations are a natural tool to help gain understanding into elliptic equations, it is natural to ask if such principles also hold at the discrete level. In this work, we prove a version of the unique continuation principle for piecewise-linear and -bilinear finite-element discretizations of the Laplacian eigenvalue problem on polygonal domains in $\mathbb{R}^2$. Namely, we show that any solution to the discretized equation $-\Delta u = \lambda u$ with vanishing Dirichlet and Neumann traces must be identically zero under certain geometric and topological assumptions on the resulting triangulation. We also provide a counterexample, showing that a nonzero \emph{inner solution} exists when the topological assumptions are not satisfied. Finally, we give an application to an eigenvalue interlacing problem, where the space of inner solutions makes an explicit appearance.
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191