BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Carlos Beltrán appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We show an alternative proof of the sharpest known lower bound for the logarithmic energy on the unit sphere $\mathbb{S}^2$. We then generalize this proof to get new lower bounds for the Green energy on the unit $n$-sphere $\mathbb{S}^n$.
In this paper, we propose a new structured Grassmannian constellation for noncoherent communications over single-input multiple-output (SIMO) Rayleigh block-fading channels. The constellation, which we call Grass-Lattice, is based on a measure preserving mapping from the unit hypercube to the Grassmannian of lines. The constellation structure allows for on-the-fly symbol generation, low-complexity decoding, and simple bit-to-symbol Gray coding. Simulation results show that Grass-Lattice has symbol and bit error rate performance close to that of a numerically optimized unstructured constellation, and is more power efficient than other structured constellations proposed in the literature and a coherent pilot-based scheme.
We consider the design of multiuser constellations for a multiple access channel (MAC) with K users, with M antennas each, that transmit simultaneously to a receiver equipped with N antennas through a Rayleigh block-fading channel, when no channel state information (CSI) is available to either the transmitter or the receiver. In full-diversity scenarios where the coherence time is at least T>= (K+1)M, the proposed constellation design criterion is based on the asymptotic expression of the multiuser pairwise error probability (PEP) derived by Brehler and Varanasi. In non-full diversity scenarios, for which the previous PEP expression is no longer valid, the proposed design criteria are based on proxies of the PEP recently proposed by Ngo and Yang. Although both the PEP expression and its bounds or proxies were previously considered intractable for optimization, in this work we derive their respective unconstrained gradients. These gradients are in turn used in the optimization of the proposed cost functions in different Riemannian manifolds representing different power constraints. In particular, in addition to the standard unitary space-time modulation (USTM) leading to optimization on the Grassmann manifold, we consider a more relaxed per-codeword power constraint leading to optimization on the so-called oblique manifold, and an average power constraint leading to optimization on the so-called trace manifold. Equipped with these theoretical tools, we design multiuser constellations for the MIMO MAC in full-diversity and non-full-diversity scenarios with state-of-the-art performance in terms of symbol error rate (SER).
We prove the existence of an open set of $n_1\times n_2 \times n_3$ tensors of rank $r$ on which a popular and efficient class of algorithms for computing tensor rank decompositions based on a reduction to a linear matrix pencil, typically followed by a generalized eigendecomposition, is arbitrarily numerically forward unstable. Our analysis shows that this problem is caused by the fact that the condition number of the tensor rank decomposition can be much larger for $n_1 \times n_2 \times 2$ tensors than for the $n_1\times n_2 \times n_3$ input tensor. Moreover, we present a lower bound for the limiting distribution of the condition number of random tensor rank decompositions of third-order tensors. The numerical experiments illustrate that for random tensor rank decompositions one should anticipate a loss of precision of a few digits.
We study expected Riesz s-energies and linear statistics of some determinantal processes on the sphere. In particular, we compute the expected Riesz and logarithmic energies of the determinantal processes given by the reproducing kernel of the space of spherical harmonics. This kernel defines the so called harmonic ensemble on the sphere. With these computations we improve previous estimates for the discrete minimal energy of configurations of points in the sphere. We prove a comparison result for Riesz 2-energies of points defined through determinantal point processes associated to isotropic kernels. As a corollary we get that the Riesz 2-energy of the harmonic ensemble is optimal among ensembles defined by isotropic kernels with the same trace. Finally, we study the variance of smooth and rough linear statistics for the harmonic ensemble and compare the results with the variance for the spherical ensemble.
We consider the problem of allocating a finite number of heat sources in the n-dimensional sphere. When only one such source -assumed to be of infinite temperature- is placed and assuming a constant cooling rate in the sphere, we prove that a (essentially) unique solution exists: the Constant Laplacian potential (CL-potential). Actually, this potential can be defined intrinsically in any CROSS (such as the real or complex projective spaces), providing a natural alternative to Riesz's potentials in manifolds lacking a standard isometric embedding into some Euclidean space. We describe an integral form of the corresponding CL-energy for the case of the sphere and prove a relation of minimizing configurations with separation distance and cap discrepancy. It follows that minimal configurations for the Riesz energy are asymptotically minimizing for the CL-energy.
We describe algorithms for computing eigenpairs (eigenvalue-eigenvector pairs) of a complex $n\times n$ matrix $A$. These algorithms are numerically stable, strongly accurate, and theoretically efficient (i.e., polynomial-time). We do not believe they outperform in practice the algorithms currently used for this computational problem. The merit of our paper is to give a positive answer to a long-standing open problem in numerical linear algebra.
Smale's 17th problem asks for an algorithm which finds an approximate zero of polynomial systems in average polynomial time (see Smale 2000). The main progress on Smale's problem is Beltrán-Pardo (2011) and Bürgisser-Cucker (2010). In this paper we will improve on both approaches and we prove an important intermediate result. Our main results are Theorem 1 on the complexity of a randomized algorithm which improves the result of Beltrán-Pardo (2011), Theorem 2 on the average of the condition number of polynomial systems which improves the estimate found in Bürgisser-Cucker (2010), and Theorem 3 on the complexity of finding a single zero of polynomial systems. This last Theorem is the main result of Bürgisser-Cucker (2010). We give a proof of it relying only on homotopy methods, thus removing the need for the elimination theory methods used in Bürgisser-Cucker (2010). We build on methods developed in Armentano et al. (2015).
We describe two algorithms for the eigenvalue, eigenvector problem which, on input a Gaussian matrix with complex entries, finish with probability 1 and in average polynomial time.
In this paper, we study the number of different interference alignment (IA) solutions in a K-user multiple-input multiple-output (MIMO) interference channel, when the alignment is performed via beamforming and no symbol extensions are allowed. We focus on the case where the number of IA equations matches the number of variables. In this situation, the number of IA solutions is finite and constant for any channel realization out of a zero-measure set and, as we prove in the paper, it is given by an integral formula that can be numerically approximated using Monte Carlo integration methods. More precisely, the number of alignment solutions is the scaled average of the determinant of a certain Hermitian matrix related to the geometry of the problem. Interestingly, while the value of this determinant at an arbitrary point can be used to check the feasibility of the IA problem, its average (properly scaled) gives the number of solutions. For single-beam systems the asymptotic growth rate of the number of solutions is analyzed and some connections with classical combinatorial problems are presented. Nonetheless, our results can be applied to arbitrary interference MIMO networks, with any number of users, antennas and streams per user.
In our previous paper [SIMAX 31 n.3 1491-1506(2010)], we studied the condition metric in the space of maximal rank matrices. Here, we show that this condition metric induces a Lipschitz-Riemann structure on that space. After investigating geodesics in such a nonsmooth structure, we show that the inverse of the smallest singular value of a matrix is a log-convex function along geodesics (Theorem 1). We also show that a similar result holds for the solution variety of linear systems (Theorem 31). Some of our intermediate results, such as Theorem 12, on the second covariant derivative or Hessian of a function with symmetries on a manifold, and Theorem 29 on piecewise self-convex functions, are of independent interest. Those results were motivated by our investigations on the com- plexity of path-following algorithms for solving polynomial systems.
We describe, for the first time, a completely rigorous homotopy (path--following) algorithm (in the Turing machine model) to find approximate zeros of systems of polynomial equations. If the coordinates of the input systems and the initial zero are rational our algorithm involves only rational computations and if the homotopy is well posed an approximate zero with integer coordinates of the target system is obtained. The total bit complexity is linear in the length of the path in the condition metric, and polynomial in the logarithm of the maximum of the condition number along the path, and in the size of the input.
Given a homotopy connecting two polynomial systems we provide a rigorous algorithm for tracking a regular homotopy path connecting an approximate zero of the start system to an approximate zero of the target system. Our method uses recent results on the complexity of homotopy continuation rooted in the alpha theory of Smale. Experimental results obtained with the implementation in the numerical algebraic geometry package of Macaulay2 demonstrate the practicality of the algorithm. In particular, we confirm the theoretical results for random linear homotopies and illustrate the plausibility of a conjecture by Shub and Smale on a good initial pair.
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 singular value of A denoted sigma_n(A). When this smallest singular value has multiplicity 1, the function A -> log (sigma_n(A)^(-2)) is a convex function with respect to the condition Riemannian structure that is t -> log (sigma_n(A(t))^(-2)) is convex, in the usual sense for any geodesic A(t). In a more abstract setting, a function alpha defined on a Riemannian manifold (M,<,>) is said to be self-convex when log alpha (gamma(t)) is convex for any geodesic in (M,<,>). Necessary and sufficient conditions for self-convexity are given when alpha is C^2. When alpha(x) = d(x,N)^(-2) where d(x,N) is the distance from x to a C^2 submanifold N of R^j we prove that alpha is self-convex when restricted to the largest open set of points x where there is a unique closest point in N to x. We also show, using this more general notion, that the square of the condition number ||A|||_F / sigma_n(A) is self-convex in projective space and the solution variety.