Smoothers for ecient multigrid methods in IGA
|
|
- Arleen Gilmore
- 5 years ago
- Views:
Transcription
1 Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third author) and J3362-N25 (second author).
2 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
3 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
4 Poisson mode probem Domain Ω R d Given function f L 2 (Ω) Find u H 1 (Ω) such that u = f u n = 0 in Ω on Ω
5 Variationa formuation Domain Ω R d Given function f L 2 (Ω) Find u H 1 (Ω) such that for a ũ H 1 (Ω). (u, ũ) H 1 (Ω) = (f, ũ) L 2 (Ω)
6 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
7 Spine spaces in one dimension Let Ω = (0, 1)
8 Spine spaces in one dimension Let Ω = (0, 1) M is the uniform subdivision of Ω = (0, 1) in n = n 0 2 subintervas T i of ength h := h 0 2 for = 0, 1, 2,...
9 Spine spaces in one dimension Let Ω = (0, 1) M is the uniform subdivision of Ω = (0, 1) in n = n 0 2 subintervas T i of ength h := h 0 2 for = 0, 1, 2,... Spine space S p,k, (Ω) is the space of a spine functions in C k (Ω), which are piecewise poynomias of degree p on each subinterva in M.
10 Spine spaces in one dimension Let Ω = (0, 1) M is the uniform subdivision of Ω = (0, 1) in n = n 0 2 subintervas T i of ength h := h 0 2 for = 0, 1, 2,... Spine space S p,k, (Ω) is the space of a spine functions in C k (Ω), which are piecewise poynomias of degree p on each subinterva in M. Maximum smoothness: S p, (Ω) := S p,p 1, (Ω)
11 Spine spaces in one dimension Let Ω = (0, 1) M is the uniform subdivision of Ω = (0, 1) in n = n 0 2 subintervas T i of ength h := h 0 2 for = 0, 1, 2,... Spine space S p,k, (Ω) is the space of a spine functions in C k (Ω), which are piecewise poynomias of degree p on each subinterva in M. Maximum smoothness: S p, (Ω) := S p,p 1, (Ω) Standard B-spine basis: φ (1) p, (x),..., φ(m ) p, (x)
12 Spine spaces in two and more dimensions Let Ω = (0, 1) d
13 Spine spaces in two and more dimensions Let Ω = (0, 1) d Tensor product B-spines: ϕ (i+m j) p, (x, y) = φ (i) p, (x)φ(j) p, (y)
14 Spine spaces in two and more dimensions Let Ω = (0, 1) d Tensor product B-spines: ϕ (i+m j) p, (x, y) = φ (i) p, (x)φ(j) p, (y) For Ω = (0, 1) d with d > 1: S p, (Ω) denotes tensor product spine space
15 Spine spaces in two and more dimensions Let Ω = (0, 1) d Tensor product B-spines: ϕ (i+m j) p, (x, y) = φ (i) p, (x)φ(j) p, (y) For Ω = (0, 1) d with d > 1: S p, (Ω) denotes tensor product spine space More genera domains: geometry mapping
16 Spine spaces in two and more dimensions Let Ω = (0, 1) d Tensor product B-spines: ϕ (i+m j) p, (x, y) = φ (i) p, (x)φ(j) p, (y) For Ω = (0, 1) d with d > 1: S p, (Ω) denotes tensor product spine space More genera domains: geometry mapping For reguar geometry mappings: mutigrid for parameter domain can be used as preconditioner
17 Discretization Variationa formuation: Find u S p, (Ω) such that for a ũ S p, (Ω) (u, ũ ) H 1 (Ω) = (f, ũ ) L 2 (Ω)
18 Discretization Variationa formuation: Find u S p, (Ω) such that for a ũ S p, (Ω) (u, ũ ) H 1 (Ω) = (f, ũ ) L 2 (Ω) Matrix-vector notation: K u = f
19 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
20 Mutigrid method One step of the mutigrid method on grid eve k appied to iterate u (0,0) = u (0) and right-hand-side f to obtain x (1) is given by: Appy ν smoothing steps u (0,m) for m = 1,..., ν. = u (0,m 1) + τ 1 K (f K u (0,m 1) )
21 Mutigrid method One step of the mutigrid method on grid eve k appied to iterate u (0,0) = u (0) and right-hand-side f to obtain x (1) is given by: Appy ν smoothing steps u (0,m) for m = 1,..., ν. = u (0,m 1) + τ Appy coarse-grid correction 1 K (f K u (0,m 1) ) Compute defect and restrict to coarser grid Sove probem on coarser grid Proongate and add resut
22 Mutigrid method One step of the mutigrid method on grid eve k appied to iterate u (0,0) = u (0) and right-hand-side f to obtain x (1) is given by: Appy ν smoothing steps u (0,m) for m = 1,..., ν. = u (0,m 1) + τ Appy coarse-grid correction 1 K (f K u (0,m 1) ) Compute defect and restrict to coarser grid Sove probem on coarser grid Proongate and add resut If reaized exacty (two-grid method): u (1) = u (0,ν) + I 1 K 1 1 I 1 (f K k u (0,ν) ) Standard arguments: convergence of two-grid method convergence of the mutigrid method (W-cyce)
23 Intergrid transfer Nested spaces: S p, 1 (Ω) S p, (Ω)
24 Intergrid transfer Nested spaces: S p, 1 (Ω) S p, (Ω) The proongation I 1 is the canonica embedding Knot insertion agorithm
25 Intergrid transfer Nested spaces: S p, 1 (Ω) S p, (Ω) The proongation I 1 is the canonica embedding Knot insertion agorithm The restriction is its transpose: I 1 = (I 1 )T
26 Hackbusch-ike convergence anaysis Convergence of two-grid method with rate q, i.e., u u(1) L q u u(0) L,
27 Hackbusch-ike convergence anaysis Convergence of two-grid method with rate q, i.e., u u(1) L q u u(0) L, in matrix notation T S ν L q,
28 Hackbusch-ike convergence anaysis Convergence of two-grid method with rate q = C AC S ν, i.e., u u(1) L q u u(0) L, in matrix notation T S ν L q, is guaranteed by L 1/2 T K 1 L 1/2 C A (approximation property) L 1/2 K S ν L 1/2 C S ν 1 (smoothing property)
29 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
30 Basis-(in)dependent smoother Richardson smoother: i.e., K := I u (0,m) = u (0,m 1) + τ (f K u (0,m 1) ),
31 Basis-independent smoother Richardson smoother (in continuous understanding): u (0,m) = u (0,m 1) + τ where K := L is the Riesz isomorphism 1 K (f K u (0,m 1) ),
32 Basis-independent smoother Richardson smoother (in continuous understanding): u (0,m) = u (0,m 1) + τ where K := L is the Riesz isomorphism 1 K (f K u (0,m 1) ), Lemma Assume that K := L and that τ is chosen such that 0 < τ 1 L 1/2 K L 1/2. Then the smoothing property is satised for C S = τ 1. Proof: Standard eigenvaue anaysis.
33 Basis-independent smoother Richardson smoother (in continuous understanding): u (0,m) = u (0,m 1) + τ where K := L is the Riesz isomorphism 1 K (f K u (0,m 1) ), Lemma Assume that K := L and that τ is chosen such that 0 < τ 1 L 1/2 K L 1/2. Then the smoothing property is satised for C S = τ 1. Proof: Standard eigenvaue anaysis. Convergence rate: q = C AC S ν = C Aτ 1 ν
34 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
35 Basis-independent smoother for one dimension Cassica anaysis: In L 2 (Ω)
36 Basis-independent smoother for one dimension Cassica anaysis: In L 2 (Ω) L = h 2 M is the propery scaed mass matrix
37 Basis-independent smoother for one dimension Cassica anaysis: In L 2 (Ω) L = h 2 M is the propery scaed mass matrix Proofs for standard smoothers (Richardson, Jacobi, Gauss-Seide) use that the mass matrix is spectray equivaent to its diagona.
38 Basis-independent smoother for one dimension Cassica anaysis: In L 2 (Ω) L = h 2 M is the propery scaed mass matrix Proofs for standard smoothers (Richardson, Jacobi, Gauss-Seide) use that the mass matrix is spectray equivaent to its diagona. This is true for B-spines, but deteriorate for increased p.
39 Basis-independent smoother for one dimension Cassica anaysis: In L 2 (Ω) L = h 2 M is the propery scaed mass matrix Proofs for standard smoothers (Richardson, Jacobi, Gauss-Seide) use that the mass matrix is spectray equivaent to its diagona. This is true for B-spines, but deteriorate for increased p. Use K := L := h 2 M (mass-richardson smoother)
40 Mass-Richardson smoother Use K := L := h 2 M
41 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) τ 1/2.
42 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) Need: robust inverse inequaity for spine space τ 1/2.
43 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) Need: robust inverse inequaity for spine space Loca Fourier anaysis suggests: Such a robust inverse inequaity hods τ 1/2.
44 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) Need: robust inverse inequaity for spine space Loca Fourier anaysis suggests: Such a robust inverse inequaity hods Numerica experiments [Hofreither, Zuehner 2014]: This approach does not work we τ 1/2.
45 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) Need: robust inverse inequaity for spine space Loca Fourier anaysis suggests: Such a robust inverse inequaity hods Numerica experiments [Hofreither, Zuehner 2014]: This approach does not work we τ 1/2. Choose u p, (x) := φ (1) p, (x) = max{0, h x} p and obtain: Such a robust inverse inequaity does not hod
46 Mass-Richardson smoother Use K := L := h 2 M Condition for choice of τ reads as: sup u R m \{0} and C S = τ 1 u K h 1 u M = sup u S p, (Ω)\{0} u H 1 (Ω) h 1 u L 2 (Ω) Need: robust inverse inequaity for spine space Loca Fourier anaysis suggests: Such a robust inverse inequaity hods Numerica experiments [Hofreither, Zuehner 2014]: This approach does not work we τ 1/2. Choose u p, (x) := φ (1) p, (x) = max{0, h x} p and obtain: Such a robust inverse inequaity does not hod sup... p τ p 2 ν p 2
47 Inverse inequaity Theorem ([T., Takacs 2015]) For a N 0 and p N, u H 1 2 3h 1 u L 2 is satised for a u S p, (Ω), where S p, (0, 1) is the space of a u S p, (0, 1) whose odd derivatives vanish at the boundary: 2+1 x 2+1 u (0) = 2+1 x 2+1 u (1) = 0 for a N 0 with < p.
48 Inverse inequaity Theorem ([T., Takacs 2015]) For a N 0 and p N, u H 1 2 3h 1 u L 2 is satised for a u S p, (Ω), where S p, (0, 1) is the space of a u S p, (0, 1) whose odd derivatives vanish at the boundary: 2+1 x 2+1 u (0) = 2+1 x 2+1 u (1) = 0 for a N 0 with < p. Note: S p, (Ω) has ony p (for p even) or p 1 (for p odd) dimensions.
49 Idea: boundary correction Theorem ([Hofreither, T., Zuehner 2015]) For a N 0 and p N, u H 1 2( )h 1 ( u 2 L 2 + H Γ, u,γ 2 H 1 ) 1/2 is satised for a u = u Γ, + u I, S p, (Ω), where u I, S (I ) p, (Ω) S p, (Ω), u Γ, (S (I ) p, (Ω)) and H Γ, is the discrete harmonic extension. Proof: need inverse inequaity on S p, (Ω) and robust approximation error estimate aso on S p, (Ω) (cf. [T., Takacs 2015]).
50 Boundary correction Reorder variabes such that: ( KΓΓ, KI T K = Γ, K I Γ, K II, ), u = ( uγ, u I, )
51 Boundary correction Reorder variabes such that: ( KΓΓ, KI T K = Γ, K I Γ, K II, ), u = ( uγ, u I, Compute the discrete Harmonic extension: ( ) u H Γ, u Γ, H 1 (Ω) = Γ, K 1 II, K = u Γ, I Γ,u KΓΓ, K T Γ, I Γ, K 1 II, K. I Γ, K )
52 Boundary correction Reorder variabes such that: ( KΓΓ, KI T K = Γ, K I Γ, K II, ), u = ( uγ, u I, Compute the discrete Harmonic extension: ( ) u H Γ, u Γ, H 1 (Ω) = Γ, K 1 II, K = u Γ, I Γ,u KΓΓ, K T Γ, I Γ, K 1 II, K. I Γ, K Generaized inverse inequaity: u K 2( ) u L ) with ( K := L := h 2 KΓΓ, KI T M + Γ, K 1 ) II, K I Γ, }{{} K :=
53 Boundary correction Reorder variabes such that: ( KΓΓ, KI T K = Γ, K I Γ, K II, ), u = ( uγ, u I, Compute the discrete Harmonic extension: ( ) u H Γ, u Γ, H 1 (Ω) = Γ, K 1 II, K = u Γ, I Γ,u KΓΓ, K T Γ, I Γ, K 1 II, K. I Γ, K Generaized inverse inequaity: with L 1/2 K L 1/2 2( ) 2 ( K := L := h 2 KΓΓ, KI T M + Γ, K 1 ) II, K I Γ, }{{} K := )
54 Boundary correction Reorder variabes such that: ( KΓΓ, KI T K = Γ, K I Γ, K II, ), u = ( uγ, u I, Compute the discrete Harmonic extension: ( ) u H Γ, u Γ, H 1 (Ω) = Γ, K 1 II, K = u Γ, I Γ,u KΓΓ, K T Γ, I Γ, K 1 II, K. I Γ, K Generaized inverse inequaity: with L 1/2 K L 1/2 2( ) 2 ( K := L := h 2 KΓΓ, KI T M + Γ, K 1 ) II, K I Γ, }{{} K := Smoothing property for 0 < τ 2 1 ( ) 2 )
55 Robust approximation error estimate Theorem ([T., Takacs 2015]) For each u H 1 (Ω), each N 0 and p N: (I Π p, )u L 2 (Ω) 2 2 h u H 1 (Ω) is satised, where Π p, is the H 1 -orthogona projection to S p, (Ω).
56 Robust approximation error estimate Theorem ([T., Takacs 2015]) For each u H 1 (Ω), each N 0 and p N: (I Π p, )u L 2 (Ω) 2 2 h u H 1 (Ω) is satised, where Π p, is the H 1 -orthogona projection to S p, (Ω). The approximation property foows using standard arguments (for the h 2 M part) u K = H Γ, u Γ, H 1 (Ω) u H 1 (Ω) = u K (for the K part)
57 Convergence theorem Theorem ([Hofreither, T., Zuehner 2015]) The two-grid method converges with rate C AC S ν if ν > C A C S smoothing steps are appied. The constants C A and C S do not depend on the grid size h and the spine degree p. The extension to the W-cyce mutigrid method is standard.
58 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
59 Matrices from tensor-product spines The mass matrix has a tensor-product structure: M = M M
60 Matrices from tensor-product spines The mass matrix has a tensor-product structure: M = M M The stiness matrix is the sum of two tensor-product matrices: K = K M + M K
61 Matrices from tensor-product spines The mass matrix has a tensor-product structure: M = M M The stiness matrix is the sum of two tensor-product matrices: K = K M + M K Tensor-product structure can be used for inverting M, but not for inverting K
62 Smoother for two dimensions K,orig = K M + M K = (h 2 M + K ) M + M (h 2 M + K ) = 2h 2 M M + K M + M K h 2 K }{{ K } h 2 K }{{ K } tensor-product ow-rank correction =: K, where K = h 1 M + K.
63 Shermann Morrisson Woodburry formua Let A R N N, B R n n, P R n N with fu rank. Then: (A PBP T ) 1 = A 1 + A 1 P(B 1 P T A 1 P) 1 P T A 1
64 Shermann Morrisson Woodburry formua Let A R N N, B R n n, P R n N with fu rank. Then: (A PBP T ) 1 = A 1 + A 1 P(B 1 P T A 1 P) 1 P T A 1 The inversion of A = K K can be done using tensor-product structure and sovers for one dimension
65 Shermann Morrisson Woodburry formua Let A R N N, B R n n, P R n N with fu rank. Then: (A PBP T ) 1 = A 1 + A 1 P(B 1 P T A 1 P) 1 P T A 1 The inversion of A = K K can be done using tensor-product structure and sovers for one dimension The rest ives ony in the boundary ayer
66 Shermann Morrisson Woodburry formua Let A R N N, B R n n, P R n N with fu rank. Then: (A PBP T ) 1 = A 1 + A 1 P(B 1 P T A 1 P) 1 P T A 1 The inversion of A = K K can be done using tensor-product structure and sovers for one dimension The rest ives ony in the boundary ayer Optima order: mutigrid sover has the same order of compexity as the mutipication with K
67 Convergence theorem Theorem ([Hofreither, T., Zuehner 2015]) The two-grid method converges with rate C AC S ν if ν > C A C S smoothing steps are appied. The constants C A and C S do not depend on the grid size h and the spine degree p. The extension to the W-cyce mutigrid method is standard.
68 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
69 One dimension p = = = p = = = ν pre + ν post = 1 + 1, τ = 0.14 Stopping criterion: Eucidean norm of the initia residua is reduced by a factor of ɛ = 10 8.
70 Two dimensions p = = p = = ν pre + ν post = 1 + 1, τ = 0.10 (for p 3) and τ = 0.11 (for p > 3) Stopping criterion: Eucidean norm of the initia residua is reduced by a factor of ɛ = 10 8.
71 Outine Preiminaries Mode probem IGA discretization Mutigrid agorithm Mutigrid framework Basis-independent smoother Smoother for one dimension Smoother for two dimensions Numerica resuts Concusions and Outook
72 Concusions (One and) two dimensions
73 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps
74 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps Optima compexity in the sense: as compex as the mutipication with K
75 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps Optima compexity in the sense: as compex as the mutipication with K Rigorous anaysis
76 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps Optima compexity in the sense: as compex as the mutipication with K Rigorous anaysis C. Hofreither, S. Takacs and W. Zuehner. A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction. Submitted (preprint: G+S-Report 33/2105), 2015.
77 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps Optima compexity in the sense: as compex as the mutipication with K Rigorous anaysis C. Hofreither, S. Takacs and W. Zuehner. A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction. Submitted (preprint: G+S-Report 33/2105), S. Takacs and T. Takacs. Approximation error estimates and inverse inequaities for B-spines of maximum smoothness. Submitted (preprint: arxiv: ), 2015.
78 Concusions (One and) two dimensions Robust convergence rates, robust number of smoothing steps Optima compexity in the sense: as compex as the mutipication with K Rigorous anaysis C. Hofreither, S. Takacs and W. Zuehner. A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction. Submitted (preprint: G+S-Report 33/2105), S. Takacs and T. Takacs. Approximation error estimates and inverse inequaities for B-spines of maximum smoothness. Submitted (preprint: arxiv: ), Thanks for your attention!
79
A Robust Multigrid Method for Isogeometric Analysis using Boundary Correction. C. Hofreither, S. Takacs, W. Zulehner. G+S Report No.
A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction C. Hofreither, S. Takacs, W. Zuehner G+S Report No. 33 Juy 2015 A Robust Mutigrid Method for Isogeometric Anaysis using Boundary
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationMultigrid Lecture Notes
Mutigrid Lecture Notes Wenqiang Feng, Steven Wise Abstract Mutigrid MG methods in numerica anaysis are a group of agorithms for soving differentia equations using a hierarchy of discretizations MG method
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationNumerische Mathematik
Numer. Math. (1999) 84: 97 119 Digita Object Identifier (DOI) 10.1007/s002119900096 Numerische Mathematik c Springer-Verag 1999 Mutigrid methods for a parameter dependent probem in prima variabes Joachim
More informationA proposed nonparametric mixture density estimation using B-spline functions
A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More information18-660: Numerical Methods for Engineering Design and Optimization
8-660: Numerica Methods for Engineering esign and Optimization in i epartment of ECE Carnegie Meon University Pittsburgh, PA 523 Side Overview Conjugate Gradient Method (Part 4) Pre-conditioning Noninear
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationThe Group Structure on a Smooth Tropical Cubic
The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More informationInvestigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l
Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationThe Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation
The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric
More informationNon-Parametric Analysis of Supernova Data and. the Dark Energy Equation of State. Department of Statistics Carnegie Mellon University
Non-Parametric Anaysis of Supernova Data and the Dark Energy Equation of State P. E. Freeman, C. R. Genovese, & L. Wasserman Department of Statistics Carnegie Meon University E-mai: pfreeman@cmu.edu INTRODUCTION
More informationB-SPLINE-BASED MONOTONE MULTIGRID METHODS
SIAM J NUMER ANAL Vo 45, No 3, pp 75 99 c 2007 Society for Industria and Appied Mathematics B-SPLINE-BASED MONOTONE MULTIGRID METHODS MARKUS HOLTZ AND ANGELA KUNOTH Abstract For the efficient numerica
More informationA two-level Schwarz preconditioner for heterogeneous problems
A two-eve Schwarz preconditioner for heterogeneous probems V. Doean 1, F. Nataf 2, R. Scheich 3 and N. Spiane 2 1 Introduction Coarse space correction is essentia to achieve agorithmic scaabiity in domain
More informationData Search Algorithms based on Quantum Walk
Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but
More informationA natural differential calculus on Lie bialgebras with dual of triangular type
Centrum voor Wiskunde en Informatica REPORTRAPPORT A natura differentia cacuus on Lie biagebras with dua of trianguar type N. van den Hijigenberg and R. Martini Department of Anaysis, Agebra and Geometry
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationc 2011 Society for Industrial and Applied Mathematics BOOTSTRAP AMG
SIAM J. SCI. COMPUT. Vo. 0, No. 0, pp. 000 000 c 2011 Society for Industria and Appied Mathematics BOOTSTRAP AMG A. BRANDT, J. BRANNICK, K. KAHL, AND I. LIVSHITS Abstract. We deveop an agebraic mutigrid
More informationORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION
J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,
More informationTHE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia
More informationB-Spline-Based Monotone Multigrid Methods. Markus Holtz, Angela Kunoth. no. 252
B-Spine-Based Monotone Mutigrid Methods Markus Hotz, Angea Kunoth no 252 Diese Arbeit ist mit Unterstützung des von der Deutschen Forschungsgemeinschaft getragenen Sonderforschungsbereiches 6 an der Universität
More informationLecture 11. Fourier transform
Lecture. Fourier transform Definition and main resuts Let f L 2 (R). The Fourier transform of a function f is a function f(α) = f(x)t iαx dx () The normaized Fourier transform of f is a function R ˆf =
More informationBALANCING REGULAR MATRIX PENCILS
BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity
More informationFFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary
More informationAnalysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case
Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationNumerical methods for elliptic partial differential equations Arnold Reusken
Numerica methods for eiptic partia differentia equations Arnod Reusken Preface This is a book on the numerica approximation of partia differentia equations. On the next page we give an overview of the
More informationPreprint BUW-SC 11/1. Bergische Universität Wuppertal. Fachbereich C Mathematik und Naturwissenschaften. Mathematik
Preprint BUW-SC 11/1 Bergische Universität Wupperta Fachbereich C Mathematik und Naturwissenschaften Mathematik A. Brandt, J. Brannick, K. Kah, and I. Livshits Bootstrap AMG Januar 2011 http://www-ai.math.uni-wupperta.de/scicomp/
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationNumerische Mathematik
First pubished in: Numer. Math. 75: 501 522 (1997) Numerische Mathematik c Springer-Verag 1997 Eectronic Edition A waveet mutieve method for i-posed probems stabiized by Tikhonov reguarization Andreas
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationarxiv: v1 [math.fa] 23 Aug 2018
An Exact Upper Bound on the L p Lebesgue Constant and The -Rényi Entropy Power Inequaity for Integer Vaued Random Variabes arxiv:808.0773v [math.fa] 3 Aug 08 Peng Xu, Mokshay Madiman, James Mebourne Abstract
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationIterative Decoding Performance Bounds for LDPC Codes on Noisy Channels
Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channes arxiv:cs/060700v1 [cs.it] 6 Ju 006 Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department University
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationAlberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain
CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D Berlin - Wilmersdorf
Konrad-Zuse-Zentrum für Informationstechnik Berin Heibronner Str. 10, D-10711 Berin - Wimersdorf Raf Kornhuber Monotone Mutigrid Methods for Eiptic Variationa Inequaities II Preprint SC 93{19 (Marz 1994)
More informationA nonlinear multigrid for imaging electrical conductivity and permittivity at low frequency
INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probems 17 (001) 39 359 www.iop.org/journas/ip PII: S066-5611(01)14169-9 A noninear mutigrid for imaging eectrica conductivity and permittivity
More informationResearch Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation
Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationParallel algebraic multilevel Schwarz preconditioners for a class of elliptic PDE systems
Parae agebraic mutieve Schwarz preconditioners for a cass of eiptic PDE systems A. Borzì Vaentina De Simone Daniea di Serafino October 3, 213 Abstract Agebraic mutieve preconditioners for agebraic probems
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS
SIAM J. NUMER. ANAL. c 1992 Society for Industria Appied Mathematics Vo. 6, No. 6, pp. 1716-1740, December 1992 011 ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS G. BEYLKIN
More informationPrimal and dual active-set methods for convex quadratic programming
Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:
More informationOn the energy distribution in Fermi Pasta Ulam lattices
Version of 3 January 01 On the energy distribution in Fermi Pasta Uam attices Ernst Hairer 1, Christian Lubich 1 Section de mathématiques, -4 rue du Lièvre, Université de Genève, CH-111 Genève 4, Switzerand.
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 7
Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine
More informationPositive real and bounded real balancing for model reduction of descriptor systems
TECHNISCHE UNIVERSITÄT BERLIN Positive rea and bounded rea baancing for mode reduction of descriptor systems Timo Reis Tatjana Styke Preprint 28/25 Preprint-Reihe des Instituts für Mathematik Technische
More informationOn the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems
Internationa Journa of Bifurcation and Chaos Vo. 25 No. 10 2015 1550131 10 pages c Word Scientific Pubishing Company DOI: 10.112/S02181271550131X On the Number of Limit Cyces for Discontinuous Generaized
More informationKatholieke Universiteit Leuven Department of Computer Science
A mutieve agorithm to compute steady states of attice Botzmann modes Giovanni Samaey Christophe Vandekerckhove Wim Vanroose Report TW 56, January 21 Kathoieke Universiteit Leuven Department of Computer
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vo. 2, No. 5, pp. 909 926 c 2000 Society for Industria and Appied Mathematics A DEFLATED VERSION OF THE CONJUGATE GRADIENT ALGORITHM Y. SAAD, M. YEUNG, J. ERHEL, AND F. GUYOMARC H
More informationApproximation algorithms for trilinear optimization with nonconvex constraints and its extensions
Approximation agorithms for triinear optimization with nonconvex constraints and its extensions Yuning Yang Schoo of Mathematics Science and LPMC Nankai University Tianjin 300071, P.R. China Emai:nk0310145@gmai.com
More informationALGORITHMIC SUMMATION OF RECIPROCALS OF PRODUCTS OF FIBONACCI NUMBERS. F. = I j. ^ = 1 ^ -, and K w = ^. 0) n=l r n «=1 -*/!
ALGORITHMIC SUMMATIO OF RECIPROCALS OF PRODUCTS OF FIBOACCI UMBERS Staney Rabinowitz MathPro Press, 2 Vine Brook Road, Westford, MA 0886 staney@tiac.net (Submitted May 997). ITRODUCTIO There is no known
More informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationAN OVERVIEW OF MULTILEVEL METHODS WITH AGGRESSIVE COARSENING AND MASSIVE POLYNOMIAL SMOOTHING
Eectronic Transactions on Numerica Anaysis Voume 44, pp 401 442, 2015 Copyright c 2015, ISSN 1068 9613 ETNA AN OVERVIEW OF MULTILEVEL METHODS WITH AGGRESSIVE COARSENING AND MASSIVE POLYNOMIAL SMOOTHING
More informationRestricted weak type on maximal linear and multilinear integral maps.
Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationEstablishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
Georgia Southern University Digita Commons@Georgia Southern Mathematica Sciences Facuty Pubications Department of Mathematica Sciences 2009 Estabishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationScott Cohen. November 10, Abstract. The method of Block Cyclic Reduction (BCR) is described in the context of
ycic Reduction Scott ohen November, 99 bstract The method of ock ycic Reduction (R) is described in the context of soving Poisson's equation with Dirichet boundary conditions The numerica instabiityof
More informationTwo-Stage Least Squares as Minimum Distance
Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationA UNIVERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIVE ALGEBRAIC MANIFOLDS
A UNIERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIE ALGEBRAIC MANIFOLDS DROR AROLIN Dedicated to M Saah Baouendi on the occasion of his 60th birthday 1 Introduction In his ceebrated
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationHigh Accuracy Split-Step Finite Difference Method for Schrödinger-KdV Equations
Commun. Theor. Phys. 70 208 43 422 Vo. 70, No. 4, October, 208 High Accuracy Spit-Step Finite Difference Method for Schrödinger-KdV Equations Feng Liao 廖锋, and Lu-Ming Zhang 张鲁明 2 Schoo of Mathematics
More informationEFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC ANALYSIS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK EFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
More informationMat 1501 lecture notes, penultimate installment
Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationOn norm equivalence between the displacement and velocity vectors for free linear dynamical systems
Kohaupt Cogent Mathematics (5 : 95699 http://dxdoiorg/8/33835595699 COMPUTATIONAL SCIENCE RESEARCH ARTICLE On norm equivaence between the dispacement and veocity vectors for free inear dynamica systems
More informationSubstructuring Preconditioners for the Bidomain Extracellular Potential Problem
Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di
More informationPath planning with PH G2 splines in R2
Path panning with PH G2 spines in R2 Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru To cite this version: Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru. Path panning with PH G2 spines
More informationMath 220B - Summer 2003 Homework 1 Solutions
Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary
More informationGlobal Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations
Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,
More informationarxiv: v1 [math.pr] 30 Nov 2017
REGULARIZATIO OF O-ORMAL MATRICES BY GAUSSIA OISE - THE BADED TOEPLITZ AD TWISTED TOEPLITZ CASES AIRBA BASAK, ELLIOT PAQUETTE, AD OFER ZEITOUI arxiv:7242v [mathpr] 3 ov 27 Abstract We consider the spectrum
More informationUniversal Consistency of Multi-Class Support Vector Classification
Universa Consistency of Muti-Cass Support Vector Cassification Tobias Gasmachers Dae Moe Institute for rtificia Inteigence IDSI, 6928 Manno-Lugano, Switzerand tobias@idsia.ch bstract Steinwart was the
More informationConvergence and quasi-optimality of adaptive finite element methods for harmonic forms
Noname manuscript No. (wi be inserted by the editor) Convergence and quasi-optimaity of adaptive finite eement methods for harmonic forms Aan Demow 1 the date of receipt and acceptance shoud be inserted
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationHomework 5 Solutions
Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More information