A New Projected Gradient Algorithm for Total Variation Based Image Denoising
|
|
- Chester Young
- 5 years ago
- Views:
Transcription
1 Noname manuscript No. (will be inserted by te editor) A New Projected Gradient Algoritm or Total Variation Based Image Denoising Ming-Jun Lai Leopold Matamba Messi te date o receipt and acceptance sould be inserted later Abstract Tis paper is concerned wit te numerical approximation o te minimizer o te continuous Rudin-Oser- Fatemi (ROF) model. A new projected gradient algoritm or computing a numerical solution o te ROF model in te discrete setting is proposed. We sow tat te new algoritm converges. Te solution o our discrete algoritm is ten used to construct a continuous piecewise linear spline wic approximates te solution o te continuous ROF model in te L norm. We sow tat i te noised image is in te space Lip(, L ()), 0 < 1, our piecewise linear approximations converge to te solution o te ROF model. Finally, we demonstrate wit numerical experiments tat te new algoritm is as good as te traditional projected gradient algoritm and Cambolle s ixed point algoritm. 1 Introduction Since te seminal work o Rudin, Oser, and Fatemi[13] total variation based model or image restoration ave received a great deal o attention. Tey are now used in image denoising, image deblurring, and image inpainting. In eac case, te problem is ormulated as a minimization o a unctional o te orm argmin Du (), subject to te constraints F (u), (1.1) were Du () is te total-variation o te unction u on, and F (u) is a suitable set o constraints satisied by u, and R is te domain o te image u. For image de-noising, te problem was ormulated as ollows: recover te true image u rom a contaminated version = u + n, were n is a wite noise wit mean 0 and standard deviation σ. Te corresponding minimization problem is argmin Du () subject to te constraints (1.) u(x)dx = (x)dx and u(x) (x) dx σ. (1.3) It was sown (c. [6,13]) tat problem (1.) is equivalent to te ollowing unconstrained minimization argmin Du () + 1 (u ) dx (1.4) u BV () λ were λ > 0 is a Lagrange multiplier, and BV () is te Banac space o unctions o bounded variation. Te existence and uniqueness o te minimizer o te above problem was establised in [1] and [6]. Ming-Jun Lai Department o Matematics University o Georgia Atens, GA 3060, USA Tel.: Fax: mjlai@mat.uga.edu Leopold Matamba Messi Department o Matematics University o Georgia Atens, GA 3060, USA Tel.: Fax: lmatamba@mat.uga.edu
2 Ming-Jun Lai, Leopold Matamba Messi To ind numerical approximations o te solution o (1.4), one as to understand ow to discretize te total variation term Du (). All te inite dierence metods proposed in te literature are based on te ollowing act (see or example [,10] or details) ( ) ( ) u u Du () = u dx = + dx, u W 1,1 (). (1.5) x 1 x As suc, te total variation is ten approximated using a combination o a quadrature ormula or te integral, and a irst order inite dierence approximation o te derivative. Several algoritms to ind approximations o te solution o problem (1.4) ave been proposed and studied in te literature. See or example [3,6,13]. Indeed, a time dependent nonlinear PDE associated wit te Euler-Lagrange equation or (1.4) was used to approximate te minimizer o (1.4) and ten solved by inite dierence metods (c. [13]) or inite elements metods (c. [7]). In te discrete setting, a primal-dual algoritm (a popular approac in optimization) based on a straigtorward discretization o te TV term and L term was introduced by Cambolle in [?] or image denoising. Recently, Duval et al [8] proposed a projected-gradient algoritm or te discrete L TV regularization to compute an approximation o te minimizer. A numerical algoritm based on completely symmetric discretization o te TV term was studied in [15]. In particular, te researcers in [15] were able to sow tat teir numerical solution approximates te minimizer o problem (1.4) in L norm. Furtermore, te convergence o te discrete solution o te inite dierence metod or time dependent nonlinear PDE was recently establised in [11]. Tese motivate us to study te convergence o te discrete solutions o Cambolle algoritm and projected gradient algoritm to te solution o problem (1.4). Ater studying te convergence issue o tese two algoritms, we realize tat one as to use a new discretization o te total variation term wic leads to a new version o te projected-gradient algoritm. Wit tis new discretization, we are able to use te solution in te discrete setting to construct a piecewise linear interpolatory spline wic can be sown to approximate te solution o te original problem (1.4). Our contributions in tis paper are te ollowing: (1) we irst prove te existence and uniqueness o te solution based on te new discrete ROF unctional (see Teorem 1); () We ormulate a projected gradient algoritm or te new minimization problem (c. Algoritm 31) and prove its convergence (c. Teorem 3); (3) We demonstrate numerically tat our new version o te projected-gradient algoritm perorms sligtly better tan te original projected-gradient algoritm in [5,8], and is at least as good as Cambolle s ixed point iterative algoritm[4], as measured by te Peak Signal to Noise Ratio (PSNR) (see Section 5); (4) We obtain an error bound in L () between te piecewise linear interpolatory spline o te discrete solution and te solution o te continuous ROF model (see Teorem ). Preliminaries and Notations In tis section, we give preliminary results, and introduce te notations tat we sall use in te paper..1 Basic notations Let be an open subset o R, we denote te indicator unction o te set by 1, x, 1 (x) = 0, x /. For a given η R, we sall denote by τ η te image o te set under te translation wit te vector η, i.e τ η := x + η : x }. For a unction u : R, we denote by τ ηu te unction wose domain is τ η and is deined by τ ηu(x) = u(x + η), x τ η. Next, we recall te deinitions o some unctional spaces tat are used in tis work. Let 1 p < be ixed. L p () is te standard Banac space o p integrable unctions } L p () := u : R : u(x) p dx <. L p loc () is te set o unctions u tat are locally p integrable, i.e } L p loc u () = : R : u(x) p dx <, K compact. K
3 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 3 Te norm o an element u L p () is given by ( 1/p u p = u(x) dx) p. It is well known tat te translation operator τ η is a bounded linear operator rom L p () into L p ( τ η ). Let > 0 be given. Te p modulus o continuity o order, o a unction u L p (), is deined by ω p(u, ) = τ ηu u p, (.1) η were η stands or te Euclidean norm o η. For u L p loc () and A a subset o wose closure Ā is a compact, te p modulus o continuity o u wit respect to A, denoted ω p(u, ) A, is given by ω p(u, ) A = ω p(u1 A, ). (.) Let 0 < 1. We denote Lip(, L ()) te subspace o L () deined by } Lip(, L ()) = u L () : ω (u, ) < >0.. A discretization o te ROF unctional. In te sequel sall denote te unit square [0, 1] o R. We subdivide into (N 1) square sub-domains o side lengt to get a uniorm quadrangulation. A triangulation is ten obtained rom by splitting eac square into two triangles using te Nortwest-Souteast diagonal as sown in Figure 1. ω i,j+1 ω i+1,j+1 T u i,j T d i,j ω i,j ω i+1,j Fig. 1: Type I triangulation o te domain wit vertexes ω i,j. T u i,j is te triangle wit vertexes ω i+1,j, ω i+1,j+1, ω i,j+1 and T d i,j is te triangle wit vertexes ω i,j, ω i+1,j, ω i,j+1. Let ω 1,1 be te lower let corner o. We denote te set o vertices o te triangulation by and deine a partition } o subordinate to V by V = ω 1,1 + Z } := ω i,j : 1 i, j N}, ( := ω i,j + [ /, /] ). Let E 1 λ (u) := Du () + (u ) dx. (.3) λ
4 4 Ming-Jun Lai, Leopold Matamba Messi Trougout te paper E λ (u) will be reerred to as te continuous ROF unctional. We are interested in devising a numerical sceme or computing an approximation o te minimizer o E λ (u). However, since Du is a measure, discretizing E λ (u) solely based on u is a delicate matter. Many researcers [5,1] ave used te closed orm o Du () given in (1.5) to designed discretization sceme o te ROF unctional E λ (u) tat combine a ormal inite dierence approximation o u wit a quadrature approximation o te integral. Our discretization o te ROF unctional E λ (u) as a similar lavor. We irst discretize te total variation as ollows J (u) := 1 i,j<n u i+1,j u i,j + u i,j+1 u i,j + 1 i,j<n u i+1,j+1 u i,j+1 + u i+1,j+1 u i+1,j, (.4) were u i,j is a suitable discretization o te unction u at te vertexes o te triangulation. For example i u is a continuous piecewise linear spline wit respect to, on can easily sow tat J (u) = Du (); tis is te reason wy we use (.4) to discretize Du (). Te resulting discrete ROF-unctional E λ, (u) is given by E λ, (u) := J (u) + 1 λ 1 i,j N ui,j i,j, u R N N (.5) were i,j := (ω i,j ) is a suitable discretization on o te datum. From tis point onwards, E λ, (u) will be reerred to as te discrete ROF unctional and te associated miniization problem sall be called te discrete ROF model. argmin E λ, (u) (.6) u R N N.3 Embedding and Projection Operators We introduce various operators tat will elp us construct a piecewise linear approximation o te solution to te continuous L ROF model. In te previous section, we gave a ormal discrete approximation o te ROF unctional. We now clariy ow tis discretization is obtained. Te sampling operator, Q : L () R N N, is deined by (Q ) i,j := 1 (x)dx. (.7) Q will also denote te projection o L () onto te linear space o piecewise constant unction wit respect to te partition i,j : 1 i, j N} o, in wic case Q is deined by Q (y) := (Q ) i,j = 1 (x)dx, or all y. (.8) Next, we denote by l (V ) te vector space R N N endowed wit te inner product u, v := N u i,j v i,j, (.9) and te corresponding norm u := u, u 1/. Te ollowing properties o Q ollows rom Jensen s Inequality. i,j=1 Lemma 1 Te sampling operator Q deined above as te ollowing properties Q, and Q, or all L (). (.10)
5 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 5 Next, we explain two metods o constructing a unction o bounded variation rom an element u l (V ). First, we create a piecewise constant unction deined on via te partition : 1 i, j N} as ollows C u(y) = u i,j, i y. (.11) We denote by C : l (V ) BV () te piecewise constant embedding o l (V ) into BV () deined by (.11). Second, we deine a piecewise linear interpolatory operator. P : l (V ) BV (). For u l (V ), P u is te continuous piecewise linear polynomial deined by P u(y) = 1 i,j N wit φ i,j te continuous piecewise linear unction suc tat u i,j φ i,j (y), (.1) φ i,j (ω i,j ) = 1, and φ i,j (ω) = 0, ω V \ ω i,j. (.13) Wit a sligt abuse o notation, P will also denote te piecewise linear interpolation o unctions o class C k ( ), k 0 over te set V ; in wic case it will be deined by P (y) = 1 i,j N (ω i,j )φ i,j (y), y, C k ( ). Lemma Suppose tat is endowed wit te triangulation. Ten or all u l (V ), tere olds C u = u l (V ) and P u 3 u l (V ), (.14) P u C u ω(u, 1) l (V ) (.15) were ω(u, m) l (V ) denotes te discrete modulus o smootness o order m N and is deined by ω(u, m) l (V ) = =( 1, ) Z, 1 i,i+ 1 N max( 1, ) m 1 j,j+ N ui+1,j+ u i,j. (.16) Proo Let u l (V ) be given. Te inequalities (.14) ollow rom te deinition o te Euclidean norm on l (V ) and te act tat te area o te port o φ i,j is at most 3. Next we prove te inequality (.15). By deinition, we ave P u C u = P u(y) C u(y) dy 1 i,j N ( ) ul,k u i,j φl,k (y) dy 1 i,j N 1 l,k N ul,k u i,j φ l,k(y)dy 1 i,j N 1 l,k N i,j ul,k u i,j ω(u, 1) l (V ), 1 i,j N 1 l N, l i =1 1 k N, k j =1 wic completes te proo. Remark 1 Let u R N N be ixed, and let ω(u, m) be te modulus o continuity o u R N N wit respect to te Euclidian norm u = u, u. It is easy to see tat ω(u, m) l (V ) ω(u, m). (.17).4 Periodic Extension operators In tis section, we construct a periodic extension o u R N N to Z. We also present te construction o a periodic extension to R o a unction deined on.
6 6 Ming-Jun Lai, Leopold Matamba Messi Te discrete periodic extension Let u R N N be ixed. We construct a discrete unction Ext(u): Z R deined in two steps as ollows: 1) First, we extend u into an element u N R N N as ollows u N i,j = u i,j i 1 i, j N, u N i+1,j i i > N, 1 j N, u i,n j+1 i 1 i N, j > N, else. u N i+1,n j+1 (.18) ) Ext(u) is te periodic extension o u N to all o Z. In te sequel, we will denote by u N te restriction o te unction Ext(u) to te set Ψ N := (i, j) Z Z: 1 i, j N}. For = ( 1, ) Z and u := (u i,j ) i,j Z a discrete unction, we deine te translation τ u o u as ollows: (τ u) i,j = u i+1,j+, i, j Z. Te restriction o τ u to a set o te orm Ψ M := (i, j): 1 i, j M} sall be denoted u M. Te continuous periodic extension Te construction o te periodic extension tat ollows was presented in [16, 15]. Let u : R be given. First, we extend u to u : [0, ] R by te ormula u (x 1, x ) = u(x 1, x )1 (x 1, x ) + u( x 1, x )1 τ1,0(x 1, x ) + u(x 1, x )1 τ0,1(x 1, x ) + u( x 1, x )1 τ1,1(x 1, x ), (.19) were 1 A denotes te indicator unction o te set A, and τ m,n = (x 1 + m, x + n): (x 1, x ) }. Te extension o u is te periodic unction Ext(u) : R R wit period [0, ] tat coincides wit u on [0, ]. Since as compact closure, it is easy to sow tat te continous periodic extension operator, Ext, maps L p () into L p loc (R ). Proposition 1 Suppose tat L (). Ten or any 0 < 1, we ave were C > 0 is a constant independent o and 1 = ω (Ext(), ) 1 Cω (, ), (.0) (m,n) Z (m,n) 1 τ (m,n). Proo Let L () be given, and η R be ixed wit η. By deinition τ η(ext()1 1) Ext()1 1 = Ext()(x + η) Ext()(x) dx 1 τ η 1 1 i,j i,j 1 1 m 1 m i =1 1 n 1 n j =1 Ext()(x + η) Ext()(x) dx τ (m,n) τ (i,j) η Ext()(x + η) Ext()(x) dx τ (i,j) ( τ η) (x + η) (x) dx = 18 τ η. τ η As a consequence, we ave ω (Ext(), ) 1 3 ω (, ). Lemma 3 For any L () and 0 < 1, tere olds were K > 0 is a constant independent o. C Q K ω (, ), (.1)
7 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 7 Proo By deinition, we ave C Q = = 1 i,j N 1 i,j N 1 i,j N 4 (x) 1 (y)dy dx 1/ ( ) 1/ 1 (x) (y) dy dx ( 4 z : z Ext()(x) Ext()(x + z) dz) dx } ( 4 z : z Ext()(x) Ext()(x + z) dz) dx } z : z } ( C Q 8πω (Ext(), ) 1. 1/ Ext()(x) Ext()(x + z) dx) dz Inequality (.1) ollows rom te latter inequality tanks to Proposition 1, and te act tat ω (Ext(), ) 1 ω (Ext(), ) 1. 1/ 1/ 3 A total variation based model or digital image denoising In te sequel we treat gray-scale images o size N N as rectangular matrices o dimension N N and denote by X := R N N te maniold o gray-scale images o size N N. 3.1 Te model and its properties To compute te total variation o elements o X, we introduce two discrete gradient operators + = ( x +, y + ) and = ( x, y ), wic are linear operators rom X into Y := X X, deined by ( x 0, i i = N or j = N +u) i,j = (3.1) u i+1,j u i,j oterwise; ( y + u) 0, i i = N or j = N i,j = (3.) u i,j+1 u i,j oterwise; and ( x 0, i i = 1 or j = 1 u) i,j = u i,j u i 1,j oterwise; ( y u) 0, i i = 1 or j = 1 i,j = u i,j u i,j 1 oterwise. (3.3) (3.4) We associate to te discrete gradient operators + and te discrete divergence operators, div + := + : Y X and div := : Y X, deined respectively by div + (p) i,j = + 0 i i = N or j = N p 1 i,j oterwise 0 i i = N or j = N p i,j oterwise 0 i i = 1 or j = N p 1 i 1,j oterwise 0 i i = N or j = 1 p i,j 1 oterwise
8 8 Ming-Jun Lai, Leopold Matamba Messi and div (p) i,j = + 0 i i = N or j = 1 p 1 i+1,j oterwise 0 i i = 1 or j = N p i,j+1 oterwise 0 i i = 1 or j = 1 p 1 i,j oterwise 0 i i = 1 or j = 1 p i,j oterwise Te total variation o a gray-scale image u X is ten deined by J(u) := 1 1 i N 1 j N ( ( + u) i,j + ( u) i,j ), (3.5) were or an element p := (p i,j ) := (p 1 i,j, p i,j) Y, pi,j is te Euclidean norm o p ij := (p 1 i,j, p i,j) in R. Consequently, te digital ROF-unctional is given by E 1 d,λ (u) := J(u) + λ u, (3.6) were X is te Euclidean norm on X = R N N, and X. Since te unctional E d,λ (u) is strictly convex and coercive, it ollows (c. [9, Proposition 1., page 35]) tat te minimization problem argmin E d,λ (u) as a unique solution. Henceort, we will reer to te problem u X argmin E d,λ (u) (3.7) u X as te digital ROF-model. Te next result, similar to te one proved in [15,16], sows tat problem (3.7) is stable under small perturbations in te datum. Lemma 4 (Stability) Let u be te minimizer o E d,λ (u) and ug te minimizer o Eg d,λ (u). Ten u u g g. (3.8). We now prove a tecnical lemma tat asserts tat te discrete ROF model is compatible wit translation o te datum Lemma 5 Let R N N be ixed and u be te minimizer o te unctional E d,λ (u). Ten, un is a minimizer o E N d,λ (u). Moreover, or any Z suc tat = 1, te minimizer o E N d,λ (u) is un,, te restriction o τu to Ψ N. E N u R N N Proo We irst sow tat u N = argmin d,λ (u). It is easy to see tat or all u RN N, we ave E N d,λ (un ) = 4E d,λ (u) Furtermore, or every u R N N, tere exists x u R N N suc tat E N N (u) E x u be te constant vector tat equals u i0,j 0 E N d,λ (xn u ). Consequently, N min E u R N N d,λ (u) = d,λ d,λ (xn u ). In eect, letting wit (i 0, j 0 ) argmin u i,j i,j, it is easy to veriy tat E N d,λ (u) 1 i,j N N min E u R N N d,λ (un ) = 4 min E u R N N d,λ (u). Hence, u N is te minimizer o E N d,λ (u) wit respect to RN N. Next let Z wit = 1. For clarity o te argument, we ix = (1, 0). An argument identical to te one above sows tat N min E u RN N d,λ (u) = min N E u RN N d,λ (un ). Let k N be ixed. Ten, or any periodic discrete unction g wit period [1, N] [1, N], we ave k E gn d,λ (un g ) = E gnk d,λ (unk g ) = min E g Nk u R Nk Nk d,λ (u), (3.9)
9 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 9 were u g is te minimizer o E gn d,λ (u), and unk g is te restriction o te periodic (wit period [1, N] [1, N]) extension o u g to [1, Nk] [1, Nk]. Moreover, or any u R N N, tere olds E gnk d,λ (unk ) = E gnk d,λ (unk ) kd(u), (3.10) were D(u) = ( + (u) 1,j + (u) 1,j ). (3.11) 1 j<n E N E N d,λ u R N N Let y R N N be suc tat y N = argmin d,λ (yn we ave (u). We want to prove tat yn ). Let us consider te minimization problem associated to Nk. Since u N k E N d,λ (yn k E N d,λ (yn ) = E Nk d,λ (ynk ) Dividing te last inequality above by k, we obtain E N d,λ (yn = E Nk d,λ (ynk ) kd(y) by equation (3.10) E Nk d,λ (unk ) kd(y) by equation (3.9) E Nk d,λ (unk, ) + k(d(u ) D(y)) by equation (3.10) ) k E N d,λ (un,) + k(d(u ) D(y)) by equation (3.9). ) E N d,λ (un Passing to te limit as k in (3.1) yields E N = u N N,, i.e E d,λ (un, ) is te minimizer o E N λ (u),,) + (D(u ) D(y)), k N, k 1. (3.1) k d,λ (yn ) E N d,λ (un, E N u R N N ). Tus, un, = argmin ). d,λ (un Our argument above works mutatis mutandis or any := ( 1, ) wit one o te components being zero. For (1, 1), (1, 1), ( 1, 1), ( 1, 1)}, te proo ollows rom te previous case by observing tat τ = τ (1,0) τ (0,). Remark We observe tat or any u R N N E λ, (u) = E / d,λ/ (u/). (3.13) Tereore, Lemma 4 and Lemma 5 remain valid or te unctional E λ, (u). Te equation (3.13) gives te relation between te discrete ROF unctional (.5) and te digital ROF unctional (3.6). 3. Primal-Dual ormulation We establis te primal-dual ormulation o te minimization problem associated to te unctional E d,λ (u). We also prove te existence and uniqueness o te minimizer. First, we establis an alternate ormula or te discrete total variation J(u) deined in (3.5). By Riesz representation Teorem, we can rewrite te discrete total variation J(u) above in te ollowing orm: J(u) = 1 = 1 = 1 ( + u) i,j, p i,j + ( u) i,j, q i,j 1 i N p i,j R, p i,j 1 q i,j R, q i,j 1 1 j N + u, p Y + 1 u, q p Y, p 1 Y q Y, q 1 div + (p) X + 1 u, div (q) p Y, p 1 u, X q Y, q 1 J(u) = p Y, p 1 q Y, q 1 u, 1 div +(p) + 1 div (q) X (3.14)
10 10 Ming-Jun Lai, Leopold Matamba Messi were or p = (p 1, p ) and q = (q 1, q ) in Y, p, q Y = 1 i,j N p 1 i,jqi,j 1 + p i,jqi,j, and p = max (p 1 i,j ij ) + (p ij ). Tereore, te minimization problem (3.7) is equivalent to te ollowing saddle-point problem argmin u X p Y, p 1 q Y, q 1 u, 1 div +(p) + 1 div (q) X + 1 λ u. (3.15) Te saddle-point problem (3.15) above is reerred to as te primal-dual ormulation o problem (3.7), wit primal variable u and dual variable p. Let and L(u : p, q) be te unctional deined on X Y Y by B Y = p Y : p 1}, L (u; p, q) := u, 1 div +(p) + 1 div (q) + 1 λ u. (3.16) We ave te ollowing result or te existence o a solution to problem (3.15). Lemma 6 Te unctional L : X Y Y R deined by (3.16) as a saddle point (ū; p, q) in te set X B Y B Y. Furtermore te component ū o any saddle point o L is te solution o te minimization problem (3.7). Proo We observe tat L(u; p, q) is a strictly convex quadratic unction in u and linear in (p, q). Moreover, bot te mappings u L(u; p, q) and (p, q) L(u; p, q) are Gâteaux-dierentiable, wit derivatives ul(u; p, q) = 1 (div +(p) + div (q)) 1 (u ), (3.17) [ λ 1 p,ql(u; p, q) = +u, 1 ] u. (3.18) Let ū = + λ (div +( p) + div ( q)), (3.19) wit ( p, q) argmin p,q B Y λ (div + (p) + div (q)) +. (3.0) We sow tat (ū; p, q) is a saddle-point or L(u; p, q) wit respect to te set X B Y B Y. To tis end, it suices to ceck tat ul(ū; p, q), u ū 0, u X p,ql(ū; p, q), (p p, q q) 0, p, q B Y. wic ollows by deinition o ū and ( p, q). In eect, by deinition o ū, we ave ul(ū; p, q) = 0. Tus, Next, it is easy to see tat or all p, q Y we ave ul(ū; p, q), u ū = 0, u X. p,ql(ū; p, q), (p p, q q) = 1 λ λ(div + p + div q) +, λ(div + (p p) + div (q q)). Now, we recall tat λ(div + ( p) + div ( q)) is te ortogonal projection o onto te closed convex set λk. Consequently, by te caracterization o te ortogonal projection, we ave λ(div + p + div q) +, λ(div + (p p) + div (q q)) 0, p, q B Y. Hence, by (3.18) p,ql(ū; p, q), (p p, q q) 0 or all p, q B Y.
11 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 11 Remark 3 I (ū; p, q) is a saddle-point o L(u; p, q) wit respect to te set X B Y B Y, ten L(ū; p, q) = in J(u) + 1 u X λ u = J(ū) + 1 ū λ = max p,q B Y λ λ 8 div +(p) + div (q) + λ. Hence, or any saddle-point (ū; p, q) o L(u; p, q), we ave ( p, q) argmin p,q B Y λ (div + (p) + div (q)) +. We ave proved te ollowing caracterization o a minimizer o te discrete ROF unctional E d,λ. Teorem 1 A vector u X is te solution o problem (3.7) i and only i (u ) is te ortogonal projection o onto te closed convex set λk := λ [div + (p) + div (q)] : p, q B Y }. (3.1) Furtermore, te solution u is given by wit ( p, q) argmin p,q B Y λ (div + (p) + div (q)) +. u = + λ (div +( p) + div ( q)), For a closed convex subset A o a Euclidean space, we denote by π A te ortogonal projection onto A. In particular or A = B Y, it can be sown tat ( ) p 1 i,j π BY (p) i,j = max(1, p i,j ), p i,j, 1 i, j N. (3.) max(1, p i,j ) Te ormula above ollows rom te observation tat B Y is te cartesian product o unit disks o R, and R can be isometrically embedded in Y, wit bot spaces endowed wit teir natural Euclidian norms. Tereore, or 1 i, j N, π BY (p) i,j is te ortogonal projection o p i,j = (p 1 i,j, p i,j) onto te unit disk o R, given by (3.). Te ollowing result is a straigtorward consequence o te caracterization o minimizers o a Gâteaux-dierentiable unction and orms te basis o our projected-gradient algoritm. Proposition I ( p, q) argmin λ (div + (p) + div (q)) +, ten or any τ > 0, tere olds p,q B Y p = π BY ( p + τ + [div + ( p) + div ( q) + /λ]), q = π BY ( q + τ [div + ( p) + div ( q) + /λ]). (3.3) Proo Let ( p, q) argmin p,q B Y λ(div + (p) + div (q) + ) be ixed. Ten, p argmin p B Y div + (p) + div ( q) + /λ, q argmin q B Y div + ( p) + div (q) + /λ). We notice tat te mappings p div + (p) + div ( q) + /λ and q div + ( p) + div (q) + /λ are Gâteauxdierentiable wit dierentials + (div + (p) + div ( q) + /λ) and (div + ( p) + div (q) + /λ), respectively. Tereore (see [9, Proposition.1, page 37]), we ave It ten ollows tat or any τ > 0 + (div + ( p) + div ( q) + /λ), p p 0, p B Y, (div + ( p) + div ( q) + /λ), q q 0, q B Y. p [ p + τ + (div + ( p) + div ( q) + /λ)], p p 0, p B Y, q [ q + τ (div + ( p) + div ( q) + /λ)], q q 0, q B Y, wic is equivalent to (3.3).
12 1 Ming-Jun Lai, Leopold Matamba Messi 3.3 A projected gradient algoritm By Teorem 1, computing te minimizer o te discrete ROF unctional E d,λ (u) is equivalent to projecting onto te set λk deined by (3.1). In view o Proposition, we propose a projected gradient algoritm or computing te projection o onto λk. We prove te convergence o te proposed algoritm by adapting te argument used in [8]. Algoritm 31 (Projected-Gradient) Let τ > 0 be ixed. For n = 0, let p 0 = q 0 = 0. Step 1: Compute u n u n = + λ [div +(p n) + div (q n)]. (3.4) Step : Compute p n+1 and q n+1 p n+1 = π BY (p n + τ + [div + (p n) + div (q n) + /λ]), q n+1 = π BY (q n + τ [div + (p n) + div (q n) + /λ]), (3.5) wit π BY given by te equation (3.). Step 3: Increment n by 1, ten go back to step 1. We ave te ollowing convergence result or Algoritm 31. Proposition 3 I 0 < τ < 1/8, ten Algoritm 31 converges. More precisely, or any ( p, q) satisying (3.0), we ave lim n (div +(p n) + div (q n)) = div + ( p) + div ( q). Proo Let ( p, q) satisying (3.3) be ixed and τ > 0 be given. Since te projection mapping π BY iner rom (3.3) tat p n+1 p + q n+1 q Id τa ( p n p + q n q ), is non-expansive, we were A : Y Y Y Y is te linear operator ( ) + div A := + + div. div + div Consequently, as long as τ > 0 is suc tat κ := Id τa = 1, we get ( ) ( ) pn+1 p pn p q n+1 q, n 0. (3.6) q n q Next, we sow tat κ = 1 or 0 < τ < 1/8. We note to begin tat A is sel-adjoint and nonnegative deinite. Indeed or p, q Y, we ave A ( ) ( ) p p, q q Y Y = + (div + (p) + div (q)), p Y + (div + (p) + div (q)), q Y = div + (p) + div (q), div + (p) X + div + (p) + div (q), div (q) X = div + (p) + div (q), div + (p) + div (q) X 0. Tus, Y Y = ker(a) F,were F is te closure i Im(A). Moreover, all o te eigenvalues o A are nonnegative and κ = max(1, 1 τ A ), were A denotes te spectral norm o A (te largest eigenvalue o A). So, or any τ suc tat 0 < τ / A, we ave Id τa = 1. Now, we compute an upper-bound or A. We recall tat by deinition A = (p,q) (0,0) A ( ) p, A q ( ) p q Y Y p + q.
13 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 13 By Caucy-Scwarz inequality and te deinition o te norm o a linear operator, we ave ( ) ( ) p p A, A q q A = Tus, it is te case tat (p,q) (0,0) p + q Y Y ( + + ) div + (p) + div (q) (p,q) (0,0) p + q ( + + ) (p,q) (0,0) ( + + ) max( div +, div ) By deinition o div +, or all u X and all p Y div + p + div + div p q + div q p + q A ( + + ) max( div +, div ). (3.7) + u, + u Y = u, div + ( + u) X div + (p), div + (p) X = p, + (div + (p)) Y. Consequently, by Caucy-Scwarz inequality and te deinition o te norm o a linear operator, we obtain + div + + and div + + div + wic sows tat + = div +. A similar argument sows tat = div. Now, or u X, + u = (u i+1,j u i,j ) + (u i,j+1 u i,j ) 1 i,j N 1 1 i,j N 1 u i+1,j + u i,j + u i,j+1 8 u. Hence, + 8. A similar argument sows tat div = 8. Tus, by (3.7) we ave tat A 16. (3.8) Tereore, or any 0 < τ < 1/8, tere olds te inequality κ = Id τ A = 1. We now sow tat te sequence (div + (p n) + div (q n)) n 0 converges by sowing tat all o its convergent subsequences ave te same limit. Let (div + (p nk ) + div (q nk )) k 0 be a convergent subsequence. Tere exists a urter subsequence, wic we denote again by n k, suc tat p nk p and q nk q. Since te projection π BY is non-expansive, it ollows rom (3.5) and (3.6) tat p nk +1 p, q nk +1 q, and ) ) (p p = ( p p q q, q q were is te Euclidean norm on Y Y induced by te Euclidean norm o Y. As a consequence o te above equality and te equation p = π BY ( p + τ + (div + ( p) + div ( q))), we ave But te latter inequality is nonsense. Tus, q = π BY ( q + τ (div + ( p) + div ( q))), ) ( p p ) = (p p ) q q ( p p q q (Id τa) q q ) ) ( p p ( p + 1 τ A p q q q q ker(a) F ) < ( p p ) ( p p i / ker(a) and 0 < τ < 1/8. q q q q ) ( p p ker(a), or equivalently, q q div + ( p p) + div ( q q) ker( + ) ker( )
14 14 Ming-Jun Lai, Leopold Matamba Messi Indeed, by deinition o te divergence operators div + ( p p) + div ( q q) = p p, + (div + ( p p) + div ( q q)) Y q q, (div + ( p p) + div ( q q)) Y = p p, 0 Y q q, 0 Y = 0. Hence, div + ( p p) + div ( q q) = 0, or equivalently, div + ( p) + div ( q) = div + ( p) + div ( q). We ave sown tat any convergent subsequence o (div + (p n) + div (q n)) n 0 converges to div + ( p) + div ( q), tus te sequence (div + (p n) + div (q n)) n 0 converges to div + ( p) + div ( q). Remark 4 One may obtain an alternating version o Algoritm 31 by modiying te second step as ollows p n+1 = π BY (p n + τ + [div + (p n) + div (q n) + /λ]), q n+1 = π BY (q n + τ [div + (p n+1 ) + div (q n) + /λ]). (3.9) Wile te proo o convergence o te alternating version o Algoritm 31 is still eluding us, te numerical experiments suggest tat one sould be able to prove a convergence result or 0 < τ < 1/4. 4 Numerical Approximation o te continuous ROF model In tis section, we study ow to approximate te minimizer o te continuous ROF model u = argmin E λ (u), (4.1) u BV () were E λ (u) is deined in (.3) Our approac consists in extending te minimizer z, o te discrete ROF unctional (.5) associated wit Q into a continuous piecewise linear polynomial P z, wit respect to te triangulation. We sall obtain an error bound in L () between P z, and u. To begin, we recall te ollowing standard act about te ROF unctional (c. [16,15]), and provide a proo or convenience. Lemma 7 Let u BV () be te minimizer o te unctional E λ (u) deined by (.3). Ten, or any v BV (), tere olds ( ) v u λ E λ (v) E λ (u ). (4.) Moreover, E λ (u ) = 1 λ ( u ). (4.3) Proo Since u is te minimizer o E λ (u) and E λ is convex, it ollows tat 0 E λ (u ). Since, E λ (u ) = Du () + (u )/λ, we iner tat ( u )/λ Du (). Tus, or any v BV (), we ave As a consequence, we ave Dv () Du () u, v u λ. (4.4) E λ (v) E λ (u ) = Dv () Du () + 1 λ u λ E λ (v) E λ (u ) 1 λ v u. ( ) v u, v u + u, v u λ + 1 λ v u To sow te equality (4.3), we coose v in (4.4) to be equal to 0 and u, respectively. Using te act tat te total variation unctional is positively 1-omogeneous, we iner tat Du () = u, u λ. Tereore E λ (u ) = u λ, u + u λ, u = 1 λ ( u ). Te next result is an upper bound o te dierence between te energy o te discrete minimizer and te energy o its piecewise linear interpolation.
15 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 15 Lemma 8 Let L () be given, and z, te minimizer o te unctional E λ, (u). Ten E λ (P z, ) E λ, (z,) C λ were C is a positive constant depending on and λ. [ ω(z,, 1) + ω (, ) ], (4.5) Proo Since J (z, ) = DP z, (), we ave ( ) λ E λ (P z, ) E λ, (z, = P z, C z, C Q P z, C z, ( ) P z, C z, + C z, + C z, C z, C Q Now, let us ind an upper bound or te dierence C z, C z, C Q. ( C z, C z, C Q = C z, (x) (x) C z, (x) C Q (x) ) dx = (C Q (x) (x)) ( C z, (x) C Q (x) (x) ) dx C Q C z, C Q. Te inequality (4.5) ollows rom Lemma togeter wit (.17), and Lemma 3 by observing tat C z, C Q 4 and C z,. Let 1 := x R : dist(x, ) < 1}. We denote by ū an extension (see [, Proposition 3.1 page 131]) o u to a unction in BV (R ) suc tat ū = 0 in R \ 1, and Dū ( ) = 0. Let ( ρ(x) = c exp 1 x 1 ) 1 x 1} (x), x R be te standard convolution kernel, and ρ ɛ} ɛ>0 te ensuing amily o molliiers. Te constant c > 0 in te deinition o ρ above is cosen suc tat ρ(x)dx = 1. For ɛ > 0 ixed, we denote by ū ɛ := ρɛ ū te ɛ-molliication o ū. Ten, it is known tat (see [,10]) Dū ɛ () Dū ( ) = Dū () = Du (). (4.6) Lemma 9 Suppose L (). Let u be te minimizer o te ROF-unctional E λ (u), and z, te minimizer o te discrete ROF-unctional E λ, (u). Ten, or any 0 < < 1 and or 1 E λ, (z,) E λ (u ) + K 1 (1 )/ + K ū ū (4.7) were ū is te /4 molliication o ū, K 1 and K are positive constants tat depend only on λ,, and. Proo Let ɛ > 0 be ixed. By deinition o z,, we ave E λ, (z,) E λ, (Q ū ɛ ) = J (Q ū ɛ ) + 1 λ C Q ū ɛ C Q. (4.8) Moreover by deinition o te operators C and Q, it is easy to establis tat C Q ū ɛ C Q ū ɛ u + ū ɛ ū ( u + ū ɛ ū ) Next, we obtain an upper bound on J (Q ū ɛ ). By deinition J (Q ū ɛ ) = 1 i,j<n + 1 i,j<n u + 6 ū ɛ ū using (4.3). (4.9) (Q ū ɛ ) i+1,j (Q ū ɛ ) i,j (Q ū ɛ ) i+1,j+1 (Q ū ɛ ) i,j+1 (Q + ū ɛ ) i,j+1 (Q ū ɛ ) i,j (Q + ū ɛ ) i+1,j+1 (Q ū ɛ ) i+1,j. (4.10)
16 16 Ming-Jun Lai, Leopold Matamba Messi Now, let 1 i, j < N be ixed. Ten, by te Taylor ormula, we ave (Q ū ɛ ) i+1,j (Q ū ɛ ) i,j = 1 ɛ ū (x 1 +, x ) ū ɛ (x 1, x ) dx = 1 ( ū ɛ (x) + ) 1 ū ɛ (x 1 + t, x )dt dx. x 1 0 Tanks to Jensen s inequality, we iner rom te latter equation tat (Q ū ɛ ) i+1,j (Q ū ɛ ) i,j 1 ū ɛ (x) x 1 dx + Bij(), 1 (4.11) were B 1 i,j() is given by B 1 i,j() = ū ɛ (x) ū ɛ x 1 x 1 (x 1 + t, x ) dtdx x 1 0 ū ɛ (x 1 + t, x ) dtdx. (4.1) Excanging te order o integration and using Caucy-Scwarz inequality, we obtain an upper bound on Bi,j() 1 as ollows: Bi,j() 1 1 ū ɛ (x) x 1 dx + + ū ɛ (x) dx. (4.13) 4 Similarly, we can sow tat (Q ū ɛ ) i,j+1 (Q ū ɛ ) i,j 1 wit B ij() deined by B i,j() = ū ɛ (x) ū ɛ x x and satisying te inequality Bi,j() 1 (x 1, x + t) dtdx ū ɛ x dx i+1,j x 1 x 1 ū ɛ (x) x dx + Bij(), (4.14) ū ɛ (x 1, x + t) dtdx, (4.15) x ū ɛ dx. (4.16) Combining (4.11) and (4.14) in equation (4.10), and using Jensen s inequality, we obtain or 1 J (Q ū ɛ ) Dū ɛ () + 1/ Dū ɛ () + ū ɛ ū ɛ x + 1 x dx ) Du () + ( Du 1/ ū ɛ () + + ū ɛ dx by (4.6). (4.17) Now, using te isometric embedding o L 1 () into te space o Radon measure, we ave or i = 1, ū ɛ x i dx = ū ɛ ( ) ρ ɛ x φ dx = ū i φ C( ) φ 1 x φ dx i φ C( ) φ 1 ρ ɛ u φ C( ) x i φ 1 Substituting te latter inequality into (4.17) gives φ K ɛ, x 1 J (Q ū ɛ ) (1 + ) Du () + K Finally, letting ɛ = /4 wit < 1 and combining (4.19) and (4.9), we get x x wit K = u ρ x i E λ, (z,) E λ (u ) + K 1 (1 )/ + K ū ū,. (4.18) ɛ. (4.19) were ū is te /4 molliication o ū, and te K is are positive constants depending only on, λ, and. We are now ready to establis an error bound or te L norm o P z, u.
17 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 17 Teorem Let u be te minimizer o te ROF unctional E λ (u) in BV (), and z, te minimizer o te discrete ROF unctional E λ, (u). Ten, or any 0 < < 1 and 0 < 1 P z, u K 3 ω (, ) + K K 5 ū ū, (4.0) were K 3, K 4, and K 5 are positive constants independent o. Proo Using (4.), we ave P z, u λ ( = λ ) E λ (P z, ) E λ (u ) Now, tanks to te Lemmas 8 and 9, i.e., (4.5) and (4.7), we ave tat ( E λ (P z, ) E λ, (z,) + E λ, (z,) E λ (u ) P z, u C(ω(z,, 1) + ω (, )) + λ(k 1 (1 )/ + K ū ū ). (4.1) Next, we notice tat by Lemmas 4 and 5, we ave ω(z,, 1) = =( 1, ) 1 i+ 1 N =1 1 j+ N =( 1, ) =1 =( 1, ) =1 =( 1, ) 1 + =1 z N,, z N, (z, ) i+1,j+ (z, ) i,j (Q ) N (Q ) N, by (3.8) (Q ) N (Q ) N ω(z,, 1) 8 (Q ) i+1,j+ (Q ) i,j. =( 1, ) 1 i+ 1 N 1 + =1 1 j+ N From te latter inequality, te deinition o Q, and Jensen s inequality, we obtain ω(z,, 1) 8 =( 1, ) 1 i+ 1 N 1 + =1 1 j+ N 8 8 =( 1, ) 1 + =1 =( 1, ) 1 + =1 1 i+ 1 N 1 j+ N (Q ) i+1,j+ (Q ) i,j 1 (x + ) (x) dx (x + ) (x) dx 8 ω (, ). ) Tus, we ave ω(z,, 1) 4 ω (, ). (4.) Finally, we combine te inequalities (4.1) and (4.), to obtain inequality (4.0) wit K 3 = 5C, K 4 = λk 1, and K 5 = λk. Corollary 1 Let 0 < β 1 be given. I Lip(β, L ()), ten P z, u 0 as 0. Proo Suppose tat Lip(β, L ()). Tere exists M β > 0 suc tat or any 0 < 1, β ω (, ) M β. As a consequence, (4.0) becomes wit = β/ P z, u K 3 M β β + K 4 ( β)/4 + K 5 ū ū. Since ū ū 0 as 0, we iner tat P z, u 0 as 0, and te proo is complete.
18 18 Ming-Jun Lai, Leopold Matamba Messi Fig. : Convergence o ALG1, ALG and ALG3 or te image in Figure 3a wit σ = 5 and λ = 1/13. In algoritms ALG1 and ALG we set τ = 0.05, wile τ = in ALG3. In tis case we used te value o λ tat gave te best PSNRs or our coices o τ. 10 ROF: λ=1/16 ALG1 ALG ALG3 Mean Square Error Iterations 5 Numerical experiments In tis section, we report te result o te numerical experiments wit Algoritm 31. We compare te perormance o te algoritm proposed above to Cambolle s ixed-point algoritm proposed in [4], and te projected-gradient algoritm proposed in [5]. Te test images used are presented in Figure 3. We sall use te ollowing abbreviations to identiy te tree algoritms under consideration ere. ALG1: Te ixed point iterative algoritm described in [4]. ALG: Te projected-gradient algoritm described in [5,8]. ALG3: Te projected-gradient algoritm presented in Algoritm 31. It sould be noted tat in our tests, we did not attempt to coose te parameters τ and λ or optimal perormance o te algoritms. Table 1 troug Table 4 below sow te capability o Algoritm 31 to remove noise or various noise level. Te inputs or all tree algoritms are obtained by adding a zero mean Gaussian noise wit standard deviation σ to te images in Figure 3. Our experiments sow tat te new projected gradient algoritm, Algoritm 31, is sligtly more eicient tan ALG. Finally in Figure, we sow te asymptotic beavior o te tree algoritms: ALG1, ALG, and ALG3. All o te tree algoritms display te same asymptotic beavior as te number o iterations goes to ininity. O course, tis agreement o asymptotic beaviors is expected since all tree algoritms in te discrete setting were proved to converge to te solution o te ROF model. Reerences 1. R. Acar and C.R. Vogel: Analysis o bounded variation penalty metods or ill-posed problems, Inverse Problems 10 (1994), L. Ambrosio, N. Fusco, and D. Pallara, Functions o bounded variation and ree discontinuity problems, Oxord University Press, USA (000). 3. T. F. Can, G. H. Golub, and P. Mulet, A nonlinear primal-dual metod or total variation-based image restoration, SIAM J. Sci. Comput. 0 (1999), no. 6, A. Cambolle, An algoritm or total variation minimization and applications. Special issue on matematics and image analysis, J. Mat. Imaging Vision 0 (004), no. 1-, pp A. Cambolle, Total variation minimization and a class o binary MRF models, in Energy minimization metods in computer vision and pattern recognition, Lecture Notes in Computer Science, Vol. 3757/005 (005), pp A. Cambolle and P.-L. Lions, Image recovery via total variation minimization and related problems, Numer. Mat. 76 (1997), no., D. C. Dobson and C. R. Vogel, Convergence o an iterative metod or total variation denoising, SIAM J. Numer. Anal. 34 (1997), V. Duval, J.-F. Aujol, and L. A. Vese, Matematical modeling o textures: Application to color image decomposition wit a projected gradient algoritm, J Mat Imaging Vis. 37 (010), I. Ekeland and R. Temam, Convex analysis and variational problems, Society or Industrial Matematics, Piladelpia, PA, USA (1987). 10. E. Giusti, Minimal suraces and unctions o bounded variation, Monograps in Matematics, 80, Birkäuser, Basel (1984).
19 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 19 Fig. 3: Te images used in te numerical experiments. Te images in te top row are o size 56 56, wile tose in te bottom row ave resolution (a) Lena (b) Peppers (c) Bank. (d) Boats 11. Q. Hong, M. J. Lai, and J. Wang, Convergence Analysis o a Finite Dierence Sceme or te Gradient Flow associated wit te ROF Model, submitted (010). 1. M. -J. Lai, B. Lucier, and J. Wang, Te convergence o a central-dierence discretization o Rudin-Oser-Fatemi model or image denoising, in Scale Space and Variational Metods in Computer Vision, Lecture Notes in Computer Science, Vol. 5567/009 (009), L.I. Rudin, S. Oser, and E. Fatemi, Nonlinear total variation based noise removal algoritms, Pysica D: Nonlinear Penomena 60 (199), no. 1-4, C. R. Vogel and M. E. Oman, Iterative metods or total variation denoising, SIAM J. Sci. Comput. 17 (1996), no. 1, J. Wang and Lucier, B., Error bounds or inite-dierence metods or Rudin-Oser-Fatemi image smooting, to appear in SIAM J. Numer. Analysis (010). 16. J. Wang, Error bound or numerical metods or te ROF image smooting model, PD Dissertation, Purdue University, IN, USA (008).
20 0 Ming-Jun Lai, Leopold Matamba Messi Table 1: Comparison o te algoritms using te image in Figure 3a. Eac algoritm is terminated wen te dierence between te Mean Square Error (MSE) at two consecutive steps is less tan σ λ τ ALG1 ALG τ ALG3 15 1/ / / / / / / / / / Table : Comparison o te Algoritms using te image in Figure 3b. Eac algoritm is terminated wen te dierence between te Mean Square Error (MSE) at two consecutive steps is less tan σ λ τ ALG1 ALG τ ALG3 15 1/ / / / / / / / / / Table 3: Comparison o te Algoritms using te image in Figure 3c. Eac algoritm is terminated wen te dierence between te Mean Square Error (MSE) at two consecutive steps is less tan σ λ τ ALG1 ALG τ ALG3 15 1/ / / / / / / / / /
21 A New Projected Gradient Algoritm or Total Variation Based Image Denoising 1 Table 4: Comparison o te Algoritms using te image in Figure 3d. Eac algoritm is terminated wen te dierence between te Mean Square Error (MSE) at two consecutive steps is less tan σ λ τ ALG1 ALG τ ALG3 15 1/ / / / / / / / / /
PIECEWISE LINEAR APPROXIMATION OF THE CONTINUOUS RUDIN-OSHER-FATEMI MODEL FOR IMAGE DENOISING
PIECEWISE LINEAR APPROXIMATION OF THE CONTINUOUS RUDIN-OSHER-FATEMI MODEL FOR IMAGE DENOISING MING-JUN LAI AND LEOPOLD MATAMBA MESSI Abstract. Tis paper is concerned wit te numerical approximation of te
More informationThe Convergence of a Central-Difference Discretization of Rudin-Osher-Fatemi Model for Image Denoising
Te Convergence of a Central-Difference Discretization of Rudin-Oser-Fatemi Model for Image Denoising Ming-Jun Lai 1, Bradley Lucier 2, and Jingyue Wang 3 1 University of Georgia, Atens GA 30602, USA mjlai@mat.uga.edu
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationConvergence of a Central Difference Discretization of ROF Model for Image Denoising
Convergence of a Central Difference Discretization of ROF Model for Image Denoising Ming-Jun Lai and Jingyue Wang Department of Matematics Te University of Georgia Atens, GA 30602 October 8, 2009 Abstract
More informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationLinearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization
Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationc 2011 Society for Industrial and Applied Mathematics Key words. total variation, variational methods, finite-difference methods, error bounds
SIAM J. NUMER. ANAL. Vol. 49, No. 2, pp. 845 868 c 211 Society for Industrial and Applied Matematics ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationOnline Appendix for Lerner Symmetry: A Modern Treatment
Online Appendix or Lerner Symmetry: A Modern Treatment Arnaud Costinot MIT Iván Werning MIT May 2018 Abstract Tis Appendix provides te proos o Teorem 1, Teorem 2, and Proposition 1. 1 Perect Competition
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationSOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL
Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationContinuity and Differentiability
Continuity and Dierentiability Tis capter requires a good understanding o its. Te concepts o continuity and dierentiability are more or less obvious etensions o te concept o its. Section - INTRODUCTION
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationDIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS
Journal o Applied Analysis Vol. 14, No. 2 2008, pp. 259 271 DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS B. BELAÏDI and A. EL FARISSI Received December 5, 2007 and,
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationMATH1901 Differential Calculus (Advanced)
MATH1901 Dierential Calculus (Advanced) Capter 3: Functions Deinitions : A B A and B are sets assigns to eac element in A eactl one element in B A is te domain o te unction B is te codomain o te unction
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM
ITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM Abstract. Te discretization of a bilaterally constrained total variation minimization problem wit conforming low order finite
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationChapter 2 Limits and Continuity. Section 2.1 Rates of Change and Limits (pp ) Section Quick Review 2.1
Section. 6. (a) N(t) t (b) days: 6 guppies week: 7 guppies (c) Nt () t t t ln ln t ln ln ln t 8. 968 Tere will be guppies ater ln 8.968 days, or ater nearly 9 days. (d) Because it suggests te number o
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationVariational Methods for Time Dependent Wave Propagation Problems
Variational Metods or Time Dependent Wave Propagation Problems Patrick Joly INRIA Rocquencourt, BP105 Le Cesnay, France (patrick.joly@inria.r) 1 Introduction Tere is an important need or numerical metods
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM
ITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM Abstract. Te discretization of a bilaterally constrained total variation minimization problem wit conforming low order finite
More informationGeneric maximum nullity of a graph
Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More informationDerivation of Some Results on the Generalized Relative Orders of Meromorphic Functions
DOI: 10.1515/awutm-2017-0004 Analele Universităţii de Vest, Timişoara Seria Matematică Inormatică LV, 1, 2017), 51 61 Derivation o Some Results on te Generalized Relative Orders o Meromorpic Functions
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationOn One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations
Pure and Applied Matematics Journal 7; 6(5: 74 ttp://wwwsciencepublisinggroupcom/j/pamj doi: 648/jpamj765 ISSN: 6979 (Print; ISSN: 698 (Online On One Justiication on te Use o Hybrids or te Solution o First
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationICES REPORT Multirate Undrained Splitting for Coupled Flow and Geomechanics in Porous Media
ICES REPORT 16-13 May 2016 Multirate Undrained Splitting or Coupled Flow and Geomecanics in Porous Media by Tameem Almani, Kundan Kumar, Mary F. Weeler Te Institute or Computational Engineering and Sciences
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationEssential Microeconomics : EQUILIBRIUM AND EFFICIENCY WITH PRODUCTION Key ideas: Walrasian equilibrium, first and second welfare theorems
Essential Microeconomics -- 5.2: EQUILIBRIUM AND EFFICIENCY WITH PRODUCTION Key ideas: Walrasian equilibrium, irst and second welare teorems A general model 2 First welare Teorem 7 Second welare teorem
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationObservability Analysis of Nonlinear Systems Using Pseudo-Linear Transformation
9t IFAC Symposium on Nonlinear Control Systems Toulouse, France, September 4-6, 2013 TC3.4 Observability Analysis o Nonlinear Systems Using Pseudo-Linear Transormation Yu Kawano Tosiyui Otsua Osaa University,
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationControllability of a one-dimensional fractional heat equation: theoretical and numerical aspects
Controllability of a one-dimensional fractional eat equation: teoretical and numerical aspects Umberto Biccari, Víctor Hernández-Santamaría To cite tis version: Umberto Biccari, Víctor Hernández-Santamaría.
More informationSemigroups of Operators
Lecture 11 Semigroups of Operators In tis Lecture we gater a few notions on one-parameter semigroups of linear operators, confining to te essential tools tat are needed in te sequel. As usual, X is a real
More informationAndrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction
HOMOGENIZATION OF FREE DISCONTINUITY PROBLEMS Andrea Braides, Anneliese Defrancesci and Enrico Vitali Introduction Following Griffit s teory, yperelastic brittle media subject to fracture can be modeled
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationAnalysis of a second order discontinuous Galerkin finite element method for the Allen-Cahn equation
Graduate Teses and Dissertations Graduate College 06 Analysis of a second order discontinuous Galerin finite element metod for te Allen-Can equation Junzao Hu Iowa State University Follow tis and additional
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationIEOR 165 Lecture 10 Distribution Estimation
IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationTHE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718XX0000-0 THE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We solve te problem of
More information