RBF-FD Approximation to Solve Poisson Equation in 3D
|
|
- Junior Moore
- 6 years ago
- Views:
Transcription
1 RBF-FD Approximation to Solve Poisson Equation in 3D Jagadeeswaran.R March 14, / 28
2 Overview Problem Setup Generalized finite difference method. Uses numerical differentiations generated by Gaussian RBF interpolation on irregular centers. Dirichlet problem for the Poisson Equation in 2D and 3D with smooth solution. RBF-QR method used when the Interpolation matrix is ill conditioned due to small Shape parameters 2 / 28
3 RBF-FD Problem Definition Dirichlet Problem for the Poisson equation: u = f on Ω u δω = g We want to find a Linear numerical differentiation formula for the Differential operator: Du(x) n w i (x)u(x i ) i=1 where w i (x) called weights, The vector w = [w 1,..., w n ] T called stencil. 3 / 28
4 RBF-FD Problem Definition... Let Discretization centers: Ξ Ω, δξ := Ξ δω, Ξ int := Ξ \ δξ Assume for each ζ Ξ int a set Ξ ζ Ξ is chosen such that ζ Ξ ζ and Ξ = ζ Ξ int Ξ ζ For each ζ Ξ int, choose a linear numerical differentiation formula for Laplace operator: u(ζ) Ξ Ξ ζ w ζ,ξ u(ξ) So we get the approximation: Ξ Ξ ζ w ζ,ξ û(ξ) = f (ζ), ζ Ξ int ; û(ξ) = g(ξ), Ξ δξ (1) 4 / 28
5 RBF-FD My Goals Add Stable Computation Algorithm from G.E. Fasshauer, M.J. McCourt, Add Support for 3 Dimensional Problems Generalize Stencil Selection for n Dimensional Domains Experiment Stable computation for n Dimensional Problems CVT Points for n Dimensions 5 / 28
6 Interpolation Problem Let a Positive(or Conditionally) definite function κ : R + R and a continuous function u : R d R Given Ξ ζ = {x 0,..., x n} R d, ζ = x 0, K j (x) = K(x x j ), K(x) = κ( x ) the RBF interpolant with a constant term is sought in the form n s(x) = a j K j (x) + c, s(x i ) = u(x i ), i = 0,..., n, j=0 n a j = 0 (2) j =0 Then the coefficients a j and c are uniquely determined by the following conditions n a j K(x i x j ) + c = u(x i ), i = 0,..., n, j=0 n a j = 0 j =0 written in matrix form as [ ] KX 1 [ ] [ ac u x 1 T = 0 0 ], K X := [e ɛ2 x i x j 2 ] n i,j=0, 1 := [1,..., 1]T (3) 6 / 28
7 RBF Differentiation Matrix Consider the approximation of a differential operator D by using RBF n w j K(x x j ) + v = DK(x x i ), i = 0,..., n, Note: This is indept ofu j=0 These weights can be found by solving the symmetric positive definite linear system [ ] KX 1 [ ] [ ] wv [Dki (x 1 T = 0 )] n i=0 0 0 (4) where, κ j (x) = K(x x j ), K(x) := κ( x ), K X := [e ɛ2 x i x j 2 ] n i,j=0, For Gaussian RBF κ(r) = e ɛr 2. Its Laplacian is: k(x) = 2ɛ 2 e ɛ2 x i x j 2 (2ɛ 2 x i x j 2 d), d := num of dimensions 7 / 28
8 Differentiation Problem The derivatives of s are good approximations of the derivatives of u if k is sufficiently smooth. An approximation of Du(x) at x 0, where D is a linear differential operator annihilating constants, may be considered in the form. Apply D to (2) Du(x 0 ) Ds(x 0 ) = n a j DK(x 0 x j ) = j=0 n w i u(x i ) where the weights w i (depending on x) exist because the coefficients a j of the interpolation function s defined as before depend linearly on the data u(x i ). This can be verified as below using (4) n j=0 [ ] [ ac T [K(x xi )] a j DK(x x j ) = n i=0 0 i=0 ] [ ] [ ac T KX 1 = 1 T 0 ] [ wv ] [ u x = 0 ] T [ wv ] = n w i u(x i ) i=0 8 / 28
9 Example Stencil Figure: typical Stencil Points in 2D Stencil used to compute Laplacian, in case of uniform points with h fill distance u(ζ) = 1 (u(ζ + (h, 0)) + u(ζ (h, 0)) + u(ζ + (0, h)) + u(ζ (0, h)) 4u(ζ)) h2 9 / 28
10 Stencil Selection algorithm in R 2 from Davydov Given ζ Ξ \ δξ, we will select a set Ξ ζ containing ζ. Let Ξ ζ = {x 0, x 1,..., x k }, ζ := x 0, where the points x 0, x 1,..., x k are ordered counterclockwise with respect to ζ. Consider the following cost function: µ(x 1,..., x k ) := k i=1 α2 i where α i denotes the angle between the rays ζx i, ζx i+1 in the counterclockwise direction. We will also need the minimum and the maximum angle: α min = min{α 1,..., α k } α max = max{α 1,..., α k } Since k i=1 α i = 2π, the expression k i=1 α2 i achieves its unique minimum for α 1 =... = α k = 2π/k, that is for the uniformly spaced directions ζx i if x 1,..., x k were chosen freely in R 2. So the goal of the algorithm below is to choose {x 1,..., x k } Ξ such that µ(x 1,..., x k is minimised while keeping the distances x i ζ as small as possible. To achieve a balance between the goals of a small µ and small distances, we introduce the restriction that x i must be among m closest points to ζ, and terminate the algorithm if the set {ζ, x 1,..., x k } satisfies α max <= uα min, where m > k and u > 1.0 are parameters to be determined empirically. 10 / 28
11 Stencil Selection algorithm in R 2 from Davydov... Input: Ξ, ζ. Output: Ξ ζ. Parameters: k (the target number of points x i ), m > k (the number of points in the local cloud) and u > 1 (the angle uniformity tolerance). Parameter values used in our numerical experiments: k = 6, m = 30, u = 3.0. I Find m nearest points {x 1,..., x m} closer to ζ in Ξ \ ζ, sorted by increasing distance to ζ, and initialise Ξ ζ := {ζ, x 1,..., x k }. If α max <= uα min, then STOP: return Ξ ζ. II For i = n + 1,..., m: 1 Compute the angles α 1,..., α k+1 formed by the extended set {x 1,..., x k+1 } = {x 1,..., x k, x i }. 2 If both angles between ζx i and its two neighbouring rays are greater than the minimum angle α min := α min(x 1,..., x k, x i ): i Find j such that α j = α. Choose p = j or p = j +1 depending on whether α j 1 < α j+1 orα j 1 >= α j+1. ii If µ(x 1,..., x k+1 x p ) < µ(x 1,..., x k ): a. Update {x 1,..., x k } = {x 1,..., x k+1 } \ {x p }. b. If α max <= uα min STOP: return the current set Ξ ζ = {ζ, x 1,..., x k }. III Observe that α max > uα min must hold for the current set Ξ ζ = {ζ, x 1,..., x k } if the algorithm has not been terminated earlier. Find j such that α j = α min (x 1,..., x k ). Choose p = j or p = j + 1 depending on whether α j 1 < α j+1 or α j 1 >= α j+1. STOP: return Ξ ζ = {ζ, x 1,..., x k } \ {x p}. 11 / 28
12 Stencil Selection algorithm in R 3 Input: Ξ, ζ. Output: Ξ ζ. Parameter values used in our numerical experiments: k = 12, m = 30, u = Cost function: Standard Deviation in Sum of Squared Distance between each other points in the Stencil (Energy). µ(x 1,..., x k ) := std({e i : i = 1,..., k}) where E i := E(x i ) = j=1,i j x i x j 2 Another option: Projected Distance: shift the current set of stencil point to make zeta as the center and normalize, i.e. {x i = (x i ζ) : i = 1,..., k}, then x i ζ compute the Squared Distance between each other points I Find m nearest points {x 1,..., x m} closer to ζ in Ξ \ ζ, sorted by increasing distance to ζ, and initialise Ξ ζ := {ζ, x 1,..., x k }. If E max <= ue min, then STOP: return Ξ ζ. II For i = n + 1,..., m: 1 Compute the angles E 1,..., E k+1 formed by the extended set {x 1,..., x k+1 } = {x 1,..., x k, x i }. 2 Find p = j: such that E j == min({e i }) 3 If µ(x 1,..., x k+1 x p ) < µ(x 1,..., x k ): a. Update {x 1,..., x k } = {x 1,..., x k+1 } \ {x p}. b. If E max <= ue min STOP: return the current set Ξ ζ = {ζ, x 1,..., x k }. III STOP: return Ξ ζ = {ζ, x 1,..., x k }. 12 / 28
13 CVT: Centroidal Voronoi Points Algorithm based iterative voronoi tessellation There is an initial assignment for the generators [Halton, random, uniform...] In each iteration, estimates is made of the volume and location of the Voronoi cells Sampling done by Monte Carlo sampling The accuracy of the resulting CVT depends in part on the number of sampling points and the number of iterations taken. Lloyd s algorithm Can be customized to generate points for any dimensions, regions Boundary points placements Algorithm called: Constrained CVT CCVTBOX software 13 / 28
14 CVT points generation 14 / 28
15 RBF-QR Since the matrix K x is extremely ill-conditioned for small ɛ, we need to use Stable computation in this case. Hilbert-Schmidt series: K(x, z) = n=1 λ nϕ n (x)ϕ n (z) When we have such an Expansion, we truncate the series to M such that K(x i, x j ) = M n=1 λ nϕ n (x i )ϕ n (x j ) and write as a finite matrix: So we get K = ΦΛΦ T 15 / 28
16 New Basis Note: λ n = since λ decays in power n, it leads to ill conditioned interpolation matrix when ɛ is small ( ) n 1 α 2 ɛ 2 α 2 +δ 2 +ɛ 2 α 2 +δ 2 +ɛ 2 We want to find a new Basis which improves the conditioning of the problems. Compute Φ = QR = Q (R 1 R 2 ) K = (QR)ΛΦ T = Q (R 1 R 2 ) ΛΦ T With more Simplifications we get a new basis: Ψ as shown next. 16 / 28
17 Finding New Basis Ψ K = ΦΛΦ T = Q (R 1 R 2 ) ( Λ1 Λ 2 ) Φ T = Q ( R 1 Λ 1 R 2 Λ 2 ) Φ T ( = QR 1 Λ 1 IN Λ 1 ) }{{} 1 R 1 1 R 2 Λ 2 Φ T }{{} ( = QR 1 Λ 1 IN Λ 1 ) ( Φ T ) }{{} 1 R R 2 Λ 2 Φ T 2 }{{} ( = QR 1 Λ 1 Φ T }{{} 1 + Λ 1 ) 1 R 1 1 R 2 Λ 2 Φ T 2 }{{} = X = Ψ So Now we have a new basis, more stable due to Λ 1 1 Λ 2 : Ψ = X 1 K = Φ T 1 + Λ 1 1 R 1 1 R 2 Λ 2 Φ T 2 18 / 28
18 Computation of RBF-FD Stencils weights by QR Method Better Conditioned basis With the new basis, [ modify the weight ] computation matrix. By left-multiplying both sides of (4) by: (QR1 Λ 1 ) [ (QR1 Λ 1 ) ] [ KX 1 1 T 0 ] [ ] [ wv = (QR1 Λ 1 ) ] [ [ ki (x 0 )] n i=0 0 ] [ ΨX h 1 T 0 ] [ ] [ wv [ ψi (x = 0 )] n i=0 0 ] (5) where, h = (QR 1 Λ 1 ) 1 1 Are We done yet? Computing h is still problematic, since it involves Λ / 28
19 Numerical Computation of RBF-FD Stencils weights by QR Method For numerical implementation Ψ is truncated: [ ΨX h 1 T 0 ] [ ] [ wv = [ ψi (x 0 )] n i=0 0 ] (6) where, Ψ X := ψ 0 (x 0 )... ψ0 (x n).. ψ n(x 0 )... ψn(x n) To improve the stability, following are done: replace h by h = (QR 1 Λ 1 ) 1 1 where Λ 1 is obtained by zero replacing the elements in Λ 1 that exceed the reciprocal of machine epsilon. The above equation is replaced by n ( ψ i (x j ) ψ i (x 0 ))w j + h i v = ψ n i (x 0 ), i = 0,..., n. w 0 = w j j=1 j=1 20 / 28
20 Numerical Computation of RBF-FD Stencils weights by QR Method ψ 0 (x 1 ) ψ 0 (x 2 )... ψ 0 (x n) h(0) ψ 1 (x 1 ) ψ1 (x 2 )... ψ1 (x n) h(1).. ψ n(x 1 ) ψn(x 2 )... ψn(x n) h(n) w 1 w 2... w n = ψ 0 (x 0 ) ψ 1 (x 0 ). ψ n(x 0 ) Also the Stencil points {x 0, x 1, x 2,..., x n} are centered to x 0 : {x 0 x 0, x 1 x 0, x 2 x 0,..., x n x 0 }. So that we have {0, x 1, x 2,..., x n}. This simplifies the computation of [ ψ i (x 0 )] n i=0 as x 0 = 0 (7) 21 / 28
21 Open Questions Points Selection? it does not work for Quasi random points (Sobol, Halton)!! Adaptive points placements, more points closer boundary, corners Boundary points placements? How many points required for the Stencil? Stencil selection criterion? 22 / 28
22 What it does? Computes Mesh Triangulations upto 5 refinements, either in Square or Circle Domain Computes Stencil for each Interior point in the domain For each ɛ compute the Stencil weights using Direct Method if the condition number of the interpolation is < Else use the Stable Method to compute the weights For the given Problem computes the Finite Difference Matrix solves the Problem using the above matrix Computes RMS Error between the actual function and approximation Plots the Error values the defined range of ɛ 23 / 28
23 RMS Error for the 2D test function Figure: RMS Error in Function approx, RMS Error in Differentiation 24 / 28
24 RMS Error for the 3D test function Figure: RMS Error in Function approx, RMS Error in Differentiation 25 / 28
25 RMS Error for the 4D test function Figure: RMS Error in Function approx, RMS Error in Differentiation 26 / 28
26 Conclusion Implemented Davydov s 2D Poisson solver method Integrated with Stable computation of Fasshauer Extended to 3D support Next planning to work on Survey of similar work done so far Work on answering some theoretical questions?? Adaptive meshless refinement Hermite interpolation Solve some real and difficult PDEs 27 / 28
27 Reference O. Davydov, D.T. Oanh, On the optimal shape parameter for Gaussian radial basis function finite difference approximation of the Poisson equation, Comput. Math. Appl. 62 (2011) G.E. Fasshauer, M.J. McCourt, Stable evaluation of Gaussian RBF interpolants, SIAM J. Sci. Comput. 34 (2) (2012) A737A CVT Points: http: //people.sc.fsu.edu/~jburkardt/m_src/cvt/cvt.html 28 / 28
MATH 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 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 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 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 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 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 informationRadial basis function partition of unity methods for PDEs
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
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 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 informationUniform Convergence of a Multilevel Energy-based Quantization Scheme
Uniform Convergence of a Multilevel Energy-based Quantization Scheme Maria Emelianenko 1 and Qiang Du 1 Pennsylvania State University, University Park, PA 16803 emeliane@math.psu.edu and qdu@math.psu.edu
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 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 informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 2: Scattered Data Interpolation in R d Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 Chapter
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 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 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 informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for 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 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 informationNumerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions
Numerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions Z. Jomaa C. Macaskill August 8, 28 Abstract We describe a 2-D finite difference algorithm for inverting
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
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 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 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 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 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 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 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 informationA Feynman-Kac Path-Integral Implementation for Poisson s Equation
for Poisson s Equation Chi-Ok Hwang and Michael Mascagni Department of Computer Science, Florida State University, 203 Love Building Tallahassee, FL 32306-4530 Abstract. This study presents a Feynman-Kac
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 informationSTABLE EVALUATION OF GAUSSIAN RBF INTERPOLANTS
STABLE EVALUATION OF GAUSSIAN RBF INTERPOLANTS GREGORY E. FASSHAUER AND MICHAEL J. MCCOURT Abstract. We provide a new way to compute and evaluate Gaussian radial basis function interpolants in a stable
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 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 informationRecent Results for Moving Least Squares Approximation
Recent Results for Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract. We describe two experiments recently conducted with the approximate moving least squares (MLS) approximation
More informationCompact Matérn Kernels and Piecewise Polynomial Splines Viewed From a Hilbert-Schmidt Perspective
Noname manuscript No (will be inserted by the editor) Compact Matérn Kernels and Piecewise Polynomial Splines Viewed From a Hilbert-Schmidt Perspective Roberto Cavoretto Gregory E Fasshauer Michael McCourt
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 informationMinimal Element Interpolation in Functions of High-Dimension
Minimal Element Interpolation in Functions of High-Dimension J. D. Jakeman, A. Narayan, D. Xiu Department of Mathematics Purdue University West Lafayette, Indiana Random Model as a function We can consider
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
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 informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
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 informationToward Approximate Moving Least Squares Approximation with Irregularly Spaced Centers
Toward Approximate Moving Least Squares Approximation with Irregularly Spaced Centers Gregory E. Fasshauer Department of Applied Mathematics, Illinois Institute of Technology, Chicago, IL 6066, U.S.A.
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 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 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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationNo USING LOCALISED QUADRATIC FUNCTIONS ON AN IRREGULAR GRID FOR PRICING HIGH-DIMENSIONAL AMERICAN OPTIONS. By S.J. Berridge, J.M.
No. 4 USING LOCALISED QUADRATIC FUNCTIONS ON AN IRREGULAR GRID FOR PRICING HIGH-DIMENSIONAL AMERICAN OPTIONS By S.J. Berridge, J.M. Schumacher February 4 ISSN 94-785 Using localised quadratic functions
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
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 informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
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 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 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 informationCompact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems
Copyright 2011 Tech Science Press SL, vol.6, no.2, pp.93-107, 2011 Compact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems T.-T. Hoang-Trieu 1, N. Mai-Duy 1 and T. Tran-Cong
More informationA fast adaptive numerical solver for nonseparable elliptic partial differential equations
J. KSIAM Vol.2. No.1. 27-39, 1998 27 A fast adaptive numerical solver for nonseparable elliptic partial differential equations June-Yub Lee Dept of Math, Ewha Womans Univ, Seoul 120-750, KOREA jylee@math.ewha.ac.kr
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
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 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 informationarxiv: v1 [math.na] 9 Sep 2014
Finite Integral ethod for Solving Poisson-type Equations on anifolds from Point Clouds with Convergence Guarantees arxiv:1409.63v1 [math.na] 9 Sep 014 Zhen Li Zuoqiang Shi Jian Sun September 10, 014 Abstract
More informationA robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme
A robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme International Conference On Preconditioning Techniques For Scientific And Industrial Applications 17-19
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationApplication of modified Leja sequences to polynomial interpolation
Special Issue for the years of the Padua points, Volume 8 5 Pages 66 74 Application of modified Leja sequences to polynomial interpolation L. P. Bos a M. Caliari a Abstract We discuss several applications
More informationLearning Vector Quantization
Learning Vector Quantization Neural Computation : Lecture 18 John A. Bullinaria, 2015 1. SOM Architecture and Algorithm 2. Vector Quantization 3. The Encoder-Decoder Model 4. Generalized Lloyd Algorithms
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
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 informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationIntroduction to Finite Di erence Methods
Introduction to Finite Di erence Methods ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: Introduction to Finite Di erence Approximation
More informationA Crash Course on Fast Multipole Method (FMM)
A Crash Course on Fast Multipole Method (FMM) Hanliang Guo Kanso Biodynamical System Lab, 2018 April What is FMM? A fast algorithm to study N body interactions (Gravitational, Coulombic interactions )
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationIntroduction to Neural Networks: Structure and Training
Introduction to Neural Networks: Structure and Training Professor Q.J. Zhang Department of Electronics Carleton University, Ottawa, Canada www.doe.carleton.ca/~qjz, qjz@doe.carleton.ca A Quick Illustration
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
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 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 informationGreen s Functions: Taking Another Look at Kernel Approximation, Radial Basis Functions and Splines
Green s Functions: Taking Another Look at Kernel Approximation, Radial Basis Functions and Splines Gregory E. Fasshauer Abstract The theories for radial basis functions (RBFs) as well as piecewise polynomial
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
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 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 informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationFINITE DIFFERENCE METHODS (II): 1D EXAMPLES IN MATLAB
1.723 - COMPUTATIONAL METHODS FOR FLOW IN POROUS MEDIA Spring 2009 FINITE DIFFERENCE METHODS (II): 1D EXAMPLES IN MATLAB Luis Cueto-Felgueroso 1. COMPUTING FINITE DIFFERENCE WEIGHTS The function fdcoefs
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationLearning Vector Quantization (LVQ)
Learning Vector Quantization (LVQ) Introduction to Neural Computation : Guest Lecture 2 John A. Bullinaria, 2007 1. The SOM Architecture and Algorithm 2. What is Vector Quantization? 3. The Encoder-Decoder
More informationThe embedded finite difference method for the Poisson equation in a domain with an irregular boundary and Dirichlet boundary conditions
The embedded finite difference method for the Poisson equation in a domain with an irregular boundary and Dirichlet boundary conditions Z. Jomaa and C. Macaskill School of Mathematics & Statistics University
More informationLeast-Squares Systems and The QR factorization
Least-Squares Systems and The QR factorization Orthogonality Least-squares systems. The Gram-Schmidt and Modified Gram-Schmidt processes. The Householder QR and the Givens QR. Orthogonality The Gram-Schmidt
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 informationPolynomial chaos expansions for sensitivity analysis
c DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for sensitivity analysis B. Sudret Chair of Risk, Safety & Uncertainty
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 5: Completely Monotone and Multiply Monotone Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
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 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 informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More information