Efficient Computation of Multidimensional Data Structures on CPU Platform

Total Page:16

File Type:pdf, Size:1020Kb

Efficient Computation of Multidimensional Data Structures on CPU Platform Eotv¨ os¨ Lor and´ University Faculty of Informatics Department of Programming Languages and Compilers Efficient computation of multidimensional data structures on CPU platform Supervisors: Departments of Computer Algebra and Prog. Lang. & Compilers Author: Attila Kov´acs Zal´an V´egv´ari Associate professor Software Information Technology BSc. Melinda T´oth Assistant lecturer ELTE-IK, PNYF Budapest, 2018. T´emabejelent˝o 1 Contents 1 Introduction 1 1.1 Fractals .................................. 1 1.1.1 Whatisafractal?......................... 1 1.1.2 Calculatingfractals. 2 1.1.3 TheEscape-TimeAlgorithm. 3 1.2 Generalization............................... 3 1.2.1 ExtendingTheEscape-TimeAlgorithm. 3 1.2.2 Extendingthespace ....................... 4 1.3 Auniversaltool.............................. 5 1.4 Scripting.................................. 5 1.5 Userinterface ............................... 6 2 User documentation 8 2.1 Thesolution................................ 8 2.2 NViewer.................................. 8 2.3 Graphicalengine ............................. 9 2.3.1 Thetransformationmatrix. 9 2.3.2 Thefunctions ........................... 10 2.3.3 Thelayers............................. 11 2.4 Goals.................................... 13 2.5 Thelanguage ............................... 13 2.5.1 Modules .............................. 14 2.5.2 Modes ............................... 14 2.5.3 Types ............................... 15 2.5.4 Parametersandstyles . 16 1 2.5.5 Queries .............................. 16 2.5.6 Definitions............................. 17 2.5.7 Nesteddefinitions. 17 2.5.8 Patterns .............................. 18 2.5.9 Branching ............................. 19 2.6 Standardlibrary ............................. 19 2.7 Toolusage................................. 23 2.7.1 Configuration ........................... 24 2.7.2 Sessions .............................. 24 2.7.3 Scripts............................... 25 2.7.4 Camera .............................. 25 2.7.5 Stylesandparameters . 26 2.7.6 Scenes ............................... 26 2.7.7 Exportingdata .......................... 27 3 Technical documentation 28 3.1 Problemdescription. 28 3.2 Enginestructure ............................. 29 3.3 Scriptstructure .............................. 30 3.3.1 Typesystem............................ 30 3.3.2 Functionstructure . 31 3.3.3 Informationflow . 31 3.4 Importingtheengine . 32 3.5 Importingthescriptparser. 32 3.6 Extendingthescriptlanguage . 33 3.7 Optimizations............................... 33 3.7.1 Functionstructure . 34 3.7.2 Constants ............................. 34 3.7.3 Memorymanagement. 34 3.8 Comparism ................................ 35 3.8.1 C++................................ 35 3.8.2 Haskell............................... 35 3.9 Examples ................................. 36 2 4 Summary 41 4.1 English................................... 41 4.2 Magyar .................................. 42 5 Results & testing 43 5.1 Testingplan................................ 43 5.2 Scripts................................... 44 Bibliography 48 6 Appendix 49 3 Abstract NViewer is a tool for visualizing N-dimensional shapes, more specifically fractals. This is not a tool to speed up something that was already known, it aims to uncover novel options for us. This tool allows us to be able to look into higher dimensions with an easy to handle camera technique. Further it allows us to discover structures unimaginable for us. NViewer provides us an easy way to transfer mathematical objects from paper to screen, in form of a scripting language. Additionally we are also providing a new type of graphical engine to render the results. Chapter 1 Introduction First, I intend to explain why this project was created and the problems it solves. The main objective is to create a tool that is easy to use not only for specific problems, but in more general cases. 1.1 Fractals Before we proceed any further, I need to clarify a number of concepts. 1.1.1 What is a fractal? A fractal is a geometrical or physical structure having an irregular or fragmented shape at all scales of measurement between a greatest and smallest scale such that certain mathematical or physical properties of the structure, as the perimeter of a curve or the flow rate in a porous medium, behave as if the dimensions of the structure (fractal dimensions) are greater than the spatial dimensions. [1] When somebody asks me what fractals are and I am to answer in a simple way, I would answer that these object represent a form of art created by mathematics and logic. 1 Figure 1.1: The mandelbrot set 1.1.2 Calculating fractals While fractals can be created in many different ways, they appear even in Nature. Plant development often shows fractal-like structures. Typical examples include plants such as broccoli, pinecone and basically most leaves. [4] In this project I will discuss only those, generated using com- puters. Some basic methods can be implemented in comenius logo in a recursive graphical representation pattern. A more complex way is to choose a set of points and run them across an iterated function system, also referred to as IFS [6]. In this case the set of points will be transformed in each step Figure 1.2: Broccoli which, in some cases, will approximate a set, which can be, to some extent, a fractal. As shown below, the nviewer project is capable of visualizing these types of fractals in an inefficient way, however my primary objective was to work with The Escape-Time Algorithm and some other similar algorithms. 2 1.1.3 The Escape-Time Algorithm The Escape-Time Algorithm is a widely used algorithm to visualize fractals. The basic idea is to use an iteration func- tion and test all points in the observed space to test if the iteration is convergent or not. If so, the point belongs to the fractal. Although this task is not entirely possible using a computer, approximations can be conducted. The most fa- mous application of this technique is the Mandelbrot Set [2]. Figure 1.3: Mandel- The convergence is tested after each iteration, as follows: if brot using low itera- the length of the current vector is greater than a specific tion settings value, it is not convergent and the iteration is interrupted. If this does not happen in a specific number of iterations, the original point (where the iteration stared), is part of the fractal at the current approximately level. For the points that are not part of the fractal, the time when they ”escape” is used for colouring, which both shows how the iterations approximates the fractal and provides an aestethical image. 1.2 Generalization Ok. Now a tool is available which is able to render 3D fractals, which is itself already useful, further, t many other objects can be created using the same approach. 1.2.1 Extending The Escape-Time Algorithm For proper functioning, the Escape-Time Algorithm requires the availability of the following features: per-pixel data management, graphical engine to render the re- sults, mathematical core to perform the necessary calculations at an acceptable speed. Once these features are available, it becomes easy to extend some of the features. The Escape-Time Algorithm provides escape-time to the points of the space, so that these can be coloured, where the colouring is carried out as a result 3 of a built-in process. Therefore, the graphical engine will not use this algorithm, but it will use more general functions that assign colours to points. This way, the colouring methods for the Escape-Time Algorithm can ban modified manually or it allows the implementation of various algorithms. At this point, we are able to create a function that colours each point in the space for rendering. It allows us to render any shape or object that can be expressed in such a function. Though this is not enough, we cannot yet apply light or any other effect to the surface. This sounds unimportant for the first time, because it seems to be only a special effect. When working in 3D, lighting is one of the most important tasks, because it allows us to visualize the surface of the object. Figure 1.4 shows a part of the mandelbox with and without lighting. If no light- ing is present, fog effect had to be used with a special colouring to make a surface patterns visible. Without these effects, only the outline of the mandelbox is visible. The second image show what it looks like, when lights are applied. There is a slight green fog that makes the distances more visible. So far this can be done using the graphical engine, but it might not satisfy all our needs, because lighting is applied in the same way to all points of the surface. This issue is continued (section 2.3.3). (a) Without lighting (b) With lighting Figure 1.4: Mandelbox with and without lighting 1.2.2 Extending the space A lot of 2D fractals have been recognized. They are easy to generate, easy to examine and at the same time they are beautiful. Additionally, three and more dimensional 4 fractals are known. Of course, besides fractals, a 4D model of the black holes has also been described [3]. This creates a need for a graphical engine that is capable of rendering N-dimensional objects. Usually, tools that operate in N-dimension, simply conduct a 2D or 3D projection of the original object. Is that enough? Projections cover some details, because only the surface will be visible. 1.3 A universal tool Most tools are designed to serve one purpose only. That is not wrong at all, because specialized tools are usually more optimized, but they have a great disadvantage: they
Recommended publications
  • Iterated Function System Quasiarcs
    CONFORMAL GEOMETRY AND DYNAMICS An Electronic Journal of the American Mathematical Society Volume 21, Pages 78–100 (February 3, 2017) http://dx.doi.org/10.1090/ecgd/305 ITERATED FUNCTION SYSTEM QUASIARCS ANNINA ISELI AND KEVIN WILDRICK Abstract. We consider a class of iterated function systems (IFSs) of contract- ing similarities of Rn, introduced by Hutchinson, for which the invariant set possesses a natural H¨older continuous parameterization by the unit interval. When such an invariant set is homeomorphic to an interval, we give necessary conditions in terms of the similarities alone for it to possess a quasisymmetric (and as a corollary, bi-H¨older) parameterization. We also give a related nec- essary condition for the invariant set of such an IFS to be homeomorphic to an interval. 1. Introduction Consider an iterated function system (IFS) of contracting similarities S = n {S1,...,SN } of R , N ≥ 2, n ≥ 1. For i =1,...,N, we will denote the scal- n n ing ratio of Si : R → R by 0 <ri < 1. In a brief remark in his influential work [18], Hutchinson introduced a class of such IFSs for which the invariant set is a Peano continuum, i.e., it possesses a continuous parameterization by the unit interval. There is a natural choice for this parameterization, which we call the Hutchinson parameterization. Definition 1.1 (Hutchinson, 1981). The pair (S,γ), where γ is the invariant set of an IFS S = {S1,...,SN } with scaling ratio list {r1,...,rN } is said to be an IFS path, if there exist points a, b ∈ Rn such that (i) S1(a)=a and SN (b)=b, (ii) Si(b)=Si+1(a), for any i ∈{1,...,N − 1}.
    [Show full text]
  • An Introduction to Apophysis © Clive Haynes MMXX
    Apophysis Fractal Generator An Introduction Clive Haynes Fractal ‘Flames’ The type of fractals generated are known as ‘Flame Fractals’ and for the curious, I append a note about their structure, gleaned from the internet, at the end of this piece. Please don’t ask me to explain it! Where to download Apophysis: go to https://sourceforge.net/projects/apophysis/ Sorry Mac users but it’s only available for Windows. To see examples of fractal images I’ve generated using Apophysis, I’ve made an Issuu e-book and here’s the URL. https://issuu.com/fotopix/docs/ordering_kaos Getting Started There’s not a defined ‘follow this method workflow’ for generating interesting fractals. It’s really a matter of considerable experimentation and the accumulation of a knowledge-base about general principles: what the numerous presets tend to do and what various options allow. Infinite combinations of variables ensure there’s also a huge serendipity factor. I’ve included a few screen-grabs to help you. The screen-grabs are detailed and you may need to enlarge them for better viewing. Once Apophysis has loaded, it will provide a Random Batch of fractal patterns. Some will be appealing whilst many others will be less favourable. To generate another set, go to File > Random Batch (shortcut Ctrl+B). Screen-grab 1 Choose a fractal pattern from the batch and it will appear in the main window (Screen-grab 1). Depending upon the complexity of the fractal and the processing power of your computer, there will be a ‘wait time’ every time you change a parameter.
    [Show full text]
  • William Broadhead Professor of History Thesis Supervisor
    Vitruvius on Architecture: A Modem Application and Stability Analysis of Classical Structures by Ana S. Escalante Submitted to the Department of Mechanical Engineering in Partial Fulfillment of the Requirements for the Degree of Bachelor of Science in Mechanical Engineering at the Lj Massachusetts Institute of Technology June 2013 0 2013 Ana S. Escalante. All rights reserved. The author hereby grants to MIT permission to reproduce and to distribute publicly paper and electronic copies of this thesis document in whole or in part in any medium now known or hereafter created Signature of Author: F , 2 I- - Department of Mechanical Engineering May 28, 2013 Certified by: William Broadhead Professor of History Thesis Supervisor Accepted by: Annette Hosoi Professor of Mechanical Engineering Undergraduate Officer 1 2 Vitruvius on Architecture: A Modem Application and Stability Analysis of Classical Structures by Ana S. Escalante Submitted to the Department of Mechanical Engineering on June 7, 2013 in Partial Fulfillment of the Requirements for the Degree of Bachelor of Science in Mechanical Engineering ABSTRACT Imperial Rome has left numerous legacies, the most well-known being its literature and monuments. Though many monuments, such as the Pantheon, are well-preserved, in cases where little physical evidence remains, historians can often use literary sources to inform reconstruction efforts. For more technical studies of Roman construction, technical literature is rare and the contemporary awareness of such literature even less known. When Vitruvius wrote De architectura,he did not intend for it to be a manual for instruction but rather a central source of general architectural knowledge. Directly aimed at architects, contractors, and other individuals involved in the design and construction of buildings, De architecturaprovides insight into contemporary technical knowledge.
    [Show full text]
  • A Mathematician's Lament
    A Mathematician’s Lament by Paul Lockhart musician wakes from a terrible nightmare. In his dream he finds himself in a society where A music education has been made mandatory. “We are helping our students become more competitive in an increasingly sound-filled world.” Educators, school systems, and the state are put in charge of this vital project. Studies are commissioned, committees are formed, and decisions are made— all without the advice or participation of a single working musician or composer. Since musicians are known to set down their ideas in the form of sheet music, these curious black dots and lines must constitute the “language of music.” It is imperative that students become fluent in this language if they are to attain any degree of musical competence; indeed, it would be ludicrous to expect a child to sing a song or play an instrument without having a thorough grounding in music notation and theory. Playing and listening to music, let alone composing an original piece, are considered very advanced topics and are generally put off until college, and more often graduate school. As for the primary and secondary schools, their mission is to train students to use this language— to jiggle symbols around according to a fixed set of rules: “Music class is where we take out our staff paper, our teacher puts some notes on the board, and we copy them or transpose them into a different key. We have to make sure to get the clefs and key signatures right, and our teacher is very picky about making sure we fill in our quarter-notes completely.
    [Show full text]
  • Fractal Dimension of Self-Affine Sets: Some Examples
    FRACTAL DIMENSION OF SELF-AFFINE SETS: SOME EXAMPLES G. A. EDGAR One of the most common mathematical ways to construct a fractal is as a \self-similar" set. A similarity in Rd is a function f : Rd ! Rd satisfying kf(x) − f(y)k = r kx − yk for some constant r. We call r the ratio of the map f. If f1; f2; ··· ; fn is a finite list of similarities, then the invariant set or attractor of the iterated function system is the compact nonempty set K satisfying K = f1[K] [ f2[K] [···[ fn[K]: The set K obtained in this way is said to be self-similar. If fi has ratio ri < 1, then there is a unique attractor K. The similarity dimension of the attractor K is the solution s of the equation n X s (1) ri = 1: i=1 This theory is due to Hausdorff [13], Moran [16], and Hutchinson [14]. The similarity dimension defined by (1) is the Hausdorff dimension of K, provided there is not \too much" overlap, as specified by Moran's open set condition. See [14], [6], [10]. REPRINT From: Measure Theory, Oberwolfach 1990, in Supplemento ai Rendiconti del Circolo Matematico di Palermo, Serie II, numero 28, anno 1992, pp. 341{358. This research was supported in part by National Science Foundation grant DMS 87-01120. Typeset by AMS-TEX. G. A. EDGAR I will be interested here in a generalization of self-similar sets, called self-affine sets. In particular, I will be interested in the computation of the Hausdorff dimension of such sets.
    [Show full text]
  • On Bounding Boxes of Iterated Function System Attractors Hsueh-Ting Chu, Chaur-Chin Chen*
    Computers & Graphics 27 (2003) 407–414 Technical section On bounding boxes of iterated function system attractors Hsueh-Ting Chu, Chaur-Chin Chen* Department of Computer Science, National Tsing Hua University, Hsinchu 300, Taiwan, ROC Abstract Before rendering 2D or 3D fractals with iterated function systems, it is necessary to calculate the bounding extent of fractals. We develop a new algorithm to compute the bounding boxwhich closely contains the entire attractor of an iterated function system. r 2003 Elsevier Science Ltd. All rights reserved. Keywords: Fractals; Iterated function system; IFS; Bounding box 1. Introduction 1.1. Iterated function systems Barnsley [1] uses iterated function systems (IFS) to Definition 1. A transform f : X-X on a metric space provide a framework for the generation of fractals. ðX; dÞ is called a contractive mapping if there is a Fractals are seen as the attractors of iterated function constant 0pso1 such that systems. Based on the framework, there are many A algorithms to generate fractal pictures [1–4]. However, dðf ðxÞ; f ðyÞÞps Á dðx; yÞ8x; y X; ð1Þ in order to generate fractals, all of these algorithms have where s is called a contractivity factor for f : to estimate the bounding boxes of fractals in advance. For instance, in the program Fractint (http://spanky. Definition 2. In a complete metric space ðX; dÞ; an triumf.ca/www/fractint/fractint.html), we have to guess iterated function system (IFS) [1] consists of a finite set the parameters of ‘‘image corners’’ before the beginning of contractive mappings w ; for i ¼ 1; 2; y; n; which is of drawing, which may not be practical.
    [Show full text]
  • Math Morphing Proximate and Evolutionary Mechanisms
    Curriculum Units by Fellows of the Yale-New Haven Teachers Institute 2009 Volume V: Evolutionary Medicine Math Morphing Proximate and Evolutionary Mechanisms Curriculum Unit 09.05.09 by Kenneth William Spinka Introduction Background Essential Questions Lesson Plans Website Student Resources Glossary Of Terms Bibliography Appendix Introduction An important theoretical development was Nikolaas Tinbergen's distinction made originally in ethology between evolutionary and proximate mechanisms; Randolph M. Nesse and George C. Williams summarize its relevance to medicine: All biological traits need two kinds of explanation: proximate and evolutionary. The proximate explanation for a disease describes what is wrong in the bodily mechanism of individuals affected Curriculum Unit 09.05.09 1 of 27 by it. An evolutionary explanation is completely different. Instead of explaining why people are different, it explains why we are all the same in ways that leave us vulnerable to disease. Why do we all have wisdom teeth, an appendix, and cells that if triggered can rampantly multiply out of control? [1] A fractal is generally "a rough or fragmented geometric shape that can be split into parts, each of which is (at least approximately) a reduced-size copy of the whole," a property called self-similarity. The term was coined by Beno?t Mandelbrot in 1975 and was derived from the Latin fractus meaning "broken" or "fractured." A mathematical fractal is based on an equation that undergoes iteration, a form of feedback based on recursion. http://www.kwsi.com/ynhti2009/image01.html A fractal often has the following features: 1. It has a fine structure at arbitrarily small scales.
    [Show full text]
  • Iterated Function System
    IARJSET ISSN (Online) 2393-8021 ISSN (Print) 2394-1588 International Advanced Research Journal in Science, Engineering and Technology ISO 3297:2007 Certified Vol. 3, Issue 8, August 2016 Iterated Function System S. C. Shrivastava Department of Applied Mathematics, Rungta College of Engineering and Technology, Bhilai C.G., India Abstract: Fractal image compression through IFS is very important for the efficient transmission and storage of digital data. Fractal is made up of the union of several copies of itself and IFS is defined by a finite number of affine transformation which characterized by Translation, scaling, shearing and rotat ion. In this paper we describe the necessary conditions to form an Iterated Function System and how fractals are generated through affine transformations. Keywords: Iterated Function System; Contraction Mapping. 1. INTRODUCTION The exploration of fractal geometry is usually traced back Metric Spaces definition: A space X with a real-valued to the publication of the book “The Fractal Geometry of function d: X × X → ℜ is called a metric space (X, d) if d Nature” [1] by the IBM mathematician Benoit B. possess the following properties: Mandelbrot. Iterated Function System is a method of constructing fractals, which consists of a set of maps that 1. d(x, y) ≥ 0 for ∀ x, y ∈ X explicitly list the similarities of the shape. Though the 2. d(x, y) = d(y, x) ∀ x, y ∈ X formal name Iterated Function Systems or IFS was coined 3. d x, y ≤ d x, z + d z, y ∀ x, y, z ∈ X . (triangle by Barnsley and Demko [2] in 1985, the basic concept is inequality).
    [Show full text]
  • An Introduction to the Mandelbrot Set
    An introduction to the Mandelbrot set Bastian Fredriksson January 2015 1 Purpose and content The purpose of this paper is to introduce the reader to the very useful subject of fractals. We will focus on the Mandelbrot set and the related Julia sets. I will show some ways of visualising these sets and how to make a program that renders them. Finally, I will explain a key exchange algorithm based on what we have learnt. 2 Introduction The Mandelbrot set and the Julia sets are sets of points in the complex plane. Julia sets were first studied by the French mathematicians Pierre Fatou and Gaston Julia in the early 20th century. However, at this point in time there were no computers, and this made it practically impossible to study the structure of the set more closely, since large amount of computational power was needed. Their discoveries was left in the dark until 1961, when a Jewish-Polish math- ematician named Benoit Mandelbrot began his research at IBM. His task was to reduce the white noise that disturbed the transmission on telephony lines [3]. It seemed like the noise came in bursts, sometimes there were a lot of distur- bance, and sometimes there was no disturbance at all. Also, if you examined a period of time with a lot of noise-problems, you could still find periods without noise [4]. Could it be possible to come up with a model that explains when there is noise or not? Mandelbrot took a quite radical approach to the problem at hand, and chose to visualise the data.
    [Show full text]
  • Iterated Function System Iterated Function System
    Iterated Function System Iterated Function System (or I.F.S. for short) is a program that explores the various aspects of the Iterated Function System. It contains three modes: 1. The Chaos Game, which implements variations of the chaos game popularized by Michael Barnsley, 2. Deterministic Mode, which creates fractals using iterated function system rules in a deterministic way, and 3. Random Mode, which creates fractals using iterated function system rules in a random way. Upon execution of I.F.S. we need to choose a mode. Let’s begin by choosing the Chaos Game. Upon selecting that mode, two windows will appear as in Figure 1. The blank one on the right is the graph window and is where the result of the execution of the chaos game will be displayed. The window on the left is the status window. It is set up to play the chaos game using three corners and a scale factor of 0.5. Thus, the next seed point will be calculated as the midpoint of the current seed point with the selected corner. It is also set up so that the corners will be graphed once they are selected, and a line is not drawn between the current seed point and the next selected corner. As a first attempt at playing the chaos game, let’s press the “Go!” button without changing any of these settings. We will first be asked to define corner #1, which can be accomplished either by typing in the xy-coordinates in the status window edit boxes, or by using the mouse to move the cross hair cursor appearing in the graph window.
    [Show full text]
  • Pointwise Regularity of Parametrized Affine Zipper Fractal Curves
    POINTWISE REGULARITY OF PARAMETERIZED AFFINE ZIPPER FRACTAL CURVES BALAZS´ BAR´ ANY,´ GERGELY KISS, AND ISTVAN´ KOLOSSVARY´ Abstract. We study the pointwise regularity of zipper fractal curves generated by affine mappings. Under the assumption of dominated splitting of index-1, we calculate the Hausdorff dimension of the level sets of the pointwise H¨olderexpo- nent for a subinterval of the spectrum. We give an equivalent characterization for the existence of regular pointwise H¨olderexponent for Lebesgue almost every point. In this case, we extend the multifractal analysis to the full spectrum. In particular, we apply our results for de Rham's curve. 1. Introduction and Statements Let us begin by recalling the general definition of fractal curves from Hutchin- son [22] and Barnsley [3]. d Definition 1.1. A system S = ff0; : : : ; fN−1g of contracting mappings of R to itself is called a zipper with vertices Z = fz0; : : : ; zN g and signature " = ("0;:::;"N−1), "i 2 f0; 1g, if the cross-condition fi(z0) = zi+"i and fi(zN ) = zi+1−"i holds for every i = 0;:::;N − 1. We call the system a self-affine zipper if the functions fi are affine contractive mappings of the form fi(x) = Aix + ti; for every i 2 f0; 1;:::;N − 1g; d×d d where Ai 2 R invertible and ti 2 R . The fractal curve generated from S is the unique non-empty compact set Γ, for which N−1 [ Γ = fi(Γ): i=0 If S is an affine zipper then we call Γ a self-affine curve.
    [Show full text]
  • Does Mathematical Beauty Pose Problems for Naturalism?1
    Does Mathematical Beauty Pose Problems for Naturalism?1 Russell W. Howell Professor of Mathematics Westmont College, Santa Barbara, CA [email protected] Numerous events occurred in 1960 whose eects could hardly have been predicted at the time: several African Americans staged a sit-in at a Greensboro lunch counter, the Soviet Union shot down Gary Powers while he was ying a U2 spy plane, the US FDA approved the use of the rst oral contraceptive, AT&T led with the Federal Communications Commission for permission to launch an experimental communications satellite, and four Presidential debates between John Kennedy and Richard Nixon aired on national television. Less well known was the publication of a paper by the physicist Eugene Wigner. Ap- pearing in Communications in Pure and Applied Mathematics, a journal certainly not widely read by the general public, it bore the mysterious title The Unreasonable Eectiveness of Mathematics in the Natural Sciences. (Wigner, 1960) Like our cultural examples of the 1960's, it has had eects beyond what most people would have imagined. Our purpose here is to tease out some strains of an important question that has emerged from Wigner's work. Wigner begins with a story about two friends who were discussing their jobs. One of them, a statistician, was working on population trends. He mentioned a paper he had produced, which contained the Gaussian distribution equation near the beginning. The statistician attempted to explain the meaning it, as well as other mathematical symbols. His friend was a bit perplexed, and was not quite sure whether the statistician was pulling his leg.
    [Show full text]