Calcium: Computing in Exact Real and Complex Fields Fredrik Johansson

Total Page:16

File Type:pdf, Size:1020Kb

Calcium: Computing in Exact Real and Complex Fields Fredrik Johansson Calcium: computing in exact real and complex fields Fredrik Johansson To cite this version: Fredrik Johansson. Calcium: computing in exact real and complex fields. ISSAC ’21, Jul 2021, Virtual Event, Russia. 10.1145/3452143.3465513. hal-02986375v2 HAL Id: hal-02986375 https://hal.inria.fr/hal-02986375v2 Submitted on 15 May 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Calcium: computing in exact real and complex fields Fredrik Johansson [email protected] Inria Bordeaux and Institut Math. Bordeaux 33400 Talence, France ABSTRACT This paper presents Calcium,1 a C library for exact computa- Calcium is a C library for real and complex numbers in a form tion in R and C. Numbers are represented as elements of fields suitable for exact algebraic and symbolic computation. Numbers Q¹a1;:::; anº where the extension numbers ak are defined symbol- ically. The system constructs fields and discovers algebraic relations are represented as elements of fields Q¹a1;:::; anº where the exten- automatically, handling algebraic and transcendental number fields sion numbers ak may be algebraic or transcendental. The system combines efficient field operations with automatic discovery and in a unified way. It is capable of deciding equality for a wide class certification of algebraic relations, resulting in a practical computa- of numbers which includes Q as a subset. We show a few basic tional model of R and C in which equality is rigorously decidable examples, here using a Python wrapper: for a large class of numbers. >>> (pi**2 - 9) / (pi + 3) 0.141593 {a-3 where a = 3.14159 [Pi]} >>> phi = (sqrt(5)+1)/2; (phi**100 - (1-phi)**100)/sqrt(5) CCS CONCEPTS 3.54225e+20 {354224848179261915075} • Computing methodologies ! Representation of exact num- >>> i**i - exp(pi / (sqrt(-2)**sqrt(2))**sqrt(2)) bers; Hybrid symbolic-numeric methods; Exact arithmetic 0 algorithms; Symbolic and algebraic manipulation. >>> log(sqrt(2)+sqrt(3)) / log(5+2*sqrt(6)) 0.500000 {1/2} >>> erf(4*atan(ca(1)/5) - atan(ca(1)/239)) + erfc(pi/4) KEYWORDS 1 Exact real arithmetic, Number fields, Transcendental numbers >>> -1e-12 < exp(pi*sqrt(163)) - 262537412640768744 < -1e-13 True ACM Reference Format: In the first example, the field is Q¹aº where a = π, and the Fredrik Johansson. 2021. Calcium: computing in exact real and complex element is a − 3. The numerical approximation (x ≈ 0:141593) fields. In Proceedings of the 2021 International Symposium on Symbolic and is computed to desired precision on demand, for example when Algebraic Computation (ISSAC ’21), July 18–23, 2021, Virtual Event, Russian printing or evaluating a numerical predicate. Examples 2–5 were ACM, New York, NY, USA,9 pages. https://doi.org/10.1145/ Federation. chosen so that the field of the result simplifies to Q. 3452143.3465513 The basic idea in this work is to combine a multivariate fraction field representation, automatic relation discovery, ideal reduction 1 INTRODUCTION (where the ideal may be incomplete), and rigorous numerical certi- fication using ball arithmetic. This approach is inspired by earlier A field K is said to be effective if its elements can be enumerated implementations of Q (including Sage [30] and Magma [3]) and f ; −; ·; /; g and the operations + = are computable. Examples include by theoretical work on transcendental fields (notably Richardson’s the rationals Q, finite fields Fq , and the algebraic numbers Q. algorithm [26]). The main improvement over previous software The fields of real and complex numbers R and C are notably non- is that we handle both algebraic and transcendental numbers in a effective, even when restricted to so-called computable numbers (a rigorous and structured way; compared to symbolic computation real number x is said to be computable if there is a program which, systems like Mathematica and Maple, our approach provides faster −n given n, outputs some xn 2 Q with jx − xn j < 2 ). The problem is field arithmetic and more reliable results. that equality is only semi-decidable: we can in general prove x , y, This paper is structured as follows. Section2 presents our high- but not x = y, as a consequence of the halting theorem. Neverthe- level strategy for exact computation, described at a general level less, as the example of Q shows, we can hope for an equality test at without reference to low-level implementation details. Section3 least for some numbers within a suitable algebraic framework. discusses the architecture of Calcium. Section4 relates our strategy to earlier work and presents some benchmark results. Permission to make digital or hard copies of all or part of this work for personal or 2 COMPUTING IN SUBFIELDS OF C classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation To simulate R and C, we may start with Q and lazily extend the on the first page. Copyrights for components of this work owned by others than the field with new numbers a as they arise in computations. A general author(s) must be honored. Abstracting with credit is permitted. To copy otherwise, or k republish, to post on servers or to redistribute to lists, requires prior specific permission way to compute in such extension fields of Q is in terms of quotient and/or a fee. Request permissions from [email protected]. rings and their fields of fractions (henceforth formal fields). ISSAC ’21, July 18–23, 2021, Virtual Event, Russian Federation © 2021 Copyright held by the owner/author(s). Publication rights licensed to ACM. 1Pronounced “kalkium” to distinguish it from the chemical element. Calcium is free and ACM ISBN 978-1-4503-8382-0/21/07...$15.00 open source (LGPL 2.1+) software. The source repository is https://github.com/fredrik- https://doi.org/10.1145/3452143.3465513 johansson/calcium and the documentation is available at http://fredrikj.net/calcium/. In the following, we assume that a1;:::; an is a finite list of • Evaluate the complex conjugation map z ! z. complex numbers. We let X1;:::; Xn denote independent formal • Evaluate numerical ordering relations (x < y, jxj < jyj, etc.). variables, we let µ : Q»X1;:::; Xn¼ ! C denote the evaluation • Exclude singularities (e.g. division by zero) and choose well- homomorphism induced by the map Xk 7! ak , and we define defined branches of multivalued functions. I := ker µ = ff 2 Q»X1;:::; Xn¼ : f ¹a1;:::; anº = 0g We therefore need a symbolic way to define extension numbers a ;:::; a so that they are explicitly (numerically) computable, and as the ideal of all algebraic relations among a ;:::; a over Q. 1 n 1 n we need to construct I from this symbolic data rather than vice Theorem 1. Assume that I is known (in the sense that an explicit versa. The following types of extensions are useful: list of generators I = hf1;:::; fmi is known). Then • Absolute algebraicp : a is a fixed algebraic constant a 2 Q, 2π i/3 5 K := Q¹a1;:::; anº Kformal := Frac¹Q»X1;:::; Xn]/Iº for example i, 2, e , or »a − a − 1 = 0; a ≈ 1:17¼. is an effective subfield of C. Such a constant can be defined canonically by its minimal polynomial over Q together with an isolating ball for a root. Proof. The isomorphism is obvious. Decidability of “=” in the • Relative algebraic: a is defined by an equation am/n = c with formal field follows from the fact that we can compute aGröb- c 2 C, or more generally P¹aº = 0 with P 2 C»X¼, together p r ner basis for I. Given formal fractions q and s with p;q;r;s 2 with an isolating complex ball for a root. Q»X1;:::; Xn¼, we can consequently decide whether ps ≡ qr mod I. • Transcendental (perhaps conjecturally): a is a symbolic con- z w Indeed, we can also decide whether q;s . 0 mod I and thereby en- stant (π, γ , etc.) or function (e , log¹zº, z , Γ¹zº, Jν ¹zº, etc.) sure that the fractions define numbers in the first place. □ evaluated at some point. • Black-box computable: a is defined by a program for numer- Some easy special cases are worth noting: ical evaluation in ball arithmetic. (We will not be able to • The trivial field K = Q (take n = 0). prove algebraic relations except self-relations like a − a = 0.) • = Q¹ ;:::; º Transcendental number fields K a1 an where the Calcium represents algebraic and transcendental extensions in numbers a1;:::; an are algebraically independent over Q. the usual way as symbolic expressions f ¹z1;:::; zp º. The argu- • Algebraic number fields K = Q¹aº Q»X]/hf ¹Xºi where ments z1;:::; zp are real or complex numbers which may belong a is an algebraic number with minimal polynomial f . to different fields, say z1 2 Q¹b1;:::;br º, z2 2 Q¹c1;:::;cs º, etc. The general case is a mixed field in which the extension numbers Each extension number defines a computable number through may be algebraic or transcendental and algebraically dependent or recursive numerical evaluation of the symbolic function and its independent in any combination. arguments in arbitrary-precision ball arithmetic. This is at least true in principle assuming that we can decide signs at discontinuities.
Recommended publications
  • Variable Planck's Constant
    Variable Planck’s Constant: Treated As A Dynamical Field And Path Integral Rand Dannenberg Ventura College, Physics and Astronomy Department, Ventura CA [email protected] [email protected] January 28, 2021 Abstract. The constant ħ is elevated to a dynamical field, coupling to other fields, and itself, through the Lagrangian density derivative terms. The spatial and temporal dependence of ħ falls directly out of the field equations themselves. Three solutions are found: a free field with a tadpole term; a standing-wave non-propagating mode; a non-oscillating non-propagating mode. The first two could be quantized. The third corresponds to a zero-momentum classical field that naturally decays spatially to a constant with no ad-hoc terms added to the Lagrangian. An attempt is made to calibrate the constants in the third solution based on experimental data. The three fields are referred to as actons. It is tentatively concluded that the acton origin coincides with a massive body, or point of infinite density, though is not mass dependent. An expression for the positional dependence of Planck’s constant is derived from a field theory in this work that matches in functional form that of one derived from considerations of Local Position Invariance violation in GR in another paper by this author. Astrophysical and Cosmological interpretations are provided. A derivation is shown for how the integrand in the path integral exponent becomes Lc/ħ(r), where Lc is the classical action. The path that makes stationary the integral in the exponent is termed the “dominant” path, and deviates from the classical path systematically due to the position dependence of ħ.
    [Show full text]
  • Computability on the Real Numbers
    4. Computability on the Real Numbers Real numbers are the basic objects in analysis. For most non-mathematicians a real number is an infinite decimal fraction, for example π = 3•14159 ... Mathematicians prefer to define the real numbers axiomatically as follows: · (R, +, , 0, 1,<) is, up to isomorphism, the only Archimedean ordered field satisfying the axiom of continuity [Die60]. The set of real numbers can also be constructed in various ways, for example by means of Dedekind cuts or by completion of the (metric space of) rational numbers. We will neglect all foundational problems and assume that the real numbers form a well-defined set R with all the properties which are proved in analysis. We will denote the R real line topology, that is, the set of all open subsets of ,byτR. In Sect. 4.1 we introduce several representations of the real numbers, three of which (and the equivalent ones) will survive as useful. We introduce a rep- n n resentation ρ of R by generalizing the definition of the main representation ρ of the set R of real numbers. In Sect. 4.2 we discuss the computable real numbers. Sect. 4.3 is devoted to computable real functions. We show that many well known functions are computable, and we show that partial sum- mation of sequences is computable and that limit operator on sequences of real numbers is computable, if a modulus of convergence is given. We also prove a computability theorem for power series. Convention 4.0.1. We still assume that Σ is a fixed finite alphabet con- taining all the symbols we will need.
    [Show full text]
  • Semiclassical Unimodular Gravity
    Preprint typeset in JHEP style - PAPER VERSION Semiclassical Unimodular Gravity Bartomeu Fiol and Jaume Garriga Departament de F´ısica Fonamental i Institut de Ci`encies del Cosmos, Universitat de Barcelona, Mart´ıi Franqu`es 1, 08028 Barcelona, Spain [email protected], [email protected] Abstract: Classically, unimodular gravity is known to be equivalent to General Rel- ativity (GR), except for the fact that the effective cosmological constant Λ has the status of an integration constant. Here, we explore various formulations of unimod- ular gravity beyond the classical limit. We first consider the non-generally covariant action formulation in which the determinant of the metric is held fixed to unity. We argue that the corresponding quantum theory is also equivalent to General Relativity for localized perturbative processes which take place in generic backgrounds of infinite volume (such as asymptotically flat spacetimes). Next, using the same action, we cal- culate semiclassical non-perturbative quantities, which we expect will be dominated by Euclidean instanton solutions. We derive the entropy/area ratio for cosmological and black hole horizons, finding agreement with GR for solutions in backgrounds of infinite volume, but disagreement for backgrounds with finite volume. In deriving the arXiv:0809.1371v3 [hep-th] 29 Jul 2010 above results, the path integral is taken over histories with fixed 4-volume. We point out that the results are different if we allow the 4-volume of the different histories to vary over a continuum range. In this ”generalized” version of unimodular gravity, one recovers the full set of Einstein’s equations in the classical limit, including the trace, so Λ is no longer an integration constant.
    [Show full text]
  • Computability and Analysis: the Legacy of Alan Turing
    Computability and analysis: the legacy of Alan Turing Jeremy Avigad Departments of Philosophy and Mathematical Sciences Carnegie Mellon University, Pittsburgh, USA [email protected] Vasco Brattka Faculty of Computer Science Universit¨at der Bundeswehr M¨unchen, Germany and Department of Mathematics and Applied Mathematics University of Cape Town, South Africa and Isaac Newton Institute for Mathematical Sciences Cambridge, United Kingdom [email protected] October 31, 2012 1 Introduction For most of its history, mathematics was algorithmic in nature. The geometric claims in Euclid’s Elements fall into two distinct categories: “problems,” which assert that a construction can be carried out to meet a given specification, and “theorems,” which assert that some property holds of a particular geometric configuration. For example, Proposition 10 of Book I reads “To bisect a given straight line.” Euclid’s “proof” gives the construction, and ends with the (Greek equivalent of) Q.E.F., for quod erat faciendum, or “that which was to be done.” arXiv:1206.3431v2 [math.LO] 30 Oct 2012 Proofs of theorems, in contrast, end with Q.E.D., for quod erat demonstran- dum, or “that which was to be shown”; but even these typically involve the construction of auxiliary geometric objects in order to verify the claim. Similarly, algebra was devoted to developing algorithms for solving equa- tions. This outlook characterized the subject from its origins in ancient Egypt and Babylon, through the ninth century work of al-Khwarizmi, to the solutions to the quadratic and cubic equations in Cardano’s Ars Magna of 1545, and to Lagrange’s study of the quintic in his R´eflexions sur la r´esolution alg´ebrique des ´equations of 1770.
    [Show full text]
  • Monodromies and Functional Determinants in the CFT Driven Quantum Cosmology
    Monodromies and functional determinants in the CFT driven quantum cosmology A.O.Barvinsky and D.V.Nesterov Theory Department, Lebedev Physics Institute, Leninsky Prospect 53, Moscow 119991, Russia Abstract We discuss the calculation of the reduced functional determinant of a special second order differential operator F = −d2/dτ 2 + g¨=g, g¨ ≡ d2g/dτ 2, with a generic function g(τ), subject to periodic boundary conditions. This implies the gauge-fixed path integral representation of this determinant and the monodromy method of its calculation. Motivations for this particular problem, coming from applications in quantum cosmology, are also briefly discussed. They include the problem of microcanonical initial conditions in cosmology driven by a conformal field theory, cosmological constant and cosmic microwave background problems. 1. Introduction Here we consider the class of problems involving the differential operator of the form d2 g¨ F = − + ; : (1.1) dτ 2 g where g = g(τ) is a rather generic function of its variable τ. From calculational viewpoint, the virtue of this operator is that g(τ) represents its explicit basis function { the solution of the homogeneous equation, F g(τ) = 0; (1.2) which immediately allows one to construct its second linearly independent solution τ dy Ψ(τ) = g(τ) (1.3) g2(y) τZ∗ and explicitly build the Green's function of F with appropriate boundary conditions. On the other hand, from physical viewpoint this operator is interesting because it describes long-wavelength per- turbations in early Universe, including the formation of observable CMB spectra [1, 2], statistical ensembles in quantum cosmology [3], etc.
    [Show full text]
  • On the Computational Complexity of Algebraic Numbers: the Hartmanis–Stearns Problem Revisited
    On the computational complexity of algebraic numbers : the Hartmanis-Stearns problem revisited Boris Adamczewski, Julien Cassaigne, Marion Le Gonidec To cite this version: Boris Adamczewski, Julien Cassaigne, Marion Le Gonidec. On the computational complexity of alge- braic numbers : the Hartmanis-Stearns problem revisited. Transactions of the American Mathematical Society, American Mathematical Society, 2020, pp.3085-3115. 10.1090/tran/7915. hal-01254293 HAL Id: hal-01254293 https://hal.archives-ouvertes.fr/hal-01254293 Submitted on 12 Jan 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. ON THE COMPUTATIONAL COMPLEXITY OF ALGEBRAIC NUMBERS: THE HARTMANIS–STEARNS PROBLEM REVISITED by Boris Adamczewski, Julien Cassaigne & Marion Le Gonidec Abstract. — We consider the complexity of integer base expansions of alge- braic irrational numbers from a computational point of view. We show that the Hartmanis–Stearns problem can be solved in a satisfactory way for the class of multistack machines. In this direction, our main result is that the base-b expansion of an algebraic irrational real number cannot be generated by a deterministic pushdown automaton. We also confirm an old claim of Cobham proving that such numbers cannot be generated by a tag machine with dilation factor larger than one.
    [Show full text]
  • This Article Was Published in an Elsevier Journal. the Attached Copy
    This article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author’s institution, sharing with colleagues and providing to institution administration. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier’s archiving and manuscript policies are encouraged to visit: http://www.elsevier.com/copyright Author's personal copy Theoretical Computer Science 394 (2008) 159–174 www.elsevier.com/locate/tcs Physical constraints on hypercomputation$ Paul Cockshotta,∗, Lewis Mackenziea, Greg Michaelsonb a Department of Computing Science, University of Glasgow, 17 Lilybank Gardens, Glasgow G12 8QQ, United Kingdom b School of Mathematical and Computer Sciences, Heriot-Watt University, Riccarton, Edinburgh EH14 4AS, United Kingdom Abstract Many attempts to transcend the fundamental limitations to computability implied by the Halting Problem for Turing Machines depend on the use of forms of hypercomputation that draw on notions of infinite or continuous, as opposed to bounded or discrete, computation. Thus, such schemes may include the deployment of actualised rather than potential infinities of physical resources, or of physical representations of real numbers to arbitrary precision. Here, we argue that such bases for hypercomputation are not materially realisable and so cannot constitute new forms of effective calculability. c 2007 Elsevier B.V. All rights reserved.
    [Show full text]
  • The Cosmological Constant Problem 1 Introduction
    The Cosmological Constant Problem and (no) solutions to it by Wolfgang G. Hollik 7/11/2017 Talk given at the Workshop Seminar Series on Vacuum Energy @ DESY. “Physics thrives on crisis. [... ] Unfortunately, we have run short on crisis lately.” (Steven Weinberg, 1989 [1]) 1 Introduction: the Cosmological Constant Besides the fact, that we are still running short on severe crises in physics also about 30 years after Weinberg’s statement, there is still yet no solution to what is called the “Cosmological Constant problem”. Neither is there any clue what the Cosmological Constant (CC) is made of and if the problem is indeed an outstanding problem. Weinberg defines and proves in his review on the Cosmological Constant [1] a “no-go” theorem. This no-go theorem is actually not a theorem on the smallness of the CC in the sense of a Vacuum Energy,1 it is rather a theorem prohibiting any kind of adjustment mechanisms that lead effectively to a universe with a flat and static (i.e. Minkowski) space-time metric in the presence of a CC. In other words: with a CC there is no Minkowski universe possible. An old crisis When Einstein first formulated his field equations for General Relativity, he was neither aware of a possible Cosmological Constant nor of an expanding universe solution to them. Originally, the proposed equations are 1 R g R = 8πG T , (1) µν − 2 µν − µν using Weinberg’s ( + ++)-metric. The matter content is given by the energy-stress tensor − Tµν, where the geometry of space-time is encoded in the Ricci tensor Rµν and the scalar µν curvature R = g Rµν of the metric gµν.
    [Show full text]
  • Transcendental Numbers and Periods
    TRANSCENDENTAL NUMBERS AND PERIODS JAMES CARLSON Contents 1. Introduction 1 1.1. Diophantine approximation I: upper bounds 2 1.2. Diophantine approximation II: lower bounds 4 1.3. Proof of the lower bound 5 2. Periods 6 2.1. Periods 6 3. Computable numbers 8 4. References 11 1. Introduction These are notes for a talk given at Colorado State University on February 12, 2015. The numbers we are familiar with fall into a hierarchy of sets of ever greater scope: N ⊂ Z ⊂ Q ⊂ Q¯ ⊂ C These are the natural numbers, the integers, the rationals, the field of algebraic numbers, and complex numbers, respectively. The field Q¯ is the set of all possible solutions of polynomial equations with rational coefficients. Such equations can be enumerated, so Q¯ , like the field of rational numbers itself, is countable. Nonetheless, the Q¯ is considerably larger than Q in the sense that it infinite dimensional rational vector space. For example, Date: 02-19-2014. 1 2 JAMES CARLSON p the numbers p, where p is prime, are linearly independent. These reults have their roots in the Greek mathematics of around 500 BC. Consider now the field of complex numbers. By the work of Georg Cantor, it is anun- countable field—a remarkable fact that rests on the Cantor’s famous diagonal argument. It follows that ”most” numbers are transcendental. That is, they satisfy no polynomial equation with rational coefficients. Consequently, there is a decomposition C = Q¯ [ (C − Q¯ ) where the first term is countable and the second is uncountable. Note also that thefirstis of measure zero while the second is of infinite measure.
    [Show full text]
  • The Representational Foundations of Computation
    The Representational Foundations of Computation Michael Rescorla Abstract: Turing computation over a non-linguistic domain presupposes a notation for the domain. Accordingly, computability theory studies notations for various non-linguistic domains. It illuminates how different ways of representing a domain support different finite mechanical procedures over that domain. Formal definitions and theorems yield a principled classification of notations based upon their computational properties. To understand computability theory, we must recognize that representation is a key target of mathematical inquiry. We must also recognize that computability theory is an intensional enterprise: it studies entities as represented in certain ways, rather than entities detached from any means of representing them. 1. A COMPUTATIONAL PERSPECTIVE ON REPRESENTATION Intuitively speaking, a function is computable when there exists a finite mechanical procedure that calculates the function’s output for each input. In the 1930s, logicians proposed rigorous mathematical formalisms for studying computable functions. The most famous formalism is the Turing machine: an abstract mathematical model of an idealized computing device with unlimited time and storage capacity. The Turing machine and other computational formalisms gave birth to computability theory: the mathematical study of computability. 2 A Turing machine operates over strings of symbols drawn from a finite alphabet. These strings comprise a formal language. In some cases, we want to study computation over the formal language itself. For example, Hilbert’s Entscheidungsproblem requests a uniform mechanical procedure that determines whether a given formula of first-order logic is valid. Items drawn from a formal language are linguistic types, whose tokens we can inscribe, concatenate, and manipulate [Parsons, 2008, pp.
    [Show full text]
  • Arxiv:2001.00578V2 [Math.HO] 3 Aug 2021
    Errata and Addenda to Mathematical Constants Steven Finch August 3, 2021 Abstract. We humbly and briefly offer corrections and supplements to Mathematical Constants (2003) and Mathematical Constants II (2019), both published by Cambridge University Press. Comments are always welcome. 1. First Volume 1.1. Pythagoras’ Constant. A geometric irrationality proof of √2 appears in [1]; the transcendence of the numbers √2 √2 i π √2 , ii , ie would follow from a proof of Schanuel’s conjecture [2]. A curious recursion in [3, 4] gives the nth digit in the binary expansion of √2. Catalan [5] proved the Wallis-like infinite product for 1/√2. More references on radical denestings include [6, 7, 8, 9]. 1.2. The Golden Mean. The cubic irrational ψ =1.3247179572... is connected to a sequence 3 ψ1 =1, ψn = 1+ ψn 1 for n 2 − ≥ which experimentally gives rise to [10] p n 1 lim (ψ ψn) 3(1 + ) =1.8168834242.... n − ψ →∞ The cubic irrational χ =1.8392867552... is mentioned elsewhere in the literature with arXiv:2001.00578v2 [math.HO] 3 Aug 2021 regard to iterative functions [11, 12, 13] (the four-numbers game is a special case of what are known as Ducci sequences), geometric constructions [14, 15] and numerical analysis [16]. Infinite radical expressions are further covered in [17, 18, 19]; more gen- eralized continued fractions appear in [20, 21]. See [22] for an interesting optimality property of the logarithmic spiral. A mean-value analog C of Viswanath’s constant 1.13198824... (the latter applies for almost every random Fibonacci sequence) was dis- covered by Rittaud [23]: C =1.2055694304..
    [Show full text]
  • The Maximum Modulus of a Trigonometric Trinomial
    The maximum modulus of a trigonometric trinomial Stefan Neuwirth Abstract Let Λ be a set of three integers and let CΛ be the space of 2π-periodic functions with spectrum in Λ endowed with the maximum modulus norm. We isolate the maximum modulus points x of trigonometric trinomials T ∈ CΛ and prove that x is unique unless |T | has an axis of symmetry. This enables us to compute the exposed and the extreme points of the unit ball of CΛ, to describe how the maximum modulus of T varies with respect to the arguments of its Fourier coefficients and to compute the norm of unimodular relative Fourier multipliers on CΛ. We obtain in particular the Sidon constant of Λ. 1 Introduction Let λ1, λ2 and λ3 be three pairwise distinct integers. Let r1, r2 and r3 be three positive real numbers. Given three real numbers t1, t2 and t3, let us consider the trigonometric trinomial i(t1+λ1x) i(t2+λ2x) i(t3+λ3x) T (x)= r1 e + r2 e + r3 e (1) for x R. The λ’s are the frequencies of the trigonometric trinomial T , the r’s are the moduli or ∈ it1 it2 it3 intensities and the t’s the arguments or phases of its Fourier coefficients r1 e , r2 e and r3 e . −1 −e iπ/3 Figure 1: The unit circle, the hypotrochoid H with equation z =4e −i2x +e ix, the segment from 1 to the unique point on H at maximum distance and the segments from e iπ/3 to the two points− on H at maximum distance.
    [Show full text]