DOCSLIB.ORG

M.Sc. in Mathematical Modelling and Scientific Computing

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
M.Sc. in Mathematical Modelling and Scientific Computing M.Sc. in Mathematical Modelling and Scientific Computing Dissertation Projects December 2019 Contents 1 Numerical Analysis Projects 3 1.1 Operator Preconditioning and Accuracy of Linear System Solutions . 3 1.2 Computing Common Zeros of Trivariate Functions . 4 1.3 Univariate Optimisation for Smooth Functions . 4 1.4 Algorithms for Least-Squares and Underdetermined Linear Systems Based on CholeskyQR . 5 1.5 Tensor and Matrix Eigenvalue Perturbation Theory . 6 1.6 Chebfun Dissertation Topics . 7 1.7 Rational Functions Dissertation Topics . 8 1.8 Low-Rank Plus Sparse Matrix Model and its Application . 9 1.9 Efficient Exemplar Selection for Representation Learning . 10 1.10 Numerical Solution of a Problem in Electrochemistry . 11 1.11 Nonconvex Optimisation . 12 2 Biological and Medical Application Projects 13 2.1 Localisation of Turing Patterning in Discrete Heterogeneous Media . 13 2.2 Distributed Delay in Reaction-Diffusion Systems . 14 2.3 Aspects of Flagellate Swimmer Dynamics and Mechanics . 16 2.4 Learning Time Dependence with Minimal Time Information . 17 2.5 Surrogate Modelling for Particle Infiltration into Tumour Spheroids . 18 2.6 Mathematical Modelling of polyelectrolyte hydrogels for Application in Regenerative Medicine . 19 3 Physical Application Projects 21 3.1 Evaporation-Driven Instabilities in Complex Fluids . 21 1 3.2 Phase Separation in Swollen Hydrogels . 22 3.3 Pattern Formation in Polymers . 23 4 Data Science 25 4.1 Computational Topology for Analysing Neuronal Branching Models and Data . 25 5 Research Interests of Academic Staff 26 2 1 Numerical Analysis Projects 1.1 Operator Preconditioning and Accuracy of Linear System Solu- tions Supervisors: Prof Yuji Nakatsukasa and Dr Carolina Urzua Torres Contact: [email protected] and [email protected] Preconditioning is widely recognised as a crucial technique for improving the convergence of iterative (usually Krylov) methods for linear systems Ax = b. A major body of research in numerical linear algebra is devoted to devising preconditioners for certain classes of matrices. One aspect that is seldom discussed is that if A is ill-conditioned, then even if a good preconditioner M is used, usually the final (relative) accuracy achievable is limited by the condition number of the original linear system O(uκ2(A)), where u is the unit roundoff. This is because applying A or M involves errors of that order. Another, more recent, line of work is operator preconditioning (e.g. [1, 3]) in the context of solving differential equations (which usually lead to a large linear system). Here the idea is to precondition the problem at the operator level to reduce the problem to that of a well-conditioned differential operator, which is then discretised to give Ax~ = ~b. This procedure hints that the effect of numerical errors could be significantly reduced, and hence the accuracy would no longer be limited by uκ2(A) but by uκ2(A~), which can be much better, potentially as small as O(u); a similar effect is described in [2]. The goal of this work is to explore operator preconditioning in terms of the solution accuracy. The project could involve a subset of the following: (i) analyzing the properties of the linear systems and solutions (forward and backward errors), (ii) roundoff error analysis, (iii) examining various equations and preconditioners, and (iv) implementing the preconditioners. The study might also lead to the design of good preconditioners at the matrix level. Exposure to basic numerical linear algebra and differential equations would be helpful for the project. References [1] R. Hiptmair. Operator preconditioning. Computers and Mathematics with Applica- tions, 52(5):699{706, 2006. [2] S. Olver and A. Townsend. A fast and well-conditioned spectral method. SIAM Rev., 55(3):462{489, 2013. [3] C. Urzua-Torres, R. Hiptmair, and C. Jerez-Hanckes. Optimal operator precondi- tioning for Galerkin boundary element methods on 3D screens. SIAM J. Numer. Anal., 2019. 3 1.2 Computing Common Zeros of Trivariate Functions Supervisor: Prof Yuji Nakatsukasa Contact: [email protected] Chebfun can deal with functions in 1D, 2D, and 3D. Its functionalities in 1D are mostly complete and highly sophisticated. In higher dimensions, and in particular 3D, the situation is sometimes different. Among the key missing operations in Chebfun3 is roots(f,g,h)—finding common zeros of three trivariate functions, which generically have zero-dimensional solutions consisting of disjoint points. The 1D version roots(f) is re- markably efficient and reliable, based on domain subdivision and colleague eigenvalue problems. The 2D version was developed in [1], which also uses these two ideas, together with resultant methods to eliminate variables. The goal of this project is to develop and implement a Chebfun3 roots(f,g,h) algorithm and command. It is envisaged that a successful algorithm will involve a mixture of mathematical analysis and numerical/computational techniques. In particular, condi- tioning and numerical stability will need to be investigated carefully, in light of the analysis in [2]. Since there is little software publicly available for trivariate rootfinding, a successful completion would be an important contribution to multivariate polynomial computing. Background in approximation theory would be very helpful. The Part C course \Ap- proximation of Functions" is highly recommended. References [1] Y. Nakatsukasa, V. Noferini, and A. Townsend. Computing the common zeros of two bivariate functions via B´ezoutresultants. Numer. Math., 129:181{209, 2015. [2] V. Noferini and A. Townsend. Numerical instability of resultant methods for multi- dimensional rootfinding. SIAM J. Numer. Anal., 54(2):719{743, 2016. 1.3 Univariate Optimisation for Smooth Functions Supervisor: Prof Yuji Nakatsukasa Contact: [email protected] Finding the (global) minimum of a univariate function minx f(x) is a classical problem in continuous optimisation. It arises naturally when a search direction is found, and one looks for an appropriate step size, e.g. in line search algorithms. For example in the currently very popular stochastic gradient descent method, usually the step size is fixed. However, an optimal step size can significantly outperform fixed and suboptimal choices. When the function is convex (or unimodal), bisection or golden section are established techniques for finding the global minimum. Another powerful approach would be to approximate the function with a polynomial and minimize the polynomial. This tech- nique is used in Chebfun (together with interval subdivision as necessary), but appears 4 not to be widely used elsewhere. We expect this approach to have a number of ad- vantages, including (i) it can easily deal with nonconvex functions, (ii) it can find local minima/maxima, and (iii) the convergence should be fast (fewer samples are needed) if f is smooth, e.g. analytic. The goals of this project include exploring these aspects, analyzing convergence, and developing efficient implementations. Extensions can be considered, such as using rational approximation or minimising bivariate functions. The student should have some prior exposure to basic optimisation [1], and ideally approximation theory [2]. The Part C courses \Continuous Optimisation" and \Ap- proximation of Functions" are recommended. References [1] J. Nocedal and S. J. Wright. Numerical Optimization. Springer New York, second edition, 1999. [2] L. N. Trefethen. Approximation Theory and Approximation Practice. SIAM, Philadelphia, 2013. 1.4 Algorithms for Least-Squares and Underdetermined Linear Sys- tems Based on CholeskyQR Supervisor: Prof Yuji Nakatsukasa Contact: [email protected] CholeskyQR|based on RT R = AT A (Cholesky factorization) and Q = AR−1|is an algorithm for computing the QR factorization of matrices, which is very attractive in speed especially in parallel computing architectures. It is unfortunately numerically unstable, and hence has not been widely used. However, the instability can be improved significantly by techniques such as repeating and shifting [1, 2]. The QR factorization is used in various problems in numerical linear algebra, including least-squares problems and underdetermined linear systems. This project aims to explore the use of CholeskyQR in these contexts: (1) We observe that when CholeskyQR is used for least-squares problems (together with a Krylov solver; a neat mixture of direct and iterative linear algebra methods), we arrive at a robust and efficient solver, achieving up to ×5 speedup. Establishing backward stability, however, remains an important open problem, which we hope to address here. Other possible directions include exploring HPC implementations and exploring its use in applications. (2) Similarly, CholeskyQR can be used for finding the minimal solution for an under- determined linear system. This line of work is at a young stage and many directions are possible: proving stability, refining the algorithm, implementation, and applications (e.g. compressed sensing). A good background in linear algebra is desirable (e.g. the part C course Numerical Linear Algebra). 5 References [1] T. Fukaya, R. Kannan, Y. Nakatsukasa, Y. Yamamoto, and Y. Yanagisawa. Shifted CholeskyQR for computing the QR factorization of ill-conditioned matrices, arXiv:1809.11085, 2018. [2] Y. Yamamoto, Y. Nakatsukasa, Y. Yanagisawa, and T. Fukaya. Roundoff error analysis of the CholeskyQR2 algorihm. Electron. Trans. Numer. Anal., 44:306{326, 2015. 1.5 Tensor and Matrix Eigenvalue Perturbation Theory
Load more

Recommended publications
  • Numerical Computation with Rational Functions
    NUMERICAL COMPUTATION WITH RATIONAL FUNCTIONS (scalars only) Nick Trefethen, University of Oxford and ENS Lyon + thanks to Silviu Filip, Abi Gopal, Stefan Güttel, Yuji Nakatsukasa, and Olivier Sète 1/26 1. Polynomial vs. rational approximations 2. Four representations of rational functions 2a. Quotient of polynomials 2b. Partial fractions 2c. Quotient of partial fractions (= barycentric) 2d. Transfer function/matrix pencil 3. The AAA algorithm with Nakatsukasa and Sète, to appear in SISC 4. Application: conformal maps with Gopal, submitted to Numer. Math. 5. Application: minimax approximation with Filip, Nakatsukasa, and Beckermann, to appear in SISC 6. Accuracy and noise 2/26 1. Polynomial vs. rational approximation Newman, 1964: approximation of |x| on [−1,1] −휋 푛 퐸푛0 ~ 0.2801. ./푛 , 퐸푛푛 ~ 8푒 Poles and zeros of r: exponentially clustered near x=0, exponentially diminishing residues. Rational approximation is nonlinear, so algorithms are nontrivial. poles There may be nonuniqueness and local minima. type (20,20) 3/26 2. Four representations of rational functions Alpert, Carpenter, Coelho, Gonnet, Greengard, Hagstrom, Koerner, Levy, Quotient of polynomials 푝(푧)/푞(푧) SK, IRF, AGH, ratdisk Pachón, Phillips, Ruttan, Sanathanen, Silantyev, Silveira, Varga, White,… 푎푘 Beylkin, Deschrijver, Dhaene, Drmač, Partial fractions ෍ 푧 − 푧 VF, exponential sums Greengard, Gustavsen, Hochman, 푘 Mohlenkamp, Monzón, Semlyen,… Berrut, Filip, Floater, Gopal, Quotient of partial fractions 푛(푧)/푑(푧) Floater-Hormann, AAA Hochman, Hormann, Ionita, Klein, Mittelmann, Nakatsukasa, Salzer, (= barycentric) Schneider, Sète, Trefethen, Werner,… 푇 −1 Antoulas, Beattie, Beckermann, Berljafa, Transfer function/matrix pencil 푐 푧퐵 − 퐴 푏 IRKA, Loewner, RKFIT Druskin, Elsworth, Gugercin, Güttel, Knizhnerman, Meerbergen, Ruhe,… Sometimes the boundaries are blurry! 4/26 2a.
    [Show full text]
  • An Extension of Matlab to Continuous Functions and Operators∗
    SIAM J. SCI. COMPUT. c 2004 Society for Industrial and Applied Mathematics Vol. 25, No. 5, pp. 1743–1770 AN EXTENSION OF MATLAB TO CONTINUOUS FUNCTIONS AND OPERATORS∗ † † ZACHARY BATTLES AND LLOYD N. TREFETHEN Abstract. An object-oriented MATLAB system is described for performing numerical linear algebra on continuous functions and operators rather than the usual discrete vectors and matrices. About eighty MATLAB functions from plot and sum to svd and cond have been overloaded so that one can work with our “chebfun” objects using almost exactly the usual MATLAB syntax. All functions live on [−1, 1] and are represented by values at sufficiently many Chebyshev points for the polynomial interpolant to be accurate to close to machine precision. Each of our overloaded operations raises questions about the proper generalization of familiar notions to the continuous context and about appropriate methods of interpolation, differentiation, integration, zerofinding, or transforms. Applications in approximation theory and numerical analysis are explored, and possible extensions for more substantial problems of scientific computing are mentioned. Key words. MATLAB, Chebyshev points, interpolation, barycentric formula, spectral methods, FFT AMS subject classifications. 41-04, 65D05 DOI. 10.1137/S1064827503430126 1. Introduction. Numerical linear algebra and functional analysis are two faces of the same subject, the study of linear mappings from one vector space to another. But it could not be said that mathematicians have settled on a language and notation that blend the discrete and continuous worlds gracefully. Numerical analysts favor a concrete, basis-dependent matrix-vector notation that may be quite foreign to the functional analysts. Sometimes the difference may seem very minor between, say, ex- pressing an inner product as (u, v)orasuTv.
    [Show full text]
  • Householder Triangularization of a Quasimatrix
    IMA Journal of Numerical Analysis (2010) 30, 887–897 doi:10.1093/imanum/drp018 Advance Access publication on August 21, 2009 Householder triangularization of a quasimatrix LLOYD N.TREFETHEN† Oxford University Computing Laboratory, Wolfson Building, Parks Road, Downloaded from Oxford OX1 3QD, UK [Received on 4 July 2008] A standard algorithm for computing the QR factorization of a matrix A is Householder triangulariza- tion. Here this idea is generalized to the situation in which A is a quasimatrix, that is, a ‘matrix’ whose http://imajna.oxfordjournals.org/ ‘columns’ are functions defined on an intervala [ , b]. Applications are mentioned to quasimatrix least squares fitting, singular value decomposition and determination of ranks, norms and condition numbers, and numerical illustrations are presented using the chebfun system. Keywords: Householder triangularization; QR factorization; chebfun; quasimatrix; singular value decom- position. 1. QR factorization and Householder triangularization at Radcliffe Science Library, Bodleian Library on April 18, 2012 Let A be an m n matrix, where m n. A (reduced) QR factorization of A is a factorization × > x x x x q q q q r r r r x x x x q q q q r r r A QR, x x x x q q q q r r , (1.1) = x x x x q q q q r x x x x q q q q x x x x q q q q where Q is m n with orthonormal columns and R is upper triangular. Here and in subsequent equations we illustrate our× matrix operations by schematic pictures in which x denotes an arbitrary entry, r denotes an entry of an upper triangular matrix, q denotes an entry of an orthonormal column and a blank denotes a zero.
    [Show full text]
  • Publications by Lloyd N. Trefethen
    Publications by Lloyd N. Trefethen I. Books I.1. Numerical Conformal Mapping, editor, 269 pages. Elsevier, 1986. I.2. Finite Difference and Spectral Methods for Ordinary and Partial Differential Equations, x+315 pages. Graduate textbook, privately published at http://comlab.ox.ac.uk/nick.trefethen/, 1996. I.3. Numerical Linear Algebra, with David Bau III, xii+361 pages. SIAM, 1997. I.4. Spectral Methods in MATLAB, xviii+165 pages. SIAM, 2000. I.5. Schwarz-Christoffel Mapping, with Tobin A. Driscoll, xvi+132 pages. Cambridge U. Press, 2002. I.6. Spectra and Pseudospectra: The Behavior of Nonnormal Matrices and Operators, with Mark Embree, xviii+606 pages. Princeton U. Press, 2005. I.7. Trefethen’s Index Cards: Forty years of Notes about People, Words, and Mathematics, xv + 368 pages. World Scientific Publishing, 2011. I.8. Approximation Theory and Approximation Practice. SIAM, to appear in 2013. In addition to these a book has been published by other authors about the international problem solving challenge I organised in 2002 (item VIII.8 below): F. Bornemann, D. Laurie, S. Wagon and D. Waldvogel, The SIAM 100-Digit Challenge, SIAM, 2004, to appear in German translation as Abenteuer Numerik: Eine Führung entlang der Herausforderung von Trefethen, Springer, 2006. II. Finite difference and spectral methods for partial differential equations II.1. Group velocity in finite difference schemes. SIAM Review 24 (1982), 113-136. II.2. Group velocity interpretation of the stability theory of Gustafsson, Kreiss, and Sundström. J. Comp. Phys. 49 (1983), 199-217. II.3. On Lp-instability and dispersion at discontinuities in finite difference schemes.
    [Show full text]
  • Efficient Shadowing of High Dimensional Chaotic Systems With
    Efficient Shadowing of High Dimensional Chaotic Systems with the Large Astrophysical N-body Problem as an Example by Wayne Hayes A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Computer Science University of Toronto c Copyright by Wayne Hayes 1995 Efficient Shadowing of High Dimensional Chaotic Systems with the Large Astrophysical N-body Problem as an Example Wayne Hayes Master of Science Department of Computer Science University of Toronto January 1995 Abstract N-body systems are chaotic. This means that numerical errors in their solution are magnified exponentially with time, perhaps making the solutions useless. Shadowing tries to show that numerical solutions of chaotic systems still have some validity. A previously published shadowing refinement algorithm is optimized to give speedups of order 60 for small problems and asymptotic speedups of O(N) on large problems. This optimized algorithm is used to shadow N-body systems with up to 25 moving particles. Preliminary results suggest that average shadow length scales roughly as 1=N, i.e., shadow lengths decrease rapidly as the number of phase-space dimensions of the system is increased. Some measures of simulation error for N-body systems are discussed that are less stringent than shadowing. Many areas of further research are discussed both for high-dimensional shadowing, and for N-body measures of error. -5 Acknowledgements I thank my supervisor, Prof. Ken Jackson, for many insightful conversations during the devel- opment of the ideas in this thesis, and for helpful comments on the text. My second reader, Dr.
    [Show full text]
  • Chebfun Guide 1St Edition for Chebfun Version 5 ! ! ! ! ! ! ! ! ! ! ! ! ! ! Edited By: Tobin A
    ! ! ! Chebfun Guide 1st Edition For Chebfun version 5 ! ! ! ! ! ! ! ! ! ! ! ! ! ! Edited by: Tobin A. Driscoll, Nicholas Hale, and Lloyd N. Trefethen ! ! ! Copyright 2014 by Tobin A. Driscoll, Nicholas Hale, and Lloyd N. Trefethen. All rights reserved. !For information, write to [email protected].! ! ! MATLAB is a registered trademark of The MathWorks, Inc. For MATLAB product information, please write to [email protected]. ! ! ! ! ! ! ! ! Dedicated to the Chebfun developers of the past, present, and future. ! Table of Contents! ! !Preface! !Part I: Functions of one variable! 1. Getting started with Chebfun Lloyd N. Trefethen! 2. Integration and differentiation Lloyd N. Trefethen! 3. Rootfinding and minima and maxima Lloyd N. Trefethen! 4. Chebfun and approximation theory Lloyd N. Trefethen! 5. Complex Chebfuns Lloyd N. Trefethen! 6. Quasimatrices and least squares Lloyd N. Trefethen! 7. Linear differential operators and equations Tobin A. Driscoll! 8. Chebfun preferences Lloyd N. Trefethen! 9. Infinite intervals, infinite function values, and singularities Lloyd N. Trefethen! 10. Nonlinear ODEs and chebgui Lloyd N. Trefethen Part II: Functions of two variables (Chebfun2) ! 11. Chebfun2: Getting started Alex Townsend! 12. Integration and differentiation Alex Townsend! 13. Rootfinding and optimisation Alex Townsend! 14. Vector calculus Alex Townsend! 15. 2D surfaces in 3D space Alex Townsend ! ! Preface! ! This guide is an introduction to the use of Chebfun, an open source software package that aims to provide “numerical computing with functions.” Chebfun extends MATLAB’s inherent facilities with vectors and matrices to functions and operators. For those already familiar with MATLAB, much of what Chebfun does will hopefully feel natural and intuitive. Conversely, for those new to !MATLAB, much of what you learn about Chebfun can be applied within native MATLAB too.
    [Show full text]
  • BIRS Proceedings 2015
    Banff International Research Station Proceedings 2015 Contents Five-day Workshop Reports 1 1 Modern Applications of Complex Variables: Modeling, Theory and Computation (15w5052) 3 2 Random Dynamical Systems and Multiplicative Ergodic Theorems Workshop (15w5059) 12 3 Mathematics of Communications: Sequences, Codes and Designs (15w5139) 24 4 Discrete Geometry and Symmetry (15w5019) 33 5 Advances in Numerical Optimal Transportation (15w5067) 41 6 Hypercontractivity and Log Sobolev inequalities in Quantum Information Theory (15w5098) 48 7 Computability in Geometry and Analysis (15w5005) 53 8 Laplacians and Heat Kernels: Theory and Applications (15w5110) 60 9 Perspectives on Parabolic Points in Holomorphic Dynamics (15w5082) 71 10 Geometric Flows: Recent Developments and Applications (15w5148) 83 11 New Perspectives for Relational Learning (15w5080) 93 12 Stochasticity and Organization of Tropical Convection (15w5023) 99 13 Groups and Geometries (15w5017) 119 14 Dispersive Hydrodynamics: the Mathematics of Dispersive Shock Waves and Applications (15w5045)129 15 Applied Probability Frontiers: Computational and Modeling Challenges (15w5160) 140 16 Advances and Challenges in ProteinRNA: Recognition, Regulation and Prediction (15w5063) 148 17 Hybrid Methods of Imaging (15w5012) 158 18 Advances in Combinatorial and Geometric Rigidity (15w5114) 175 19 Lifting Problems and Galois Theory (15w5035) 184 20 New Trends in Nonlinear Elliptic Equations (15w5004) 193 21 The Use of Linear Algebraic Groups in Geometry and Number Theory (15w5016) 202 i ii 22
    [Show full text]
  • January 2014
    LONDONLONDON MATHEMATICALMATHEMATICAL SOCIETYSOCIETY NEWSLETTER No. 432 January 2014 Society MeetingsSociety RETIRING MEMBERS OF COUNCIL andMeetings Events AND COMMITTEE CHAIRS and Events GRAEME SEGAL (President) 2014 After serving for two years, Dr introduced electronic voting for Graeme Segal, FRS, handed over elections, which has seen turnout Friday 28 February the LMS Presidency at the AGM more than double. Mary Cartwright in November 2013. Prior to being The past two years have been Lecture, York President, he had served for a challenging period for math- page 27 two years as a Member-at-Large ematical sciences research policy Monday 31 March of Council. Dr Segal is held in in the UK. Dr Segal has led the Northern Regional high esteem across the whole Society in dealing with a range of Meeting, Durham mathematics community and important issues: the introduc- 1 his standing, as well as his sig- tion of open access publishing Tuesday 8 April nificant efforts, have been enor- to which the Society made sig- Society Meeting mously valuable to the Society in nificant policy submissions, the at BMC, London meeting its objectives. issues surrounding impact and page 22 During his tenure as LMS the 2014 REF, the consequences 14–17 April President he has continued the of changes in science funding and Invited Lectures, work of previous presidents university student tuition fees, University of to increase member engage- and the EPSRC Shaping Capabil- East Anglia ment with the Society and to ity agenda. Dr Segal also worked page 28 modernise the way the Society closely with the Education Friday 25 April operates.
    [Show full text]
  • New Directions in Numerical Computation
    COMMUNICATION New Directions in Numerical Computation Tobin A. Driscoll, Endre Süli, and Alex Townsend, Editors In August 2015 a distinguished collec- use of smooth func- tion of numerical analysts gathered at tions in two and higher Oxford to celebrate the sixtieth birth- dimensions was, until day of Nick Trefethen FRS and consider recently, limited to rela- the future of numerical analysis. Some tively simple domains. of the plenary speakers provided short Fundamental research essays for Notices. The full collection is on infinitely smooth online.1 two-dimensional inter- In a 1992 essay, “The Definition of polants may lead to Numerical Analysis”,2 Trefethen writes interesting new ap- of the field, “[O]ur central mission is proaches, yet the number of scientists to compute quantities that are typically working on them ap- uncomputable, from an analytical point Photo courtesy of Burchard Kaup. Nick Trefethen pears to be surprisingly of view, and to do it at lightning speed.” Jean-Paul Berrut small. These essays explore a few of the particulars of that mission. Bengt Fornberg Jean-Paul Berrut The Method—Not the Machine It used to be said that improvements in sci- A Puzzling Question about Numerical entific computing capabilities originate about Analysis equally from advances in algorithms and in hard- Why are so many academic mathematicians con- ware. During the last decade or two, the focus tent with piecewise smooth approximations to has shifted to building heroic-scale supercom- solutions to functional equations when those so- puting facilities. Maybe this is partly because lutions are known a priori to be smooth? Chebfun processing characteristics are quantifiable and impressively and beautifully demonstrates the easy to show in lists, with national and world records falling incessantly.
    [Show full text]
  • Short of Proof: How Many Digits Are Nonetheless Correct?
    Short of Proof: How Many Digits are Nonetheless Correct? Folkmar Bornemann Warm-Up Quiz >> exp(pi*sqrt(67)/6) - sqrt(5280) ans = 6.12120487630818e-08 How many digits are correct? No tools allowed . ICMS 2018, NOTRE DAME,JULY 26 FOLKMAR BORNEMANN 1 CAVEAT At its highest level, numerical analysis is a mixture of science, art and bar-room brawl. —T.W. Körner: The Pleasures of Counting ICMS 2018, NOTRE DAME,JULY 26 FOLKMAR BORNEMANN 2 THE TITLE-QUESTION IS A WEE BIT OLDFASHIONED, ISN’TIT? ultimately, numerical analysis is about efficient algorithms/software . so, answering the title-question is left to the discretion of the educated user? prevalent attitudes towards the accuracy of numerical output level 0: scared! (errors all over the place: roundoff, approximation, . ) level 1: far more accurate than ever needed in practice (lots of plots) level 2: as good as it gets in the realm of noisy data (a.k.a. backward stability) level 3: adaptivity steered by accuracy requirement (“user tolerance”) in short 0: pessimistic; 1: optimistic; 2: blame the data; 3: gimme steering wheel ICMS 2018, NOTRE DAME,JULY 26 FOLKMAR BORNEMANN 3 Prologue Setting the Stage The mathematical proposition has been given the stamp of incontestability: “dispute other things; this is immovable.” —L. Wittgenstein: On Certainty ICMS 2018, NOTRE DAME,JULY 26 FOLKMAR BORNEMANN 4 POINT OF ENTRY I: THE SIAM 100-DIGIT CHALLENGE prehistory: Problem Solving Squad @ Stanford George Pólya George Forsythe Bob Floyd Don Knuth Nick Trefethen ! ! ! ! the challenge (Nick Trefethen, SIAM News & Science, January ’02) 10 problems, the answer to each is a real number sought to 10 sig.
    [Show full text]
  • Master of Technology (M.Tech.) in Modeling and Simulation
    Centre for Modeling and Simulation university of pune • pune 411 007 • india cms.unipune.ac.in +91.(20).2560.1448 • +91.(20).2569.0842 offi[email protected] Master of Technology (M.Tech.) in Modeling and Simulation Approved by the University of Pune Board of Studies: Scientific Computing July 12, 2007 • Version 3.0 Minor Corrections: July 19, 2011 2 About This Document This document outlines an academic programme, the Master of Technology (M.Tech.) in Modeling and Simulation, designed by the people at the Centre for Modeling and Simulation, University of Pune. This academic programme has been approved by the University of Pune. Feedback Requested. The utility of modeling and simulation as a methodology is extensive, and the community that uses it, both from academics and from industry, is very diverse. We would thus highly appreciate your feedback and suggestions on any aspect of this programme { We believe that your feedback will help us make this programme better. Citing This Document This is a public document available at http://cms.unipune.ac.in/reports/. Complete citation for this document is as follows: Mihir Arjunwadkar, Abhay Parvate, Sukratu Barve, P. M. Gade, D. G. Kanhere, Sheelan Chowdhary, Deepak Bankar, and Anjali Kshirsagar Master of Technology (M.Tech.) in Modeling and Simulation. Public Document CMS-PD-20120305 of the Centre for Modeling and Simulation, University of Pune, Pune 411 007, INDIA (2012); available at http://cms.unipune.ac.in/reports/. Acknowledgments Many people have contributed to the creation of this this programme and this document. Their time and effort, contributions and insights, and comments and suggestions have all gone a long way towards bringing this programme to its present form.
    [Show full text]
  • De Morgan Association Newsletter
    UCL DEPARTMENT OF MATHEMATICS De Morgan Association Newsletter from the Department of Mathematics UCL Issue - 20 - 2012 Editor - Michael O'Neill DE MORGAN ASSOCIATION DINNER Wednesday 6 June 2012 The annual dinner of the De Morgan Association was held on Wednesday 6 June 2012 in the Jeremy Bentham Room at UCL, preceded by sherry in the Haldane Room. This year the Guest of Honour was Professor Jeremy Gray, Professor of the History of Mathematics, Open University and Teaching Fellow at UCL. Professor Gray’s address highlighted the role played by UCL, or the University of London as it was when it was founded in 1826, in opening up university education to gifted students, whatever their background. It is hard to believe that the creation of the University of London, with its universal religious tolerance, was opposed by the universities of Oxford and Cambridge as well as the Church of England and even Parliament. This edition of the Newsletter is issue number 20 and it marks 10 years since I took over as Editor from Patricia Rothman, who had edited the Newsletter since its introduction 10 years earlier when David Larman was Head of Department. I therefore thought that it was now an opportune time for me to pass the torch to a new generation, and I am delighted that my successor as Editor will be Professor Ted Johnson. Professor Jeremy Gray speaking at the The past 10 years have been eventful years for the De Morgan Association Dinner Department of Mathematics, in common with all university departments in the country, and it has been interesting to observe how the changing face of the Department has been reflected in the articles contributed to the Newsletter.
    [Show full text]