{"title":"Computing the Conformal Barycenter","authors":"J. Cantarella, Henrik Schumacher","doi":"10.1137/21M1449282","DOIUrl":"https://doi.org/10.1137/21M1449282","url":null,"abstract":"The conformal barycenter of a point cloud on the sphere at infinity of the Poincare ball model of hyperbolic space is a hyperbolic analogue of the geometric median of a point cloud in Euclidean space. It was defined by Douady and Earle as part of a construction of a conformally natural way to extend homeomorphisms of the circle to homeomorphisms of the disk, and it plays a central role in Millson and Kapovich's model of the configuration space of cyclic linkages with fixed edgelengths. \u0000In this paper we consider the problem of computing the conformal barycenter. Abikoff and Ye have given an iterative algorithm for measures on $mathbb{S}^1$ which is guaranteed to converge. We analyze Riemannian versions of Newton's method computed in the intrinsic geometry of the Poincare ball model. We give Newton-Kantorovich (NK) conditions under which we show that Newton's method with fixed step size is guaranteed to converge quadratically to the conformal barycenter for measures on any $mathbb{S}^d$ (including infinite-dimensional spheres). For measures given by $n$ atoms on a finite dimensional sphere which obey the NK conditions, we give an explicit linear bound on the computation time required to approximate the conformal barycenter to fixed error. We prove that our NK conditions hold for all but exponentially few $n$ atom measures. For all measures with a unique conformal barycenter we show that a regularized Newton's method with line search will always converge (eventually superlinearly) to the conformal barycenter. Though we do not have hard time bounds for this algorithm, experiments show that it is extremely efficient in practice and in particular much faster than the Abikoff-Ye iteration.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"70 1","pages":"503-530"},"PeriodicalIF":1.2,"publicationDate":"2020-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87427417","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Carlos Am'endola, Kathlén Kohn, Philipp Reichenbach, A. Seigal
{"title":"Invariant Theory and Scaling Algorithms for Maximum Likelihood Estimation","authors":"Carlos Am'endola, Kathlén Kohn, Philipp Reichenbach, A. Seigal","doi":"10.1137/20M1328932","DOIUrl":"https://doi.org/10.1137/20M1328932","url":null,"abstract":"We show that maximum likelihood estimation in statistics is equivalent to finding the capacity in invariant theory, in two statistical settings: log-linear models and Gaussian transformation families.The former includes the classical independence model while the latter includes matrix normal models and Gaussian graphical models given by transitive directed acyclic graphs. We use stability under group actions to characterize boundedness of the likelihood, and existence and uniqueness of the maximum likelihood estimate. Our approach reveals promising consequences of the interplay between invariant theory and statistics. In particular, existing scaling algorithms from statistics can be used in invariant theory, and vice versa.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"38 1","pages":"304-337"},"PeriodicalIF":1.2,"publicationDate":"2020-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87697676","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Abstract Vergleichsstellensätze for Preordered Semifields and Semirings I","authors":"T. Fritz","doi":"10.1137/22M1498413","DOIUrl":"https://doi.org/10.1137/22M1498413","url":null,"abstract":"Real algebra is usually thought of as the study of certain kinds of preorders on fields and rings. Among its core themes are the separation theorems known as Positivstellens\"atze. However, there is a nascent subfield of real algebra which studies preordered semirings and semifields, which is motivated by applications to probability, graph theory and theoretical computer science, among others. Here, we contribute to this subfield by developing a number of foundational results for it, with two abstract Vergleichsstellens\"atze being our main theorems. Our first Vergleichsstellensatz states that every semifield preorder is the intersection of its total extensions. We apply this to derive our second main result, a Vergleichsstellensatz for certain non-Archimedean preordered semirings in which the homomorphisms to the tropical reals play an important role. We show how this result recovers the existing Vergleichsstellensatz of Strassen and (through the latter) the classical Positivstellensatz of Krivine--Kadison--Dubois.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"48 1","pages":"505-547"},"PeriodicalIF":1.2,"publicationDate":"2020-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90008740","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multi-Trek Separation in Linear Structural Equation Models","authors":"Elina Robeva, Jean-Baptiste Seby","doi":"10.1137/20M1316470","DOIUrl":"https://doi.org/10.1137/20M1316470","url":null,"abstract":"Building on the theory of causal discovery from observational data, we study interactions between multiple (sets of) random variables in a linear structural equation model with non-Gaussian error terms. We give a correspondence between structure in the higher order cumulants and combinatorial structure in the causal graph. It has previously been shown that low rank of the covariance matrix corresponds to trek separation in the graph. Generalizing this criterion to multiple sets of vertices, we characterize when determinants of subtensors of the higher order cumulant tensors vanish. This criterion applies when hidden variables are present as well. For instance, it allows us to identify the presence of a hidden common cause of k of the observed variables.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"157 1","pages":"278-303"},"PeriodicalIF":1.2,"publicationDate":"2020-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76699715","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Distributions of Functions Related to Parametric Integer Optimization","authors":"Timm Oertel, Joseph Paat, R. Weismantel","doi":"10.1137/19M1275954","DOIUrl":"https://doi.org/10.1137/19M1275954","url":null,"abstract":"We consider the asymptotic distribution of the IP sparsity function, which measures the minimal support of optimal IP solutions, and the IP to LP distance function, which measures the distance between optimal IP and LP solutions. To this end, we create a framework for studying the asymptotic distribution of general functions related to integer optimization. While there has been a significant amount of research focused around the extreme values that these functions can attain, little is known about their typical values. Each of these functions is defined for a fixed constraint matrix and objective vector while the right hand sides are treated as input. We show that the typical values of these functions are smaller than the known worst case bounds by providing a spectrum of probability-like results that govern their overall asymptotic distributions.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"23 1","pages":"422-440"},"PeriodicalIF":1.2,"publicationDate":"2019-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87267415","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Representations of Partial Leaf Sets in Phylogenetic Tree Space","authors":"Gillian Grindstaff, Megan Owen","doi":"10.1137/18m1235855","DOIUrl":"https://doi.org/10.1137/18m1235855","url":null,"abstract":"The metric space of phylogenetic trees defined by Billera, Holmes, and Vogtmann [Adv. in Appl. Math., 27 (2001), pp. 733--767], which we refer to as BHV space, provides a natural geometric setting ...","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"58 1","pages":"691-720"},"PeriodicalIF":1.2,"publicationDate":"2019-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83655833","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Whitney Numbers of Combinatorial Geometries and Higher-Weight Dowling Lattices","authors":"A. Ravagnani","doi":"10.1137/20m1382635","DOIUrl":"https://doi.org/10.1137/20m1382635","url":null,"abstract":"We study the Whitney numbers of the first kind of combinatorial geometries. The first part of the paper is devoted to general results relating the Mobius functions of nested atomistic lattices, extending some classical theorems in combinatorics. We then specialize our results to restriction geometries, i.e., to sublattices $mathcal{L}(A)$ of the lattice of subspaces of an $mathbb{F}_q$-linear space, say $X$, generated by a set of projective points $A subseteq X$. In this context, we introduce the notion of subspace distribution, and show that partial knowledge of the latter is equivalent to partial knowledge of the Whitney numbers of $mathcal{L}(A)$. This refines a classical result by Dowling. \u0000The most interesting applications of our results are to be seen in the theory of higher-weight Dowling lattices (HWDLs), to which we dovote the second and most substantive part of the paper. These combinatorial geometries were introduced by Dowling in 1971 in connection with fundamental problems in coding theory, and further studied, among others, by Zaslavsky, Bonin, Kung, Brini, and Games. To date, still very little is known about these lattices. In particular, the techniques to compute their Whitney numbers have not been discovered yet. In this paper, we bring forward the theory of HWDLs, computing their Whitney numbers for new infinite families of parameters. Moreover, we show that the second Whitney numbers of HWDLs are polynomials in the underlying field size $q$, whose coefficients are expressions involving the Bernoulli numbers. This reveals a new link between combinatorics, coding theory, and number theory. We also study the asymptotics of the Whitney numbers of HWDLs as the field size grows, giving upper bounds and exact estimates in some cases. In passing, we obtain new results on the density functions of error-correcting codes.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"66 1","pages":"156-189"},"PeriodicalIF":1.2,"publicationDate":"2019-09-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"72788425","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"p-Adic Integral Geometry","authors":"Avinash Kulkarni, A. Lerário","doi":"10.1137/19m1284737","DOIUrl":"https://doi.org/10.1137/19m1284737","url":null,"abstract":"We prove a $p$-adic version of the Integral Geometry Formula for averaging the intersection of two $p$-adic projective algebraic sets. We apply this result to give bounds on the number of points in the modulo $p^m$ reduction of a projective set (reproving a result by Oesterl'e) and to the study of random $p$-adic polynomial systems of equations.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"1 1","pages":"28-59"},"PeriodicalIF":1.2,"publicationDate":"2019-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83706808","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convex Algebraic Geometry of Curvature Operators","authors":"R. G. Bettiol, Mario Kummer, R. Mendes","doi":"10.1137/20M1350777","DOIUrl":"https://doi.org/10.1137/20M1350777","url":null,"abstract":"We study the structure of the set of algebraic curvature operators satisfying a sectional curvature bound under the light of the emerging field of Convex Algebraic Geometry. More precisely, we determine in which dimensions $n$ this convex semialgebraic set is a spectrahedron or a spectrahedral shadow; in particular, for $ngeq5$, these give new counter-examples to the Helton--Nie Conjecture. Moreover, efficient algorithms are provided if $n=4$ to test membership in such a set. For $ngeq5$, algorithms using semidefinite programming are obtained from hierarchies of inner approximations by spectrahedral shadows and outer relaxations by spectrahedra.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"10 1","pages":"200-228"},"PeriodicalIF":1.2,"publicationDate":"2019-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"72970617","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Algebraic Boundary of the Sonc-Cone","authors":"Jens Forsgård, T. Wolff","doi":"10.1137/20m1325484","DOIUrl":"https://doi.org/10.1137/20m1325484","url":null,"abstract":"In this article, we explore a connection between nonnegativity, the theory of A-discriminants, and tropical geometry. We show that the algebraic strata of the boundary of the sonc cone are parametrized by families of tropical hypersurfaces. Each strata is contained in a rational variety called a positive discriminant. As an application, we characterization generic support sets for which the sonc cone is equal to the sparse nonnegativity cone, and we give a complete description of the semi-algebraic stratification of the boundary of the sonc cone in the univariate case.","PeriodicalId":48489,"journal":{"name":"SIAM Journal on Applied Algebra and Geometry","volume":"1 1","pages":"468-502"},"PeriodicalIF":1.2,"publicationDate":"2019-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75364295","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}