-
Detection of manatee vocalisations using the Audio Spectrogram Transformer
Authors:
Stefano Schiappacasse,
Taco de Wolff,
Yann Henaut,
Regina Cervera,
Aviva Charles,
Felipe Tobar
Abstract:
The Antillean manatee (\emph{Trichechus manatus}) is an endangered herbivorous aquatic mammal whose role as an ecological balancer and umbrella species underscores the importance of its conservation. An innovative approach to monitor manatee populations is passive acoustic monitoring (PAM), where vocalisations are extracted from submarine audio. We propose a novel end-to-end approach to detect man…
▽ More
The Antillean manatee (\emph{Trichechus manatus}) is an endangered herbivorous aquatic mammal whose role as an ecological balancer and umbrella species underscores the importance of its conservation. An innovative approach to monitor manatee populations is passive acoustic monitoring (PAM), where vocalisations are extracted from submarine audio. We propose a novel end-to-end approach to detect manatee vocalisations building on the Audio Spectrogram Transformer (AST). In a transfer learning spirit, we fine-tune AST to detect manatee calls by redesigning its filterbanks and adapting a real-world dataset containing partial positive labels. Our experimental evaluation reveals the two key features of the proposed model: i) it performs on par with the state of the art without requiring hand-tuned denoising or detection stages, and ii) it can successfully identify missed vocalisations in the training dataset, thus reducing the workload of expert bioacoustic labellers. This work is a preliminary relevant step to develop novel, user-friendly tools for the conservation of the different species of manatees.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
Revealing hidden physical nonclassicality with nonnegative polynomials
Authors:
Ties-A. Ohst,
Benjamin Yadin,
Birte Ostermann,
Timo de Wolff,
Otfried Gühne,
Hai-Chau Nguyen
Abstract:
Understanding quantum phenomena which go beyond classical concepts is a focus of modern quantum physics. Here, we show how the theory of nonnegative polynomials emerging around Hilbert's 17th problem, can be used to optimally exploit data capturing the nonclassical nature of light. Specifically, we show that nonnegative polynomials can reveal nonclassicality in data even when it is hidden from sta…
▽ More
Understanding quantum phenomena which go beyond classical concepts is a focus of modern quantum physics. Here, we show how the theory of nonnegative polynomials emerging around Hilbert's 17th problem, can be used to optimally exploit data capturing the nonclassical nature of light. Specifically, we show that nonnegative polynomials can reveal nonclassicality in data even when it is hidden from standard detection methods up to now. Moreover, the abstract language of nonnegative polynomials also leads to a unified mathematical approach to nonclassicality for light and spin systems, allowing us to map methods for one to the other. Conversely, the physical problems arising also inspire several mathematical insights into characterisation of nonnegative polynomials.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
A Speed-up for Helsgaun's TSP Heuristic by Relaxing the Positive Gain Criterion
Authors:
Sabrina C. L. Ammann,
Birte Ostermann,
Sebastian Stiller,
Timo de Wolff
Abstract:
The Traveling Salesman Problem (TSP) is one of the most extensively researched and widely applied combinatorial optimization problems. It is NP-hard even in the symmetric and metric case. Building upon elaborate research, state-of-the-art exact solvers such as CONCORDE can solve TSP instances with several ten thousand vertices. A key ingredient for these integer programming approaches are fast heu…
▽ More
The Traveling Salesman Problem (TSP) is one of the most extensively researched and widely applied combinatorial optimization problems. It is NP-hard even in the symmetric and metric case. Building upon elaborate research, state-of-the-art exact solvers such as CONCORDE can solve TSP instances with several ten thousand vertices. A key ingredient for these integer programming approaches are fast heuristics to find a good initial solution, in particular the Lin-Kernighan-Helsgaun (LKH) heuristic. For instances with few hundred vertices heuristics like LKH often find an optimal solution. In this work we develop variations of LKH that perform significantly better on large instances. LKH repeatedly improves an initially random tour by exchanging edges along alternating circles. Thereby, it respects several criteria designed to quickly find alternating circles that give a feasible improvement of the tour. Among those criteria, the positive gain criterion stayed mostly untouched in previous research. It requires that, while constructing an alternating circle, the total gain has to be positive after each pair of edges. We relax this criterion carefully leading to improvement steps hitherto undiscovered by LKH. We confirm this improvement experimentally via extensive simulations on various benchmark libraries for TSP. Our computational study shows that for large instances our method is on average 13% faster than the latest version of LKH.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Realistic Runtime Analysis for Quantum Simplex Computation
Authors:
Sabrina Ammann,
Maximilian Hess,
Debora Ramacciotti,
Sándor P. Fekete,
Paulina L. A. Goedicke,
David Gross,
Andreea Lefterovici,
Tobias J. Osborne,
Michael Perk,
Antonio Rotundo,
S. E. Skelton,
Sebastian Stiller,
Timo de Wolff
Abstract:
In recent years, strong expectations have been raised for the possible power of quantum computing for solving difficult optimization problems, based on theoretical, asymptotic worst-case bounds. Can we expect this to have consequences for Linear and Integer Programming when solving instances of practically relevant size, a fundamental goal of Mathematical Programming, Operations Research and Algor…
▽ More
In recent years, strong expectations have been raised for the possible power of quantum computing for solving difficult optimization problems, based on theoretical, asymptotic worst-case bounds. Can we expect this to have consequences for Linear and Integer Programming when solving instances of practically relevant size, a fundamental goal of Mathematical Programming, Operations Research and Algorithm Engineering? Answering this question faces a crucial impediment: The lack of sufficiently large quantum platforms prevents performing real-world tests for comparison with classical methods.
In this paper, we present a quantum analog for classical runtime analysis when solving real-world instances of important optimization problems. To this end, we measure the expected practical performance of quantum computers by analyzing the expected gate complexity of a quantum algorithm. The lack of practical quantum platforms for experimental comparison is addressed by hybrid benchmarking, in which the algorithm is performed on a classical system, logging the expected cost of the various subroutines that are employed by the quantum versions. In particular, we provide an analysis of quantum methods for Linear Programming, for which recent work has provided asymptotic speedup through quantum subroutines for the Simplex method. We show that a practical quantum advantage for realistic problem sizes would require quantum gate operation times that are considerably below current physical limitations.
△ Less
Submitted 16 November, 2023;
originally announced November 2023.
-
Initial Application of SONC to Lyapunov Stability of Dynamical Systems
Authors:
Janin Heuer,
Timo de Wolff
Abstract:
Certifying the stability of dynamical systems is a central and challenging task in control theory and systems analysis. To tackle these problems we present an algorithmic approach to finding polynomial Lyapunov functions. Our method relies on sums of nonnegative circuit functions (SONC), a certificate of nonnegativity of real polynomials. We show that both the problem of verifying as well as the m…
▽ More
Certifying the stability of dynamical systems is a central and challenging task in control theory and systems analysis. To tackle these problems we present an algorithmic approach to finding polynomial Lyapunov functions. Our method relies on sums of nonnegative circuit functions (SONC), a certificate of nonnegativity of real polynomials. We show that both the problem of verifying as well as the more difficult task of finding Lyapunov functions can be carried out via relative entropy programming when using SONC certificates. This approach is analogue yet independent to finding Lyapunov functions via sums of squares (SOS) certificates and semidefinite programming. Furthermore, we explore whether using the related, recently introduced DSONC certificate is advantageous compared to SONC for this type of problem. We implemented our results, and present examples to show their applicability.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
A Generalized Muirhead Inequality and Symmetric Sums of Nonnegative Circuits
Authors:
Janin Heuer,
Ngoc Mai Tran,
Timo de Wolff
Abstract:
Circuit polynomials are a certificate of nonnegativity for real polynomials, which can be derived via a generalization of the classical inequality of arithmetic and geometric means. In this article, we show that similarly nonnegativity of symmetric real polynomials can be certified via a generalization of the classical Muirhead inequality. Moreover, we show that a nonnegative symmetric polynomial…
▽ More
Circuit polynomials are a certificate of nonnegativity for real polynomials, which can be derived via a generalization of the classical inequality of arithmetic and geometric means. In this article, we show that similarly nonnegativity of symmetric real polynomials can be certified via a generalization of the classical Muirhead inequality. Moreover, we show that a nonnegative symmetric polynomial admits a decomposition into sums of nonnegative circuit polynomials if and only if it satisfies said generalized Muirhead condition. The latter re-proves a result by Moustrou, Naumann, Riener, Theobald, and Verdure for the case of the symmetric group in a shortened and more elementary way.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
Computationally-efficient initialisation of GPs: The generalised variogram method
Authors:
Felipe Tobar,
Elsa Cazelles,
Taco de Wolff
Abstract:
We present a computationally-efficient strategy to initialise the hyperparameters of a Gaussian process (GP) avoiding the computation of the likelihood function. Our strategy can be used as a pretraining stage to find initial conditions for maximum-likelihood (ML) training, or as a standalone method to compute hyperparameters values to be plugged in directly into the GP model. Motivated by the fac…
▽ More
We present a computationally-efficient strategy to initialise the hyperparameters of a Gaussian process (GP) avoiding the computation of the likelihood function. Our strategy can be used as a pretraining stage to find initial conditions for maximum-likelihood (ML) training, or as a standalone method to compute hyperparameters values to be plugged in directly into the GP model. Motivated by the fact that training a GP via ML is equivalent (on average) to minimising the KL-divergence between the true and learnt model, we set to explore different metrics/divergences among GPs that are computationally inexpensive and provide hyperparameter values that are close to those found via ML. In practice, we identify the GP hyperparameters by projecting the empirical covariance or (Fourier) power spectrum onto a parametric family, thus proposing and studying various measures of discrepancy operating on the temporal and frequency domains. Our contribution extends the variogram method developed by the geostatistics literature and, accordingly, it is referred to as the generalised variogram method (GVM). In addition to the theoretical presentation of GVM, we provide experimental validation in terms of accuracy, consistency with ML and computational complexity for different kernels using synthetic and real-world data.
△ Less
Submitted 26 April, 2023; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Learning Variational Models with Unrolling and Bilevel Optimization
Authors:
Christoph Brauer,
Niklas Breustedt,
Timo de Wolff,
Dirk A. Lorenz
Abstract:
In this paper we consider the problem of learning variational models in the context of supervised learning via risk minimization. Our goal is to provide a deeper understanding of the two approaches of learning of variational models via bilevel optimization and via algorithm unrolling. The former considers the variational model as a lower level optimization problem below the risk minimization probl…
▽ More
In this paper we consider the problem of learning variational models in the context of supervised learning via risk minimization. Our goal is to provide a deeper understanding of the two approaches of learning of variational models via bilevel optimization and via algorithm unrolling. The former considers the variational model as a lower level optimization problem below the risk minimization problem, while the latter replaces the lower level optimization problem by an algorithm that solves said problem approximately. Both approaches are used in practice, but unrolling is much simpler from a computational point of view. To analyze and compare the two approaches, we consider a simple toy model, and compute all risks and the respective estimators explicitly. We show that unrolling can be better than the bilevel optimization approach, but also that the performance of unrolling can depend significantly on further parameters, sometimes in unexpected ways: While the stepsize of the unrolled algorithm matters a lot (and learning the stepsize gives a significant improvement), the number of unrolled iterations plays a minor role.
△ Less
Submitted 6 September, 2023; v1 submitted 26 September, 2022;
originally announced September 2022.
-
Parameter region for multistationarity in $n-$site phosphorylation networks
Authors:
Elisenda Feliu,
Nidhi Kaihnsa,
Timo de Wolff,
Oğuzhan Yürük
Abstract:
Multisite phosphorylation is a signaling mechanism well known to give rise to multiple steady states, a property termed multistationarity. When phosphorylation occurs in a sequential and distributive manner, we obtain a family of networks indexed by the number of phosphorylation sites $n$. This work addresses the problem of understanding the parameter region where this family of networks displays…
▽ More
Multisite phosphorylation is a signaling mechanism well known to give rise to multiple steady states, a property termed multistationarity. When phosphorylation occurs in a sequential and distributive manner, we obtain a family of networks indexed by the number of phosphorylation sites $n$. This work addresses the problem of understanding the parameter region where this family of networks displays multistationarity, by focusing on the projection of this region onto the set of kinetic parameters. The problem is phrased in the context of real algebraic geometry and reduced to studying whether a polynomial, defined as the determinant of a parametric matrix of size three, attains negative values over the positive orthant. The coefficients of the polynomial are functions of the kinetic parameters. For any $n$, we provide sufficient conditions for the polynomial to be positive and hence, preclude multistationarity, and also sufficient conditions for it to attain negative values and hence, enable multistationarity. These conditions are derived by exploiting the structure of the polynomial, its Newton polytope, and employing circuit polynomials. A relevant consequence of our results is that the set of kinetic parameters that enable or preclude multistationarity are both connected for all $n$.
△ Less
Submitted 5 March, 2023; v1 submitted 17 June, 2022;
originally announced June 2022.
-
The Duality of SONC: Advances in Circuit-based Certificates
Authors:
Janin Heuer,
Timo de Wolff
Abstract:
The cone of sums of nonnegative circuits (SONCs) is a subset of the cone of nonnegative polynomials / exponential sums, which has been studied extensively in recent years. In this article, we construct a subset of the SONC cone which we call the DSONC cone. The DSONC cone can be seen as an extension of the dual SONC cone; membership can be tested via linear programming. We show that the DSONC cone…
▽ More
The cone of sums of nonnegative circuits (SONCs) is a subset of the cone of nonnegative polynomials / exponential sums, which has been studied extensively in recent years. In this article, we construct a subset of the SONC cone which we call the DSONC cone. The DSONC cone can be seen as an extension of the dual SONC cone; membership can be tested via linear programming. We show that the DSONC cone is a proper, full-dimensional cone, we provide a description of its extreme rays, and collect several properties that parallel those of the SONC cone. Moreover, we show that functions in the DSONC cone cannot have real zeros, which yields that DSONC cone does not intersect the boundary of the SONC cone. Furthermore, we discuss the intersection of the DSONC cone with the SOS and SDSOS cones. Finally, we show that circuit functions in the boundary of the DSONC cone are determined by points of equilibria, which hence are the analogues to singular points in the primal SONC cone, and relate the DSONC cone to tropical geometry.
△ Less
Submitted 8 April, 2022;
originally announced April 2022.
-
Towards Optimally Weighted Physics-Informed Neural Networks in Ocean Modelling
Authors:
Taco de Wolff,
Hugo Carrillo,
Luis Martí,
Nayat Sanchez-Pi
Abstract:
The carbon pump of the world's ocean plays a vital role in the biosphere and climate of the earth, urging improved understanding of the functions and influences of the ocean for climate change analyses. State-of-the-art techniques are required to develop models that can capture the complexity of ocean currents and temperature flows. This work explores the benefits of using physics-informed neural…
▽ More
The carbon pump of the world's ocean plays a vital role in the biosphere and climate of the earth, urging improved understanding of the functions and influences of the ocean for climate change analyses. State-of-the-art techniques are required to develop models that can capture the complexity of ocean currents and temperature flows. This work explores the benefits of using physics-informed neural networks (PINNs) for solving partial differential equations related to ocean modeling; such as the Burgers, wave, and advection-diffusion equations. We explore the trade-offs of using data vs. physical models in PINNs for solving partial differential equations. PINNs account for the deviation from physical laws in order to improve learning and generalization. We observed how the relative weight between the data and physical model in the loss function influence training results, where small data sets benefit more from the added physics information.
△ Less
Submitted 16 June, 2021;
originally announced June 2021.
-
Re-entrant tensegrity: A three-periodic, chiral, tensegrity structure that is auxetic
Authors:
Mathias Oster,
Marcelo A. Dias,
Timo de Wolff,
Myfanwy E. Evans
Abstract:
We present a three-periodic, chiral, tensegrity structure and demonstrate that it is auxetic. Our tensegrity structure is constructed using the chiral symmetry Π+ cylinder packing, transforming the cylinders themselves to the elastic elements and cylinder contacts to incompressible rods. The resulting structure displays local re-entrant geometry at its vertices, and is shown to be auxetic when mod…
▽ More
We present a three-periodic, chiral, tensegrity structure and demonstrate that it is auxetic. Our tensegrity structure is constructed using the chiral symmetry Π+ cylinder packing, transforming the cylinders themselves to the elastic elements and cylinder contacts to incompressible rods. The resulting structure displays local re-entrant geometry at its vertices, and is shown to be auxetic when modelled as an equilibrium configuration of spatial constraints subject to a quasi-static deformation. When the tensegrity structure is subsequently modelled as a lattice material with elastic elements, the auxetic behaviour is again confirmed through finite element modelling as well as demonstrated with a 3d printed example. The cubic symmetry of the structure means that the behaviour is independent of the chosen stretching direction and the auxetic behaviour is observed in both perpendicular directions. This structure could be the simplest three-dimensional analogue to the two-dimensional re-entrant honeycomb. This, alongside the chirality of the structure, make it an interesting design target for multifunctional materials.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Assessment of event-triggered policies of nonpharmaceutical interventions based on epidemiological indicators
Authors:
Carla Castillo-Laborde,
Taco de Wolff,
Pedro Gajardo,
Rodrigo Lecaros,
Gerard Olivar,
Hector Ramirez C
Abstract:
Nonpharmaceutical interventions (NPI) such as banning public events or instituting lockdowns have been widely applied around the world to control the current COVID-19 pandemic. Typically, this type of intervention is imposed when an epidemiological indicator in a given population exceeds a certain threshold. Then, the nonpharmaceutical intervention is lifted when the levels of the indicator used h…
▽ More
Nonpharmaceutical interventions (NPI) such as banning public events or instituting lockdowns have been widely applied around the world to control the current COVID-19 pandemic. Typically, this type of intervention is imposed when an epidemiological indicator in a given population exceeds a certain threshold. Then, the nonpharmaceutical intervention is lifted when the levels of the indicator used have decreased sufficiently. What is the best indicator to use? In this paper, we propose a mathematical framework to try to answer this question. More specifically, the proposed framework permits to assess and compare different event-triggered controls based on epidemiological indicators. Our methodology consists of considering some outcomes that are consequences of the nonpharmaceutical interventions that a decision maker aims to make as low as possible. The peak demand for intensive care units (ICU) and the total number of days in lockdown are examples of such outcomes. If an epidemiological indicator is used to trigger the interventions, there is naturally a trade-off between the outcomes that can be seen as a curve parameterized by the trigger threshold to be used. The computation of these curves for a group of indicators then allows the selection of the best indicator the curve of which dominates the curves of the other indicators. This methodology is illustrated using indicators in the context of COVID-19 using deterministic compartmental models in discrete-time, although the framework can be adapted for a larger class of models.
△ Less
Submitted 30 October, 2020;
originally announced November 2020.
-
Computing the Real Isolated Points of an Algebraic Hypersurface
Authors:
Huu Phuoc Le,
Mohab Safey El Din,
Timo de Wolff
Abstract:
Let $\mathbb{R}$ be the field of real numbers. We consider the problem of computing the real isolated points of a real algebraic set in $\mathbb{R}^n$ given as the vanishing set of a polynomial system. This problem plays an important role for studying rigidity properties of mechanism in material designs. In this paper, we design an algorithm which solves this problem. It is based on the computatio…
▽ More
Let $\mathbb{R}$ be the field of real numbers. We consider the problem of computing the real isolated points of a real algebraic set in $\mathbb{R}^n$ given as the vanishing set of a polynomial system. This problem plays an important role for studying rigidity properties of mechanism in material designs. In this paper, we design an algorithm which solves this problem. It is based on the computations of critical points as well as roadmaps for answering connectivity queries in real algebraic sets. This leads to a probabilistic algorithm of complexity $(nd)^{O(n\log(n))}$ for computing the real isolated points of real algebraic hypersurfaces of degree $d$. It allows us to solve in practice instances which are out of reach of the state-of-the-art.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Evaluation of Pool-based Testing Approaches to Enable Population-wide Screening for COVID-19
Authors:
Timo de Wolff,
Dirk Pflüger,
Michael Rehme,
Janin Heuer,
Martin-Immanuel Bittner
Abstract:
Background: Rapid testing for an infection is paramount during a pandemic to prevent continued viral spread and excess morbidity and mortality. This study aimed to determine whether alternative testing strategies based on sample pooling can increase the speed and throughput of screening for SARS-CoV-2.
Methods: A mathematical modelling approach was chosen to simulate six different testing strate…
▽ More
Background: Rapid testing for an infection is paramount during a pandemic to prevent continued viral spread and excess morbidity and mortality. This study aimed to determine whether alternative testing strategies based on sample pooling can increase the speed and throughput of screening for SARS-CoV-2.
Methods: A mathematical modelling approach was chosen to simulate six different testing strategies based on key input parameters (infection rate, test characteristics, population size, testing capacity etc.). The situations in five countries (US, DE, UK, IT and SG) currently experiencing COVID-19 outbreaks were simulated to reflect a broad variety of population sizes and testing capacities. The primary study outcome measurements that were finalised prior to any data collection were time and number of tests required; number of cases identified; and number of false positives.
Findings: The performance of all tested methods depends on the input parameters, i.e. the specific circumstances of a screening campaign. To screen one tenth of each country's population at an infection rate of 1% - e.g. when prioritising frontline medical staff and public workers -, realistic optimised testing strategies enable such a campaign to be completed in ca. 29 days in the US, 71 in the UK, 25 in Singapore, 17 in Italy and 10 in Germany (ca. eight times faster compared to individual testing). When infection rates are considerably lower, or when employing an optimal, yet logistically more complex pooling method, the gains are more pronounced. Pool-based approaches also reduces the number of false positive diagnoses by 50%.
Interpretation: The results of this study provide a clear rationale for adoption of pool-based testing strategies to increase speed and throughput of testing for SARS-CoV-2. The current individual testing approach unnecessarily wastes valuable time and resources.
△ Less
Submitted 8 October, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
Global Optimization via the Dual SONC Cone and Linear Programming
Authors:
Mareike Dressler,
Janin Heuer,
Helen Naumann,
Timo de Wolff
Abstract:
Using the dual cone of sums of nonnegative circuits (SONC), we provide a relaxation of the global optimization problem to minimize an exponential sum and, as a special case, a multivariate real polynomial. Our approach builds on two key observations. First, that the dual SONC cone is contained in the primal one. Hence, containment in this cone is a certificate of nonnegativity. Second, we show tha…
▽ More
Using the dual cone of sums of nonnegative circuits (SONC), we provide a relaxation of the global optimization problem to minimize an exponential sum and, as a special case, a multivariate real polynomial. Our approach builds on two key observations. First, that the dual SONC cone is contained in the primal one. Hence, containment in this cone is a certificate of nonnegativity. Second, we show that membership in the dual cone can be verified by a linear program. We implement the algorithm and present initial experimental results comparing our method to existing approaches.
△ Less
Submitted 22 October, 2020; v1 submitted 21 February, 2020;
originally announced February 2020.
-
Gaussian process imputation of multiple financial series
Authors:
Taco de Wolff,
Alejandro Cuevas,
Felipe Tobar
Abstract:
In Financial Signal Processing, multiple time series such as financial indicators, stock prices and exchange rates are strongly coupled due to their dependence on the latent state of the market and therefore they are required to be jointly analysed. We focus on learning the relationships among financial time series by modelling them through a multi-output Gaussian process (MOGP) with expressive co…
▽ More
In Financial Signal Processing, multiple time series such as financial indicators, stock prices and exchange rates are strongly coupled due to their dependence on the latent state of the market and therefore they are required to be jointly analysed. We focus on learning the relationships among financial time series by modelling them through a multi-output Gaussian process (MOGP) with expressive covariance functions. Learning these market dependencies among financial series is crucial for the imputation and prediction of financial observations. The proposed model is validated experimentally on two real-world financial datasets for which their correlations across channels are analysed. We compare our model against other MOGPs and the independent Gaussian process on real financial data.
△ Less
Submitted 11 February, 2020;
originally announced February 2020.
-
MOGPTK: The Multi-Output Gaussian Process Toolkit
Authors:
Taco de Wolff,
Alejandro Cuevas,
Felipe Tobar
Abstract:
We present MOGPTK, a Python package for multi-channel data modelling using Gaussian processes (GP). The aim of this toolkit is to make multi-output GP (MOGP) models accessible to researchers, data scientists, and practitioners alike. MOGPTK uses a Python front-end, relies on the GPflow suite and is built on a TensorFlow back-end, thus enabling GPU-accelerated training. The toolkit facilitates impl…
▽ More
We present MOGPTK, a Python package for multi-channel data modelling using Gaussian processes (GP). The aim of this toolkit is to make multi-output GP (MOGP) models accessible to researchers, data scientists, and practitioners alike. MOGPTK uses a Python front-end, relies on the GPflow suite and is built on a TensorFlow back-end, thus enabling GPU-accelerated training. The toolkit facilitates implementing the entire pipeline of GP modelling, including data loading, parameter initialization, model learning, parameter interpretation, up to data imputation and extrapolation. MOGPTK implements the main multi-output covariance kernels from literature, as well as spectral-based parameter initialization strategies. The source code, tutorials and examples in the form of Jupyter notebooks, together with the API documentation, can be found at http://github.com/GAMES-UChile/mogptk
△ Less
Submitted 9 February, 2020;
originally announced February 2020.
-
The kinetic space of multistationarity in dual phosphorylation
Authors:
E. Feliu,
N. Kaihnsa,
T. de Wolff,
O. Yürük
Abstract:
Multistationarity in molecular systems underlies switch-like responses in cellular decision making. Determining whether and when a system displays multistationarity is in general a difficult problem. In this work we completely determine the set of kinetic parameters that enable multistationarity in a ubiquitous motif involved in cell signaling, namely a dual phosphorylation cycle. In addition we s…
▽ More
Multistationarity in molecular systems underlies switch-like responses in cellular decision making. Determining whether and when a system displays multistationarity is in general a difficult problem. In this work we completely determine the set of kinetic parameters that enable multistationarity in a ubiquitous motif involved in cell signaling, namely a dual phosphorylation cycle. In addition we show that the regions of multistationarity and monostationarity are both path connected.
We model the dynamics of the concentrations of the proteins over time by means of a parametrized polynomial ordinary differential equation (ODE) system arising from the mass-action assumption. Since this system has three linear first integrals defined by the total amounts of the substrate and the two enzymes, we study for what parameter values the ODE system has at least two positive steady states after suitably choosing the total amounts. We employ a suite of techniques from (real) algebraic geometry, which in particular concern the study of the signs of a multivariate polynomial over the positive orthant and sums of nonnegative circuit polynomials.
△ Less
Submitted 14 April, 2020; v1 submitted 22 January, 2020;
originally announced January 2020.
-
A Polyhedral Homotopy Algorithm For Real Zeros
Authors:
Alperen A. Ergür,
Timo de Wolff
Abstract:
We design a homotopy continuation algorithm, that is based on numerically tracking Viro's patchworking method, for finding real zeros of sparse polynomial systems. The algorithm is targeted for polynomial systems with coefficients satisfying certain concavity conditions. It operates entirely over the real numbers and tracks the optimal number of solution paths. In more technical terms; we design a…
▽ More
We design a homotopy continuation algorithm, that is based on numerically tracking Viro's patchworking method, for finding real zeros of sparse polynomial systems. The algorithm is targeted for polynomial systems with coefficients satisfying certain concavity conditions. It operates entirely over the real numbers and tracks the optimal number of solution paths. In more technical terms; we design an algorithm that correctly counts and finds the real zeros of polynomial systems that are located in the unbounded components of the complement of the underlying A-discriminant amoeba.
△ Less
Submitted 31 July, 2022; v1 submitted 4 October, 2019;
originally announced October 2019.
-
Initial Steps in the Classification of Maximal Mediated Sets
Authors:
Jacob Hartzer,
Olivia Röhrig,
Timo de Wolff,
Oğuzhan Yürük
Abstract:
Maximal mediated sets (MMS), introduced by Reznick, are distinguished subsets of lattice points in integral polytopes with even vertices. MMS of Newton polytopes of AGI-forms and nonnegative circuit polynomials determine whether these polynomials are sums of squares. In this article, we take initial steps in classifying MMS both theoretically and practically. Theoretically, we show that MMS of sim…
▽ More
Maximal mediated sets (MMS), introduced by Reznick, are distinguished subsets of lattice points in integral polytopes with even vertices. MMS of Newton polytopes of AGI-forms and nonnegative circuit polynomials determine whether these polynomials are sums of squares. In this article, we take initial steps in classifying MMS both theoretically and practically. Theoretically, we show that MMS of simplices are isomorphic if and only if the simplices generate the same lattice up to permutations. Furthermore, we generalize a result of Iliman and the third author. Practically, we fully characterize the MMS for all simplices of sufficiently small dimensions and maximal 1-norms. In particular, we experimentally prove a conjecture by Reznick for 2 dimensional simplices up to maximal 1-norm 150 and provide indications on the distribution of the density of MMS.
△ Less
Submitted 10 July, 2020; v1 submitted 1 October, 2019;
originally announced October 2019.
-
The algebraic boundary of the sonc cone
Authors:
Jens Forsgård,
Timo de Wolff
Abstract:
In this article, we explore the connections between nonnegativity, the theory of $A$-discriminants, and tropical geometry. For an integral support set $A \subset \mathbb{Z}^n$, we cover the boundary of the sonc-cone by semi-algebraic sets that are parametrized by families of tropical hypersurfaces. As an application, we characterization generic support sets for which the sonc-cone is equal to the…
▽ More
In this article, we explore the connections between nonnegativity, the theory of $A$-discriminants, and tropical geometry. For an integral support set $A \subset \mathbb{Z}^n$, we cover the boundary of the sonc-cone by semi-algebraic sets that are parametrized by families of tropical hypersurfaces. As an application, we characterization generic support sets for which the sonc-cone is equal to the sparse nonnegativity cone, and we describe a semi-algebraic stratification of the boundary of the sonc-cone in the univariate case.
△ Less
Submitted 20 August, 2021; v1 submitted 12 May, 2019;
originally announced May 2019.
-
New Dependencies of Hierarchies in Polynomial Optimization
Authors:
Adam Kurpisz,
Timo de Wolff
Abstract:
We compare four key hierarchies for solving Constrained Polynomial Optimization Problems (CPOP): Sum of Squares (SOS), Sum of Diagonally Dominant Polynomials (SDSOS), Sum of Nonnegative Circuits (SONC), and the Sherali Adams (SA) hierarchies. We prove a collection of dependencies among these hierarchies both for general CPOPs and for optimization problems on the Boolean hypercube. Key results incl…
▽ More
We compare four key hierarchies for solving Constrained Polynomial Optimization Problems (CPOP): Sum of Squares (SOS), Sum of Diagonally Dominant Polynomials (SDSOS), Sum of Nonnegative Circuits (SONC), and the Sherali Adams (SA) hierarchies. We prove a collection of dependencies among these hierarchies both for general CPOPs and for optimization problems on the Boolean hypercube. Key results include for the general case that the SONC and SOS hierarchy are polynomially incomparable, while SDSOS is contained in SONC. A direct consequence is the non-existence of a Putinar-like Positivstellensatz for SDSOS. On the Boolean hypercube, we show as a main result that Schmüdgen-like versions of the hierarchies SDSOS*, SONC*, and SA* are polynomially equivalent. Moreover, we show that SA* is contained in any Schmüdgen-like hierarchy that provides a O(n) degree bound.
△ Less
Submitted 12 March, 2019;
originally announced March 2019.
-
Exact Optimization via Sums of Nonnegative Circuits and Sums of AM/GM Exponentials
Authors:
Victor Magron,
Henning Seidler,
Timo de Wolff
Abstract:
We provide two hybrid numeric-symbolic optimization algorithms, computing exact sums of nonnegative circuits (SONC) and sums of arithmetic-geometric-exponentials (SAGE) decompositions. Moreover, we provide a hybrid numeric-symbolic decision algorithm for polynomials lying in the interior of the SAGE cone. Each framework, inspired by previous contributions of Parrilo and Peyrl, is a rounding-projec…
▽ More
We provide two hybrid numeric-symbolic optimization algorithms, computing exact sums of nonnegative circuits (SONC) and sums of arithmetic-geometric-exponentials (SAGE) decompositions. Moreover, we provide a hybrid numeric-symbolic decision algorithm for polynomials lying in the interior of the SAGE cone. Each framework, inspired by previous contributions of Parrilo and Peyrl, is a rounding-projection procedure.
For a polynomial lying in the interior of the SAGE cone, we prove that the decision algorithm terminates within a number of arithmetic operations, which is polynomial in the number of terms of the input, and linear in the distance to the boundary of the cone. We also provide experimental comparisons regarding the implementation of the two optimization algorithms.
△ Less
Submitted 20 August, 2021; v1 submitted 6 February, 2019;
originally announced February 2019.
-
Computing Elimination Ideals and Discriminants of Likelihood Equations
Authors:
Xiaoxian Tang,
Timo De Wolff,
Rukai Zhao
Abstract:
We develop a probabilistic algorithm for computing elimination ideals of likelihood equations, which is for larger models by far more efficient than directly computing Groebner bases or the interpolation method proposed in the first author's previous work. The efficiency is improved by a theoretical result showing that the sum of data variables appears in most coefficients of the generator polynom…
▽ More
We develop a probabilistic algorithm for computing elimination ideals of likelihood equations, which is for larger models by far more efficient than directly computing Groebner bases or the interpolation method proposed in the first author's previous work. The efficiency is improved by a theoretical result showing that the sum of data variables appears in most coefficients of the generator polynomial of elimination ideal. Furthermore, applying the known structures of Newton polytopes of discriminants, we can also efficiently deduce discriminants of the elimination ideals. For instance, the discriminants of 3 by 3 matrix model and one Jukes-Cantor model in phylogenetics (with sizes over 30 GB and 8 GB text files, respectively) can be computed by our methods.
△ Less
Submitted 12 October, 2018;
originally announced October 2018.
-
An Experimental Comparison of SONC and SOS Certificates for Unconstrained Optimization
Authors:
Henning Seidler,
Timo de Wolff
Abstract:
Finding the minimum of a multivariate real polynomial is a well-known hard problem with various applications. We present a polynomial time algorithm to approximate such lower bounds via sums of nonnegative circuit polynomials (SONC). As a main result, we carry out the first large-scale comparison of SONC, using this algorithm and different geometric programming (GP) solvers, with the classical sum…
▽ More
Finding the minimum of a multivariate real polynomial is a well-known hard problem with various applications. We present a polynomial time algorithm to approximate such lower bounds via sums of nonnegative circuit polynomials (SONC). As a main result, we carry out the first large-scale comparison of SONC, using this algorithm and different geometric programming (GP) solvers, with the classical sums of squares (SOS) approach, using several of the most common semidefinite programming (SDP) solvers. SONC yields bounds competitive to SOS in several cases, but using significantly less time and memory. In particular, SONC/GP can handle much larger problem instances than SOS/SDP.
△ Less
Submitted 25 August, 2018;
originally announced August 2018.
-
Optimization over the Boolean Hypercube via Sums of Nonnegative Circuit Polynomials
Authors:
Mareike Dressler,
Adam Kurpisz,
Timo de Wolff
Abstract:
Various key problems from theoretical computer science can be expressed as polynomial optimization problems over the boolean hypercube. One particularly successful way to prove complexity bounds for these types of problems are based on sums of squares (SOS) as nonnegativity certificates. In this article, we initiate the analysis of optimization problems over the boolean hypercube via a recent, alt…
▽ More
Various key problems from theoretical computer science can be expressed as polynomial optimization problems over the boolean hypercube. One particularly successful way to prove complexity bounds for these types of problems are based on sums of squares (SOS) as nonnegativity certificates. In this article, we initiate the analysis of optimization problems over the boolean hypercube via a recent, alternative certificate called sums of nonnegative circuit polynomials (SONC). We show that key results for SOS based certificates remain valid: First, for polynomials, which are nonnegative over the $n$-variate boolean hypercube with constraints of degree $d$ there exists a SONC certificate of degree at most $n+d$. Second, if there exists a degree $d$ SONC certificate for nonnegativity of a polynomial over the boolean hypercube, then there also exists a short degree $d$ SONC certificate, that includes at most $n^{O(d)}$ nonnegative circuit polynomials.
△ Less
Submitted 27 February, 2018;
originally announced February 2018.
-
Nondegenerate multistationarity in small reaction networks
Authors:
Anne Shiu,
Timo de Wolff
Abstract:
Much attention has been focused in recent years on the following algebraic problem arising from applications: which chemical reaction networks, when taken with mass-action kinetics, admit multiple positive steady states? The interest behind this question is in steady states that are stable. As a step toward this difficult question, here we address the question of multiple nondegenerate positive st…
▽ More
Much attention has been focused in recent years on the following algebraic problem arising from applications: which chemical reaction networks, when taken with mass-action kinetics, admit multiple positive steady states? The interest behind this question is in steady states that are stable. As a step toward this difficult question, here we address the question of multiple nondegenerate positive steady states. Mathematically, this asks whether certain families of parametrized, real, sparse polynomial systems ever admit multiple positive real roots that are simple. Our main results settle this problem for certain types of small networks, and our techniques point the way forward for larger networks.
△ Less
Submitted 5 June, 2018; v1 submitted 31 January, 2018;
originally announced February 2018.
-
The Lattice of Amoebas
Authors:
Jens Forsgård,
Timo de Wolff
Abstract:
We study amoebas of exponential sums as functions of the support set $A$. To any amoeba, we associate a set of approximating sections of amoebas, which we call caissons. We show that a bounded modular lattice of subspaces of a certain vector space induces a lattice structure on the set of caissons. Our results unifies the theories of lopsided amoebas and amoebas of exponential sums. As an applicat…
▽ More
We study amoebas of exponential sums as functions of the support set $A$. To any amoeba, we associate a set of approximating sections of amoebas, which we call caissons. We show that a bounded modular lattice of subspaces of a certain vector space induces a lattice structure on the set of caissons. Our results unifies the theories of lopsided amoebas and amoebas of exponential sums. As an application, we show that our theory of caissons yields improved certificates for existence of certain components of the complement of an amoeba.
△ Less
Submitted 10 November, 2018; v1 submitted 7 November, 2017;
originally announced November 2017.
-
Lopsided Approximation of Amoebas
Authors:
Jens Forsgård,
Laura Felicia Matusevich,
Nathan Mehlhop,
Timo de Wolff
Abstract:
The amoeba of a Laurent polynomial is the image of the corresponding hypersurface under the coordinatewise log absolute value map. In this article, we demonstrate that a theoretical amoeba approximation method due to Purbhoo can be used efficiently in practice. To do this, we resolve the main bottleneck in Purbhoo's method by exploiting relations between cyclic resultants. We use the same approach…
▽ More
The amoeba of a Laurent polynomial is the image of the corresponding hypersurface under the coordinatewise log absolute value map. In this article, we demonstrate that a theoretical amoeba approximation method due to Purbhoo can be used efficiently in practice. To do this, we resolve the main bottleneck in Purbhoo's method by exploiting relations between cyclic resultants. We use the same approach to give an approximation of the Log preimage of the amoeba of a Laurent polynomial using semi-algebraic sets. We also provide a SINGULAR/SAGE implementation of these algorithms, which shows a significant speedup when our specialized cyclic resultant computation is used, versus a general purpose resultant algorithm.
△ Less
Submitted 9 September, 2017; v1 submitted 30 August, 2016;
originally announced August 2016.
-
A Positivstellensatz for Sums of Nonnegative Circuit Polynomials
Authors:
Mareike Dressler,
Sadik Iliman,
Timo de Wolff
Abstract:
Recently, the second and the third author developed sums of nonnegative circuit polynomials (SONC) as a new certificate of nonnegativity for real polynomials, which is independent of sums of squares.
In this article we show that the SONC cone is full-dimensional in the cone of nonnegative polynomials. We establish a Positivstellensatz which guarantees that every polynomial which is positive on a…
▽ More
Recently, the second and the third author developed sums of nonnegative circuit polynomials (SONC) as a new certificate of nonnegativity for real polynomials, which is independent of sums of squares.
In this article we show that the SONC cone is full-dimensional in the cone of nonnegative polynomials. We establish a Positivstellensatz which guarantees that every polynomial which is positive on a given compact, semi-algebraic set can be represented by the constraints of the set and SONC polynomials. Based on this Positivstellensatz we provide a hierarchy of lower bounds converging against the minimum of a polynomial on a given compact set $K$. Moreover, we show that these new bounds can be computed efficiently via interior point methods using results about relative entropy functions.
△ Less
Submitted 17 March, 2017; v1 submitted 20 July, 2016;
originally announced July 2016.
-
An Approach to Constrained Polynomial Optimization via Nonnegative Circuit Polynomials and Geometric Programming
Authors:
Mareike Dressler,
Sadik Iliman,
Timo de Wolff
Abstract:
In this article we combine two developments in polynomial optimization. On the one hand, we consider nonnegativity certificates based on sums of nonnegative circuit polynomials, which were recently introduced by the second and the third author. On the other hand, we investigate geometric programming methods for constrained polynomial optimization problems, which were recently developed by Ghasemi…
▽ More
In this article we combine two developments in polynomial optimization. On the one hand, we consider nonnegativity certificates based on sums of nonnegative circuit polynomials, which were recently introduced by the second and the third author. On the other hand, we investigate geometric programming methods for constrained polynomial optimization problems, which were recently developed by Ghasemi and Marshall. We show that the combination of both results yields a new method to solve certain classes of constrained polynomial optimization problems. We test the new method experimentally and compare it to semidefinite programming in various examples.
△ Less
Submitted 5 June, 2018; v1 submitted 19 February, 2016;
originally announced February 2016.
-
Imaginary projections of polynomials
Authors:
Thorsten Jörgens,
Thorsten Theobald,
Timo de Wolff
Abstract:
We introduce the imaginary projection of a multivariate polynomial $f \in \mathbb{C}[\mathbf{z}]$ as the projection of the variety of $f$ onto its imaginary part, $\mathcal{I}(f) \ = \ \{\text{Im}(\mathbf{z}) \, : \, \mathbf{z} \in \mathcal{V}(f) \}$. Since a polynomial $f$ is stable if and only if $\mathcal{I}(f) \cap \mathbb{R}_{>0}^n \ = \ \emptyset$, the notion offers a novel geometric view un…
▽ More
We introduce the imaginary projection of a multivariate polynomial $f \in \mathbb{C}[\mathbf{z}]$ as the projection of the variety of $f$ onto its imaginary part, $\mathcal{I}(f) \ = \ \{\text{Im}(\mathbf{z}) \, : \, \mathbf{z} \in \mathcal{V}(f) \}$. Since a polynomial $f$ is stable if and only if $\mathcal{I}(f) \cap \mathbb{R}_{>0}^n \ = \ \emptyset$, the notion offers a novel geometric view underlying stability questions of polynomials.
We show that the connected components of the complement of the closure of the imaginary projections are convex, thus opening a central connection to the theory of amoebas and coamoebas. Building upon this, the paper establishes structural properties of the components of the complement, such as lower bounds on their maximal number, proves a complete classification of the imaginary projections of quadratic polynomials and characterizes the limit directions for polynomials of arbitrary degree.
△ Less
Submitted 31 May, 2018; v1 submitted 5 February, 2016;
originally announced February 2016.
-
Intersections of Amoebas
Authors:
Martina Juhnke-Kubitzke,
Timo de Wolff
Abstract:
Amoebas are projections of complex algebraic varieties in the algebraic torus under a Log-absolute value map, which have connections to various mathematical subjects. While amoebas of hypersurfaces have been intensively studied in recent years, the non-hypersurface case is barely understood so far.
We investigate intersections of amoebas of $n$ hypersurfaces in $(\mathbb{C}^*)^n$, which are cano…
▽ More
Amoebas are projections of complex algebraic varieties in the algebraic torus under a Log-absolute value map, which have connections to various mathematical subjects. While amoebas of hypersurfaces have been intensively studied in recent years, the non-hypersurface case is barely understood so far.
We investigate intersections of amoebas of $n$ hypersurfaces in $(\mathbb{C}^*)^n$, which are canonical supersets of amoebas given by non-hypersurface varieties. Our main results are amoeba analogs of Bernstein's Theorem and Bézout's Theorem providing an upper bound for the number of connected components of such intersections. Moreover, we show that the \emph{order map} for hypersurface amoebas can be generalized in a natural way to intersections of amoebas. In particular, analogous to the case of amoebas of hypersurfaces, the restriction of this generalized order map to a single connected component is still $1$-to-$1$.
△ Less
Submitted 4 February, 2017; v1 submitted 28 October, 2015;
originally announced October 2015.
-
A Sharp Upper Bound for the Complexity of Labeled Oriented Trees
Authors:
Moritz Christmann,
Timo de Wolff
Abstract:
A labeled oriented graph (LOG) is an oriented graph with a labeling function from the edge set into the vertex set. The complexity of a LOG is the minimal cardinality of an initial set $S$ of vertices such that every vertex can be reached successively from $S$ only using edges with labels in $S$ or already visited vertices. We give a constructive proof of a conjecture by Rosebrock stating that for…
▽ More
A labeled oriented graph (LOG) is an oriented graph with a labeling function from the edge set into the vertex set. The complexity of a LOG is the minimal cardinality of an initial set $S$ of vertices such that every vertex can be reached successively from $S$ only using edges with labels in $S$ or already visited vertices. We give a constructive proof of a conjecture by Rosebrock stating that for an interior reduced, connected LOG with $m$ vertices the complexity is at most $(m+1) / 2$ and show that this bound is sharp.
Due to results of Howie labeled oriented trees (LOTs) yield crucial candidates for counterexamples of the Whitehead Conjecture stating that every subcomplex of an aspherical 2-complex is aspherical. We explicitly describe the structure of LOTs of maximal complexity $(m+1)/2$. We conclude that the 2-complexes associated to these LOTs are always aspherical excluding them from the list of possible counterexamples.
△ Less
Submitted 23 December, 2014;
originally announced December 2014.
-
Norms of Roots of Trinomials
Authors:
Thorsten Theobald,
Timo de Wolff
Abstract:
The behavior of norms of roots of univariate trinomials $z^{s+t} + p z^t + q \in \mathbb{C}[z]$ for fixed support $A = \{0,t,s+t\} \subset \mathbb{N}$ with respect to the choice of coefficients $p,q \in \mathbb{C}$ is a classical late 19th and early 20th century problem. Although algebraically characterized by P.\ Bohl in 1908, the geometry and topology of the corresponding parameter space of coef…
▽ More
The behavior of norms of roots of univariate trinomials $z^{s+t} + p z^t + q \in \mathbb{C}[z]$ for fixed support $A = \{0,t,s+t\} \subset \mathbb{N}$ with respect to the choice of coefficients $p,q \in \mathbb{C}$ is a classical late 19th and early 20th century problem. Although algebraically characterized by P.\ Bohl in 1908, the geometry and topology of the corresponding parameter space of coefficients had yet to be revealed. Assuming $s$ and $t$ to be coprime we provide such a characterization for the space of trinomials by reinterpreting the problem in terms of amoeba theory. The roots of given norm are parameterized in terms of a hypotrochoid curve along a $\mathbb{C}$-slice of the space of trinomials, with multiple roots of this norm appearing exactly on the singularities. As a main result, we show that the set of all trinomials with support $A$ and certain roots of identical norm, as well as its complement can be deformation retracted to the torus knot $K(s+t,s)$, and thus are connected but not simply connected. An exception is the case where the $t$-th smallest norm coincides with the $(t+1)$-st smallest norm. Here, the complement has a different topology since it has fundamental group $\mathbb{Z}^2$.
△ Less
Submitted 24 October, 2015; v1 submitted 24 November, 2014;
originally announced November 2014.
-
Lower Bounds for Polynomials with Simplex Newton Polytopes Based on Geometric Programming
Authors:
Sadik Iliman,
Timo de Wolff
Abstract:
In this article, we propose a geometric programming method in order to compute lower bounds for real polynomials. We provide new sufficient conditions for polynomials to be nonnegative as well as to have a sum of binomial squares representation. These criteria rely on the coefficients and the support of a polynomial and generalize all previous ones by Lasserre, Ghasemi, Marshall, Fidalgo and Kovac…
▽ More
In this article, we propose a geometric programming method in order to compute lower bounds for real polynomials. We provide new sufficient conditions for polynomials to be nonnegative as well as to have a sum of binomial squares representation. These criteria rely on the coefficients and the support of a polynomial and generalize all previous ones by Lasserre, Ghasemi, Marshall, Fidalgo and Kovacec to polynomials with arbitrary simplex Newton polytopes.
This generalization yields a geometric programming approach for computing lower bounds for polynomials that significantly extends the geometric programming method proposed by Ghasemi and Marshall. Furthermore, it shows that geometric programming is strongly related to nonnegativity certificates based on sums of nonnegative circuit polynomials, which were recently introduced by the authors.
△ Less
Submitted 24 February, 2016; v1 submitted 25 February, 2014;
originally announced February 2014.
-
Amoebas, Nonnegative Polynomials and Sums of Squares Supported on Circuits
Authors:
Sadik Iliman,
Timo de Wolff
Abstract:
We completely characterize sections of the cones of nonnegative polynomials, convex polynomials and sums of squares with polynomials supported on circuits, a genuine class of sparse polynomials. In particular, nonnegativity is characterized by an invariant, which can be immediately derived from the initial polynomial. Furthermore, nonnegativity of such polynomials $f$ coincides with solidness of t…
▽ More
We completely characterize sections of the cones of nonnegative polynomials, convex polynomials and sums of squares with polynomials supported on circuits, a genuine class of sparse polynomials. In particular, nonnegativity is characterized by an invariant, which can be immediately derived from the initial polynomial. Furthermore, nonnegativity of such polynomials $f$ coincides with solidness of the amoeba of $f$, i.e., the Log-absolute-value image of the algebraic variety $\mathcal{V}(f) \subset (\mathbb{C}^*)^n$ of $f$.
These results generalize earlier works both in amoeba theory and real algebraic geometry by Fidalgo, Kovacec, Reznick, Theobald and de Wolff and solve an open problem by Reznick. They establish the first direct connection between amoeba theory and nonnegativity of real polynomials. Additionally, these statements yield a completely new class of nonnegativity certificates independent from sums of squares certificates.
△ Less
Submitted 24 October, 2015; v1 submitted 3 February, 2014;
originally announced February 2014.
-
The Boundary of Amoebas
Authors:
Franziska Schroeter,
Timo de Wolff
Abstract:
The computation of amoebas has been a challenging open problem for the last dozen years. The most natural approach, namely to compute an amoeba via its boundary, has not been practical so far since only a superset of the boundary, the contour, is understood in theory and computable in practice.
We define and characterize the extended boundary of an amoeba, which is sensitive to some degeneration…
▽ More
The computation of amoebas has been a challenging open problem for the last dozen years. The most natural approach, namely to compute an amoeba via its boundary, has not been practical so far since only a superset of the boundary, the contour, is understood in theory and computable in practice.
We define and characterize the extended boundary of an amoeba, which is sensitive to some degenerations that the topological boundary does not detect. Our description of the extended boundary also allows us to distinguish between the contour and the boundary. This gives rise not only to new structural results in amoeba theory, but in particular allows us to compute hypersurface amoebas via their boundary in any dimension. In dimension two this can be done using Gröbner bases alone.
We introduce the concept of amoeba bases, which are sufficient for understanding the amoeba of an ideal. We show that our characterization of the boundary is essential for the computation of these amoeba bases and we illustrate the potential of this concept by constructing amoeba bases for linear systems of equations.
△ Less
Submitted 11 June, 2016; v1 submitted 28 October, 2013;
originally announced October 2013.
-
Low Dimensional Test Sets for Nonnegativity of Even Symmetric Forms
Authors:
Sadik Iliman,
Timo de Wolff
Abstract:
An important theorem by Timofte states that nonnegativity of real $n$-variate symmetric polynomials of degree $d$ can be decided at test sets given by all points with at most $\lfloor\frac{d}{2}\rfloor$ distinct components. However, if the degree is sufficiently larger than the number of variables, then the theorem obviously does not provide nontrivial information. Our approach is to look at…
▽ More
An important theorem by Timofte states that nonnegativity of real $n$-variate symmetric polynomials of degree $d$ can be decided at test sets given by all points with at most $\lfloor\frac{d}{2}\rfloor$ distinct components. However, if the degree is sufficiently larger than the number of variables, then the theorem obviously does not provide nontrivial information. Our approach is to look at $(m + 1)$-dimensional subspaces of even symmetric forms of degree 4d, at which nonnegativity can be checked at $(m - 1)$-points, i.e., points with at most $m - 1 \in \N$ distinct components, where $m$ is independent of the degree of the forms and better than Timofte's bound. Furthermore, for fixed $k \in \N$, we tackle problems concerning the maximum dimension of such subspaces, at which nonnegativity can be checked at all $k$-points, as well as the geometrical and topological structure of the set of all forms whose nonnegativity can be decided at all $k$-points.
△ Less
Submitted 18 March, 2013;
originally announced March 2013.
-
Separating inequalities for nonnegative polynomials that are not sums of squares
Authors:
Sadik Iliman,
Timo de Wolff
Abstract:
Ternary sextics and quaternary quartics are the smallest cases where there exist nonnegative polynomials that are not sums of squares (SOS). A complete classification of the difference between these cones was given by G. Blekherman via analyzing the extreme rays of the corresponding dual cones. However, an exact computational approach in order to build separating extreme rays for nonnegative polyn…
▽ More
Ternary sextics and quaternary quartics are the smallest cases where there exist nonnegative polynomials that are not sums of squares (SOS). A complete classification of the difference between these cones was given by G. Blekherman via analyzing the extreme rays of the corresponding dual cones. However, an exact computational approach in order to build separating extreme rays for nonnegative polynomials that are not sums of squares is a widely open problem. We provide a method substantially simplifying this computation for certain classes of polynomials on the boundary of the PSD cones. In particular, our method yields separating extreme rays for every nonnegative ternary sextic with at least seven zeros. As an application to further instances, we compute a rational certificate proving that the Motzkin polynomial is not SOS.
△ Less
Submitted 1 August, 2012; v1 submitted 19 January, 2012;
originally announced January 2012.
-
Amoebas of genus at most one
Authors:
Thorsten Theobald,
Timo de Wolff
Abstract:
The amoeba of a Laurent polynomial $f \in \C[z_1^{\pm 1},\ldots,z_n^{\pm 1}]$ is the image of its zero set $\mathcal{V}(f)$ under the log-absolute-value map. Understanding the space of amoebas (i.e., the decomposition of the space of all polynomials, say, with given support or Newton polytope, with regard to the existing complement components) is a widely open problem.
In this paper we investiga…
▽ More
The amoeba of a Laurent polynomial $f \in \C[z_1^{\pm 1},\ldots,z_n^{\pm 1}]$ is the image of its zero set $\mathcal{V}(f)$ under the log-absolute-value map. Understanding the space of amoebas (i.e., the decomposition of the space of all polynomials, say, with given support or Newton polytope, with regard to the existing complement components) is a widely open problem.
In this paper we investigate the class of polynomials $f$ whose Newton polytope $\New(f)$ is a simplex and whose support $A$ contains exactly one point in the interior of $\New(f)$. Amoebas of polynomials in this class may have at most one bounded complement component. We provide various results on the space of these amoebas. In particular, we give upper and lower bounds in terms of the coefficients of $f$ for the existence of this complement component and show that the upper bound becomes sharp under some extremal condition. We establish connections from our bounds to Purbhoo's lopsidedness criterion and to the theory of $A$-discriminants.
Finally, we provide a complete classification of the space of amoebas for the case that the exponent of the inner monomial is the barycenter of the simplex Newton polytope. In particular, we show that the set of all polynomials with amoebas of genus 1 is path-connected in the corresponding space of amoebas, which proves a special case of the question on connectivity (for general Newton polytopes) stated by H. Rullgård.
△ Less
Submitted 23 April, 2013; v1 submitted 11 August, 2011;
originally announced August 2011.
-
Approximating amoebas and coamoebas by sums of squares
Authors:
Thorsten Theobald,
Timo de Wolff
Abstract:
Amoebas and coamoebas are the logarithmic images of algebraic varieties and the images of algebraic varieties under the arg-map, respectively. We present new techniques for computational problems on amoebas and coamoebas, thus establishing new connections between (co-)amoebas, semialgebraic and convex algebraic geometry and semidefinite programming.
Our approach is based on formulating the membe…
▽ More
Amoebas and coamoebas are the logarithmic images of algebraic varieties and the images of algebraic varieties under the arg-map, respectively. We present new techniques for computational problems on amoebas and coamoebas, thus establishing new connections between (co-)amoebas, semialgebraic and convex algebraic geometry and semidefinite programming.
Our approach is based on formulating the membership problem in amoebas (respectively coamoebas) as a suitable real algebraic feasibility problem. Using the real Nullstellensatz, this allows to tackle the problem by sums of squares techniques and semidefinite programming. Our method yields polynomial identities as certificates of non-containment of a point in an amoeba or coamoeba. As the main theoretical result, we establish some degree bounds on the polynomial certificates. Moreover, we provide some actual computations of amoebas based on the sums of squares approach.
△ Less
Submitted 13 June, 2013; v1 submitted 21 January, 2011;
originally announced January 2011.
-
Polytopes with Special Simplices
Authors:
Timo de Wolff
Abstract:
For a polytope P a simplex S with vertex set V(S) is called a special simplex if every facet of P contains all but exactly one vertex of S. For such polytopes P with face complex F(P) containing a special simplex the subcomplex F(P) / V(S) of all faces not containing vertices of S is the boundary of a polytope Q - the basis polytope of P. If additionally the dimension of the affine basis space of…
▽ More
For a polytope P a simplex S with vertex set V(S) is called a special simplex if every facet of P contains all but exactly one vertex of S. For such polytopes P with face complex F(P) containing a special simplex the subcomplex F(P) / V(S) of all faces not containing vertices of S is the boundary of a polytope Q - the basis polytope of P. If additionally the dimension of the affine basis space of F(P) / V(S) equals dim(Q), we call P meek; otherwise we call P wild. We give a full combinatorial classification and techniques for geometric construction of the class of meek polytopes with special simplices. We show that every wild polytope P' with special simplex can be constructed out of a particular meek one P by intersecting P with particular hyperplanes. It is non-trivial to find all these hyperplanes for an arbitrary basis polytope; we give an exact description for 2-basis polytopes. Furthermore we show that the f-vector of each wild polytope with special simplex is component wise bounded above by the f-vector of a particular meek one which can be computed explicitly. Finally, we discuss the n-cube as a non-trivial example of a wild polytope with special simplex and prove that its basis polytope is the zonotope given by the Minkowski sum of the (n-1)-cube and the vector (1,...,1). Polytopes with special simplex have applications on Ehrhart theory, toric rings and were just used by Francisco Santos to construct a counter-example disproving the Hirsch conjecture.
△ Less
Submitted 30 September, 2010;
originally announced September 2010.