Efficient Solvers for Stochastic Finite Element Saddle Point Problems
|
|
- Elaine Todd
- 5 years ago
- Views:
Transcription
1 Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 1/3
2 Collaborators: Elisabeth Ullmann, University of Freiberg, Germany. Andrew Gordon, University of Manchester, UK. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 2/3
3 Deterministic PDEs Galerkin mixed finite element approximation for two field PDE problems on D R 2,3, leads to variational problems of the form find u h (x) V h and p h (x) W h such that: a (u h, v) + b(p h, v) = g (v), v V h b(w, u h ) = f (w), w W h and saddle point systems A BT B 0 u p = g f for which the study of efficient solvers is mature. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 3/3
4 PDEs with Random Data Stochastic Galerkin mixed finite element approximation for two field PDE problems with random data, on D Γ, with Γ R M, leads to find u hd (x, ξ) V h S d and p hd (x, ξ) W h S d such that: â (u hd, v) + ˆb(p hd, v) = ĝ (v), v V h S d ˆb(w, uhd ) = ˆf (w), w W h S d and saddle point systems  ˆB T ˆB 0 u p = g f which are much larger and more complicated to solve. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 4/3
5 Model Problem: Darcy Flow T p = u, in D u = 0, in D p = g, on D D T p n = 0, on D N = D\ D D Solution variables: p = p(x) hydraulic head (pressure) u = u(x) velocity field Data: T = T(x) transmissivity coefficients g = g(x) boundary data D R d spatial domain Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 5/3
6 Outline Stochastic Galerkin methods for: T (x, ω) 1 u(x, ω) p (x, ω) = 0, u(x, ω) = 0 Solving single very large saddle point system Stochastic collocation methods for: T (x, ω) 1 u(x, ω) p (x, ω) = 0, u(x, ω) = 0 Solving many small deterministic saddle point systems (see A. Gordon s poster) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 6/3
7 Stochastic PDE Problem Let T(x, ω) : D Ω R be a correlated random field. Approximate T 1 (x, ω) by a function T 1 M (x, ξ) involving a finite number of random variables ξ = (ξ 1,..., ξ M ) taking values in Γ = M i=1 Γ i R M. Find p(x, ξ), u(x, ξ) such that P -a.s. T 1 M (x, ξ)u(x, ξ) p (x, ξ) = 0, u(x, ξ) = 0 in D Γ, p(x, ξ) = g(x) on D D Γ, u(x, ξ) n = 0 on D N Γ. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 7/3
8 Finite-dimensional Noise Assumption Popular choices are: Truncated Karhunen-Loève expansion (linear) T 1 M (x, ξ) = µ(x) + M k=1 t k (x) ξ k, t k (x) = λ k c k (x) {λ k, c k (x)}, k = 1, 2,... are the eigenvalues and eigenfunctions of covariance {ξ 1, ξ 2,...} are uncorrelated, with mean zero and unit variance Truncated Polynomial Chaos expansion (nonlinear) T 1 M (x, ξ(ω)) = µ(x) + N k=1 t k (x) ψ k (ξ), t k (x) = T 1 ψ k (ξ) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 8/3
9 Example - KL expansion Consider the covariance function C T 1(x, y) = exp ( x y 1 ) D = [0, 1] [0, 1] and Gaussian random variables. Realisations of T 1 M with M = 5, 20, 50, with µ(x) = Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 9/3
10 Polynomial Chaos If the random variables are independent then the joint density function has the form: ρ(ξ) = i ρ i (ξ i ) and E [g(ξ)] = g(ξ) = Γ ρ(y)g(y) dy. Orthonormal polynomials in ξ are constructed via ψ i (ξ) = M s=1 ψ is (ξ s ) where ψ is (ξ s ) is a univariate polynomial of degree i s, orthonormal w.r.t to, = E [ ]. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 10/3
11 Weak Formulation V = L 2 ρ(γ, H(div; D)) W = L 2 ρ(γ, L 2 (D)) We seek u(x, ξ) V and p(x, ξ) W such that: ( T 1 M u, v ) + (p, v) = (g, v n) ΓD, (w, u) = 0, v(x, ξ) V and w(x, ξ) W. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 11/3
12 Stochastic Galerkin Equations Find u hd (x, ξ) V h S d and p hd (x, ξ) W h S d satisfying: ( ) T 1 M u hd, v + (p hd, v) = (g, v n) ΓD, (w, u hd ) = 0, v V h S d and w W h S d. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 12/3
13 Stochastic Galerkin Equations Find u hd (x, ξ) V h S d and p hd (x, ξ) W h S d satisfying: ( ) T 1 M u hd, v + (p hd, v) = (g, v n) ΓD, (w, u hd ) = 0, v V h S d and w W h S d. V h H(div; D), W h L 2 (D) are a deterministic inf-sup stable pairing S d L 2 (Γ) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 12/3
14 Total degree polynomials For S d we choose M-variate polynomials of total degree d. Use orthonormal polynomial chaos basis S d = span { } ψ 1 (ξ),..., ψ Nξ (ξ), N ξ = (M + d)! M!d! M = 1 M = 5 M = 10 d = 1 d = 2 d = 3 d = 1 d = 2 d = 3 d = 1 d = 2 d = Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 13/3
15 Saddle point systems V h = span {ϕ i (x)} N u i=1, W h = span {φ j (x)} N p j=1, S d = span {ψ k (ξ)} N ξ k=1 Â ˆB T u = g ˆB 0 p f (never assembled) of dimension N x N ξ N x N ξ where N x = N u + N p and Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 14/3
16 Saddle point systems V h = span {ϕ i (x)} N u i=1, W h = span {φ j (x)} N p j=1, S d = span {ψ k (ξ)} N ξ k=1  ˆB T u = g ˆB 0 p f (never assembled) of dimension N x N ξ N x N ξ where N x = N u + N p and  ˆB T ˆB 0 = I A 0 + N k=1 G k A k I B T I B 0 [A 0 ] ij = D µ(x) ϕ i ϕ j dx, [A k ] ij = D t k (x) ϕ i ϕ j dx. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 14/3
17 Linear (KL) Case T 1 M (x, ξ) = µ(x) + M k=1 t k (x)ξ k N = M (the number of random variables) G k matrices have at most two non-zeros entries per row. [G k ] rs = ξ k ψ r (ξ)ψ s (ξ), k = 1 : M N k=1 G k is sparse so  is block sparse. Matrix-vector products are cheap  O(M(N u N ξ )) ˆB, ˆB T O(N p N ξ ) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 15/3
18 Matrix structure M = 2, d = 2 (left) and M = 4, d = 2 (right) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 16/3
19 Nonlinear (PC) Case T 1 M (x, ξ) = µ(x) + N k=1 t k (x)ψ k N = (M+2d)! M!(2d)! >> N ξ G k matrices have non-trivial sparsity pattern [G k ] rs = ψ k (ξ)ψ r (ξ)ψ s (ξ), k = 1 : N M k=1 G k is dense so  is block dense. Matrix-vector products are expensive  O(N(N u N ξ )) O(N(N u N 2 ξ )) ˆB, ˆB T O(N p N ξ ) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 17/3
20 Matrix structure Want to perform as few matrix vector products with this as possible! Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 18/3
21 (Minres) preconditioning strategies There are two well-known classes of block-diagonal preconditioners. Schur-Complement Preconditioning  0 0 ˆB 1 ˆBT Augmented Preconditioning  + γ 1 ˆBT Ŵ 1 ˆB 0 0 γŵ Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 19/3
22 (Minres) preconditioning strategies There are two well-known classes of block-diagonal preconditioners. Schur-Complement Preconditioning  0 0 ˆB 1 ˆBT Augmented Preconditioning  + γ 1 ˆBT Ŵ 1 ˆB 0 0 γŵ Both require cheap, robust approximations to  Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 20/3
23 Schur-Complement Preconditioning Approximate  L diag(a 0) where L is s.p.d. A practical preconditioner is: P = L diag(a 0) 0 0 ˆB (L diag(a0 )) 1 ˆBT = L diag(a 0) 0 0 L ( Bdiag(A 0 ) 1 B T ). Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 21/3
24 Schur-Complement Preconditioning Approximate  L diag(a 0) where L is s.p.d. A practical preconditioner is: P = L diag(a 0) 0 0 ˆB (L diag(a0 )) 1 ˆBT = L diag(a 0) 0 0 L ( Bdiag(A 0 ) 1 B T ). For the model problem, Bdiag(A 0 ) 1 B T µ and optimal deterministic elliptic PDE solvers (eg. AMG) can be applied. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 21/3
25 Schur-Complement Preconditioning Approximate  L diag(a 0) where L is s.p.d. A practical preconditioner is: P = L diag(a 0) 0 0 ˆB (L diag(a0 )) 1 ˆBT = L diag(a 0) 0 0 L ( Bdiag(A 0 ) 1 B T ). For the model problem, Bdiag(A 0 ) 1 B T µ and optimal deterministic elliptic PDE solvers (eg. AMG) can be applied. Simplest choice: L = I. Approximating P 1 r requires: One solve with a diagonal matrix N ξ V-cycles of multigrid (e.g. AMG) on a deterministic matrix of dimension N p. Cost is O(N ξ (N u + N p )). Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 21/3
26 Convergence analysis P = I diag(a 0) 0 0 I amg(bdiag(a 0 ) 1 B T ) Theorem Let 0 < ˆν 1... ˆν N be the eigenvalues of diag(â) 1 Â. The eigenvalues of the preconditioned saddle-point matrix lie in the union of the bounded intervals, [ ( ) 1 ˆν 1 ˆν Θ2, 1 ( )] ˆν N ˆν 2N 2 + 4θ2 [ ˆν 1, 1 ( )] ˆν N + ˆν 2N 2 + 4Θ2. When L = I, this preconditioner is mean-based. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 22/3
27 Linear (KL) case Theorem E.g. if lowest-order Raviart-Thomas elements are used, the eigenvalues of diag(â) 1 Â lie in the bounded interval [ 1 2 c dτ, c dτ ] where τ = 3σ 2µ M k=1 λk t k (x), σ and µ are the standard deviation and mean of T 1 (x, ω), and c d is a constant (possibly) depending on d. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 23/3
28 Linear (KL) case Theorem E.g. if lowest-order Raviart-Thomas elements are used, the eigenvalues of diag(â) 1 Â lie in the bounded interval [ 1 2 c dτ, c dτ ] where τ = 3σ 2µ M k=1 λk t k (x), σ and µ are the standard deviation and mean of T 1 (x, ω), and c d is a constant (possibly) depending on d. When T 1 M σ µ is a KL expansion this is OK because is small for well-posed problem. c d = O(1) (bounded rvs), c d = O( d) (Gaussian rvs) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 23/3
29 Example - KL case Let D = [0, 1] [0, 1], f = 0 with mixed bcs. We choose a Bessel covariance function for the random input with µ(x) = 1 and λ = 1 M = 6 random variables capture 98% of the variance. u n =0 p=1 p=0 u n =0 Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 24/3
30 Mean of numerical solution Pressure (left), Flux (right) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 25/3
31 Variance of numerical solution Pressure (left), y component (middle) and x component (right) of the Flux x 10 3 x Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 26/3
32 Preconditioned minres iterations d N ξ (N u + N p ) 344,064 1,032,192 2,580,480 5,677,056 σ µ σ µ σ µ = 0.1 total iters N V 1,260 3,864 10,080 22,176 total solve time 14.0s 45.35s s s = 0.2 total iters N V 1,540 4,956 13,020 29,106 total solve time 17.18s 58.51s s = 0.3 total iters N V 1,848 6,216 16,800 39,732 total solve time 20.66s 72.97s s Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 27/3
33 Nonlinear (PC case) Theorem E.g. if lowest-order Raviart-Thomas elements are used, the eigenvalues of diag(â) 1 Â lie in the bounded interval [ 1 2 c dτ, c dτ ] where, for lognormal diffusion coefficients in particular, τ depends nonlinearly on σ G c d = O (exp(dm)) The Schur-complement preconditioner is too weak for large d and σ G requiring an excessive number of matrix-vector products. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 28/3
34 Preconditioned minres iterations Previous test problem. T M (x, ξ) = exp(g M (x, ξ)) where underlying Gaussian field G(x, ω) has Bessel covariance function. M = 5. d N ξ (N u + N p ) 31, , , ,184 N σ G = 0.2 total iters total solve time 2s 22s 245s 2979s σ G = 0.4 total iters total solve time 3s 35s 534s 6321s σ G = 0.6 total iters total solve time 3s 53s 1369s 13794s σ G = 0.8 total iters total solve time 4s 84s 2604s 29808s Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 29/3
35 Augmented Preconditioning Consider the ideal preconditioner P = Â + γ 1 ˆBT Ŵ 1 ˆB 0 0 γŵ Theorem Let Ŵ = I M, where w T (I M) w = w 2 L 2 (Γ,L 2. The eigenvalues of the (D)) preconditioned saddle-point matrix are bounded and lie in ( 1, ˆβ 2 t min γ + ˆβ 2 t min ] {1} where ˆβ is the inf-sup constant and T M t min a.e. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 30/3
36 Augmented Preconditioning Suppose we approximate  by L A 0 where L is any s.p.d N ξ N ξ matrix and replace Ŵ with L 1 M then, P = L ( A 0 + γ 1 B T M 1 B ) 0 0 L 1 γm Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 31/3
37 Augmented Preconditioning Suppose we approximate  by L A 0 where L is any s.p.d N ξ N ξ matrix and replace Ŵ with L 1 M then, P = L ( A 0 + γ 1 B T M 1 B ) 0 0 L 1 γm For this model problem, v T ( A 0 + γ 1 B T M 1 B ) v = ( T 1 M v, v ) + γ 1 ( v, v). The matrix in red is a discretisation of a weighted H(div) operator on the velocity space. Optimal deterministic solvers (geometric multigrid) can be exploited. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 31/3
38 Augmented Preconditioning Simplest choice: L = I ( mean-based preconditioner) P = L ( A 0 + γ 1 B T M 1 B ) 0 0 L γm Approximating P 1 r requires: N ξ solves with M and N p multiplications with L. N ξ V-cycles of multigrid (e.g. Arnold-Falk-Winther) on a deterministic matrix of dimension N u and N u solves with L. Even if L is dense, still cheaper than a matrix-vector product with saddle point matrix! Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 32/3
39 Convergence analysis P = L ( A 0 + γ 1 B T M 1 B ) 0 0 L γm Spectral bounds for preconditioned system matrix depend on choice of L and the efficiency of the approximation  L A 0 and γ. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 33/3
40 Convergence analysis P = L ( A 0 + γ 1 B T M 1 B ) 0 0 L γm Spectral bounds for preconditioned system matrix depend on choice of L and the efficiency of the approximation  L A 0 and γ. However, asymptotically, as γ 0, we can show all N ξ N p negative eigenvalues cluster at 1. N ξ N p positive eigenvalues cluster at +1 N ξ (N u N p ) positive eigenvalues lie in [ˆν 1, ˆν n ] where ˆν 1 v T Âv v T (L A 0 ) v ˆν n Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 33/3
41 Preconditioned minres iterations Previous test problem. Now T M (x, ξ) = exp(g M (x, ξ)) where underlying Gaussian field G(x, ω) has Bessel covariance function. M = 5. Choose L = I and γ = d N ξ (N u + N p ) 31, , , ,184 N σ G = 0.2 total iters total solve time 3s 9s 65s 399s σ G = 0.4 total iters total solve time 3s 12s 97s 668s σ G = 0.6 total iters total solve time 4s 17s 140s 1063s σ G = 0.8 total iters total solve time 5s 19s 203s 1658s Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 34/3
42 Improved choice of L Could use Pitsianis-Van Loan so-called best Kronecker product approximation : L = argmin{h R N ξ N ξ : Â H A 0 F } L is the symmetric and positive definite matrix L = I + N k=1 tr(a T k A 0) tr(a T 0 A 0) G k Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 35/3
43 Preconditioned minres iterations Previous test problem. Now T M (x, ξ) = exp(g M (x, ξ)) where underlying Gaussian field G(x, ω) has Bessel covariance function. M = 5. d N ξ (N u + N p ) 31, , , ,184 N σ G = 0.2 total iters total solve time 3s 8s 54s 328s σ G = 0.4 total iters total solve time 3s 11s 73s 429s σ G = 0.6 total iters total solve time 4s 12s 86s 576s σ G = 0.8 total iters total solve time 4s 14s 106s 723s Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 36/3
44 Average time in seconds per minres iteration and (av. mat-vec time + av. preconditioner time) d=2 d=3 d=4 L = I 0.58 ( ) 4.00 ( ) ( ) Improved L 0.59 ( ) 4.08 ( ) ( ) Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 37/3
45 Summary If random inputs are KL expansions: Mean-based preconditioners are cheap and robust within the range of statistical parameters where the problem is well-posed. Schur-complement preconditioners are OK as matrix-vector products are cheap. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 38/3
46 Summary If random inputs are KL expansions: Mean-based preconditioners are cheap and robust within the range of statistical parameters where the problem is well-posed. Schur-complement preconditioners are OK as matrix-vector products are cheap. If the random inputs are PC expansions Mean-based preconditioners are still cheap but are inadequate as σ and d increase. Augmented preconditioners are promising because fewer matrix-vector products are required Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 38/3
47 Summary If random inputs are KL expansions: Mean-based preconditioners are cheap and robust within the range of statistical parameters where the problem is well-posed. Schur-complement preconditioners are OK as matrix-vector products are cheap. If the random inputs are PC expansions Mean-based preconditioners are still cheap but are inadequate as σ and d increase. Augmented preconditioners are promising because fewer matrix-vector products are required But: Also need to look at stopping criteria for iteration. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 38/3
48 References C.E.Powell, E. Ullmann, Preconditioning stochastic Galerkin saddle point systems, Submitted. H. Elman, D. Furnival, C.E. Powell, H(div) preconditioning for a mixed finite element formulation of the stochastic diffusion equation. Math. Comp O.G. Ernst, C.E. Powell, D. Silvester and E. Ullmann, Efficient solvers for a Linear Stochastic Galerkin Mixed Formulation of the Steady-State Diffusion Equation, SIAM J. Sci. Comput C.E. Powell, H. Elman, Block-diagonal preconditioning for spectral stochastic finite element systems, IMA J. Numer Anal. 2008/9. Efficient Solvers for Stochastic Finite Element Saddle Point Problems p. 39/3
A Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationH(div) Preconditioning for a Mixed Finite Element Formulation of the Stochastic Diffusion Problem 1
University of Maryland Department of Computer Science CS-TR-4918 University of Maryland Institute for Advanced Computer Studies UMIACS-TR-2008-15 H(div) Preconditioning for a Mixed Finite Element Formulation
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationBlock-diagonal preconditioning for spectral stochastic finite-element systems
IMA Journal of Numerical Analysis (2009) 29, 350 375 doi:0.093/imanum/drn04 Advance Access publication on April 4, 2008 Block-diagonal preconditioning for spectral stochastic finite-element systems CATHERINE
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 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 informationarxiv: v2 [math.na] 8 Apr 2017
A LOW-RANK MULTIGRID METHOD FOR THE STOCHASTIC STEADY-STATE DIFFUSION PROBLEM HOWARD C. ELMAN AND TENGFEI SU arxiv:1612.05496v2 [math.na] 8 Apr 2017 Abstract. We study a multigrid method for solving large
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA
International Journal for Uncertainty Quantification, 11):19 33, 2011 ASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA Howard C. Elman, 1, Christopher W.
More informationarxiv: v1 [math.na] 18 Dec 2017
A Bramble-Pasciak conjugate gradient method for discrete Stokes problems with lognormal random viscosity arxiv:1712.06472v1 [math.na] 18 Dec 2017 Christopher Müller 1,2, Sebastian Ullmann 1,2, and Jens
More informationUncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media
Uncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media M.Köppel, C.Rohde Institute for Applied Analysis and Numerical Simulation Inria, Nov 15th, 2016 Porous Media Examples: sponge,
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationSolving the stochastic steady-state diffusion problem using multigrid
IMA Journal of Numerical Analysis (2007) 27, 675 688 doi:10.1093/imanum/drm006 Advance Access publication on April 9, 2007 Solving the stochastic steady-state diffusion problem using multigrid HOWARD ELMAN
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 informationCBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods. Crowder, Adam J and Powell, Catherine E
CBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods Crowder, Adam J and Powell, Catherine E 2017 MIMS EPrint: 2017.18 Manchester Institute for Mathematical Sciences
More informationNumerical Solution for Random Forced SPDE via Galerkin Finite Element Method
Journal of mathematics and computer science 9 (014), 71-8 umerical Solution for Random Forced SE via Galerkin Finite Element ethod Rasoul aseri, A. alek 1 epartment of Applied athematics, Faculty of athematical
More informationCollocation based high dimensional model representation for stochastic partial differential equations
Collocation based high dimensional model representation for stochastic partial differential equations S Adhikari 1 1 Swansea University, UK ECCM 2010: IV European Conference on Computational Mechanics,
More informationIterative solvers for saddle point algebraic linear systems: tools of the trade. V. Simoncini
Iterative solvers for saddle point algebraic linear systems: tools of the trade V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The problem
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 1 / 23 Lecture outline
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
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 informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationAdaptive Collocation with Kernel Density Estimation
Examples of with Kernel Density Estimation Howard C. Elman Department of Computer Science University of Maryland at College Park Christopher W. Miller Applied Mathematics and Scientific Computing Program
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationNumerical Approximation of Stochastic Elliptic Partial Differential Equations
Numerical Approximation of Stochastic Elliptic Partial Differential Equations Hermann G. Matthies, Andreas Keese Institut für Wissenschaftliches Rechnen Technische Universität Braunschweig wire@tu-bs.de
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 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 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 informationANALYSIS OF AUGMENTED LAGRANGIAN-BASED PRECONDITIONERS FOR THE STEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
ANALYSIS OF AUGMENTED LAGRANGIAN-BASED PRECONDITIONERS FOR THE STEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS MICHELE BENZI AND ZHEN WANG Abstract. We analyze a class of modified augmented Lagrangian-based
More informationQuantifying Uncertainty: Modern Computational Representation of Probability and Applications
Quantifying Uncertainty: Modern Computational Representation of Probability and Applications Hermann G. Matthies with Andreas Keese Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de
More informationScalable kernel methods and their use in black-box optimization
with derivatives Scalable kernel methods and their use in black-box optimization David Eriksson Center for Applied Mathematics Cornell University dme65@cornell.edu November 9, 2018 1 2 3 4 1/37 with derivatives
More informationSolution Strategies for Stochastic Finite Element Discretizations
Solution Strategies for Stochastic Finite Element Discretizations Von der Faultät für Mathemati und Informati der Technischen Universität Bergaademie Freiberg genehmigte DISSERTATION zur Erlangung des
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More informationEFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS
EFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS DONGBIN XIU AND JIE SHEN Abstract. We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We
More informationChebyshev semi-iteration in Preconditioning
Report no. 08/14 Chebyshev semi-iteration in Preconditioning Andrew J. Wathen Oxford University Computing Laboratory Tyrone Rees Oxford University Computing Laboratory Dedicated to Victor Pereyra on his
More informationSparse polynomial chaos expansions in engineering applications
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Sparse polynomial chaos expansions in engineering applications B. Sudret G. Blatman (EDF R&D,
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationWhat Should We Do with Radioactive Waste?
What Should We Do with Radioactive Waste? Andrew Cliffe School of Mathematical Sciences University of Nottingham Simulation of Flow in Porous Media and Applications in Waste Management and CO 2 Sequestration
More informationEfficient Solvers for the Navier Stokes Equations in Rotation Form
Efficient Solvers for the Navier Stokes Equations in Rotation Form Computer Research Institute Seminar Purdue University March 4, 2005 Michele Benzi Emory University Atlanta, GA Thanks to: NSF (MPS/Computational
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 informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationSampling and Low-Rank Tensor Approximations
Sampling and Low-Rank Tensor Approximations Hermann G. Matthies Alexander Litvinenko, Tarek A. El-Moshely +, Brunswick, Germany + MIT, Cambridge, MA, USA wire@tu-bs.de http://www.wire.tu-bs.de $Id: 2_Sydney-MCQMC.tex,v.3
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationarxiv: v2 [math.na] 8 Sep 2017
arxiv:1704.06339v [math.na] 8 Sep 017 A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations Juan Galvis O. Andrés Cuervo September 3, 018 Abstract We use a Monte
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationNonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed
Nonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed Eveline Rosseel Herbert De Gersem Stefan Vandewalle Report TW 541, July 29 Katholieke
More informationUncertainty analysis of large-scale systems using domain decomposition
Center for Turbulence Research Annual Research Briefs 2007 143 Uncertainty analysis of large-scale systems using domain decomposition By D. Ghosh, C. Farhat AND P. Avery 1. Motivation and objectives A
More informationOn the interplay between discretization and preconditioning of Krylov subspace methods
On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
More informationc 2011 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 33, No. 5, pp. 2761 2784 c 2011 Society for Industrial and Applied Mathematics ANALYSIS OF AUGMENTED LAGRANGIAN-BASED PRECONDITIONERS FOR THE STEADY INCOMPRESSIBLE NAVIER STOKES
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationMultigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients
Multigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Random fields
More informationIncomplete Cholesky preconditioners that exploit the low-rank property
anapov@ulb.ac.be ; http://homepages.ulb.ac.be/ anapov/ 1 / 35 Incomplete Cholesky preconditioners that exploit the low-rank property (theory and practice) Artem Napov Service de Métrologie Nucléaire, Université
More informationPARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 569 PARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD Riki Honda
More informationAMS Mathematics Subject Classification : 65F10,65F50. Key words and phrases: ILUS factorization, preconditioning, Schur complement, 1.
J. Appl. Math. & Computing Vol. 15(2004), No. 1, pp. 299-312 BILUS: A BLOCK VERSION OF ILUS FACTORIZATION DAVOD KHOJASTEH SALKUYEH AND FAEZEH TOUTOUNIAN Abstract. ILUS factorization has many desirable
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 informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
More informationNon-Intrusive Solution of Stochastic and Parametric Equations
Non-Intrusive Solution of Stochastic and Parametric Equations Hermann G. Matthies a Loïc Giraldi b, Alexander Litvinenko c, Dishi Liu d, and Anthony Nouy b a,, Brunswick, Germany b École Centrale de Nantes,
More informationQuasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient
Quasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient F. Nobile, L. Tamellini, R. Tempone, F. Tesei CSQI - MATHICSE, EPFL, Switzerland Dipartimento di
More informationUNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA
UNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationSTOCHASTIC SAMPLING METHODS
STOCHASTIC SAMPLING METHODS APPROXIMATING QUANTITIES OF INTEREST USING SAMPLING METHODS Recall that quantities of interest often require the evaluation of stochastic integrals of functions of the solutions
More informationFEM and sparse linear system solving
FEM & sparse linear system solving, Lecture 9, Nov 19, 2017 1/36 Lecture 9, Nov 17, 2017: Krylov space methods http://people.inf.ethz.ch/arbenz/fem17 Peter Arbenz Computer Science Department, ETH Zürich
More informationLOW-RANK TENSOR METHODS for SIMULATION, OPTIMIZATION and UNCERTAINTY QUANTIFICATION of PARAMETRIC PDEs
LOW-RANK TENSOR METHODS for SIMULATION, OPTIMIZATION and UNCERTAINTY QUANTIFICATION of PARAMETRIC PDEs Peter Benner Joint work with Sergey Dolgov (U Bath), Martin Stoll (TU Chemnitz) and Akwum Onwunta
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationSpectral Properties of Saddle Point Linear Systems and Relations to Iterative Solvers Part I: Spectral Properties. V. Simoncini
Spectral Properties of Saddle Point Linear Systems and Relations to Iterative Solvers Part I: Spectral Properties V. Simoncini Dipartimento di Matematica, Università di ologna valeria@dm.unibo.it 1 Outline
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
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 informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
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 informationPreconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners
Preconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, Georgia, USA
More informationUncertainty Quantification: Does it need efficient linear algebra?
USA trip October 2014 p. 1/62 Uncertainty Quantification: Does it need efficient linear algebra? David Silvester University of Manchester, UK Catherine Powell University of Manchester,UK Yes. USA trip
More informationA Chebyshev-based two-stage iterative method as an alternative to the direct solution of linear systems
A Chebyshev-based two-stage iterative method as an alternative to the direct solution of linear systems Mario Arioli m.arioli@rl.ac.uk CCLRC-Rutherford Appleton Laboratory with Daniel Ruiz (E.N.S.E.E.I.H.T)
More information4.6 Iterative Solvers for Linear Systems
4.6 Iterative Solvers for Linear Systems Why use iterative methods? Virtually all direct methods for solving Ax = b require O(n 3 ) floating point operations. In practical applications the matrix A often
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
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 informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationBook of Abstracts - Extract th Annual Meeting. March 23-27, 2015 Lecce, Italy. jahrestagung.gamm-ev.de
GESELLSCHAFT für ANGEWANDTE MATHEMATIK und MECHANIK e.v. INTERNATIONAL ASSOCIATION of APPLIED MATHEMATICS and MECHANICS 86 th Annual Meeting of the International Association of Applied Mathematics and
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 informationApproximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity
Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Marcus Sarkis Jointly work with Juan Galvis USC09 Marcus Sarkis (WPI) SPDE without Uniform Ellipticity USC09
More informationDynamic response of structures with uncertain properties
Dynamic response of structures with uncertain properties S. Adhikari 1 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Bay Campus, Fabian Way, Swansea, SA1 8EN, UK International
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationComputational methods for large-scale matrix equations: recent advances and applications. V. Simoncini
Computational methods for large-scale matrix equations: recent advances and applications V. Simoncini Dipartimento di Matematica Alma Mater Studiorum - Università di Bologna valeria.simoncini@unibo.it
More information