An adaptive RBF-based Semi-Lagrangian scheme for HJB equations
|
|
- Lillian Gordon
- 5 years ago
- Views:
Transcription
1 An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal control and related fields, Linz, 22 Nov 2016 joint work with G. Ferretti, O. Junge, A. Schrieber UNIVERSITÀ DEGLI STUDI ROMA TRE Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, UNIVERSITÀ 22 Nov 2016 DEGLI 1 / STUDI 41
2 Numerical Methods in Optimal Control: algorithms, analysis and applications Rome (INdAM), June 19 23, Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
3 Outline 1 Optimal Control via Dynamic Programming 2 Semi-Lagrangian approximation Construction of the SL scheme Convergence of the SL scheme Lagrange vs Least squares space reconstruction operators 3 Shepard/SL scheme The single-level case Full-grid multilevel case Adaptive multilevel case 4 Numerical examples 1D examples 2D example 5 Implementation and efficiency issues 6 Conclusions Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
4 Optimal Control and dynamic programming Controlled dynamical system { ẏ(s) = f (y(s), u(s)), y(t 0 ) = x, the objective of Optimal Control theory is to find and optimal pair (y, u ) minimizing the Cost functional (e.g., of the infinite horizon form) J x (u) = 0 g(y x (s; u), u(s))e λs ds. In Dynamic Prograqmming (DP) techniques, the starting point is to introduce the value function: v(x) = inf u U J x(u) which represents the best possible performance of the system when starting from the state x. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
5 Optimal Control and dynamic programming Controlled dynamical system { ẏ(s) = f (y(s), u(s)), y(t 0 ) = x, the objective of Optimal Control theory is to find and optimal pair (y, u ) minimizing the Cost functional (e.g., of the infinite horizon form) J x (u) = 0 g(y x (s; u), u(s))e λs ds. In Dynamic Prograqmming (DP) techniques, the starting point is to introduce the value function: v(x) = inf u U J x(u) which represents the best possible performance of the system when starting from the state x. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
6 From Dynamic Programming to Hamilton Jacobi Bellman equation DP Principle gives a characterization of value function via the functional equation (which holds for any τ > 0): { τ } v(x) = inf g(y x (s; u), u(s))e λs ds + e λτ v(y x (τ; u)) u U 0 Assuming v differentiable and passing to the limit for τ 0, we obtain the Hamilton Jacobi Bellman (HJB) equation: λv(x) + sup { f (x, u) Dv(x) g(x, u)} = 0. u U In general, in control probelms, v / C 1, and therefore a suitable weak concept of solution (usually, in the viscosity sense) must be used. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
7 Construction of a Semi-Lagrangian scheme for the HJB equation (1) Time discretization Once fixed a time step t, the system evolution is approximated by a one-step scheme, and the integral by a numerical quadrature, e.g., Euler scheme and rectangle quadrature: { } t v(x) = inf u U 0 g(y x (s; u), u(s))e λs ds + e λ t v(y x ( t; u)) v t { (x) = min t g(x, u) + e λ t v t (x + t f (x, u)) }. u U Clearly, more accurate choices are possible. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
8 Construction of a Semi-Lagrangian scheme for the HJB equation (1) Time discretization Once fixed a time step t, the system evolution is approximated by a one-step scheme, and the integral by a numerical quadrature, e.g., Euler scheme and rectangle quadrature: { } t v(x) = inf u U 0 g(y x (s; u), u(s))e λs ds + e λ t v(y x ( t; u)) v t { (x) = min t g(x, u) + e λ t v t (x + t f (x, u)) }. u U Clearly, more accurate choices are possible. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
9 Construction of a Semi-Lagrangian scheme for the HJB equation (2) Space discretization A space grid of nodes {x 1,..., x N } with space step x is set, the semi-discrete equation is collocated at the nodes, and the computation of v is replaced by an interpolation I : v t (x) = min u U { t g(x, u) + e λ t v t (x + t f (x, u)) }. { v j = min t g(xj, u) + e λ t I [V ](x j + t f (x j, u)) }. u U In a similar way, time-dependent HJB equations can also be handled Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
10 Construction of a Semi-Lagrangian scheme for the HJB equation (2) Space discretization A space grid of nodes {x 1,..., x N } with space step x is set, the semi-discrete equation is collocated at the nodes, and the computation of v is replaced by an interpolation I : v t (x) = min u U { t g(x, u) + e λ t v t (x + t f (x, u)) }. { v j = min t g(xj, u) + e λ t I [V ](x j + t f (x j, u)) }. u U In a similar way, time-dependent HJB equations can also be handled Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
11 Accuracy, stability and convergence issues Accuracy: consistency rate follow from the rate of convergence of the elementary building blocks (one-step scheme, quadrature formula, space reconstruction) Monotonicity: this is the most usual stability requirement, sometimes relaxed as quasi-monotonicity: a monotonicity defect of o( t) is allowed Convergence: so far, the most general convergence theory (Barles Souganidis theorem) requires consistency, L stability and monotonicity (or at least quasi-monotonicity) Convergence estimates: for SL schemes and Lipschitz solutions, the typical convergence estimate reads v j v(x j ) C 1 t p + C 2 x t Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
12 Accuracy, stability and convergence issues Accuracy: consistency rate follow from the rate of convergence of the elementary building blocks (one-step scheme, quadrature formula, space reconstruction) Monotonicity: this is the most usual stability requirement, sometimes relaxed as quasi-monotonicity: a monotonicity defect of o( t) is allowed Convergence: so far, the most general convergence theory (Barles Souganidis theorem) requires consistency, L stability and monotonicity (or at least quasi-monotonicity) Convergence estimates: for SL schemes and Lipschitz solutions, the typical convergence estimate reads v j v(x j ) C 1 t p + C 2 x t Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
13 Accuracy, stability and convergence issues Accuracy: consistency rate follow from the rate of convergence of the elementary building blocks (one-step scheme, quadrature formula, space reconstruction) Monotonicity: this is the most usual stability requirement, sometimes relaxed as quasi-monotonicity: a monotonicity defect of o( t) is allowed Convergence: so far, the most general convergence theory (Barles Souganidis theorem) requires consistency, L stability and monotonicity (or at least quasi-monotonicity) Convergence estimates: for SL schemes and Lipschitz solutions, the typical convergence estimate reads v j v(x j ) C 1 t p + C 2 x t Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
14 Lagrange interpolation The basic choice for the space reconstruction operator is the P 1 /Q 1 (piecewise linear/multilinear) Lagrange interpolation I 1 [ ]. In a single space dimension and on a generic interval I j = [x j, x j+1 ], it reads: I 1 [V ](x) = x x j+1 v j + x x j v j+1 x j x j+1 x j+1 x j = x j+1 x v j + x x j x x v j+1 = ψ j (x)v j + ψ j+1 (x)v j+1 For x I j, this is a convex combination of v j and v j+1. This ensures monotonicity and L stability. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
15 Shepard s method A second space reconstruction will be used here, namely Shepard s method. Rather than an interpolation, this is a weighted least squares technique, based on compactly supported, positive Radial Basis Functions (RBFs). Let X = {x 1,, x N } Ω R d be the set of grid nodes Let w(x, y) : Ω Ω R suitable weight functions, such that w(x, y) = 0 if x y 2 is large, that is, w(x, y) = Φ ρ (x y) with Φ ρ : R d R ρ plays the role of a scale parameter for w (e.g., the radius of its support) Shepard s approximation for f I [f ](x) = N i=1 w(x i, x) N j=1 w(x j, x) f (x i) Radial Basis Functions are usually considered as more suitable for local refinements this feature can be exploited for adaptive methods Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
16 Shepard s method A second space reconstruction will be used here, namely Shepard s method. Rather than an interpolation, this is a weighted least squares technique, based on compactly supported, positive Radial Basis Functions (RBFs). Let X = {x 1,, x N } Ω R d be the set of grid nodes Let w(x, y) : Ω Ω R suitable weight functions, such that w(x, y) = 0 if x y 2 is large, that is, w(x, y) = Φ ρ (x y) with Φ ρ : R d R ρ plays the role of a scale parameter for w (e.g., the radius of its support) Shepard s approximation for f I [f ](x) = N i=1 w(x i, x) N j=1 w(x j, x) f (x i) Radial Basis Functions are usually considered as more suitable for local refinements this feature can be exploited for adaptive methods Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
17 Shepard s method A second space reconstruction will be used here, namely Shepard s method. Rather than an interpolation, this is a weighted least squares technique, based on compactly supported, positive Radial Basis Functions (RBFs). Let X = {x 1,, x N } Ω R d be the set of grid nodes Let w(x, y) : Ω Ω R suitable weight functions, such that w(x, y) = 0 if x y 2 is large, that is, w(x, y) = Φ ρ (x y) with Φ ρ : R d R ρ plays the role of a scale parameter for w (e.g., the radius of its support) Shepard s approximation for f I [f ](x) = N i=1 w(x i, x) N j=1 w(x j, x) f (x i) Radial Basis Functions are usually considered as more suitable for local refinements this feature can be exploited for adaptive methods Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
18 Shepard s method A second space reconstruction will be used here, namely Shepard s method. Rather than an interpolation, this is a weighted least squares technique, based on compactly supported, positive Radial Basis Functions (RBFs). Let X = {x 1,, x N } Ω R d be the set of grid nodes Let w(x, y) : Ω Ω R suitable weight functions, such that w(x, y) = 0 if x y 2 is large, that is, w(x, y) = Φ ρ (x y) with Φ ρ : R d R ρ plays the role of a scale parameter for w (e.g., the radius of its support) Shepard s approximation for f I [f ](x) = N i=1 w(x i, x) N j=1 w(x j, x) f (x i) Radial Basis Functions are usually considered as more suitable for local refinements this feature can be exploited for adaptive methods Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
19 Single-level Shepard/SL scheme Shepard s approximation for f I [f ](x) = N i=1 w(x i, x) N j=1 w(x j, x) f (x i) = N ψ i (x)f (x i ) Basis functions ψ( ) are positive and with unitary sum, therefore the reconstruction is monotone By Barles Souganidis theorem, the Shepard implementation of the SL scheme is therefore monotone and convergent The reconstruction error is again O( x) on Lipschitz functions, and as a consequence we obtain again i=1 v j v(x j ) C 1 t p + C 2 x t Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
20 Some remarks on least squares RBF reconstructions (1) In presence of singularities, Shepard s method cuts the smaller scales, thus causing a relevant error with respect to the true value of the function Shepard s reconstruction is exact on constants, but is not in general an interpolation, and may have nonzero residuals on the nodes themselves The residual is concentrated in the neighbourhood of singularities in the first derivatives, and this can indicate the need for refinement f (x, y) = 9 (x 5) 2 Function f and reconstruction Residual at nodes Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
21 Some remarks on least squares RBF reconstructions (1) In presence of singularities, Shepard s method cuts the smaller scales, thus causing a relevant error with respect to the true value of the function Shepard s reconstruction is exact on constants, but is not in general an interpolation, and may have nonzero residuals on the nodes themselves The residual is concentrated in the neighbourhood of singularities in the first derivatives, and this can indicate the need for refinement f (x, y) = 9 (x 5) 2 Function f and reconstruction Residual at nodes Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
22 Some remarks on least squares RBF reconstructions (1) In presence of singularities, Shepard s method cuts the smaller scales, thus causing a relevant error with respect to the true value of the function Shepard s reconstruction is exact on constants, but is not in general an interpolation, and may have nonzero residuals on the nodes themselves The residual is concentrated in the neighbourhood of singularities in the first derivatives, and this can indicate the need for refinement f (x, y) = 9 (x 5) 2 Function f and reconstruction Residual at nodes Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
23 Some remarks on least squares RBF reconstructions (2) Constants: Shepard s method has no numerical error, both at the nodes and between nodes Linear functions: On structured uniform grids, Shepard s method is exact at the nodes, and has a local error expressed as f (x) I [f ](x) C 1 (ρ)f (x) x with C 1 (ρ) decreasing with the scale parameter ρ Quadratic functions: There is a nonzero residual at the nodes, and this introduces a further term in the error: f (x) I [f ](x) C 1 (ρ)f (x) x + C 2 (ρ)f (x) x 2 with C 2 (ρ) increasing with the scale parameter ρ Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
24 Some remarks on least squares RBF reconstructions (2) Constants: Shepard s method has no numerical error, both at the nodes and between nodes Linear functions: On structured uniform grids, Shepard s method is exact at the nodes, and has a local error expressed as f (x) I [f ](x) C 1 (ρ)f (x) x with C 1 (ρ) decreasing with the scale parameter ρ Quadratic functions: There is a nonzero residual at the nodes, and this introduces a further term in the error: f (x) I [f ](x) C 1 (ρ)f (x) x + C 2 (ρ)f (x) x 2 with C 2 (ρ) increasing with the scale parameter ρ Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
25 Some remarks on least squares RBF reconstructions (2) Constants: Shepard s method has no numerical error, both at the nodes and between nodes Linear functions: On structured uniform grids, Shepard s method is exact at the nodes, and has a local error expressed as f (x) I [f ](x) C 1 (ρ)f (x) x with C 1 (ρ) decreasing with the scale parameter ρ Quadratic functions: There is a nonzero residual at the nodes, and this introduces a further term in the error: f (x) I [f ](x) C 1 (ρ)f (x) x + C 2 (ρ)f (x) x 2 with C 2 (ρ) increasing with the scale parameter ρ Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
26 Shepard/SL scheme on the eikonal equation We apply the Shepard/SL scheme to the problem vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Shepard s reconstruction smoothes out the singularity Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
27 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
28 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
29 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
30 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
31 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
32 Remarks Efficiency issues the Shepard/SL scheme is slower than the P 1 /SL scheme At nodes close to the boundary, numerical solution is less accurate Numerical solution is smoothed out at switching points/curves Why choosing adaptivity? A better resolution of singularities may be expected Refinements are avoided in smooth regions (where the scheme is already accurate enough) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
33 Full-grid multilevel approximation Let X 1, X 2,..., X n be a sequence of space grids for the levels 1, 2,..., n, and let w j (x, y) = δ d j Φ((x y)/δ j ) be weight functions with supports of radius δ j (in our case, δ j 2). At each level, the multilevel approximation is constructed working on the residual of the previous levels, each time on a finer grid. Let e 0 = f, and let s j := s ej 1,X j be the approximation of the residual e j 1 at the level j and e j its residual. Then, the multilevel reconstruction reads: e 0 = f f 1 = s 1, e 1 = f s 1, f 2 = s 1 + s 2,. e j 1 = f (s 1 + s s j 1 ), f j = s 1 + s s j. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
34 Full-grid multilevel approximation Let X 1, X 2,..., X n be a sequence of space grids for the levels 1, 2,..., n, and let w j (x, y) = δ d j Φ((x y)/δ j ) be weight functions with supports of radius δ j (in our case, δ j 2). At each level, the multilevel approximation is constructed working on the residual of the previous levels, each time on a finer grid. Let e 0 = f, and let s j := s ej 1,X j be the approximation of the residual e j 1 at the level j and e j its residual. Then, the multilevel reconstruction reads: e 0 = f f 1 = s 1, e 1 = f s 1, f 2 = s 1 + s 2,. e j 1 = f (s 1 + s s j 1 ), f j = s 1 + s s j. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
35 Full-grid multilevel approximation Let X 1, X 2,..., X n be a sequence of space grids for the levels 1, 2,..., n, and let w j (x, y) = δ d j Φ((x y)/δ j ) be weight functions with supports of radius δ j (in our case, δ j 2). At each level, the multilevel approximation is constructed working on the residual of the previous levels, each time on a finer grid. Let e 0 = f, and let s j := s ej 1,X j be the approximation of the residual e j 1 at the level j and e j its residual. Then, the multilevel reconstruction reads: e 0 = f f 1 = s 1, e 1 = f s 1, f 2 = s 1 + s 2,. e j 1 = f (s 1 + s s j 1 ), f j = s 1 + s s j. Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
36 Nested multilevel grids Multilevel 1D grid binary tree Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
37 Adaptive multilevel approximation Space reconstruction is no longer monotone if implemented in multilevel form A convergence result can be proved for full-grid, N-levels approximation, without exploiting Barles Souganidis theorem. The error estimate is of the form v j v(x j ) C 1 t p + C 2 x 2 N 1 t In the adaptive code: At a certain node x j, a new level is created only if e j > τ (with τ a suitable refinement threshold) The threshold τ must discriminate between smooth and nonsmooth regions. For example, defining τ = x 3/2 sets the threshold at an intermediate order of magnitude (Convergence analysis for the adaptive case in progress) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
38 Adaptive multilevel approximation Space reconstruction is no longer monotone if implemented in multilevel form A convergence result can be proved for full-grid, N-levels approximation, without exploiting Barles Souganidis theorem. The error estimate is of the form v j v(x j ) C 1 t p + C 2 x 2 N 1 t In the adaptive code: At a certain node x j, a new level is created only if e j > τ (with τ a suitable refinement threshold) The threshold τ must discriminate between smooth and nonsmooth regions. For example, defining τ = x 3/2 sets the threshold at an intermediate order of magnitude (Convergence analysis for the adaptive case in progress) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
39 Adaptive multilevel approximation Space reconstruction is no longer monotone if implemented in multilevel form A convergence result can be proved for full-grid, N-levels approximation, without exploiting Barles Souganidis theorem. The error estimate is of the form v j v(x j ) C 1 t p + C 2 x 2 N 1 t In the adaptive code: At a certain node x j, a new level is created only if e j > τ (with τ a suitable refinement threshold) The threshold τ must discriminate between smooth and nonsmooth regions. For example, defining τ = x 3/2 sets the threshold at an intermediate order of magnitude (Convergence analysis for the adaptive case in progress) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
40 Adaptive multilevel approximation Space reconstruction is no longer monotone if implemented in multilevel form A convergence result can be proved for full-grid, N-levels approximation, without exploiting Barles Souganidis theorem. The error estimate is of the form v j v(x j ) C 1 t p + C 2 x 2 N 1 t In the adaptive code: At a certain node x j, a new level is created only if e j > τ (with τ a suitable refinement threshold) The threshold τ must discriminate between smooth and nonsmooth regions. For example, defining τ = x 3/2 sets the threshold at an intermediate order of magnitude (Convergence analysis for the adaptive case in progress) Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
41 Pseudo-code (1D case) Require: nodes, τ, n_steps, maximum number of levels N_lev. Ensure: Value function, value, control, feedback. 1: for n = 1,..., n_steps do 2: for each node x do 3: build the first level approximation of value and feedback 4: end for 5: for i = 1,..., N_lev 1 do 6: for each node x at level i do 7: if value(x) SHEPARD(x) > τ then 8: generate 2 nodes at the level i + 1 9: end if 10: end for 11: for each new node x do 12: compute the new approximation of value and feedback 13: compute the space reconstruction of value via the function SHEPARD 14: save in value the new residual 15: end for 16: end for 17: end for Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
42 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
43 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
44 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
45 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
46 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
47 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
48 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
49 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
50 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
51 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
52 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
53 A time-dependent example 1D vx 2 {v t + 2 = 0 v(x, 0) = max(1 (x 2) 2, 0) Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
54 A stationary example 1D Maximization problem: f (x, u) = u µx + mx ρ n ρ + x ρ g(x, u) = au σ kβx 2 with U = [0, 0.4], λ = 0.1, Ω = [0, 4], and constants a = 2, σ = β = k = 1/2, m = n = 1, ρ = 2, µ = Single level Multilevel Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
55 A stationary example 2D (1) This test is a reformulation in terms of an infinite horizon problem of a chemotherapy control problem. The dynamics is given by the compartment model: { ẋ 1 (t) = a 1 x 1 (t) + 2(1 u(t))a 2 x 2 (t) ẋ 2 (t) = a 1 x 1 (t) a 2 x 2 (t), in which the two state compartments denote two successive phases of cellular growth, with the drug u acting only on the second phase. The running cost is g(x, u) = u + r 1 ẋ 1 + r 2 ẋ 2 = u + r 1 ( a 1 x 1 + 2(1 u)a 2 x 2 ) + r 2 (a 1 x 1 a 2 x 2 ) thus aiming at a reduction of the tumour mass, but avoiding an excessive dose of drug. Here, the parameters have been set as U = [0, 1], λ = 0.3, Ω = [0, 1] 2, and, following the literature, a 1 = 0.197, a 2 = 0.356, r 1 = 6.94 and r 2 = Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
56 A stationary example 2D (2) Single level solution Single level solution Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
57 A stationary example 2D (3) Adaptive 3-level solution Adaptive 3-level grid Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
58 Data structure for the multilevel grid (1) Sparse memorization of the multilevel grid: each row of the matrix contains a list of the nodes at the corresponding level Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
59 Data structure for the multilevel grid (2) Full memorization of the multilevel grid: each row of the matrix contains the full grid at the corresponding level Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
60 Efficiency issues For the implementation in form of a multiple (level-by-level) list: Memorization is very compact, but execution times are high Difficult to keep memory of the nonzero basis functions redundant computation For a full memorization of the multilevel grid: Execution times are faster, but large memory occupation Possible intermediate solutions: Node list, with implementation of a dependency tree of the multilevel grid (nonzero basis functions are found via a tree-based splitting of the computational domain) Full grid with a sparse matrix memorization Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
61 Efficiency issues For the implementation in form of a multiple (level-by-level) list: Memorization is very compact, but execution times are high Difficult to keep memory of the nonzero basis functions redundant computation For a full memorization of the multilevel grid: Execution times are faster, but large memory occupation Possible intermediate solutions: Node list, with implementation of a dependency tree of the multilevel grid (nonzero basis functions are found via a tree-based splitting of the computational domain) Full grid with a sparse matrix memorization Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
62 Efficiency issues For the implementation in form of a multiple (level-by-level) list: Memorization is very compact, but execution times are high Difficult to keep memory of the nonzero basis functions redundant computation For a full memorization of the multilevel grid: Execution times are faster, but large memory occupation Possible intermediate solutions: Node list, with implementation of a dependency tree of the multilevel grid (nonzero basis functions are found via a tree-based splitting of the computational domain) Full grid with a sparse matrix memorization Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
63 Conclusions Pros: Cons: Better approximation in nonsmooth regions Good robustness of the refinement indicator At the moment, the adaptive code is less efficient than the single-level, Lagrange-based scheme Complicate dependence on the choice of the basis and of the refinement threshold Future tasks: Improve the sparse memorization to speed up the reconstruction procedure Fully unstructured refinement Convergence analysis for the adaptive case Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
64 Conclusions Pros: Cons: Better approximation in nonsmooth regions Good robustness of the refinement indicator At the moment, the adaptive code is less efficient than the single-level, Lagrange-based scheme Complicate dependence on the choice of the basis and of the refinement threshold Future tasks: Improve the sparse memorization to speed up the reconstruction procedure Fully unstructured refinement Convergence analysis for the adaptive case Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
65 Conclusions Pros: Cons: Better approximation in nonsmooth regions Good robustness of the refinement indicator At the moment, the adaptive code is less efficient than the single-level, Lagrange-based scheme Complicate dependence on the choice of the basis and of the refinement threshold Future tasks: Improve the sparse memorization to speed up the reconstruction procedure Fully unstructured refinement Convergence analysis for the adaptive case Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
66 References M. Falcone, R. Ferretti, Semi-Lagrangian Approximation Schemes for Linear and Hamilton-Jacobi Equations, SIAM (2014). G. Ferretti, R. Ferretti, O. Junge, A. Schreiber, A multilevel radial basis function scheme for the HJB equation, submitted for IFAC O. Junge, A. Schreiber, Dynamic Programming using Radial Basis Functions, Discrete and Continuous Dynamical Systems, A, 35 (2015). U. Ledzewicz H. Schättler, Optimal bang-bang controls for a two-compartment model in cancer chemotherapy, Journal of Optimization Theory and Applications, 114 (2002), H. Wendland, Local polynomial reproduction and moving least squares approximation, IMA J. Num. Anal., 21 (2001), H. Wendland, Multiscale analysis in Sobolev spaces on bounded domains, Numerische Mathematik, 116 (2010), H. Wendland, Scattered data approximation, Cambridge University Press, Cambridge UK (2005). Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
67 Thank you for your attention! Roberto Ferretti (Roma Tre) Adaptive multilevel Shepard for HJB SSCMSE Linz, 22 Nov / 41
The HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationNumerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle
Numerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle Ph.D. course in OPTIMAL CONTROL Emiliano Cristiani (IAC CNR) e.cristiani@iac.cnr.it
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationpacman A classical arcade videogame powered by Hamilton-Jacobi equations Simone Cacace Universita degli Studi Roma Tre
HJ pacman A classical arcade videogame powered by Hamilton-Jacobi equations Simone Cacace Universita degli Studi Roma Tre Control of State Constrained Dynamical Systems 25-29 September 2017, Padova Main
More informationDynamic programming using radial basis functions and Shepard approximations
Dynamic programming using radial basis functions and Shepard approximations Oliver Junge, Alex Schreiber Fakultät für Mathematik Technische Universität München Workshop on algorithms for dynamical systems
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationMean field games and related models
Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:
More informationSuboptimal feedback control of PDEs by solving HJB equations on adaptive sparse grids
Wegelerstraße 6 53115 Bonn Germany phone +49 228 73-3427 fax +49 228 73-7527 www.ins.uni-bonn.de Jochen Garcke and Axel Kröner Suboptimal feedback control of PDEs by solving HJB equations on adaptive sparse
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 informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
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 informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
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 informationA Monotone Scheme for Hamilton-Jacobi Equations via the Nonstandard Finite Difference Method
A Monotone Scheme for Hamilton-Jacobi Equations via the Nonstandard Finite Difference Method Roumen Anguelov, Jean M-S Lubuma and Froduald Minani Department of Mathematics and Applied Mathematics University
More informationConvergence of a large time-step scheme for Mean Curvature Motion
Outline Convergence of a large time-step scheme for Mean Curvature Motion M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma in collaboration with E. Carlini (SAPIENZA) and R. Ferretti
More informationOn the approximation of the principal eigenvalue for a class of nonlinear elliptic operators
On the approximation of the principal eigenvalue for a class of nonlinear elliptic operators Fabio Camilli ("Sapienza" Università di Roma) joint work with I.Birindelli ("Sapienza") I.Capuzzo Dolcetta ("Sapienza")
More informationOptimal Control. Lecture 18. Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline Hamilton-Jacobi-Bellman (HJB) Equation Iterative solution of HJB Equation
More informationarxiv: v1 [math.na] 29 Jul 2018
AN EFFICIENT DP ALGORITHM ON A TREE-STRUCTURE FOR FINITE HORIZON OPTIMAL CONTROL PROBLEMS ALESSANDRO ALLA, MAURIZIO FALCONE, LUCA SALUZZI arxiv:87.8v [math.na] 29 Jul 28 Abstract. The classical Dynamic
More informationA Semi-Lagrangian discretization of non linear fokker Planck equations
A Semi-Lagrangian discretization of non linear fokker Planck equations E. Carlini Università Sapienza di Roma joint works with F.J. Silva RICAM, Linz November 2-25, 206 / 3 Outline A Semi-Lagrangian scheme
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 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 informationHJ equations. Reachability analysis. Optimal control problems
HJ equations. Reachability analysis. Optimal control problems Hasnaa Zidani 1 1 ENSTA Paris-Tech & INRIA-Saclay Graz, 8-11 September 2014 H. Zidani (ENSTA & Inria) HJ equations. Reachability analysis -
More informationDynamic programming with radial basis functions and Shepard s method
Technische Universität München Zentrum Mathematik Wissenschaftliches Rechnen Dynamic programming with radial basis functions and Shepard s method Alex Joachim Ernst Schreiber Vollständiger Abdruck der
More informationZur Konstruktion robust stabilisierender Regler mit Hilfe von mengenorientierten numerischen Verfahren
Zur Konstruktion robust stabilisierender Regler mit Hilfe von mengenorientierten numerischen Verfahren Oliver Junge Zentrum Mathematik Technische Universität München mit Lars Grüne, Bayreuth Outline optimal
More informationValue Function and Optimal Trajectories for some State Constrained Control Problems
Value Function and Optimal Trajectories for some State Constrained Control Problems Hasnaa Zidani ENSTA ParisTech, Univ. of Paris-saclay "Control of State Constrained Dynamical Systems" Università di Padova,
More informationOptimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations
Martino Bardi Italo Capuzzo-Dolcetta Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations Birkhauser Boston Basel Berlin Contents Preface Basic notations xi xv Chapter I. Outline
More informationarxiv: v1 [math.na] 28 Dec 2018
ERROR ESTIMATES FOR A TREE STRUCTURE ALGORITHM SOLVING FINITE HORIZON CONTROL PROBLEMS LUCA SALUZZI, ALESSANDRO ALLA, MAURIZIO FALCONE, arxiv:1812.11194v1 [math.na] 28 Dec 2018 Abstract. In the Dynamic
More informationA Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory
A Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory M. Falcone and S. Finzi Vita Abstract We introduce and analyze a new scheme for the approximation of the two-layer model proposed
More informationThuong Nguyen. SADCO Internal Review Metting
Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)
More informationLipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator
Lipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator Thi Tuyen Nguyen Ph.D student of University of Rennes 1 Joint work with: Prof. E. Chasseigne(University of
More informationPreconditioned space-time boundary element methods for the heat equation
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 Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationLecture 3: Hamilton-Jacobi-Bellman Equations. Distributional Macroeconomics. Benjamin Moll. Part II of ECON Harvard University, Spring
Lecture 3: Hamilton-Jacobi-Bellman Equations Distributional Macroeconomics Part II of ECON 2149 Benjamin Moll Harvard University, Spring 2018 1 Outline 1. Hamilton-Jacobi-Bellman equations in deterministic
More informationComputation of Local ISS Lyapunov Function Via Linear Programming
Computation of Local ISS Lyapunov Function Via Linear Programming Huijuan Li Joint work with Robert Baier, Lars Grüne, Sigurdur F. Hafstein and Fabian Wirth Institut für Mathematik, Universität Würzburg
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
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 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 informationChapter 5. Pontryagin s Minimum Principle (Constrained OCP)
Chapter 5 Pontryagin s Minimum Principle (Constrained OCP) 1 Pontryagin s Minimum Principle Plant: (5-1) u () t U PI: (5-2) Boundary condition: The goal is to find Optimal Control. 2 Pontryagin s Minimum
More informationThe Hopf - Lax solution for state dependent Hamilton - Jacobi equations
The Hopf - Lax solution for state dependent Hamilton - Jacobi equations Italo Capuzzo Dolcetta Dipartimento di Matematica, Università di Roma - La Sapienza 1 Introduction Consider the Hamilton - Jacobi
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationarxiv: v2 [math.na] 20 Jan 2016
Piecewise Constant Policy Approximations to Hamilton-Jacobi-Bellman Equations C. Reisinger and P.A. Forsyth arxiv:1503.05864v2 [math.na] 20 Jan 2016 January 21, 2016 Abstract An advantageous feature of
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 informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationConvergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics
Convergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics Shiming Yang and Matthias K. Gobbert Abstract. The convergence order of finite elements is related to the polynomial order
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 informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
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 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 informationThe Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations
The Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations Iain Smears INRIA Paris Linz, November 2016 joint work with Endre Süli, University of Oxford Overview
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationHigh-order filtered schemes for time-dependent second order HJB equations
High-order filtered schemes for time-dependent second order HJB equations Olivier Bokanowski, Athena Picarelli, and Christoph Reisinger November 15, 2016 Abstract In this paper, we present and analyse
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationConvex computation of the region of attraction for polynomial control systems
Convex computation of the region of attraction for polynomial control systems Didier Henrion LAAS-CNRS Toulouse & CTU Prague Milan Korda EPFL Lausanne Region of Attraction (ROA) ẋ = f (x,u), x(t) X, u(t)
More informationOptimal Control. Macroeconomics II SMU. Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112
Optimal Control Ömer Özak SMU Macroeconomics II Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112 Review of the Theory of Optimal Control Section 1 Review of the Theory of Optimal Control Ömer
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationMean-Field Games with non-convex Hamiltonian
Mean-Field Games with non-convex Hamiltonian Martino Bardi Dipartimento di Matematica "Tullio Levi-Civita" Università di Padova Workshop "Optimal Control and MFGs" Pavia, September 19 21, 2018 Martino
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationCONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS AND FREE BOUNDARY PROBLEMS
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 879 895 c 26 Society for Industrial and Applied Mathematics CONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS
More informationSuboptimal feedback control of PDEs by solving HJB equations on adaptive sparse grids
Suboptimal feedback control of PDEs by solving HJB equations on adaptive sparse grids Jochen Garcke, Axel Kröner To cite this version: Jochen Garcke, Axel Kröner. Suboptimal feedback control of PDEs by
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationNeighboring feasible trajectories in infinite dimension
Neighboring feasible trajectories in infinite dimension Marco Mazzola Université Pierre et Marie Curie (Paris 6) H. Frankowska and E. M. Marchini Control of State Constrained Dynamical Systems Padova,
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. COTROL OPTIM. Vol. 43, o. 4, pp. 1222 1233 c 2004 Society for Industrial and Applied Mathematics OZERO-SUM STOCHASTIC DIFFERETIAL GAMES WITH DISCOTIUOUS FEEDBACK PAOLA MAUCCI Abstract. The existence
More informationAccepted Manuscript. A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
Accepted Manuscript A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang, Hongkai Zhao PII: S-()- DOI:./j.jcp... Reference: YJCPH To
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationMulti-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form
Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form Yipeng Yang * Under the supervision of Dr. Michael Taksar Department of Mathematics University of Missouri-Columbia Oct
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 informationGENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS. Jong Kyu Kim, Salahuddin, and Won Hee Lim
Korean J. Math. 25 (2017), No. 4, pp. 469 481 https://doi.org/10.11568/kjm.2017.25.4.469 GENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS Jong Kyu Kim, Salahuddin, and Won Hee Lim Abstract. In this
More informationCalculating the domain of attraction: Zubov s method and extensions
Calculating the domain of attraction: Zubov s method and extensions Fabio Camilli 1 Lars Grüne 2 Fabian Wirth 3 1 University of L Aquila, Italy 2 University of Bayreuth, Germany 3 Hamilton Institute, NUI
More informationComputational Issues in Nonlinear Dynamics and Control
Computational Issues in Nonlinear Dynamics and Control Arthur J. Krener ajkrener@ucdavis.edu Supported by AFOSR and NSF Typical Problems Numerical Computation of Invariant Manifolds Typical Problems Numerical
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationRegion of attraction approximations for polynomial dynamical systems
Region of attraction approximations for polynomial dynamical systems Milan Korda EPFL Lausanne Didier Henrion LAAS-CNRS Toulouse & CTU Prague Colin N. Jones EPFL Lausanne Region of Attraction (ROA) ẋ(t)
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 informationINFINITE TIME HORIZON OPTIMAL CONTROL OF THE SEMILINEAR HEAT EQUATION
Nonlinear Funct. Anal. & Appl., Vol. 7, No. (22), pp. 69 83 INFINITE TIME HORIZON OPTIMAL CONTROL OF THE SEMILINEAR HEAT EQUATION Mihai Sîrbu Abstract. We consider here the infinite horizon control problem
More informationPseudo-Poincaré Inequalities and Applications to Sobolev Inequalities
Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Laurent Saloff-Coste Abstract Most smoothing procedures are via averaging. Pseudo-Poincaré inequalities give a basic L p -norm control
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationOn the Bellman equation for control problems with exit times and unbounded cost functionals 1
On the Bellman equation for control problems with exit times and unbounded cost functionals 1 Michael Malisoff Department of Mathematics, Hill Center-Busch Campus Rutgers University, 11 Frelinghuysen Road
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
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 informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationRapporto di Ricerca CS
UNIVERSITÀ CA FOSCARI DI VENEZIA Dipartimento di Informatica Technical Report Series in Computer Science Rapporto di Ricerca CS-7-7 Ottobre 7 Flavio Sartoretto Meshless Petrov-Galerkin methods Dipartimento
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
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 informationGLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS
GLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS GUY BARLES, PIERRE CARDALIAGUET, OLIVIER LEY & RÉGIS MONNEAU Abstract. We are interested in nonlocal Eikonal Equations arising in the study
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 informationQUADRATIC MAJORIZATION 1. INTRODUCTION
QUADRATIC MAJORIZATION JAN DE LEEUW 1. INTRODUCTION Majorization methods are used extensively to solve complicated multivariate optimizaton problems. We refer to de Leeuw [1994]; Heiser [1995]; Lange et
More informationConvergence of a Generalized Fast Marching Method for an Eikonal equation with a Velocity Changing Sign
Convergence of a Generalized Fast Marching Method for an Eikonal equation with a Velocity Changing Sign E. Carlini, M. Falcone, N. Forcadel, R. Monneau December 28, 2007 Abstract We present a new Fast
More informationNumerical method for HJB equations. Optimal control problems and differential games (lecture 3/3)
Numerical method for HJB equations. Optimal control problems and differential games (lecture 3/3) Maurizio Falcone (La Sapienza) & Hasnaa Zidani (ENSTA) ANOC, 23 27 April 2012 M. Falcone & H. Zidani ()
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationarxiv: v2 [math.ap] 28 Nov 2016
ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game
More informationThermodynamic limit and phase transitions in non-cooperative games: some mean-field examples
Thermodynamic limit and phase transitions in non-cooperative games: some mean-field examples Conference on the occasion of Giovanni Colombo and Franco Ra... https://events.math.unipd.it/oscgc18/ Optimization,
More information