H 2 -matrices with adaptive bases
|
|
- Arleen Carol Jennings
- 5 years ago
- Views:
Transcription
1 1 H 2 -matrices with adaptive bases Steffen Börm MPI für Mathematik in den Naturwissenschaften Inselstraße 22 26, Leipzig
2 Problem 2 Goal: Treat certain large dense matrices in linear complexity. Model problem: Discretization of integral operators of the form K[u](x) = κ(x, y)u(y) dy. Ω Discretization: Galerkin method with FE basis (ϕ i ) n i=1 leads to a matrix K R n n with K ij = ϕ i, K[ϕ j ] L 2 = ϕ i (x)κ(x, y)ϕ j (y) dy dx. (1 i, j n) Ω Ω Problem: K is dense. Idea: Use properties of the kernel function to approximate K. Examples: Wavelets, multipole expansion, interpolation.
3 Interpolation 3 Idea: We consider a subdomain τ σ Ω Ω and replace the kernel function κ by an interpolant κ τ,σ (x, y) := k k κ(x τ ν, x σ µ)l τ ν(x)l σ µ(y). (k n) ν=1 µ=1 Result: Local matrix in factorized form K τ,σ ij := ϕ i (x)κ(x, y)ϕ j (y) dy dx = τ k σ k ν=1 µ=1 κ(x τ ν, x σ µ) } {{ } =:Sνµ τ,σ = (V τ S τ,σ (V σ ) ) ij τ ϕ i (x)l τ ν(x) dx τ } {{ } =:Viν τ σ ϕ i (x) κ τ,σ (x, y)ϕ j (y) dy dx ϕ j (y)l σ µ(y) dy σ } {{ } =:Vjµ σ n k n k
4 Local interpolation 4 Problem: Typical kernel functions are not globally smooth. Idea: They are locally smooth outside of the diagonal x = y. Examples: 1/ x y, log x y. Error bound: For m-th order tensor interpolation on axis-parallel boxes B τ τ and B σ σ, we get ( κ τ,σ κ L (B τ B σ ) C apx (m) dist(bτ, B σ ) m ) diam(b τ B σ. ) Rank: m interpolation points per coordinate k = m d. Admissibility: Uniform exponential convergence requires diam(b τ B σ ) 2η dist(b τ, B σ ).
5 Block structure 5 Approach: Split Ω Ω into admissible subdomains and a small remainder. Cluster tree: Hierarchy of subdomains of Ω Nearfield: Small blocks, stored in standard format. Farfield: Admissible blocks, stored in V τ S τ,σ (V σ ) format. Partition: Collection of farand nearfield blocks. Construction of cluster tree and block partition can be accomplished by general algorithms in O(n) operations.
6 Nested bases 6 Goal: Treat the cluster basis (V τ ) τ T efficiently. Idea: For a cluster τ and τ sons(τ) we have L τ ν = k ν =1 L τ ν(x τ ν )Lτ. (1 ν k) ν Nested bases: Consider i τ τ. Setting T τ ν ν := Lτ ν(x τ ν ), we find k V τ iν = τ ϕ i (x)l τ ν(x) dx = ν =1 T τ ν ν τ ϕ i (x)l τ ν (x) dx = (V τ T τ ) iν. Consequence: Store V τ only for leaves and use T τ for all other clusters. Complexity: Storage O(nk), matrixvector O(nk).
7 Orthogonalization 7 Problem: Efficiency depends on the rank k. Can we reduce this rank by eliminating redundant expansion functions? Idea: Orthogonalize columns of V τ. Gram-Schmidt: Find Ṽ τ = V τ Z τ such that (Ṽ τ ) Ṽ τ = I holds. Nested structure: Using P τ with V τ = Ṽ τ P τ, we find V τ = V τ 1 T τ 1 = Ṽ τ 1 P τ 1 T τ 1 = Ṽ τ 1 V τ 2 T τ 2 Ṽ τ 2 P τ 2 T τ 2 Ṽ τ 2 = Ṽ τ V τ Z τ P τ 1 T τ 1, P τ 2 T τ 2 so it is sufficient to orthogonalize a (2k) k matrix. The resulting cluster basis is nested due to Ṽ τ = V τ Z τ = Ṽ τ 1 (P τ 1 T τ 1 Z τ ) = Ṽ τ 1 T τ 1. Ṽ τ 2 (P τ 2 T τ 2 Z τ ) Ṽ τ 2 τ T 2 Result: Complexity O(nk 2 ).
8 Orthogonalization: Simple example 8 Example: Unit sphere, piecewise constant trial and test functions. Interpolation Orthogonalized n Build MVM Mem/n Expl Impl MVM Mem/n H 2 -matrix constructed by cubic interpolation and subsequent orthogonalization, relative error Hardware: SunFire 6800 with 900 MHz UltraSparc IIIcu processors Software: HLib, see
9 Recompression algorithm: Problem 9 Problem: Interpolation does not take into account that a subset of polynomials may be sufficient to approximate the kernel function. Idea: Use local singular value decompositions to find optimal cluster bases Take kernel and geometry into account. Important: Since the cluster bases are nested, the father clusters depend on the son clusters. Consequence: If τ σ is a farfield block and τ sons(τ), we have K τ σ = (V τ S τ,σ (V σ ) ) τ σ = V τ T τ S τ,σ (V σ ), i.e., the choice of V τ its other ancestors. influences all farfield blocks connected to τ or
10 Recompression algorithm: Local problems and recursion 10 Block row: R τ := {σ : τ + τ : τ + σ is admissible}. Goal: Find an orthogonal matrix V τ R τ kτ that maximizes σ R τ (V τ ) K τ σ 2 F. Optimal solution: Use eigenvectors of corresponding Gram matrix. Nested structure: Ensured by bottom-up recursion and projections. Result: Algorithm with complexity O(nk 2 ), truncation error can be controlled directly.
11 Recompression: Simple example 11 Example: Unit sphere, piecewise constant trial and test functions. Interpolation Recompressed n Build MVM Mem/n Build MVM Mem/n H 2 -matrix constructed by cubic interpolation. Recompression with relative error bound 10 3, maximal rank 16. Hardware: SunFire 6800 with 900 MHz UltraSparc IIIcu processors. Software: HLib, see
12 Recompression: Netgen example 12 Example: Grid Crank shaft of the Netgen package. n Build Mem/n Near/n MVM Rel. error 0.05% 0.07% 0.08% DMem/n
13 Recompression: Netgen geometries and eddy current model 13 Z Z b(v, φ) = Γ Z φ(y), v(x) x Φ(x, y), n(x) dy dx Z φ(y),n(x) x Φ(x, y), v(x) dy dx Γ Γ Γ n Direct MVM Mem/DoF Recompressed MVM Mem/DoF Rel. error
14 Summary 14 H 2 -matrices: Nested structure leads to complexity O(nk) for storage requirements and matrix-vector multiplication. Black box: The approximation is constructed by using only kernel evaluations κ(x τ ν, x σ µ) and works for all asymptotically smooth kernel functions. Recompression: Use O(nk 2 ) algorithm to reduce the rank k, find a quasi-optimal cluster basis adapted to geometry, discretization and kernel function. Work in progress: H 2 -matrix arithmetics, HCA.
15 H 2 -matrix arithmetics 15 Goal: Perform matrix addition and multiplication efficiently. Addition and scaling: For an admissible block τ σ, we have (A + λb) τ σ = V τ (S τ,σ A + λsτ,σ B )(V σ ), so H 2 -matrices form a vector space. Addition and scaling can be performed in O(nk) operations. ( + λ ) Multiplication: Defined recursively. A 11 A 12 B 11 B 12 = A 11B 11 + A 12 B 21 A 11 B 12 + A 12 B 22 A 21 A 22 B 21 B 22 A 21 B 11 + A 22 B 21 A 21 B 12 + A 22 B 22 Use truncation if product does not match prescribed format.
16 H 2 -matrix multiplication I 16 := Situation: Let A = V τ S τ,ϱ A (V ϱ ), B not admissible. Compute best approximation of AB in the form C = V τ S τ,σ C (V σ ). Orthogonality: If V τ and V σ are orthogonal: S τ,σ C = (V τ ) ABV σ = (V τ ) V τ S τ,ϱ A (V ϱ ) BV σ = S τ,ϱ A (V ϱ ) BV σ. }{{} =: S ϱ,σ B Problem: Computing (V ϱ ) BV σ directly is too expensive. Idea: Prepare S ϱ,σ B := (V ϱ ) BV σ in advance. Matrix forward transformation: Since the cluster bases are nested, all matrices can be prepared in O(nk3 ) operations. S ϱ,σ B
17 H 2 -matrix multiplication II 17 := Situation: Let A = V τ S τ,ϱ A (V ϱ ) and B = V ϱ S ϱ,σ B (V σ ). Advantage: The product AB already has the factorized form AB = V τ S τ,ϱ A (V ϱ ) V ϱ S ϱ,σ B (V σ ) = V τ (S τ,ϱ A Sϱ,σ B )(V σ ). Problem: Converting AB directly to the format of C will lead to inacceptable complexity. S τ,σ C Idea: Store result in multiplication is complete. := S τ,ϱ A Sϱ,σ B and distribute after the Matrix backward transformation: Since the cluster bases are nested, all matrices can be distributed in O(nk3 ) operations. S τ,σ C
18 H 2 -matrix multiplication III 18 := Situation: Let B = V ϱ S ϱ,σ B (V σ ). Approach: Split B and proceed by recursion: B ϱ σ = V ϱ T ϱ S ϱ,σ B (T σ ) (V σ ). }{{} =: bs ϱ,σ B Complexity: Only O(Cspn) 2 such splitting operations are required Complexity O(Cspnk 2 3 ). Result: The best approximation C of AB in a prescribed H 2 -matrix format can be computed in O(C 2 spnk 3 ) operations.
19 H 2 -matrix multiplication: Fixed structure 19 Example: Discretized double layer potential on the unit sphere K, approximate KK. n Build C sp H 2 -MVM H 2 -MMM H-MMM Error < Important: Here, the result matrix is projected to the same cluster bases and the same block structure as K. Question: Can we reduce the error by using adapted cluster bases?
20 H 2 -matrix multiplication: Adapted structure 20 Modification: Choose cluster bases for result matrix adaptively to improve the precision. n C sp Old MMM Error New basis New MMM Error Advantage: Precision significantly improved, computation time only slightly increased. Current work: Improve the efficiency of the construction of adapted cluster bases.
21 Hybrid Cross Approximation 21 Problem: For each admissble block V τ S τ,σ (V σ ) of the H 2 -matrix, the kernel function κ has to be evaluated in m 2d points. Very time-consuming if high precision is required. Approach: Apply adaptive cross approximation to the coefficient matrices S τ,σ to find rank k representation S τ,σ X τ,σ (Y τ,σ ). Result: Only m d k m 2d kernel evaluations required to approximate coefficient matrix S τ,σ. Important: ACA is only applied to coefficient matrices. Construction still linear in n, no problems with quadrature. m Std 24.1s 64.7s 160.8s 399.2s HCA 22.3s 49.2s 90.0s 222.7s
22 Conclusion 22 H 2 -matrices: Data-sparse approximation of certain dense matrices with linear complexity in the number of degrees of freedom. Adaptive cluster bases: Find quasi-optimal expansion systems for arbitrary matrices. In the case of integral operators, point evaluations of the kernel function are sufficient to construct efficient and reliable H 2 -matrix approximations. H 2 -matrix arithmetics: Approximate A + λb and AB in O(nk 2 ) and O(nk 3 ) operations. Next goal: Approximate A 1 and LU. HCA: Reduce the number of kernel evaluations required to construct the H 2 -matrix. Software: HLib package for hierarchical matrices and H 2 -matrices, including applications to elliptic PDEs and integral equations.
23 Interpolation vs. Multipole 23 Construction: Pointwise evaluation and integrals of polynomials. Tν t ν = L t ν(x t ν ), S t,s νµ = κ(x t ν, x s µ), Viν t = ϕ i (x)l t ν(x) dx. Interpolation Multipole - Rank k = p d + Rank k = p d 1 - Low frequencies only + Low and high frequencies + General kernels - Special kernels only + Anisotropic clusters - Spherical clusters only Γ Goal: Reduce rank for interpolation without sacrificing precision.
Fast evaluation of boundary integral operators arising from an eddy current problem. MPI MIS Leipzig
1 Fast evaluation of boundary integral operators arising from an eddy current problem Steffen Börm MPI MIS Leipzig Jörg Ostrowski ETH Zürich Induction heating 2 Induction heating: Time harmonic eddy current
More informationfür Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Approximation of integral operators by H 2 -matrices with adaptive bases (revised version: December 2004) by Steffen Börm Preprint
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 informationAPPROXIMATING GAUSSIAN PROCESSES
1 / 23 APPROXIMATING GAUSSIAN PROCESSES WITH H 2 -MATRICES Steffen Börm 1 Jochen Garcke 2 1 Christian-Albrechts-Universität zu Kiel 2 Universität Bonn and Fraunhofer SCAI 2 / 23 OUTLINE 1 GAUSSIAN PROCESSES
More informationHierarchical Matrices. Jon Cockayne April 18, 2017
Hierarchical Matrices Jon Cockayne April 18, 2017 1 Sources Introduction to Hierarchical Matrices with Applications [Börm et al., 2003] 2 Sources Introduction to Hierarchical Matrices with Applications
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationfür Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Introduction to Hierarchical Matrices with Applications by Steffen Börm, Lars Grasedyck, and Wolfgang Hackbusch Preprint no.: 18 2002
More informationHierarchical matrix techniques for low and high frequency Helmholtz problems
Hierarchical matrix techniques for low and high frequency Helmholtz problems Lehel Banjai Wolfgang Hackbusch 12th December 2006 Abstract In this paper, we discuss the application of hierarchical matrix
More informationMath 671: Tensor Train decomposition methods
Math 671: Eduardo Corona 1 1 University of Michigan at Ann Arbor December 8, 2016 Table of Contents 1 Preliminaries and goal 2 Unfolding matrices for tensorized arrays The Tensor Train decomposition 3
More informationFast matrix algebra for dense matrices with rank-deficient off-diagonal blocks
CHAPTER 2 Fast matrix algebra for dense matrices with rank-deficient off-diagonal blocks Chapter summary: The chapter describes techniques for rapidly performing algebraic operations on dense matrices
More informationMax-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Hierarchical matrices and the High-Frequency Fast Multipole Method for the Helmholtz Equation with Decay revised version: March 014)
More informationTHE H 2 -matrix is a general mathematical framework [1],
Accuracy Directly Controlled Fast Direct Solutions of General H -Matrices and ts Application to Electrically Large ntegral-equation-based Electromagnetic Analysis Miaomiao Ma, Student Member, EEE, and
More informationMULTI-LAYER HIERARCHICAL STRUCTURES AND FACTORIZATIONS
MULTI-LAYER HIERARCHICAL STRUCTURES AND FACTORIZATIONS JIANLIN XIA Abstract. We propose multi-layer hierarchically semiseparable MHS structures for the fast factorizations of dense matrices arising from
More informationHybrid Cross Approximation for the Electric Field Integral Equation
Progress In Electromagnetics Research M, Vol. 75, 79 90, 2018 Hybrid Cross Approximation for the Electric Field Integral Equation Priscillia Daquin, Ronan Perrussel, and Jean-René Poirier * Abstract The
More informationAn H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages
PIERS ONLINE, VOL. 6, NO. 7, 2010 679 An H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages Haixin Liu and Dan Jiao School of Electrical
More informationTheory and implementation of H-matrix based iterative and direct solvers for Helmholtz and elastodynamic oscillatory kernels
Theory and implementation of H-matrix based iterative and direct solvers for Helmholtz and elastodynamic oscillatory kernels Stéphanie Chaillat, Luca Desiderio, Patrick Ciarlet To cite this version: Stéphanie
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationCOMPARED to other computational electromagnetic
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 58, NO. 12, DECEMBER 2010 3697 Existence of -Matrix Representations of the Inverse Finite-Element Matrix of Electrodynamic Problems and -Based
More informationLecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:
tifica Lecture 1: Computation of Karhunen-Loeve Expansion: Alexander Litvinenko http://sri-uq.kaust.edu.sa/ Stochastic PDEs We consider div(κ(x, ω) u) = f (x, ω) in G, u = 0 on G, with stochastic coefficients
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationAn Introduction to Hierachical (H ) Rank and TT Rank of Tensors with Examples
An Introduction to Hierachical (H ) Rank and TT Rank of Tensors with Examples Lars Grasedyck and Wolfgang Hackbusch Bericht Nr. 329 August 2011 Key words: MSC: hierarchical Tucker tensor rank tensor approximation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationTeil 2: Hierarchische Matrizen
Teil 2: Hierarchische Matrizen Mario Bebendorf, Universität Leipzig bebendorf@math.uni-leipzig.de Bebendorf, Mario (Univ. Leipzig) Kompaktkurs lineare GS und H-Matrizen 24. 2. April 2 1 / 1 1 1 2 1 2 2
More informationMaster Thesis Literature Review: Efficiency Improvement for Panel Codes. TU Delft MARIN
Master Thesis Literature Review: Efficiency Improvement for Panel Codes TU Delft MARIN Elisa Ang Yun Mei Student Number: 4420888 1-27-2015 CONTENTS 1 Introduction... 1 1.1 Problem Statement... 1 1.2 Current
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationContents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information
Contents Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information xi xiv xvii xix 1 Preliminaries 1 1.0 Introduction.............................
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 informationExperiences with Model Reduction and Interpolation
Experiences with Model Reduction and Interpolation Paul Constantine Stanford University, Sandia National Laboratories Qiqi Wang (MIT) David Gleich (Purdue) Emory University March 7, 2012 Joe Ruthruff (SNL)
More information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
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 671: Tensor Train decomposition methods II
Math 671: Tensor Train decomposition methods II Eduardo Corona 1 1 University of Michigan at Ann Arbor December 13, 2016 Table of Contents 1 What we ve talked about so far: 2 The Tensor Train decomposition
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationAn Adaptive Hierarchical Matrix on Point Iterative Poisson Solver
Malaysian Journal of Mathematical Sciences 10(3): 369 382 (2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal An Adaptive Hierarchical Matrix on Point
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationA New Scheme for the Tensor Representation
J Fourier Anal Appl (2009) 15: 706 722 DOI 10.1007/s00041-009-9094-9 A New Scheme for the Tensor Representation W. Hackbusch S. Kühn Received: 18 December 2008 / Revised: 29 June 2009 / Published online:
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
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 informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationSolving an Elliptic PDE Eigenvalue Problem via Automated Multi-Level Substructuring and Hierarchical Matrices
Solving an Elliptic PDE Eigenvalue Problem via Automated Multi-Level Substructuring and Hierarchical Matrices Peter Gerds and Lars Grasedyck Bericht Nr. 30 März 2014 Key words: automated multi-level substructuring,
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationfür Mathematik in den Naturwissenschaften Leipzig
ŠܹÈÐ Ò ¹ÁÒ Ø ØÙØ für Mathematik in den Naturwissenschaften Leipzig Nonlinear multigrid for the solution of large scale Riccati equations in low-rank and H-matrix format (revised version: May 2007) by
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationLecture 1: Introduction to low-rank tensor representation/approximation. Center for Uncertainty Quantification. Alexander Litvinenko
tifica Lecture 1: Introduction to low-rank tensor representation/approximation Alexander Litvinenko http://sri-uq.kaust.edu.sa/ KAUST Figure : KAUST campus, 5 years old, approx. 7000 people (include 1400
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationTensor-Product Representation of Operators and Functions (7 introductory lectures) Boris N. Khoromskij
1 Everything should be made as simple as possible, but not simpler. A. Einstein (1879-1955) Tensor-Product Representation of Operators and Functions (7 introductory lectures) Boris N. Khoromskij University
More informationFast Direct Methods for Gaussian Processes
Fast Direct Methods for Gaussian Processes Mike O Neil Departments of Mathematics New York University oneil@cims.nyu.edu December 12, 2015 1 Collaborators This is joint work with: Siva Ambikasaran Dan
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit VII Sparse Matrix
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 1: Direct Methods Dianne P. O Leary c 2008
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More information1. Structured representation of high-order tensors revisited. 2. Multi-linear algebra (MLA) with Kronecker-product data.
Lect. 4. Toward MLA in tensor-product formats B. Khoromskij, Leipzig 2007(L4) 1 Contents of Lecture 4 1. Structured representation of high-order tensors revisited. - Tucker model. - Canonical (PARAFAC)
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
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 informationhomogeneous 71 hyperplane 10 hyperplane 34 hyperplane 69 identity map 171 identity map 186 identity map 206 identity matrix 110 identity matrix 45
address 12 adjoint matrix 118 alternating 112 alternating 203 angle 159 angle 33 angle 60 area 120 associative 180 augmented matrix 11 axes 5 Axiom of Choice 153 basis 178 basis 210 basis 74 basis test
More informationThe Fast Multipole Method and other Fast Summation Techniques
The Fast Multipole Method and other Fast Summation Techniques Gunnar Martinsson The University of Colorado at Boulder (The factor of 1/2π is suppressed.) Problem definition: Consider the task of evaluation
More informationSOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS. Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA
1 SOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA 2 OUTLINE Sparse matrix storage format Basic factorization
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Steffen Börm Compiled July 12, 2018, 12:01 All rights reserved. Contents 1. Introduction 5 2. Finite difference methods 7 2.1. Potential equation.............................
More informationApplication of the LLL Algorithm in Sphere Decoding
Application of the LLL Algorithm in Sphere Decoding Sanzheng Qiao Department of Computing and Software McMaster University August 20, 2008 Outline 1 Introduction Application Integer Least Squares 2 Sphere
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationSparse BLAS-3 Reduction
Sparse BLAS-3 Reduction to Banded Upper Triangular (Spar3Bnd) Gary Howell, HPC/OIT NC State University gary howell@ncsu.edu Sparse BLAS-3 Reduction p.1/27 Acknowledgements James Demmel, Gene Golub, Franc
More informationA fast randomized algorithm for computing a Hierarchically Semi-Separable representation of a matrix
A fast randomized algorithm for computing a Hierarchically Semi-Separable representation of a matrix P.G. Martinsson, Department of Applied Mathematics, University of Colorado at Boulder Abstract: Randomized
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationConjugate Gradient Method
Conjugate Gradient Method direct and indirect methods positive definite linear systems Krylov sequence spectral analysis of Krylov sequence preconditioning Prof. S. Boyd, EE364b, Stanford University Three
More informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationFast Nyström Methods for Parabolic Boundary Integral Equations
Fast Nyström Methods for Parabolic Boundary Integral Equations Johannes Tausch Abstract Time dependence in parabolic boundary integral operators appears in form of an integral over the previous time evolution
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationCourse Notes: Week 1
Course Notes: Week 1 Math 270C: Applied Numerical Linear Algebra 1 Lecture 1: Introduction (3/28/11) We will focus on iterative methods for solving linear systems of equations (and some discussion of eigenvalues
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationHierarchical Matrices in Computations of Electron Dynamics
Hierarchical Matrices in Computations of Electron Dynamics Othmar Koch Christopher Ede Gerald Jordan Armin Scrinzi July 20, 2009 Abstract We discuss the approximation of the meanfield terms appearing in
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationLecture 8: Fast Linear Solvers (Part 7)
Lecture 8: Fast Linear Solvers (Part 7) 1 Modified Gram-Schmidt Process with Reorthogonalization Test Reorthogonalization If Av k 2 + δ v k+1 2 = Av k 2 to working precision. δ = 10 3 2 Householder Arnoldi
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 19: Computing the SVD; Sparse Linear Systems Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More information