On the regularizing power of multigrid-type algorithms
|
|
- Elisabeth Sharp
- 6 years ago
- Views:
Transcription
1 On the regularizing power of multigrid-type algorithms Marco Donatelli Dipartimento di Fisica e Matematica Università dell Insubria
2 Outline 1 Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods 2 Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost 3 Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization 4 Conclusions
3 Restoration of blurred and noisy images Outline 1 Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods 2 Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost 3 Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization 4 Conclusions
4 Restoration of blurred and noisy images The model problem Image restoration with Boundary Conditions Using Boundary Conditions (BCs), the restored image f is obtained solving: (in some way...) g = blurred image ξ = noise (random vector) Af = g + ξ A = two-level matrix depending on PSF and BCs BCs Dirichlet periodic Neumann (reflective) anti-reflective A Toeplitz circulant DCT III DST I + low-rank
5 Restoration of blurred and noisy images Properties of the PSF Generating function of PSF Gaussian PSF PSF Generating function z(x) The ill-conditioned subspace is mainly constituted by the high frequencies.
6 Restoration of blurred and noisy images Iterative regularization methods Iterative regularization methods Several iterative methods (e.g., Landweber, CGNE,...) have regularization properties. Reason They firstly reduce the error in the low frequencies (well-conditioned subspace). When they arrive to reduce the algebraic error in the high frequencies then the restoration error increases because of the noise.
7 Multigrid regularization Outline 1 Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods 2 Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost 3 Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization 4 Conclusions
8 Multigrid regularization Multigrid Methods Multigrid structure Multigrid Idea Project the system in a subspace, solve the resulting system in this subspace and interpolate the solution in order to improve the previous approximation. Multigrid components The Multigrid combines two iterative methods: Smoother: a classic iterative method, Coarse Grid Correction: projection, solution of the restricted problem, interpolation. At the lower level(s) it works on the error equation!
9 Multigrid regularization Multigrid Methods The algorithm Two-Grid Methods (TGM) The j-th iteration for the system Ax = b: (1) x = Smooth(A,x (j),b, ν) (2) r 1 = P(b A x) (3) A 1 = PAP H (4) e 1 = A 1 1 r 1 (5) x (j+1) = x (j) + P H e 1 Multigrid (MGM): the step (4) becomes a recursive application of the algorithm.
10 Multigrid regularization Multigrid Methods The Algebraic Multigrid (AMG) The AMG uses information on the coefficient matrix and no geometric information on the problem. Different classic smoothers have a similar behavior: in the initial iterations they are not able to reduce effectively the error in the subspace generated by the eigenvectors associated to small eigenvalues (ill-conditioned subspace) To obtain a fast solver the projector is chosen in order to project the error equation in such subspace. For instance, for Toeplitz and algebra of matrices, see [Aricò, Donatelli, Serra Capizzano, SIMAX, Vol pp ].
11 Multigrid regularization Multigrid Methods Algebraic analysis of the Geometric Multigrid The MGM is an optimal solver for elliptic PDEs For elliptic PDEs the Poisson s problem ill-conditioned subspace is made by low frequencies (complemen tary with respect to the gaussian blur) π 3.5 Projector P = 1 16 P T (full weighting) (linear interpolation)
12 Multigrid regularization Multigrid Methods Multigrid for structured matrices Preserve the structure In order to apply recursively the MGM it is necessary to keep the same structure at each level (Toeplitz,... ). For every structure arising from the proposed BCs there are projectors that preserve the same structure. Projector At i-th level P i = K Ni A Ni (p), where K Ni R N i 4 N i is the cutting matrix that preserves the structure at the lower level. p(x, y) is the generating function of the projector, which selects the subspace where to project the linear system.
13 Multigrid regularization Multigrid Methods Multigrid, structured matrices and images The cutting matrix K ni in 1D circulant Toeplitz&DST I DCT III [ ] [ ] [ ] Low-pass filter Project in the low frequencies reduces the noise effects: 2D p(x, y) = (1 + cos(x))(1 + cos(y)) ց Full weighting ր Linear interpolation
14 Multigrid regularization Iterative Multigrid regularization Image restoration and Multigrid In the image deblurring the ill-conditioned subspace is related to high frequencies, while the well-conditioned subspace is generated by low frequencies. In order to obtain a fast convergence the algebraic multigrid projects in the high frequencies where the noise lives = noise explosion already at the first iteration (it requires Tikhonov regularization [NLAA in press]). In this case the low-pass filter projects in the well-conditioned subspace (low frequencies) = it is slowly convergent but it can be a good iterative regularizer.
15 Multigrid regularization Iterative Multigrid regularization Iterative multigrid regularization The Multigrid as an iterative regularization method If we have an iterative regularization method we can improve its regularizing properties and/or accelerate its convergence using it as smoother in a Multigrid algorithm. Regularization The regularization properties of the smoother are preserved since it is combined with a low-pass filter.
16 Multigrid regularization Iterative Multigrid regularization Two-Level (TL) regularization Idea Project into the low frequencies and then apply an iterative regularization method. TL as a specialization of TGM Smoother: iterative regularization method Projector: low-pass filter TL Algorithm 1 No smoothing at step (1): x = x (j) 2 Step (4): e 1 = A 1 1 r 1 Smooth(A 1,e 1,r 1, 1)
17 Multigrid regularization Iterative Multigrid regularization Multigrid regularization (applying recursively the TL) V-cycle Using a greater number of recursive calls (e.g. W-cycle), the algorithm works more in the well-conditioned subspace but it is more difficult to define an early stopping criterium.
18 Multigrid regularization Computational Cost Computational Cost Precomputing Assumptions: images n n and PSF m m with m n. Computation and memorization of the PSFs (coefficient matrices) for each coarse grids (A i+1 = P T i A i P i ) O(log(n)) V-cycle The problem at level j has size n j n j with n j = n/2 j. Cost of projection at level j: ց 7 4 n2 j, ր 7 8 n2 j flops. Let W(n) the computational cost of one smoother iteration W(n) = cn 2 + O(n), c 1.
19 Multigrid regularization Computational Cost Computational Cost V-cycle (PSF m m with m n) The total cost at level j is about c j = W(n j ) n2 j flops. The total cost of one V-cycle iteration is: log 21 2 (n) 1 8 n2 + c j < 7 2 n2 + 4 ( n ) 3 W W(n). j=1
20 Multigrid regularization Computational Cost Computational cost for more complex cycles Increasing the number of recursive calls Only one recursive call (projection) V scheme. Increasing the number of recursive calls (consecutive projections) the algorithm works more in the low frequencies subspace, however it is more difficult to define a stopping criterium (the method accelerates). γ = number of recursive calls, the cost of one iteration is C(γ,n) C(γ,n) ( cγ ) n 2 4 γ { O(n 2 log(n)), γ = 4 O(n 2log 4 (γ) ), γ > W(n), γ = 1 W(n), γ = 2 3W(n), γ = 3
21 Numerical experiments Outline 1 Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods 2 Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost 3 Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization 4 Conclusions
22 Numerical experiments An airplane An airplane Periodic BCs Gaussian PSF (A spd) noise = 1% Original Image Inner part Observed image Restored with MGM
23 Numerical experiments An airplane Restoration error e j = f f (j) 2 / f 2 restoration error at the j-th iteration. Minimum restoration error Method min j) j=1,... arg min j) j=1,... CG Richardson TL(CG) TL(Rich) MGM(Rich, 1) MGM(Rich, 2) CGNE RichNE
24 Numerical experiments An airplane Noise = 10% For CG and Richardson it is better to resort to normal equations. Minimum restoration error Method min j) j=1,... arg min j) j=1,... CGNE RichNE TL(CGNE) TL(RichNE) MGM(RichNE,1) MGM(RichNE,2) MGM(Rich,1) MGM(Rich,2)
25 Numerical experiments An astronomic example with nonnegativity constraint Saturn Periodic BCs and dark background (exact) Gaussian PSF + noise = 5% Original Image PSF Blurred + SNR = 20
26 Numerical experiments An astronomic example with nonnegativity constraint Nonnegativity constraint Does Multigrid preserve the nonnegativity? The projector (full weighting and linear interpolation) preserves the nonnegativity. Also with a projected smoother, the multigrid does not preserve the nonnegativity (cause: direct resolution at the coarser level). Preserving the nonnegativity with multigrid Project the computed solution at each iteration. It is better if we use a projected smoother too.
27 Numerical experiments An astronomic example with nonnegativity constraint Minimum restoration error Method min j) j=1,... argmin j) j=1,... CG Rich MGM(Rich +,1) MGM(Rich,1) MGM(Rich +,1) MGM(Rich +,2) RichNE CGNE MGM(CGNE +,2) MGM(CGNE +,3)
28 Numerical experiments An astronomic example with nonnegativity constraint Graph of the restoration error
29 Numerical experiments An astronomic example with nonnegativity constraint Restored images Rich + (9 iter.). MGM(Rich +,2) + (12 iter.) CGNE + (885 iter.).
30 Numerical experiments An astronomic example with nonnegativity constraint Restored images MGM(CGNE +,2) + (109 iter.) CGNE + (885 iter.).
31 Numerical experiments Strengthen the projector Strengthen the projector With multigrid we can strengthen the filtering properties of the projector instead of to use a smoother for normal equations: p(x, y) = (1 + cos(x)) α (1 + cos(y)) α, α N +. Computational cost increases The computational cost of one W-Cycle iteration become C(α,n) q(α)n 2, q(x) polynomial of degree 2 Linear system conditioning: let f j be the generating function of A j and β j the order of its unique zero in [π, π], then β j = β 0 + 2jα (cond. numb. = O(n β j j )).
32 Numerical experiments Strengthen the projector Strengthen the projector Projection and generalized re-blurring The projector performs a double re-blurring since A i+1 = P T i A i P i. Minimum restoration error with MGM(Rich +,2) + varying α α min j) j=1,... argmin j) j=1,
33 Numerical experiments Strengthen the projector Strengthen the projector MGM(Rich +,2) + (α = 1 in 12 iter.) MGM(Rich +,2) + (α = 4 in 38 iter.) MGM(CGNE +,2) + (109 iter.)
34 Numerical experiments Direct multigrid regularization The γ regularization Varying γ, the proposed multigrid is a direct (one step) regularization method with regularization parameter γ. The airplane example with noise = 1% e 1 = restoration error after one iteration of MGM(Rich,γ). The CGNE reaches a minimum error equal to after 2500 iterations! γ e
35 Numerical experiments Direct multigrid regularization The γ regularization (the airplane with noise = 1%) MGM(Rich,1) MGM(Rich,2) MGM(Rich,4) MGM(Rich,6) Restoration error norm vs. number of iterations The computational cost increases with γ but not so much (e.g. γ = 8 O(N 1.5 ) where N = n 2 ).
36 Conclusions Outline 1 Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods 2 Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost 3 Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization 4 Conclusions
37 Conclusions Summarizing... Multigrid regularization method It is a general framework which can be used to improve the regularization properties of an iterative regularizing method. It leads to a smaller relative error and a more flat error curve with respect to the smoother applied alone. It is fast and it obtains a good restored image without resort to normal equations. It can be combined with several techniques (e.g., nonnegativity constraints) and it can lead to several generalizations (e.g., γ regularization). Reference M. Donatelli and S. Serra Capizzano, On the regularizing power of multigrid-type algorithms, SIAM J. Sci. Comput., to appear.
38 Conclusions Future work Theoretical Theoretical analysis of the regularization property. Application PSF strictly non-symmetric. Combination with techniques for edge enhancing (Wavelet, Total Variation,... ). Numerical/Simulating A complete experimentation with all the proposed BCs (multigrid methods already exist for the arising matrices).
Image deblurring by multigrid methods. Department of Physics and Mathematics University of Insubria
Image deblurring by multigrid methods Marco Donatelli Stefano Serra-Capizzano Department of Physics and Mathematics University of Insubria Outline 1 Restoration of blurred and noisy images The model problem
More informationOn the regularizing power of multigrid-type algorithms. Marco Donatelli. Stefano Serra Capizzano
On the regularizing power of multigrid-type algorithms Marco Donatelli Stefano Serra Capizzano Università dell Insubria, Dip. Fisica e Matematica - Sede di Como Via Valleggio 11-22100 Como, Italy (donatelli@uninsubria.it)
More informationRegularization by multigrid-type algorithms
Regularization by multigrid-type algorithms Marco Donatelli Department of Physics and Mathematics University of Insubria Joint work with S. Serra-Capizzano Outline 1 Restoration of blurred and noisy images
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 20: Sparse Linear Systems; Direct Methods vs. Iterative Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 26
More informationIntroduction to Multigrid and its Parallelization
Introduction to Multigrid and its Parallelization! Thomas D. Economon Lecture 14a May 28, 2014 Announcements 2 HW 1 & 2 have been returned. Any questions? Final projects are due June 11, 5 pm. If you are
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationPARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS. Ioana Chiorean
5 Kragujevac J. Math. 25 (2003) 5 18. PARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS Ioana Chiorean Babeş-Bolyai University, Department of Mathematics, Cluj-Napoca, Romania (Received May 28,
More informationMultigrid Methods for Markov Chains
Multigrid Methods for Markov Chains Hans De Sterck Department of Applied Mathematics, University of Waterloo collaborators Killian Miller Department of Applied Mathematics, University of Waterloo, Canada
More informationAn edge-preserving multilevel method for deblurring, denoising, and segmentation
An edge-preserving multilevel method for deblurring, denoising, and segmentation S. Morigi 1, L. Reichel 2, and F. Sgallari 1 1 Dept. of Mathematics-CIRAM, University of Bologna, Bologna, Italy {morigi,sgallari}@dm.unibo.it
More informationFast Iterative Solvers for Markov Chains, with Application to Google's PageRank. Hans De Sterck
Fast Iterative Solvers for Markov Chains, with Application to Google's PageRank Hans De Sterck Department of Applied Mathematics University of Waterloo, Ontario, Canada joint work with Steve McCormick,
More informationAdaptive Multiple-Frame Image Super- Resolution Based on U-Curve
Adaptive Multiple-Frame Image Super- Resolution Based on U-Curve IEEE Transaction on Image Processing, Vol. 19, No. 12, 2010 Qiangqiang Yuan, Liangpei Zhang, Huanfeng Shen, and Pingxiang Li Presented by
More informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
More informationPROJECT REPORT. Symbol based Multigrid method. Hreinn Juliusson, Johanna Brodin, Tianhao Zhang Project in Computational Science: Report
Symbol based Multigrid method Hreinn Juliusson, Johanna Brodin, Tianhao Zhang Project in Computational Science: Report February 20, 2018 PROJECT REPORT Department of Information Technology Abstract Discretization
More informationA Hybrid Geometric+Algebraic Multigrid Method with Semi-Iterative Smoothers
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 013; 00:1 18 Published online in Wiley InterScience www.interscience.wiley.com). A Hybrid Geometric+Algebraic Multigrid Method with
More informationSmoothers. < interactive example > Partial Differential Equations Numerical Methods for PDEs Sparse Linear Systems
Smoothers Partial Differential Equations Disappointing convergence rates observed for stationary iterative methods are asymptotic Much better progress may be made initially before eventually settling into
More informationChemnitz Scientific Computing Preprints
Roman Unger Obstacle Description with Radial Basis Functions for Contact Problems in Elasticity CSC/09-01 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN
More informationAdaptive Smoothed Aggregation (αsa) Multigrid
SIAM REVIEW Vol. 47,No. 2,pp. 317 346 c 2005 Society for Industrial and Applied Mathematics Adaptive Smoothed Aggregation (αsa) Multigrid M. Brezina R. Falgout S. MacLachlan T. Manteuffel S. McCormick
More informationPRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING. Nam-Yong Lee. and Bradley J. Lucier
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX PRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING Nam-Yong Lee Department
More informationWhat is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear.
AMSC 600/CMSC 760 Fall 2007 Solution of Sparse Linear Systems Multigrid, Part 1 Dianne P. O Leary c 2006, 2007 What is Multigrid? Originally, multigrid algorithms were proposed as an iterative method to
More informationThe Pre-Image Problem in Kernel Methods
The Pre-Image Problem in Kernel Methods James Kwok Ivor Tsang Department of Computer Science Hong Kong University of Science and Technology Hong Kong The Pre-Image Problem in Kernel Methods ICML-2003 1
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More informationWhy Use the GPU? How to Exploit? New Hardware Features. Sparse Matrix Solvers on the GPU: Conjugate Gradients and Multigrid. Semiconductor trends
Imagine stream processor; Bill Dally, Stanford Connection Machine CM; Thinking Machines Sparse Matrix Solvers on the GPU: Conjugate Gradients and Multigrid Jeffrey Bolz Eitan Grinspun Caltech Ian Farmer
More informationPRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING. Nam-Yong Lee. Bradley J. Lucier. (Communicated by Hao-Min Zhou)
Inverse Problems and Imaging Volume 10, No. 1, 2016, 195 225 doi:10.3934/ipi.2016.10.195 PRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING Nam-Yong Lee Department of
More informationADAPTIVE SMOOTHED AGGREGATION (αsa)
ADAPTIVE SMOOTHED AGGREGATION (αsa) M. BREZINA, R. FALGOUT, S. MACLACHLAN, T. MANTEUFFEL, S. MCCORMICK, AND J. RUGE Abstract. Substantial effort has been focused over the last two decades on developing
More informationSuppression of Missing Data Artifacts. for Deblurring Images Corrupted by
Suppression of Missing Data Artifacts for Deblurring Images Corrupted by Random Valued Noise Nam-Yong Lee Department of Applied Mathematics Institute of Basic Sciences, Inje University Gimhae, Gyeongnam
More informationLocally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling
Locally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling Moritz Baecher May 15, 29 1 Introduction Edge-preserving smoothing and super-resolution are classic and important
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationEfficient Finite Element Geometric Multigrid Solvers for Unstructured Grids on GPUs
Efficient Finite Element Geometric Multigrid Solvers for Unstructured Grids on GPUs Markus Geveler, Dirk Ribbrock, Dominik Göddeke, Peter Zajac, Stefan Turek Institut für Angewandte Mathematik TU Dortmund,
More information.. Lecture 2. learning and regularization. from interpolation to approximation.
.. Lecture. learning and regularization. from interpolation to approximation. Stéphane Canu and Cheng Soon Ong stephane.canu@insarouen.fr asi.insarouen.fr~scanu. RSISE ANU NICTA, Canberra INSA, Rouen RSISE,
More informationSolution of Ambrosio-Tortorelli Model for Image Segmentation by Generalized Relaxation Method
Solution of Ambrosio-Tortorelli Model for Image Segmentation by Generalized Relaxation Method Pasqua D Ambra, Gaetano Tartaglione Institute for High-Performance Computing and Networking (ICAR), CNR, Via
More informationLecture 15: More Iterative Ideas
Lecture 15: More Iterative Ideas David Bindel 15 Mar 2010 Logistics HW 2 due! Some notes on HW 2. Where we are / where we re going More iterative ideas. Intro to HW 3. More HW 2 notes See solution code!
More informationTowards a complete FEM-based simulation toolkit on GPUs: Geometric Multigrid solvers
Towards a complete FEM-based simulation toolkit on GPUs: Geometric Multigrid solvers Markus Geveler, Dirk Ribbrock, Dominik Göddeke, Peter Zajac, Stefan Turek Institut für Angewandte Mathematik TU Dortmund,
More informationLaboratorio di Problemi Inversi Esercitazione 3: regolarizzazione iterativa, metodo di Landweber
Laboratorio di Problemi Inversi Esercitazione 3: regolarizzazione iterativa, metodo di Landweber Luca Calatroni Dipartimento di Matematica, Universitá degli studi di Genova May 18, 2016. Luca Calatroni
More informationTHE DEVELOPMENT OF THE POTENTIAL AND ACADMIC PROGRAMMES OF WROCLAW UNIVERISTY OF TECH- NOLOGY ITERATIVE LINEAR SOLVERS
ITERATIVE LIEAR SOLVERS. Objectives The goals of the laboratory workshop are as follows: to learn basic properties of iterative methods for solving linear least squares problems, to study the properties
More informationAdvanced Computer Graphics
G22.2274 001, Fall 2010 Advanced Computer Graphics Project details and tools 1 Projects Details of each project are on the website under Projects Please review all the projects and come see me if you would
More informationsmooth coefficients H. Köstler, U. Rüde
A robust multigrid solver for the optical flow problem with non- smooth coefficients H. Köstler, U. Rüde Overview Optical Flow Problem Data term and various regularizers A Robust Multigrid Solver Galerkin
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationPhD Student. Associate Professor, Co-Director, Center for Computational Earth and Environmental Science. Abdulrahman Manea.
Abdulrahman Manea PhD Student Hamdi Tchelepi Associate Professor, Co-Director, Center for Computational Earth and Environmental Science Energy Resources Engineering Department School of Earth Sciences
More informationADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N.
ADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N. Dartmouth, MA USA Abstract: The significant progress in ultrasonic NDE systems has now
More informationEfficient Imaging Algorithms on Many-Core Platforms
Efficient Imaging Algorithms on Many-Core Platforms H. Köstler Dagstuhl, 22.11.2011 Contents Imaging Applications HDR Compression performance of PDE-based models Image Denoising performance of patch-based
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationContents. Implementing the QR factorization The algebraic eigenvalue problem. Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationA Comparison of Algebraic Multigrid Preconditioners using Graphics Processing Units and Multi-Core Central Processing Units
A Comparison of Algebraic Multigrid Preconditioners using Graphics Processing Units and Multi-Core Central Processing Units Markus Wagner, Karl Rupp,2, Josef Weinbub Institute for Microelectronics, TU
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More information2D vector fields 3. Contents. Line Integral Convolution (LIC) Image based flow visualization Vector field topology. Fast LIC Oriented LIC
2D vector fields 3 Scientific Visualization (Part 8) PD Dr.-Ing. Peter Hastreiter Contents Line Integral Convolution (LIC) Fast LIC Oriented LIC Image based flow visualization Vector field topology 2 Applied
More informationApplication of the topological gradient method to color image restoration and classification
Application of the topological gradient method to color image restoration and classification Didier Auroux a Lamia Jaafar Belaid b a Institut de Mathématiques de Toulouse, Université Paul Sabatier Toulouse
More informationMultigrid solvers M. M. Sussman sussmanm@math.pitt.edu Office Hours: 11:10AM-12:10PM, Thack 622 May 12 June 19, 2014 1 / 43 Multigrid Geometrical multigrid Introduction Details of GMG Summary Algebraic
More informationA hybrid GMRES and TV-norm based method for image restoration
A hybrid GMRES and TV-norm based method for image restoration D. Calvetti a, B. Lewis b and L. Reichel c a Department of Mathematics, Case Western Reserve University, Cleveland, OH 44106 b Rocketcalc,
More informationIterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms
Iterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms By:- Nitin Kamra Indian Institute of Technology, Delhi Advisor:- Prof. Ulrich Reude 1. Introduction to Linear
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationClustering: Classic Methods and Modern Views
Clustering: Classic Methods and Modern Views Marina Meilă University of Washington mmp@stat.washington.edu June 22, 2015 Lorentz Center Workshop on Clusters, Games and Axioms Outline Paradigms for clustering
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationBlur Space Iterative De-blurring
Blur Space Iterative De-blurring RADU CIPRIAN BILCU 1, MEJDI TRIMECHE 2, SAKARI ALENIUS 3, MARKKU VEHVILAINEN 4 1,2,3,4 Multimedia Technologies Laboratory, Nokia Research Center Visiokatu 1, FIN-33720,
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn School of Mathematical Sciences Tel Aviv University Michael S. Floater Department of Informatics University of
More informationData-driven modeling: A low-rank approximation problem
1 / 34 Data-driven modeling: A low-rank approximation problem Ivan Markovsky Vrije Universiteit Brussel 2 / 34 Outline Setup: data-driven modeling Problems: system identification, machine learning,...
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationDistributed NVAMG. Design and Implementation of a Scalable Algebraic Multigrid Framework for a Cluster of GPUs
Distributed NVAMG Design and Implementation of a Scalable Algebraic Multigrid Framework for a Cluster of GPUs Istvan Reguly (istvan.reguly at oerc.ox.ac.uk) Oxford e-research Centre NVIDIA Summer Internship
More informationLagrangian methods for the regularization of discrete ill-posed problems. G. Landi
Lagrangian methods for the regularization of discrete ill-posed problems G. Landi Abstract In many science and engineering applications, the discretization of linear illposed problems gives rise to large
More informationPartition of unity algorithm for two-dimensional interpolation using compactly supported radial basis functions
Communications in Applied and Industrial Mathematics, ISSN 2038-0909, e-431 DOI: 10.1685/journal.caim.431 Partition of unity algorithm for two-dimensional interpolation using compactly supported radial
More informationparameter. Systems with this structure arise, for example, in image deblurring. Blurred, noisy image data is modeled by d ij = Z Z k(x i? x; y j? y)u(
Two-Level Preconditioners for Ill-Conditioned Linear Systems with Semidenite Regularization Kyle L. Riley a Curtis R. Vogel b a Department of Mathematical Sciences, Montana State University, Bozeman, MT
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationAN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES
AN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES Nader Moayeri and Konstantinos Konstantinides Hewlett-Packard Laboratories 1501 Page Mill Road Palo Alto, CA 94304-1120 moayeri,konstant@hpl.hp.com
More informationA Lagrange method based L-curve for image restoration
Journal of Physics: Conference Series OPEN ACCESS A Lagrange method based L-curve for image restoration To cite this article: G Landi 2013 J. Phys.: Conf. Ser. 464 012011 View the article online for updates
More informationMultigrid. James Demmel. Poisson s equation in 1D: T = Graph and stencil. 2 u/ x 2 = f(x)
Multigrid James Demmel www.cs.berkeley.edu/~demmel/ma221_spr16 Poisson s equation in 1D: 2 u/ x 2 = f(x) T = 2-1 -1 2-1 -1 2-1 -1 2-1 -1 2 Graph and stencil -1 2-1 1 2D Poisson s equation Similar to the
More informationAmgX 2.0: Scaling toward CORAL Joe Eaton, November 19, 2015
AmgX 2.0: Scaling toward CORAL Joe Eaton, November 19, 2015 Agenda Introduction to AmgX Current Capabilities Scaling V2.0 Roadmap for the future 2 AmgX Fast, scalable linear solvers, emphasis on iterative
More informationHighly Parallel Multigrid Solvers for Multicore and Manycore Processors
Highly Parallel Multigrid Solvers for Multicore and Manycore Processors Oleg Bessonov (B) Institute for Problems in Mechanics of the Russian Academy of Sciences, 101, Vernadsky Avenue, 119526 Moscow, Russia
More informationData-driven modeling: A low-rank approximation problem
1 / 31 Data-driven modeling: A low-rank approximation problem Ivan Markovsky Vrije Universiteit Brussel 2 / 31 Outline Setup: data-driven modeling Problems: system identification, machine learning,...
More informationSegmentation of MR Images of a Beating Heart
Segmentation of MR Images of a Beating Heart Avinash Ravichandran Abstract Heart Arrhythmia is currently treated using invasive procedures. In order use to non invasive procedures accurate imaging modalities
More informationAPPROXIMATING PDE s IN L 1
APPROXIMATING PDE s IN L 1 Veselin Dobrev Jean-Luc Guermond Bojan Popov Department of Mathematics Texas A&M University NONLINEAR APPROXIMATION TECHNIQUES USING L 1 Texas A&M May 16-18, 2008 Outline 1 Outline
More informationACCELERATING CFD AND RESERVOIR SIMULATIONS WITH ALGEBRAIC MULTI GRID Chris Gottbrath, Nov 2016
ACCELERATING CFD AND RESERVOIR SIMULATIONS WITH ALGEBRAIC MULTI GRID Chris Gottbrath, Nov 2016 Challenges What is Algebraic Multi-Grid (AMG)? AGENDA Why use AMG? When to use AMG? NVIDIA AmgX Results 2
More informationPROGRAMMING OF MULTIGRID METHODS
PROGRAMMING OF MULTIGRID METHODS LONG CHEN In this note, we explain the implementation detail of multigrid methods. We will use the approach by space decomposition and subspace correction method; see Chapter:
More informationEfficient O(N log N) algorithms for scattered data interpolation
Efficient O(N log N) algorithms for scattered data interpolation Nail Gumerov University of Maryland Institute for Advanced Computer Studies Joint work with Ramani Duraiswami February Fourier Talks 2007
More informationBlind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization
Blind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization Amandeep Kaur CEM Kapurthala, Punjab Vinay Chopra DAVIET Jalandhar Punjab ABSTRACT Image Restoration is a field of
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationImage Deconvolution.
Image Deconvolution. Mathematics of Imaging. HW3 Jihwan Kim Abstract This homework is to implement image deconvolution methods, especially focused on a ExpectationMaximization(EM) algorithm. Most of this
More informationSuperresolution Using Preconditioned Conjugate Gradient Method
Superresolution Using Preconditioned Conjugate Gradient Method Changjiang Yang, Ramani Duraiswami and Larry Davis Computer Vision Laboratory University of Maryland, College Park, MD 20742 E-mail: {yangcj,ramani,lsd}@umiacs.umd.edu
More informationvector space retrieval many slides courtesy James Amherst
vector space retrieval many slides courtesy James Allan@umass Amherst 1 what is a retrieval model? Model is an idealization or abstraction of an actual process Mathematical models are used to study the
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn Michael S. Floater Kai Hormann Abstract. We present a new four-point subdivision scheme that generates C 2 curves.
More informationA Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems
A Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems Panos Parpas May 9, 2017 Abstract Composite optimization models consist of the minimization of the sum of a smooth
More information(Sparse) Linear Solvers
(Sparse) Linear Solvers Ax = B Why? Many geometry processing applications boil down to: solve one or more linear systems Parameterization Editing Reconstruction Fairing Morphing 2 Don t you just invert
More informationAn algebraic multigrid method for isotropic linear elasticity problems on anisotropic meshes
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed. Engng. 2010; 26:534 553 Published online 10 July 2008 in Wiley InterScience (www.interscience.wiley.com)..1140
More informationSolving basis pursuit: infeasible-point subgradient algorithm, computational comparison, and improvements
Solving basis pursuit: infeasible-point subgradient algorithm, computational comparison, and improvements Andreas Tillmann ( joint work with D. Lorenz and M. Pfetsch ) Technische Universität Braunschweig
More informationBlind Image Deblurring Using Dark Channel Prior
Blind Image Deblurring Using Dark Channel Prior Jinshan Pan 1,2,3, Deqing Sun 2,4, Hanspeter Pfister 2, and Ming-Hsuan Yang 3 1 Dalian University of Technology 2 Harvard University 3 UC Merced 4 NVIDIA
More informationSubdivision Curves and Surfaces
Subdivision Surfaces or How to Generate a Smooth Mesh?? Subdivision Curves and Surfaces Subdivision given polyline(2d)/mesh(3d) recursively modify & add vertices to achieve smooth curve/surface Each iteration
More informationA node-based agglomeration AMG solver for linear elasticity in thin bodies
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:1 6 [Version: 2002/09/18 v1.02] A node-based agglomeration AMG solver for linear elasticity in thin bodies Prasad
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationEfficient Linear System Solvers for Mesh Processing
Efficient Linear System Solvers for Mesh Processing Mario Botsch, David Bommes, and Leif Kobbelt Computer Graphics Group, RWTH Aachen Technical University Abstract. The use of polygonal mesh representations
More informationBindel, Spring 2015 Numerical Analysis (CS 4220) Figure 1: A blurred mystery photo taken at the Ithaca SPCA. Proj 2: Where Are My Glasses?
Figure 1: A blurred mystery photo taken at the Ithaca SPCA. Proj 2: Where Are My Glasses? 1 Introduction The image in Figure 1 is a blurred version of a picture that I took at the local SPCA. Your mission
More informationResilient geometric finite-element multigrid algorithms using minimised checkpointing
Resilient geometric finite-element multigrid algorithms using minimised checkpointing Dominik Göddeke, Mirco Altenbernd, Dirk Ribbrock Institut für Angewandte Mathematik (LS3) Fakultät für Mathematik TU
More informationNatural image modeling using complex wavelets
Natural image modeling using complex wavelets André Jalobeanu Bayesian Learning group, NASA Ames - Moffett Field, CA Laure Blanc-Féraud, Josiane Zerubia Ariana, CNRS / INRIA / UNSA - Sophia Antipolis,
More informationImplementation of the Non-Local Bayes (NL-Bayes) Image Denoising Algorithm
2014/07/01 v0.5 IPOL article class Published in Image Processing On Line on 2013 06 17. Submitted on 2012 05 22, accepted on 2013 02 14. ISSN 2105 1232 c 2013 IPOL & the authors CC BY NC SA This article
More informationAn algebraic multi-grid implementation in FEniCS for solving 3d fracture problems using XFEM
An algebraic multi-grid implementation in FEniCS for solving 3d fracture problems using XFEM Axel Gerstenberger, Garth N. Wells, Chris Richardson, David Bernstein, and Ray Tuminaro FEniCS 14, Paris, France
More informationRecent Developments in Model-based Derivative-free Optimization
Recent Developments in Model-based Derivative-free Optimization Seppo Pulkkinen April 23, 2010 Introduction Problem definition The problem we are considering is a nonlinear optimization problem with constraints:
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
More informationSuper-Resolution. Many slides from Miki Elad Technion Yosi Rubner RTC and more
Super-Resolution Many slides from Mii Elad Technion Yosi Rubner RTC and more 1 Example - Video 53 images, ratio 1:4 2 Example Surveillance 40 images ratio 1:4 3 Example Enhance Mosaics 4 5 Super-Resolution
More informationImage Restoration under Significant Additive Noise
1 Image Restoration under Significant Additive Noise Wenyi Zhao 1 and Art Pope Sarnoff Corporation 01 Washington Road, Princeton, NJ 08540, USA email: { wzhao, apope }@ieee.org Tel: 408-53-178 Abstract
More information