Radial basis function partition of unity methods for PDEs
|
|
- Miranda Hubbard
- 6 years ago
- Views:
Transcription
1 Radial basis function partition of unity methods for PDEs Elisabeth Larsson, Scientific Computing, Uppsala University Credit goes to a number of collaborators Alfa Ali Alison Lina Victor Igor Heryudono Safdari Ramage von Sydow Shcherbakov Tominec Localized Kernel-Based Meshless Methods for Partial Differential Equations, ICERM, Aug 8, 2017 E. Larsson, Aug 8, 2017 (1 : 28)
2 RBF partition of unity methods for PDEs E. Larsson, Aug 8, 2017 (2 : 28)
3 Motivation for Global RBF approximation + Ease of implementation in any dimension. + Flexibility with respect to geometry. + Potentially spectral convergence rates. Computationally expensive for large problems. Local RBF approximations on patches are blended into a global solution using a partition of unity. Provides spectral or high-order convergence. Solves the computational cost issues. Allows for local adaptivity. Wendland (2002), Fasshauer (2007), Cavoretto, De Rossi, Perracchione et al., Larsson, Heryudono et al. E. Larsson, Aug 8, 2017 (3 : 28)
4 The RBF partition of unity method Global approximation P ũ(x) = w j (x)ũ j (x) j=1 PU weight functions Generate weight functions from compactly supported C 2 Wendland functions ψ(ρ) = (4ρ + 1)(1 ρ) 4 + using Shepard s method w i (x) = ψ i (x) M j=1 ψ j (x). Ω j Cover Each x Ω must be in the interior of at least one Ω j. Patches that do not contain unique points are pruned. E. Larsson, Aug 8, 2017 (4 : 28)
5 Differentiating approximations Applying an operator globally M ũ = w i ũ i + 2 w i ũ i + w i ũ i i=1 Local differentiation matrices Let u i be the vector of nodal values in patch Ω i, then u i = Aλ i, where A ij = φ j (x i ) Lu i = A L A 1 u i, where A L ij = Lφ j (x i ). The global differentiation matrix Local contributions are added into the global matrix. 300 E. Larsson, Aug 8, 2017 (5 : 28) nz = 10801
6 An collocation method Choices & Implications Nodes and evaluation points coincide. Square matrix, iterative solver available (Heryudono, Larsson, Ramage, von Sydow 2015). Global node set. Solutions ũ i (x k ) = ũ j (x k ) for x k in overlap regions. Patches are cut by the domain boundary. Potentially strange shapes and lowered local order. E. Larsson, Aug 8, 2017 (6 : 28)
7 An least squares method Choices & Implications Each patch has an identical node layout. Computational cost for setup is drastically reduced. Evaluation nodes are uniform. Easy to generate both local and global high quality node sets. Patches have nodes outside the domain. Good for local order, but requires denser evaluation points. E. Larsson, Aug 8, 2017 (7 : 28)
8 The interpolation error E α = D α (I (u) u) = M j=1 β α The weight functions For C k weight functions and α k D α w j L (Ω j ) C α H α j ( ) α D β w j D α β (I (u j ) u j ) β, H j = diam(ω j ). The local RBF interpolants (Gaussians) Define the local fill distance h j (Rieger, Zwicknagl 2010) D α (I (u j ) u j ) L ( Ω j ) c α,jh m j d 2 α j u j N ( Ω j ), D α (I (u j ) u j ) L ( Ω j ) eγ α,j log(h j )/ h j u j N ( Ωj ). E. Larsson, Aug 8, 2017 (8 : 28)
9 interpolation error estimates Algebraic estimate for H j /h j = c E α L (Ω) K max 1 j M C j H m j d 2 α j u N ( Ωj ) K Maximum # of Ω j overlapping at one point m j Related to the local # of points Ω j Ω j Ω Spectral estimate for fixed partitions E α L (Ω) K max 1 j M Ce γ j log(h j )/ h j u N ( Ωj ) Implications Bad patch reduces global order. Two refinement modes. Guidelines for adaptive refinement. E. Larsson, Aug 8, 2017 (9 : 28)
10 Error estimate for PDE approximation The PDE estimate ũ u L (Ω) C P E L + C P L,X L + Y,X (C M δ M + E L ), where C P is a well-posedness constant and C M δ M is a small multiple of the machine precision. Implications Interpolation error E L provides convergence rate. Norm of inverse/pseudoinverse can be large. Matrix norm better with oversampling. Finite precision accuracy limit involves matrix norm. Follows strategies from Schaback (2007) and Schaback (2016) E. Larsson, Aug 8, 2017 (10 : 28)
11 Does require stable methods? In order to achieve convergence we have two options Refine patches such that diameter H decreases. Increase node numbers such that N j increases. In both cases, theory assumes ε fixed. The effect of patch refinement H = 1, ε = 4 H = 0.5, ε = 4 H = 0.25, ε = H H H The RBF QR method: Stable as ε 0 for N 1 Effectively a change to a stable basis. Fornberg, Piret (2007), Fornberg, Larsson, Flyer (2011), Larsson, Lehto, Heryudono, Fornberg (2013) E. Larsson, Aug 8, 2017 (11 : 28)
12 Effects on the local matrices Local contribution to a global Laplacian L j = (Wj A j + 2Wj A j + W j A j )A 1 j. Typically: A j ill-conditioned, L j better conditioned. for accuracy Stable for small RBF shape parameters ε Change of basis à = AQR1 T D1 T Same result in theory à L à 1 = A L A 1 More accurate in practice Relative error in A j A 1 j without N=10 N=20 N= ε E. Larsson, Aug 8, 2017 (12 : 28)
13 Poisson test problems in 2-D Domain Ω = [ 2, 2] 2. Uniform nodes in the collocation case. u R (x, y) = 1 25x 2 +25y 2 +1 u T (x, y) = sin(2(x 0.1) 2 ) cos((x 0.3) 2 )+sin 2 ((y 0.5) 2 ) E. Larsson, Aug 8, 2017 (13 : 28)
14 Error results with and without RBF QR Error Least squares Fixed shape ε = 0.5 or scaled such that εh = c Left: 5 5 patches Spectral mode RBF QR Direct Scaled Points per dimension Right: 55 points per patch Error Algebraic mode p= 7 RBF QR Direct Scaled p= Patches per dimension With RBF QR better results for H/h large. Scaled approach good until saturation. E. Larsson, Aug 8, 2017 (14 : 28)
15 as a function of patch size Runge Trig Error H Error H Collocation (dashed lines) and Least Squares (solid lines). Points per patch n = 28, 55, 91. Theoretical rates p = 4, 7, 10. Numerical rates p 3.9, 6.9, 9.8. E. Larsson, Aug 8, 2017 (15 : 28)
16 Spectral convergence for fixed patches Error Runge /h Error Trig /h Collocation (dashed lines) and Least Squares (solid lines). LS-RBF-PU is significantly more accurate due to the constant number of nodes per patch. E. Larsson, Aug 8, 2017 (16 : 28)
17 and large scale problems The global error estimate ũ u L (Ω) C P E L + C P L,X L + Y,X (C M δ M + E L ) The dark horse is the stability matrix norm The stability norm is related to conditioning. In the collocation case, L 1 X,X grows with N. How does it behave with least squares? E. Larsson, Aug 8, 2017 (17 : 28)
18 The stability matrix norm: fill distance Fixed patch size with patches. Oversampling M/N = 1.1, 1.2, 1.5 Stability norm /h Error /h Collocation (dashed) and Least Squares (solid) Stability norm grows exponentially as h decreases. Oversampling reduces the norm. Collocation is less robust. E. Larsson, Aug 8, 2017 (18 : 28)
19 Stability norm: Patch size Fixed number of points per patch n = 28, 55, 91 Results as a function of patch diameter H Stability(H) Error(H) Stability norm H Error H Collocation (dashed) and LS (solid) The norm does not grow for LS- (!) E. Larsson, Aug 8, 2017 (19 : 28)
20 Oversampling provides stability (at a cost). For robustness in h, M/N needs to grow with N. E. Larsson, Aug 8, 2017 (20 : 28) Stability norm: Oversampling Fixed number of local points n = 55. Patches: 5 5, 10 10, Stability norm Stability(M/N) β Collocation (dashed) and LS (solid)
21 Poisson test problems in 3-D ( ) u(x) = sin π(x1 0.5)x 3 log(x 2 +3) 1 0 E. Larsson, Aug 8, 2017 (21 : 28)
22 in 3-D Error n = 20, 35, 56, 84, 120, 165 points per patch Sphere Pepper Spectral p = 1.7 p = 4.1 p = 5.4 p = 6.6 p = 6.1 p = H Error p = 1.5 p = 3.6 p = 3.9 p = 6.2 p = 5.6 p = H Error /h order is a bit better than expected. Refinement modes work as expected. E. Larsson, Aug 8, 2017 (22 : 28)
23 Computational time comparison in 2-D Time Collocation LS Error Experiments with fixed n = 28, 55, 91. Time Error For a fixed problem size, RBF-PU-LS involves more work, but yields a smaller error. Overall, the LS approach is the winner. Strategy: Points per patch adjusted to tolerance. E. Larsson, Aug 8, 2017 (23 : 28)
24 Solution using iterative methods So far implemented for the collocation approach. First step: Re-ordering of the unknowns to improve structure. Patches: Preceded and followed by a neighbour. Nodes xk : Define home patch Ωj such that wj wi (xk ). Within patch: Sub-order according to node memberships. Heryudono, Larsson, Ramage, and von Sydow (2015) E. Larsson, Aug 8, 2017 (24 : 28)
25 Matrix structures with snake order I The patch order fills in gaps in the band. I The node order minimizes the central band width. I Chosen preconditioner, ILU(0) of central band. Structure 2-D Structure 3-D nz = E. Larsson, Aug 8, 2017 (25 : 28) 300
26 Results for the iterative method contd. Results in 2-D with Halton nodes N # it no prec # it ILU(0) Time gain Results in 3-D with Halton nodes N # it no prec # it ILU(0) Time gain Memory gain not to be forgotten. E. Larsson, Aug 8, 2017 (26 : 28)
27 Current project: With Alfa Heryudono, Danilo Stocchino and Stefano De Marchi Box structure helps book keeping. Overlaps complicate things. Evaluation points placed using node placing algorithm. Fornberg & Flyer (2015) E. Larsson, Aug 8, 2017 (27 : 28)
28 Results is a flexible tool for solving PDEs. By using template patches, LS becomes computationally efficient. LS is numerically stable for large scale problems. Things to do Change patch type to reduce unnecessary overlap. Adaptive algorithms based on LS. Time-stability for LS. For collocation and hyperbolic PDEs on the sphere, see poster by Igor Tominec et al. E. Larsson, Aug 8, 2017 (28 : 28)
Radial basis function approximation
Radial basis function approximation PhD student course in Approximation Theory Elisabeth Larsson 2017-09-18 E. Larsson, 2017-09-18 (1 : 77) Global RBF approximation Stable evaluation methods RBF partition
More informationA scaling perspective on accuracy and convergence in RBF approximations
A scaling perspective on accuracy and convergence in RBF approximations Elisabeth Larsson With: Bengt Fornberg, Erik Lehto, and Natasha Flyer Division of Scientific Computing Department of Information
More informationPreconditioning for radial basis function partition of unity methods
Journal of Scientific Computing manuscript No. (will be inserted by the editor) Preconditioning for radial basis function partition of unity methods Alfa Heryudono Elisabeth Larsson Alison Ramage Lina
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationRBF-Partition of Unity Method: an overview of recent results
RBF-Partition of Unity Method: an overview of recent results Alessandra De Rossi a In collaboration with Roberto Cavoretto a, Emma Perracchione b a Department of Mathematics G. Peano University of Torino
More informationRadial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance
1 / 27 Radial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance S. Milovanović, L. von Sydow Uppsala University Department of Information Technology Division
More informationRadial basis functions topics in slides
Radial basis functions topics in 40 +1 slides Stefano De Marchi Department of Mathematics Tullio Levi-Civita University of Padova Napoli, 22nd March 2018 Outline 1 From splines to RBF 2 Error estimates,
More informationSTABLE CALCULATION OF GAUSSIAN-BASED RBF-FD STENCILS
STABLE CALCULATION OF GAUSSIAN-BASED RBF-FD STENCILS BENGT FORNBERG, ERIK LEHTO, AND COLLIN POWELL Abstract. Traditional finite difference (FD) methods are designed to be exact for low degree polynomials.
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 40: Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationSpace-Time Localized Radial Basis Function Collocation Methods for PDEs
Space-Time Localized Radial Basis Function Collocation Methods for PDEs Alfa Heryudono Department of Mathematics University of Massachusetts Dartmouth ICERM Topical Workshop Providence August 27 This research
More informationRBF-FD Approximation to Solve Poisson Equation in 3D
RBF-FD Approximation to Solve Poisson Equation in 3D Jagadeeswaran.R March 14, 2014 1 / 28 Overview Problem Setup Generalized finite difference method. Uses numerical differentiations generated by Gaussian
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 informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 3: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu
More informationNumerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method To cite this article: R
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationConsistency Estimates for gfd Methods and Selection of Sets of Influence
Consistency Estimates for gfd Methods and Selection of Sets of Influence Oleg Davydov University of Giessen, Germany Localized Kernel-Based Meshless Methods for PDEs ICERM / Brown University 7 11 August
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationFast and stable rational RBF-based Partition of Unity interpolation
Fast and stable rational RBF-based Partition of Unity interpolation S. De Marchi, A. Martínez, E. Perracchione Dipartimento di Matematica Tullio Levi-Civita, Università di Padova, Italia Abstract We perform
More informationA new stable basis for radial basis function interpolation
A new stable basis for radial basis function interpolation Stefano De Marchi and Gabriele Santin Department of Mathematics University of Padua (Italy) Abstract It is well-known that radial basis function
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 informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationRBF Collocation Methods and Pseudospectral Methods
RBF Collocation Methods and Pseudospectral Methods G. E. Fasshauer Draft: November 19, 24 Abstract We show how the collocation framework that is prevalent in the radial basis function literature can be
More informationTheoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions
Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions Elisabeth Larsson Bengt Fornberg June 0, 003 Abstract Multivariate interpolation of smooth
More informationA flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds
A flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds Andriy Sokolov Dmitri Kuzmin, Oleg Davydov and Stefan Turek Institut für Angewandte Mathematik (LS3) TU Dortmund
More informationACCURACY AND STABILITY OF GLOBAL RADIAL BASIS FUNCTION METHODS FOR THE NUMERICAL SOLUTION OF PARTIAL DIFFERENTIAL EQUATIONS. by Rodrigo B.
ACCURACY AND STABILITY OF GLOBAL RADIAL BASIS FUNCTION METHODS FOR THE NUMERICAL SOLUTION OF PARTIAL DIFFERENTIAL EQUATIONS by Rodrigo B. Platte A dissertation submitted to the Faculty of the University
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationNumerical solution of surface PDEs with Radial Basis Functions
Numerical solution of surface PDEs with Radial Basis Functions Andriy Sokolov Institut für Angewandte Mathematik (LS3) TU Dortmund andriy.sokolov@math.tu-dortmund.de TU Dortmund June 1, 2017 Outline 1
More informationGreedy Kernel Techniques with Applications to Machine Learning
Greedy Kernel Techniques with Applications to Machine Learning Robert Schaback Jochen Werner Göttingen University Institut für Numerische und Angewandte Mathematik http://www.num.math.uni-goettingen.de/schaback
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
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 informationRBF kernel method and its applications to clinical data
Proceedings of Kernel-based Methods and Function Approximation 216, Volume 9 216 Pages 13 18 RBF kernel method and its applications to clinical data Emma Perracchione a Ilaria Stura b Abstract In this
More informationStable Parameterization Schemes for Gaussians
Stable Parameterization Schemes for Gaussians Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado Denver ICOSAHOM 2014 Salt Lake City June 24, 2014 michael.mccourt@ucdenver.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationLab 1: Iterative Methods for Solving Linear Systems
Lab 1: Iterative Methods for Solving Linear Systems January 22, 2017 Introduction Many real world applications require the solution to very large and sparse linear systems where direct methods such as
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationAll Well Posed Problems have Uniformly Stable and Convergent Discretizations
Noname manuscript No. (will be inserted by the editor) All Well Posed Problems have Uniformly Stable and Convergent Discretizations Robert Schaback Received: date / Accepted: date Abstract This paper considers
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationSOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS
SOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS NATASHA FLYER 1. Introduction Mathematical modeling of space and climate phenomena generally requires
More informationUsing Gaussian eigenfunctions to solve boundary value problems
Advances in Applied Mathematics and Mechanics Adv Appl Math Mech, Vol xx, No x, pp xx-xx DOI: xxxx xxx 200x Using Gaussian eigenfunctions to solve boundary value problems Michael McCourt Center for Applied
More informationNo. 579 November 2017
No. 579 November 27 Numerical study of the RBF-FD level set based method for partial differential equations on evolving-in-time surfaces A. Sokolov, O. Davydov, S. Turek ISSN: 29-767 Numerical study of
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationNumerical analysis of heat conduction problems on irregular domains by means of a collocation meshless method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on irregular domains by means of a collocation meshless method To cite this article: R Zamolo and
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 39: Non-Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationCorrespondence should be addressed to Qibai Huang,
Mathematical Problems in Engineering Volume 2012, Article ID 632072, 20 pages doi:10.1155/2012/632072 Research Article An Improved Collocation Meshless Method Based on the Variable Shaped Radial Basis
More informationMultivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness
Multivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness Guohui Song John Riddle Gregory E. Fasshauer Fred J. Hickernell Abstract In this paper, we consider multivariate
More informationRecent developments in the dual receiprocity method using compactly supported radial basis functions
Recent developments in the dual receiprocity method using compactly supported radial basis functions C.S. Chen, M.A. Golberg and R.A. Schaback 3 Department of Mathematical Sciences University of Nevada,
More informationBreaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points
Breaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points Michael Griebel Christian Rieger Peter Zaspel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. UMER. AAL. Vol. 43, o. 2, pp. 75 766 c 25 Society for Industrial and Applied Mathematics POLYOMIALS AD POTETIAL THEORY FOR GAUSSIA RADIAL BASIS FUCTIO ITERPOLATIO RODRIGO B. PLATTE AD TOBI A. DRISCOLL
More informationThe automatic solution of PDEs using a global spectral method
The automatic solution of PDEs using a global spectral method 2014 CBMS-NSF Conference, 29th June 2014 Alex Townsend PhD student University of Oxford (with Sheehan Olver) Supervised by Nick Trefethen,
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationStability and Lebesgue constants in RBF interpolation
constants in RBF Stefano 1 1 Dept. of Computer Science, University of Verona http://www.sci.univr.it/~demarchi Göttingen, 20 September 2008 Good Outline Good Good Stability is very important in numerical
More informationDOUBLY STOCHASTIC RADIAL BASIS FUNCTION METHODS
DOUBLY STOCHASTIC RADIAL BASIS FUNCTION METHODS FENGLIAN YANG, LIANG YAN, AND LEEVAN LING Abstract. We propose a doubly stochastic radial basis function (DSRBF) method for function recoveries. Instead
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationA RADIAL BASIS FUNCTION (RBF) COMPACT FINITE DIFFERENCE (FD) SCHEME FOR REACTION-DIFFUSION EQUATIONS ON SURFACES
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 A RADIAL BASIS FUNCTION (RBF) COMPACT FINITE DIFFERENCE (FD) SCHEME FOR REACTION-DIFFUSION EQUATIONS ON SURFACES ERIK LEHTO, VARUN SHANKAR, AND GRADY
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 informationOversampling for the partition of unity parallel finite element algorithm
International conference of Computational Mathematics and its application Urumqi, China, July 25-27, 203 Oversampling for the partition of unity parallel finite element algorithm HAIBIAO ZHENG School of
More informationA Hybrid Meshless-Collocation/Spectral-Element Method for Elliptic Partial Differential Equations
International Journal of Computer Mathematics Vol., No., Month x, A Hybrid Meshless-Collocation/Spectral-Element Method for Elliptic Partial Differential Equations Christopher D. Blakely University of
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590
More informationMultiscale RBF collocation for solving PDEs on spheres
Multiscale RBF collocation for solving PDEs on spheres Q. T. Le Gia I. H. Sloan H. Wendland Received: date / Accepted: date Abstract In this paper, we discuss multiscale radial basis function collocation
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationA orthonormal basis for Radial Basis Function approximation
A orthonormal basis for Radial Basis Function approximation 9th ISAAC Congress Krakow, August 5-9, 2013 Gabriele Santin, joint work with S. De Marchi Department of Mathematics. Doctoral School in Mathematical
More informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
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 informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationA Numerical Study of the Accuracy and Stability of Symmetric and Asymmetric RBF Collocation methods for Hyperbolic PDEs
A Numerical Study of the Accuracy and Stability of Symmetric and Asymmetric RBF Collocation methods for Hyperbolic PDEs Scott A. Sarra Marshall University May 9, 27 Abstract Differentiation matrices associated
More informationdifferential problems
A numerical study of 2D integrated RBFNs incorporating Cartesian grids for solving 2D elliptic differential problems N. Mai-Duy and T. Tran-Cong Computational Engineering and Science Research Centre Faculty
More informationNational Taiwan University
National Taiwan University Meshless Methods for Scientific Computing (Advisor: C.S. Chen, D.L. oung) Final Project Department: Mechanical Engineering Student: Kai-Nung Cheng SID: D9956 Date: Jan. 8 The
More informationA Radial Basis Function (RBF) Compact Finite Difference (FD) Scheme for Reaction-Diffusion Equations on Surfaces
Boise State University ScholarWorks Mathematics Faculty Publications and Presentations Department of Mathematics 1-1-2017 A Radial Basis Function (RBF) Compact Finite Difference (FD) Scheme for Reaction-Diffusion
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationA RADIAL BASIS FUNCTION METHOD FOR THE SHALLOW WATER EQUATIONS ON A SPHERE
A RADIAL BASIS FUNCTION METHOD FOR THE SHALLOW WATER EQUATIONS ON A SPHERE NATASHA FLYER AND GRADY B. WRIGHT Abstract. The paper derives the first known numerical shallow water model on the sphere using
More informationMeshfree Exponential Integrators
Meshfree joint work with A. Ostermann (Innsbruck) M. Caliari (Verona) Leopold Franzens Universität Innsbruck Innovative Integrators 3 October 2 Meshfree Problem class: Goal: Time-dependent PDE s with dominating
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationOn Multivariate Newton Interpolation at Discrete Leja Points
On Multivariate Newton Interpolation at Discrete Leja Points L. Bos 1, S. De Marchi 2, A. Sommariva 2, M. Vianello 2 September 25, 2011 Abstract The basic LU factorization with row pivoting, applied to
More informationSparse Linear Systems. Iterative Methods for Sparse Linear Systems. Motivation for Studying Sparse Linear Systems. Partial Differential Equations
Sparse Linear Systems Iterative Methods for Sparse Linear Systems Matrix Computations and Applications, Lecture C11 Fredrik Bengzon, Robert Söderlund We consider the problem of solving the linear system
More informationSolving linear equations with Gaussian Elimination (I)
Term Projects Solving linear equations with Gaussian Elimination The QR Algorithm for Symmetric Eigenvalue Problem The QR Algorithm for The SVD Quasi-Newton Methods Solving linear equations with Gaussian
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationStability of radial basis function methods for convection problems on the circle and sphere
Stability of radial basis function methods for convection problems on the circle and sphere Jordan M. Martel and Rodrigo B. Platte Abstract This paper investigates the stability of the Radial Basis Function
More informationComputing the common zeros of two bivariate functions via Bézout resultants Colorado State University, 26th September 2013
Work supported by supported by EPSRC grant EP/P505666/1. Computing the common zeros of two bivariate functions via Bézout resultants Colorado State University, 26th September 2013 Alex Townsend PhD student
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
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 informationStudia Geophysica et Geodaetica. Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids
Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids Journal: Studia Geophysica et Geodaetica Manuscript ID: SGEG-0-00 Manuscript Type: Original Article Date Submitted
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationA New Trust Region Algorithm Using Radial Basis Function Models
A New Trust Region Algorithm Using Radial Basis Function Models Seppo Pulkkinen University of Turku Department of Mathematics July 14, 2010 Outline 1 Introduction 2 Background Taylor series approximations
More informationA Meshfree Method for Solving the Monge Ampère Equation Klaus Böhmer 1, Robert Schaback 2
INTRODUCTION A Meshfree Method for Solving the Monge Ampère Equation Klaus Böhmer, Robert Schaback 2 Abstract: This paper solves the two-dimensional Dirichlet problem for the Monge-Ampère equation by a
More informationMeshless algorithm for partial differential equations on open and singular surfaces
Hong Kong Baptist University HKBU Institutional Repository Open Access Theses and Dissertations Electronic Theses and Dissertations 3-11-216 Meshless algorithm for partial differential equations on open
More informationLecture Notes for LG s Diff. Analysis
Lecture Notes for LG s Diff. Analysis trans. Paul Gallagher Feb. 18, 2015 1 Schauder Estimate Recall the following proposition: Proposition 1.1 (Baby Schauder). If 0 < λ [a ij ] C α B, Lu = 0, then a ij
More informationA Cartesian-Grid Discretisation Scheme Based on Local Integrated RBFNs for Two-Dimensional Elliptic Problems
Copyright 2009 Tech Science Press CMES, vol.51, no.3, pp.213-238, 2009 A Cartesian-Grid Discretisation Scheme Based on Local Integrated RBFNs for Two-Dimensional Elliptic Problems N. Mai-Duy 1 and T. Tran-Cong
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 information