Computation of operators in wavelet coordinates
|
|
- Matthew Moris Cooper
- 5 years ago
- Views:
Transcription
1 Computation of operators in wavelet coordinates Tsogtgerel Gantumur and Rob Stevenson Department of Mathematics Utrecht University Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.1/22
2 Overview Linear operator equation Lu = g with L : H H Riesz basis Ψ = {ψ λ } of H, e.g. u = λ u λψ λ Infinite dimensional matrix-vector system Lu = g, with u = (u λ ) λ and L : l 2 l 2 Convergent iterations such as u (i+1) = u (i) + α[g Lu (i) ] We can approximate Lu (i) by a finitely supported vector How cheap can we compute this approximation? The answer will depend on L and Ψ Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.2/22
3 Linear operator equations Let Ω be an n-dimensional domain or smooth manifold H t H t (Ω) be a subspace, and H t be its dual space Consider the problem of finding u from Lu = g where L : H t H t is a self-adjoint elliptic operator of order 2t and g H t is a linear functional Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.3/22
4 Differential operators Partial differential operators of order 2t v, Lu = α v, a αβ β u, α, β t Example: The reaction-diffusion equation (t = 1) v, Lu = v u + κ 2 vu, Ω Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.4/22
5 Singular integral operators Boundary integral operators [Lu](x) = K(x, y)u(y)dω y with the kernel K(x, y) singular at x = y Ω Example: The single layer operator for the Laplace BVP in 3-d domain (t = 1 2 ) K(x, y) = 1 4π x y Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.5/22
6 Multiresolution analysis S 0 S 1... H t and S0 S 1... H t dim S j, dim S j = O(2 jn ) (dyadic refinements) S j contains all piecewise pols of degree d 1 S j contains all piecewise pols of degree d 1 S j is globally C r -smooth Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.6/22
7 Wavelet bases Ψ = {ψ λ : λ Λ} is a Riesz basis for H t each v H t has a unique expansion v = λ Λ v λ ψ λ s.t. c v v H t C v For every index λ Λ, there is a number λ IN 0 called the level of ψ λ span{ψ λ : λ j} = S j ψ λ, v = 0 for any v S λ 1 diam(supp ψ λ ) = O(2 λ ) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.7/22
8 Typical wavelets ψ λ ψ µ x ψ λ is a piecewise polynomial of degree d 1 x k ψ λ (x)dx = 0 for k < d ( d vanishing moments) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.8/22
9 Galerkin methods Wavelet basis Ψ j := {ψ λ : λ j} of S j Stiffness L (j) = Lψ λ, ψ µ λ, µ j load g (j) = g, ψ λ λ j Linear equation in IR N j (N j := dim S j ) L (j) u (j) = g (j) L (j) : IR N j IR N j SPD and g (j) IR N j u (j) = λ [u (j)] λ ψ λ approximates the solution of Lu = g Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.9/22
10 Galerkin approximation If u H s for some s [t, d] N j = dim S j = O(2 jn ) ε (j) O(N s t n j ) ε (j) := u (j) u H t O(2 j(s t) ) Solve L (j) u (j) = g (j) with CG complexity O(N j ) Similar estimates for FEM Better convergence? Adaptive methods? Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.10/22
11 Nonlinear approximation Given u = (u λ ) λ l 2 Approximate u using N coeffs u λ λ Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.11/22
12 Nonlinear approximation Given u = (u λ ) λ l 2 Approximate u using N coeffs u λ λ (arranged) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.11/22
13 Nonlinear approximation Given u = (u λ ) λ l 2 Approximate u using N coeffs u λ } {{ } N biggest u N best approximation of u with #supp u N N Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.11/22
14 Nonlinear approximation Given u = (u λ ) λ l 2 Approximate u using N coeffs [u N u] λ N zeroes u N best approximation of u with #supp u N N Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.11/22
15 Nonlinear vs. linear approximation If u B s τ,τ with 1 τ = s t n for some s < d ε N = u N u O(N s t n ) If u H s for some s d, uniform refinement ε (j) = u (j) u O(N s t n j ) B s τ,τ is bigger than H s Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.12/22
16 Besov vs. Sobolev regularity s d 1 τ = s t n B s τ,τ t τ [Dahlke, DeVore]: u B d τ,τ with 1 τ = d t n "often" Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.13/22
17 Equivalent problem in l 2 [Cohen, Dahmen, DeVore 02] Wavelet basis Ψ = {ψ λ : λ Λ} Stiffness L = Lψ λ, ψ µ λ,µ and load g = g, ψ λ λ Linear equation in l 2 (Λ) Lu = g L : l 2 (Λ) l 2 (Λ) SPD and g l 2 (Λ) u = λ u λψ λ is the solution of Lu = g Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.14/22
18 Richardson iterations in l 2 u (0) = 0 u (i+1) = u (i) + α[g Lu (i) ] i = 0, 1,... g and Lu (i) are infinitely supported Approximate them by finitely supported sequences Algorithm SOLVE[N, L, g] u [N] (N operations) #supp u [N] O(N) and ε [N] = u [N] u 0 as N ε [N] speed of convergence? Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.15/22
19 Complexity of SOLVE Matrix L is called q -computable, when for each N one can construct an infinite matrix L N s.t. for any q < q, L N L O(N q ) having in each column O(N) non-zero entries whose computation takes O(N) operations [CDD 02]: Suppose that u N u O(N s ) [s < d t n ] L is q -computable with q > s then for suitable g, u [N] = SOLVE[N, L, g] satisfies u [N] u O(N s ) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.16/22
20 Computability [L λ,µ ] λ Λ µ-th column L λ,µ Approximate by N entries? λ Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.17/22
21 Computability [L λ,µ ] λ Λ µ-th column arranged by modulus L λ,µ N biggest entries? λ (arranged) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.17/22
22 Computability [L λ,µ ] λ Λ µ-th column L λ,µ }{{} N Compute the N biggest entries Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.17/22
23 Computability The µ-th column of the difference [L N L] λ,µ } {{ } N Need to locate the biggest entries a priori Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.17/22
24 Compressibility L is called q -compressible, when L is q -computable assuming each entry of L is available at unit cost [CDD 01], [Stevenson 04]: Suppose {ψ λ } are piecewise polynomial wavelets that are sufficiently smooth and have sufficiently many vanishing moments L is either differential or singular integral operator then L is q -compressible for some q d t n (> s) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.18/22
25 Computability Distribute computational works over the entries W Require: shaded area = O(N) N Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.19/22
26 Computability Distribute computational works over the entries N θ p(x) = N θ x ϱ when θ ϱ < 1 N N 0 N θ x ϱ dx = 1 1 ϱ N 1+θ ϱ = O(N) Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.20/22
27 Computability [T.G., Stevenson 04], [T.G. 04]: Suppose {ψ λ } are piecewise polynomial wavelets that are sufficiently smooth and have sufficiently many vanishing moments L is either differential or singular integral operator then L is q -computable for some q d t n (> s) So the adaptive wavelet method has the optimal convergence rate and optimal computational complexity Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.21/22
28 References A. Cohen, W. Dahmen, and R. DeVore. Adaptive wavelet methods II - Beyond the elliptic case. Found. Comput. Math., 2(3): , R.P. Stevenson. On the compressibility of operators in wavelet coordinates. SIAM J. Math. Anal., 35(5): , T. Gantumur and R.P. Stevenson. Computation of differential operators in wavelet coordinates. Technical Report 1306, Utrecht University, August Submitted. T. Gantumur. Computation of singular integral operators in wavelet coordinates. Technical Report, Utrecht University, To appear. This work was supported by the Netherlands Organization for Scientific Research and by the EC-IHP project Breaking Complexity. Tsogtgerel Gantumur - Computation of operators in wavelet coordinates - Sixth Minisimposium TIANA. Amsterdam. Oct p.22/22
Adaptive Wavelet Algorithms
Adaptive Wavelet Algorithms for solving operator equations Tsogtgerel Gantumur August 2006 Contents Notations and acronyms vii 1 Introduction 1 1.1 Background.............................. 1 1.2 Thesis
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationSpace-Time Adaptive Wavelet Methods for Parabolic Evolution Problems
Space-Time Adaptive Wavelet Methods for Parabolic Evolution Problems Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam joint work with Christoph Schwab (ETH, Zürich), Nabi
More informationBesov regularity for the solution of the Stokes system in polyhedral cones
the Stokes system Department of Mathematics and Computer Science Philipps-University Marburg Summer School New Trends and Directions in Harmonic Analysis, Fractional Operator Theory, and Image Analysis
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationBesov regularity for operator equations on patchwise smooth manifolds
on patchwise smooth manifolds Markus Weimar Philipps-University Marburg Joint work with Stephan Dahlke (PU Marburg) Mecklenburger Workshop Approximationsmethoden und schnelle Algorithmen Hasenwinkel, March
More informationAn optimal adaptive finite element method. Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam
An optimal adaptive finite element method Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam Contents model problem + (A)FEM newest vertex bisection convergence
More informationHow much should we rely on Besov spaces as a framework for the mathematical study of images?
How much should we rely on Besov spaces as a framework for the mathematical study of images? C. Sinan Güntürk Princeton University, Program in Applied and Computational Mathematics Abstract Relations between
More informationQuarkonial frames of wavelet type - Stability, approximation and compression properties
Quarkonial frames of wavelet type - Stability, approximation and compression properties Stephan Dahlke 1 Peter Oswald 2 Thorsten Raasch 3 ESI Workshop Wavelet methods in scientific computing Vienna, November
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
More informationAn Adaptive Wavelet Viscosity Method for Hyperbolic Conservation Laws
An Adaptive Wavelet Viscosity Method for Hyperbolic Conservation Laws Daniel Castaño Díez, 1 Max Gunzburger, 2 Angela Kunoth 3 1 EMBL Heidelberg, Meyerhofstraße 1, 69117 Heidelberg, Germany 2 School of
More informationDISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS
DISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS D. Černá, V. Finěk Department of Mathematics and Didactics of Mathematics, Technical University in Liberec Abstract Wavelets and a discrete
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More informationBesov Regularity and Approximation of a Certain Class of Random Fields
Besov Regularity and Approximation of a Certain Class of Random Fields Nicolas Döhring (TU Kaiserslautern) joint work with F. Lindner, R. Schilling (TU Dresden), K. Ritter (TU Kaiserslautern), T. Raasch
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationAdaptive Wavelet BEM for Boundary Integral Equations: Theory and Numerical Experiments
Adaptive Wavelet BEM for Boundary Integral Equations: Theory and Numerical Experiments S. Dahlke, H. Harbrecht, M. Utzinger, M. Weimar October 7, 06 Abstract In this paper, we are concerned with the numerical
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationrecent developments of approximation theory and greedy algorithms
recent developments of approximation theory and greedy algorithms Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Reduced Order Modeling in
More informationThe adaptive tensor product wavelet scheme: Sparse matrices and the application to singularly perturbed problems
IMA Journal of Numerical Analysis () Page of 8 doi:.93/imanum/drn The adaptive tensor product wavelet scheme: Sparse matrices and the application to singularly perturbed problems NABI CHEGINI AND ROB STEVENSON
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationBesov regularity of solutions of the p-laplace equation
Besov regularity of solutions of the p-laplace equation Benjamin Scharf Technische Universität München, Department of Mathematics, Applied Numerical Analysis benjamin.scharf@ma.tum.de joint work with Lars
More informationLecture 3: Initial and boundary value problems in 1 dimension
Lecture 3: Initial and boundary value problems in 1 dimension Gantumur Tsogtgerel Assistant professor of Mathematics Math 319: Introduction to PDEs McGill University, Montréal January 10, 2011 Euler s
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationAdaptive Wavelet Methods for Elliptic SPDEs
Adaptive Wavelet Methods for Elliptic SPDEs Klaus Ritter Computational Stochastics TU Kaiserslautern 1/1 Introduction In most papers on approximation of SPDEs non-adaptive discretization discretization
More informationNEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS
NEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS YOUNGMI HUR AND AMOS RON Abstract. The classical Haar wavelet system of L 2 (R n ) is commonly considered to be very local in space. We introduce and study
More informationAdaptive Wavelet Methods on Unbounded Domains
Adaptive Wavelet Methods on Unbounded Domains Sebastian Kestler und Karsten Urban Preprint Series: 2-4 Fakultät für Mathematik und Wirtschaftswissenschaften UNIVERSITÄT ULM ADAPTIVE WAVELET METHODS ON
More informationADAPTIVE WAVELET SCHEMES FOR PARABOLIC PROBLEMS: SPARSE MATRICES AND NUMERICAL RESULTS
ADAPTVE WAVELET SCHEMES FOR PARABOLC PROBLEMS: SPARSE MATRCES AND NUMERCAL RESULTS NAB CHEGN AND ROB STEVENSON Abstract. A simultaneous space-time variational formulation of a parabolic evolution problem
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationHow smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron
How smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron A closed subspace V of L 2 := L 2 (IR d ) is called PSI (principal shift-invariant) if it is the
More informationDFG-Schwerpunktprogramm 1324
DFG-Schwerpunktprogramm 1324 Extraktion quantifizierbarer Information aus komplexen Systemen Adaptive Wavelet Methods for Elliptic Stochastic Partial Differential Equations P. A. Cioica, S. Dahlke, N.
More informationNonlinear and Adaptive Frame Approximation Schemes for Elliptic PDEs: Theory and Numerical Experiments
Nonlinear and Adaptive Frame Approximation Schemes for Elliptic PDEs: Theory and Numerical Experiments Stephan Dahlke, Massimo Fornasier, Miriam Primbs, Thorsten Raasch, and Manuel Werner Abstract This
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationNonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames
Nonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames Martin Ehler September 26, 2008 Abstract In the present paper, we study nonlinear approximation properties of multivariate
More informationOn the Numerical Evaluation of Fractional Sobolev Norms. Carsten Burstedde. no. 268
On the Numerical Evaluation of Fractional Sobolev Norms Carsten Burstedde no. 268 Diese Arbeit ist mit Unterstützung des von der Deutschen Forschungsgemeinschaft getragenen Sonderforschungsbereiches 611
More informationINSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD
INSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD () Instanton (definition) (2) ADHM construction (3) Compactification. Instantons.. Notation. Throughout this talk, we will use the following notation:
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationWavelet bases on a triangle
Wavelet bases on a triangle N. Ajmi, A. Jouini and P.G. Lemarié-Rieusset Abdellatif.jouini@fst.rnu.tn neyla.ajmi@laposte.net Département de Mathématiques Faculté des Sciences de Tunis 1060 Campus-Tunis,
More informationRate optimal adaptive FEM with inexact solver for strongly monotone operators
ASC Report No. 19/2016 Rate optimal adaptive FEM with inexact solver for strongly monotone operators G. Gantner, A. Haberl, D. Praetorius, B. Stiftner Institute for Analysis and Scientific Computing Vienna
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationUNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES. Tight compactly supported wavelet frames of arbitrarily high smoothness
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Tight compactly supported wavelet frames of arbitrarily high smoothness Karlheinz Gröchenig Amos Ron Department of Mathematics U-9 University
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods The Connection to Green s Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationMultilevel Preconditioning and Adaptive Sparse Solution of Inverse Problems
Multilevel and Adaptive Sparse of Inverse Problems Fachbereich Mathematik und Informatik Philipps Universität Marburg Workshop Sparsity and Computation, Bonn, 7. 11.6.2010 (joint work with M. Fornasier
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationSparse Wavelet Representation of Differential Operators with Piecewise Polynomial Coefficients
Article Sparse Wavelet Representation of Differential Operators with Piecewise Polynomial Coefficients Dana Černá * and Václav Finěk Department of Mathematics and Didactics of Mathematics, Technical University
More informationFrames. Hongkai Xiong 熊红凯 Department of Electronic Engineering Shanghai Jiao Tong University
Frames Hongkai Xiong 熊红凯 http://ivm.sjtu.edu.cn Department of Electronic Engineering Shanghai Jiao Tong University 2/39 Frames 1 2 3 Frames and Riesz Bases Translation-Invariant Dyadic Wavelet Transform
More information7.3 The Jacobi and Gauss-Siedel Iterative Techniques. Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP.
7.3 The Jacobi and Gauss-Siedel Iterative Techniques Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP. 7.3 The Jacobi and Gauss-Siedel Iterative Techniques
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationECE 901 Lecture 16: Wavelet Approximation Theory
ECE 91 Lecture 16: Wavelet Approximation Theory R. Nowak 5/17/29 1 Introduction In Lecture 4 and 15, we investigated the problem of denoising a smooth signal in additive white noise. In Lecture 4, we considered
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 27, No. 4, pp. 1347 1370 c 2006 Society for Industrial and Applied Mathematics WAVELET GALERKIN SCHEMES FOR BOUNDARY INTEGRAL EQUATIONS IMPLEMENTATION AND QUADRATURE HELMUT HARBRECHT
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
More informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig H 2 -matrix approximation of integral operators by interpolation by Wolfgang Hackbusch and Steen Borm Preprint no.: 04 200 H 2 -Matrix
More informationSMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN
SMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN Abstract. In [B. Han and Z. Shen, SIAM J. Math. Anal., 38 (006), 530 556], a family of univariate short support Riesz
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationOn the Computations of Eigenvalues of the Fourth-order Sturm Liouville Problems
Int. J. Open Problems Compt. Math., Vol. 4, No. 3, September 2011 ISSN 1998-6262; Copyright c ICSRS Publication, 2011 www.i-csrs.org On the Computations of Eigenvalues of the Fourth-order Sturm Liouville
More informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
More informationS.Novikov. Singular Solitons and Spectral Theory
S.Novikov Singular Solitons and Spectral Theory Moscow, August 2014 Collaborators: P.Grinevich References: Novikov s Homepage www.mi.ras.ru/ snovikov click Publications, items 175,176,182, 184. New Results
More informationAlmost diagonal matrices and Besov-type spaces based on wavelet expansions
arxiv:1408.4581v1 [math.fa] 20 Aug 2014 Almost diagonal matrices and Besov-type spaces based on wavelet expansions Markus Weimar August 21, 2014 Abstract This paper is concerned with problems in the context
More informationFOUNDATIONS OF COMPUTATIONAL MATHEMATICS The Journal of the Society for the Foundations of Computational Mathematics
Found. Comput. Math. (2002) 2:203 245 2002 SFoCM DOI: 10.1007/s102080010027 FOUNDATIONS OF COMPUTATIONAL MATHEMATICS The Journal of the Society for the Foundations of Computational Mathematics Adaptive
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationIII. Iterative Monte Carlo Methods for Linear Equations
III. Iterative Monte Carlo Methods for Linear Equations In general, Monte Carlo numerical algorithms may be divided into two classes direct algorithms and iterative algorithms. The direct algorithms provide
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationNumerical Solution of a Coefficient Identification Problem in the Poisson equation
Swiss Federal Institute of Technology Zurich Seminar for Applied Mathematics Department of Mathematics Bachelor Thesis Spring semester 2014 Thomas Haener Numerical Solution of a Coefficient Identification
More informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationQualifying Exams I, 2014 Spring
Qualifying Exams I, 2014 Spring 1. (Algebra) Let k = F q be a finite field with q elements. Count the number of monic irreducible polynomials of degree 12 over k. 2. (Algebraic Geometry) (a) Show that
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationLearning gradients: prescriptive models
Department of Statistical Science Institute for Genome Sciences & Policy Department of Computer Science Duke University May 11, 2007 Relevant papers Learning Coordinate Covariances via Gradients. Sayan
More informationOn the fast algorithm for multiplication of functions in the wavelet bases
Published in Proceedings of the International Conference Wavelets and Applications, Toulouse, 1992; Y. Meyer and S. Roques, edt., Editions Frontieres, 1993 On the fast algorithm for multiplication of functions
More informationReduced Modeling in Data Assimilation
Reduced Modeling in Data Assimilation Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Challenges in high dimensional analysis and computation
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationResolving the White Noise Paradox in the Regularisation of Inverse Problems
1 / 32 Resolving the White Noise Paradox in the Regularisation of Inverse Problems Hanne Kekkonen joint work with Matti Lassas and Samuli Siltanen Department of Mathematics and Statistics University of
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More information