View: session overviewtalk overview
Abstract. Subtropical real root finding is a heuristic method designed to solve large multivariate polynomials. This method considers the polynomial as a collection of exponent vectors associated with sign information on the coefficients. It employs linear programming to heuristically identify points where the values of the polynomial have opposite signs. In successful cases, the intermediate value theorem is utilized to compute a zero of the polynomial represented as a tuple of real algebraic numbers. By implementing this approach using the Reduce computer algebra system and theGurobi LP solver, we have achieved success in solving hundreds of problems derived from established mathematical models in chemistry and systems biology. These problems involved polynomials with up to 800,000 monomials in 10 variables and degrees up to 12, with an incompleteness failure rate of approximately 10 percent.
The method has been adapted to satisfiability modulo theories solving (SMT), where the focus went from one single polynomial equation to several simultaneous strict polynomial inequalities. We found that such problems represent more than 40 percent of the QF_NRA section of the SMT-LIB library of benchmarks. Technically, input problems in the theory QF_NRA are reduced problems in QF_LRA. Then SMT itself can be used in place of LP solving, which admits linear problems that are not necessarily pure conjunctions of constraints. Systematic experiments on the SMT-LIB have shown that the method is not a sufficiently strong decision procedure by itself but a valuable heuristic to use within a portfolio of techniques.
The talk gives an overview of the method and its application, along with theoretical results regarding its incompleteness.
Abstract. Rewriting is a formalism widely used in computer science and mathematical logic. The classical formalism has been extended, in the context of functional languages, with an order over the rules and, in the context of rewrite based languages, with the negation over patterns. We have proposed a concise and clear algorithm computing the difference over patterns which can be used to define generic encodings of constructor term rewriting systems with negation and order into classical term rewriting systems. As a direct consequence, established methods used for term rewriting systems can be applied to analyze the properties of the extended systems. The approach can also be seen as a generic compiler which targets any language providing basic pattern matching primitives. The formalism provides also a new method for deciding if a set of patterns subsumes a given pattern and thus, for checking the completeness of a set of patterns, the presence of useless patterns, or the absence of some patterns. The latter is particularly useful when one wants to verify that some patterns are absent from the result of a particular transformation and several extensions have been proposed to statically verify the absence of specified patterns from the reachable terms of constructor based term rewriting systems.
11:40 | Uniform SAmplINg with BOLTZmann PRESENTER: Martin Pépin ABSTRACT. Usain Boltz is a fast Python library for the uniform random generation of tree-like structures. It allows the user to specify both (1) the data structure they wish to sample, using simple combinators similar to those of context-free grammars, and (2) their memory representation. The underlying algorithms are optimised Boltzmann samplers allowing to get approximate-size uniform random generation in linear time. Experimental results show that Usain Boltz matches the performance of the experimental Arbogen package for OCaml, and out-performs the Boltzmann brain Haskell library, while being easier to integrate into existing scientific tools such as Sagemath. |
12:00 | Towards Learning Infinite SMT Models PRESENTER: Mikolas Janota ABSTRACT. This short paper proposes to learn models of Satisfiability Modulo Theories (SMT) formulas during solving. Specifically, we focus on infinite models for problems in the logic of linear arithmetic with uninterpreted functions (UFLIA). The constructed models are piecewise linear. Such models are useful for satisfiable problems but also provide an alternative driver for model based quantifier instantiation (MBQI). |
12:20 | On the Variants of Subset Sum: Projected and Unbounded PRESENTER: Mahesh Rajasree ABSTRACT. In this work, we focus on variants of Subset Sum. We first define a new variant called the Unique Projection Subset Sum (u-PSSUM) problem: Given (a_1, \hdots, a_n) \in Z_{\ge 0}^n, such that for every t, \sum_{i \in [n]} c_i a_i=t, for c_i \in \{0,1\} has unique solution if exists, u-PSSUM asks for a vector \Vec{x} = (x_1, \dots, x_n), such that there exists an S \subseteq [n], where \sum_{i \in S} x_i = t, and ||\Vec{a} - \Vec{x}||_p is minimized, where ||\cdot||_p denotes the \ell_p norm. We present a deterministic O(nt)-time algorithm and a randomized \widetilde{O}(n+t)-time algorithm for PSSUM, in \ell_1-norm. The second one is already a known variant called Unbounded Subset Sum (UBSSUM) problem, which takes a tuple of non-negative integers (a_1, \dots, a_n, t), as an input, and asks whether there exists non-negative integers \beta_1, \dots, \beta_n, such that \sum_{i=1}^n \beta_i a_i = t. We present the first polynomial time reductions from UBSSUM to Closest Vector Problem (CVP) in \ell_1 and \ell_{\infty} norms. We also give new algorithms for two variants of UBSSUM problem. |
Abstract. TLA+ is a language for the formal description of (mainly) concurrent and distributed algorithms and systems, based on Zermelo-Fraenkel set theory and linear-time temporal logic. Verifying properties of TLA+ specifications is supported by model checking (TLC and Apalache) and by theorem proving (TLAPS). This tutorial will provide an introduction to the language and its accompanying tools, using the problem of distributed termination detection as a running example.
15:50 | Hybrid Intervals and Symbolic Block Matrices ABSTRACT. Structured matrices with symbolic sizes appear frequently in the literature, especially in the description of algorithms for linear algebra. Recent work has treated these symbolic structured matrices themselves as computational objects, showing how to add matrices with blocks of different symbolic sizes in a general way while avoiding a combinatorial explosion of cases. The present article introduces the concept of hybrid intervals, which may have negative multiplicity. Various operations on hybrid intervals have compact and elegant formulations that do not require cases to handle different orders of the end points. This makes them useful to represent symbolic block matrix structures and to express arithmetic on symbolic block matrices compactly. The present article shows how hybrid intervals can be used to formulate symbolic block matrix addition and multiplication in a compact and uniform way. |
16:10 | Symbolic and numeric computation of symmetries for a class of Schrodinger Equations ABSTRACT. An important and challenging computational problem is to identify and include the missing compatibility (integrability) conditions for general systems of partial differential equations. The inclusion of such missing conditions is executed by the application of differential elimination algorithms. Differential equations arising during modeling generally contain both exactly known coefficients and coefficients known approximately from data. We focus on our recent work on approximate differential elimination methods and in particular their application to the determination of approximate symmetries. We illustrate this with applications to a class of Schrodinger equations. |
16:30 | Experimental comparison of real root isolation implementations PRESENTER: Zafeirakis Zafeirakopoulos ABSTRACT. Real solving of univariate polynomials is a fundamental algebraic question with important applications. This short paper focuses on the comparison of state-of-the-art free software implementations for isolating real roots of univariate polynomials. The experimental comparison is done on an extensive dataset of polynomials of different types of roots, as well as varying coefficient bitsize and degree. |
16:50 | On systems of linear homogeneous algebraic equations over truncated formal series PRESENTER: Sergei Abramov ABSTRACT. Let $K$ be a field, $K[[x]]$ the ring of formal power series, and $K((x))$ the field of formal Laurent series over $K$. Let $A$ be a nonzero polynomial matrix of degree $l$ which is an $l$-truncation of a matrix $M\in K[[x]]^{n\times n}$ which is not known in full. We establish here that it is impossible for the system $My=0$ to have a nonzero solution for every prolongation $M$ of a singular matrix $A$, and that for any truncated matrix $A$ there exists a prolongation giving a nonsingular matrix belonging to $K[[x]]^{n\times n}$. An algorithm is proposed which, given $A$, checks if the system $My=0$ has no non-zero solution for any prolongation $M$, and if this system, depending on $M$, may or may not have a nonzero solution in $K((x))^n$. |
17:10 | New Techniques in Numerical Analysis for Artificial Intelligence ABSTRACT. Artificial intelligence (AI) has revolutionized numerous industries by providing intelligent solutions to complex problems. Numerical analysis plays a pivotal role in enhancing the performance and reliability of AI algorithms. This research paper explores new implementation techniques in numerical analysis for artificial intelligence. We investigate various areas where numerical analysis methods can be leveraged to improve AI algorithms, such as optimization, machine learning, and deep learning. Additionally, we discuss advancements in numerical techniques, including numerical optimization algorithms, numerical stability, and uncertainty quantification. The results demonstrate the effectiveness of these techniques in enhancing the performance and robustness of AI systems. |