Multigrid Acceleration of the Horn-Schunck Algorithm for the Optical Flow Problem
|
|
- Sandra Chapman
- 6 years ago
- Views:
Transcription
1 Multigrid Acceleration of the Horn-Schunck Algorithm for the Optical Flow Problem El Mostafa Kalmoun Ulrich Ruede Institute of Computer Science 10 Friedrich Alexander University Germany
2 Overview Introduction and Related Works Optical Flow Constraints Regularization Horn Schunck Algorithm Multigrid Scheme A Simple Illustration Variational Multigrid Experimental Results Slide 2 of 23
3 Introduction and Related Works Frame at time t y x Slide 3 of 23
4 Introduction and Related Works Frame at time t Frame at time t + dt y dy dx y+dy x x+dx Slide 3 of 23
5 Introduction and Related Works Frame at time t Frame at time t + dt y dy dx y+dy x x+dx The optical flow at the pixel (x,y) is the 2D-velocity vector (u,v) = ( dx dt, dy dt ) Slide 3 of 23
6 Introduction and Related Works Frame at time t Frame at time t + dt y dy dx y+dy x x+dx The optical flow at the pixel (x,y) is the 2D-velocity vector (u,v) = ( dx dt, dy dt ) An approximation of the 2D-motion Slide 3 of 23
7 Introduction and Related Works Frame at time t Frame at time t + dt y dy dx y+dy x x+dx The optical flow at the pixel (x,y) is the 2D-velocity vector (u,v) = ( dx dt, dy dt ) An approximation of the 2D-motion Optical Flow Motion Slide 3 of 23
8 Introduction and Related Works Term originates with James Gibson in 1979 Quadratic smoothness Horn Schunck (1981) Registration technique with local constraints Lucas Kanade (1981) Oriented smoothness Nagel Enkelmann ( ) Multigrid relaxation Terzopolous (1986) Performance evaluation of popular algorithms Barron Fleet Beauchemin (1994) General anisotropic smoothness Weickert (1996) Optimal control framework Borzi Ito Kunisch (2002) Slide 4 of 23
9 Optical Flow Constraints I(x,y,t) : The image intensity of the pixel (x,y) at time t. I x,i y,i t : Spatial and temporal derivatives of I. Slide 5 of 23
10 Optical Flow Constraints I(x,y,t) : The image intensity of the pixel (x,y) at time t. I x,i y,i t : Spatial and temporal derivatives of I. Assumption : Objects keep the same intensity over time Slide 5 of 23
11 Optical Flow Constraints I(x,y,t) : The image intensity of the pixel (x,y) at time t. I x,i y,i t : Spatial and temporal derivatives of I. Assumption : Objects keep the same intensity over time I(x,y,t) = I(x + dx,y + dy,t + dt) Slide 5 of 23
12 Optical Flow Constraints I(x,y,t) : The image intensity of the pixel (x,y) at time t. I x,i y,i t : Spatial and temporal derivatives of I. Assumption : Objects keep the same intensity over time I(x,y,t) = I(x + dx,y + dy,t + dt) Taylor Expansion = I x u + I y v + I t = 0 The optical flow constraint equation (OFCE) Slide 5 of 23
13 Optical Flow Constraints An equivalent form of the (OFCE) : I. w = I t DI.( w,1) = 0 where I = (I x,i y ), DI = (I x,i y,i t ) and w = (u,v) Slide 6 of 23
14 Optical Flow Constraints An equivalent form of the (OFCE) : u I. w = I t DI.( w,1) = 0 where I = (I x,i y ), DI = (I x,i y,i t ) and w = (u,v) v Slide 6 of 23
15 Optical Flow Constraints An equivalent form of the (OFCE) : u I. w = I t DI.( w,1) = 0 where I = (I x,i y ), DI = (I x,i y,i t ) I and w = (u,v) v Slide 6 of 23
16 Optical Flow Constraints An equivalent form of the (OFCE) : I. w = I t DI.( w,1) = 0 I t / I u where I = (I x,i y ), DI = (I x,i y,i t ) I and w = (u,v) v Slide 6 of 23
17 Optical Flow Constraints An equivalent form of the (OFCE) : I. w = I t DI.( w,1) = 0 I t / I u where I = (I x,i y ), DI = (I x,i y,i t ) I and w = (u,v) v OFC line Slide 6 of 23
18 Optical Flow Constraints An equivalent form of the (OFCE) : I. w = I t DI.( w,1) = 0 I t / I u where I = (I x,i y ), DI = (I x,i y,i t ) w 1 I and w = (u,v) v OFC line Slide 6 of 23
19 Optical Flow Constraints An equivalent form of the (OFCE) : where I. w = I t DI.( w,1) = 0 I t / I w 2 u I = (I x,i y ), DI = (I x,i y,i t ) w 1 I and w = (u,v) v OFC line Slide 6 of 23
20 Optical Flow Constraints An equivalent form of the (OFCE) : where I. w = I t DI.( w,1) = 0 I t / I w 2 u I = (I x,i y ), DI = (I x,i y,i t ) w 1 I and w = (u,v) v OFC line We can only calculate the normal component of the velocity w and not the tangent flow Slide 6 of 23
21 Optical Flow Constraints An equivalent form of the (OFCE) : where I. w = I t DI.( w,1) = 0 I t / I w 2 u I = (I x,i y ), DI = (I x,i y,i t ) w 1 I and w = (u,v) v OFC line We can only calculate the normal component of the velocity w and not the tangent flow The aperture problem. Slide 6 of 23
22 Regularization The (OFCE) is replaced by min (u,v) x y E(u,v)dxdy (P) Slide 7 of 23
23 Regularization The (OFCE) is replaced by where min (u,v) x y E(u,v)dxdy (P) E = E d + αe r, E d (u,v) = (I x u + I y v + I t ) 2 (the data term) E r is the regularization term Slide 7 of 23
24 Regularization The (OFCE) is replaced by where min (u,v) x y E(u,v)dxdy (P) E = E d + αe r, E d (u,v) = (I x u + I y v + I t ) 2 (the data term) E r is the regularization term and α is a positive scalar for adjustment between E d and E r Slide 7 of 23
25 Regularization The (OFCE) is replaced by where min (u,v) x y E(u,v)dxdy (P) E = E d + αe r, E d (u,v) = (I x u + I y v + I t ) 2 (the data term) E r is the regularization term and α is a positive scalar for adjustment between E d and E r The method for solving (P) will depend on the choice of E r. Slide 7 of 23
26 Horn-Schunck Algorithm A standard choice of E r is the isotropic stabilizer : E r (u,v) = u 2 + v 2 Slide 8 of 23
27 Horn-Schunck Algorithm A standard choice of E r is the isotropic stabilizer : From calculus of variations, we get E r (u,v) = u 2 + v 2 α u I x (I x u + I y v + I t ) = 0 α v I y (I x u + I y v + I t ) = 0 Slide 8 of 23
28 Horn-Schunck Algorithm A standard choice of E r is the isotropic stabilizer : From calculus of variations, we get E r (u,v) = u 2 + v 2 α u I x (I x u + I y v + I t ) = 0 α v I y (I x u + I y v + I t ) = 0 We discretize the Laplacien by the standard 5 point stencil and use w = w w where w is the average of the neighbors. Slide 8 of 23
29 Horn-Schunck Algorithm Spatial and Temporal Image Derivatives Masks Slide 9 of 23
30 Horn-Schunck Algorithm Spatial and Temporal Image Derivatives Masks M x = 1 4 ( ) ( 1 1 M y = ) M t = 1 4 ( ) Slide 9 of 23
31 Horn-Schunck Algorithm Spatial and Temporal Image Derivatives Masks M x = 1 4 ( ) ( 1 1 M y = ) M t = 1 4 ( I x = M x (I 1 + I 2 ) I y = M y (I 1 + I 2 ) I t = M t (I 2 I 1 ) ) Coupled Gauss-Seidel relaxation Slide 9 of 23
32 Horn-Schunck Algorithm Spatial and Temporal Image Derivatives Masks M x = 1 4 ( ) ( 1 1 M y = ) M t = 1 4 ( I x = M x (I 1 + I 2 ) I y = M y (I 1 + I 2 ) I t = M t (I 2 I 1 ) ) Coupled Gauss-Seidel relaxation u k+1 = ū k I x I x ū k + I y v k + I t α + I 2 x + I 2 y v k+1 = v k I y I x ū k + I y v k + I t α + I 2 x + I 2 y Slide 9 of 23
33 Multigrid Scheme Multigrid Components Slide 10 of 23
34 Multigrid Scheme Multigrid Components Vertex-centered grid and standard coarsening Slide 10 of 23
35 Multigrid Scheme Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Slide 10 of 23
36 Multigrid Scheme Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Slide 10 of 23
37 Multigrid Scheme Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Discretization coarse grid approximation (DCA approach) Slide 10 of 23
38 A Simple Illustration Intensity value : I(x,y,t) = x + y +t Slide 11 of 23
39 A Simple Illustration Intensity value : I(x,y,t) = x + y +t Slide 11 of 23
40 A Simple Illustration Intensity value : I(x,y,t) = x + y +t Corresponding system of PDE s : α u = u + v + 1 α v = u + v + 1 Slide 11 of 23
41 A Simple Illustration Intensity value : I(x,y,t) = x + y +t Corresponding system of PDE s : α u = u + v + 1 α v = u + v + 1 Slide 11 of 23
42 A Simple Illustration Intensity value : I(x,y,t) = x + y +t Corresponding system of PDE s : α u = u + v + 1 α v = u + v + 1 Method α = 1 and u 0 v 0 Convergence rate Horn-Schunck V(1,0) V(1,1) V(2,1) V(3,3) Slide 11 of 23
43 Variational Multigrid Write the system of the two PDE s as ( u L v ) = ( f g ) Slide 12 of 23
44 Variational Multigrid Write the system of the two PDE s as ( u L v ) = ( f g ) where f = I x I t, g = I y I t ( α 0 L = L d + L r, L r = 0 α and L d is a 2x2 block diagonal matrix with entries ) ( I 2 x I x I y I x I y I 2 y ) Slide 12 of 23
45 Variational Multigrid Write the system of the two PDE s as ( u L v ) = ( f g ) where f = I x I t, g = I y I t ( α 0 L = L d + L r, L r = 0 α and L d is a 2x2 block diagonal matrix with entries ) ( I 2 x I x I y I x I y I 2 y ) L d (positive semi-definite) Slide 12 of 23
46 Variational Multigrid Write the system of the two PDE s as ( u L v ) = ( f g ) where f = I x I t, g = I y I t ( α 0 L = L d + L r, L r = 0 α and L d is a 2x2 block diagonal matrix with entries ) ( I 2 x I x I y I x I y I 2 y ) L d (positive semi-definite) + L r (positive definite) Slide 12 of 23
47 Variational Multigrid Write the system of the two PDE s as ( u L v ) = ( f g ) where f = I x I t, g = I y I t ( α 0 L = L d + L r, L r = 0 α and L d is a 2x2 block diagonal matrix with entries ) ( I 2 x I x I y I x I y I 2 y ) L d (positive semi-definite) + L r (positive definite) = L (positive definite). Slide 12 of 23
48 Variational Multigrid Variational Minimization Form Lw = F w = argmin a(z) Ω Slide 13 of 23
49 Variational Multigrid Variational Minimization Form Lw = F w = argmin a(z) Ω a(z) = 1 2 (Lz,z) (F,z) Slide 13 of 23
50 Variational Multigrid Variational Minimization Form Lw = F w = argmin a(z) Ω a(z) = 1 2 (Lz,z) (F,z) Galerkin Approach Let I h H : ΩH Ω h be a full rank linear mapping. Slide 13 of 23
51 Variational Multigrid Variational Minimization Form Lw = F w = argmin a(z) Ω a(z) = 1 2 (Lz,z) (F,z) Galerkin Approach Let I h H : ΩH Ω h be a full rank linear mapping. An optimal coarse grid correction I h H w H of w h is characterized by ((I h H) T L h I h H)w H = (I h H) T (F L h w h ) Slide 13 of 23
52 Variational Multigrid The CGO is chosen then as follows L H = I H h L h I h H and I H h = (I h H) T Slide 14 of 23
53 Variational Multigrid The CGO is chosen then as follows L H = I H h L h I h H and I H h = (I h H) T For our system, we get L H = = ( R 0 0 R )( L 1 h L 2 h L 2 h ( RL 1 h P RLh 2 P RLh 2 P RL3 h P ) L 3 h )( P 0 0 P ) Slide 14 of 23
54 Variational Multigrid Multigrid Components Slide 15 of 23
55 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Slide 15 of 23
56 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) Slide 15 of 23
57 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme Slide 15 of 23
58 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES Slide 15 of 23
59 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES COMPUTE SPATIAL and TEMPORAL FRAMES DERIVATIVES Slide 15 of 23
60 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES COMPUTE SPATIAL and TEMPORAL FRAMES DERIVATIVES CALCULATE the FINEST GRID OPERATOR and the RHS Slide 15 of 23
61 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES COMPUTE SPATIAL and TEMPORAL FRAMES DERIVATIVES CALCULATE the FINEST GRID OPERATOR and the RHS SET UP the COARSE GRID OPERATORS HIEARACHY with GALERKIN Slide 15 of 23
62 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES COMPUTE SPATIAL and TEMPORAL FRAMES DERIVATIVES CALCULATE the FINEST GRID OPERATOR and the RHS SOLVE with MULTIGRID V-CYCLING SET UP the COARSE GRID OPERATORS HIEARACHY with GALERKIN Slide 15 of 23
63 Variational Multigrid Multigrid Components Vertex-centered grid and standard coarsening Coupled lexicographic point Gauss-Seidel smoother Full weighting and bilinear interpolation Galerkin coarse grid approximation (GCA approach) The General Scheme READ FRAMES COMPUTE SPATIAL and TEMPORAL FRAMES DERIVATIVES CALCULATE the FINEST GRID OPERATOR and the RHS WRITE the OPTICAL FLOW FIELD SOLVE with MULTIGRID V-CYCLING SET UP the COARSE GRID OPERATORS HIEARACHY with GALERKIN Slide 15 of 23
64 Experimental Results Slide 16 of 23
65 Experimental Results A rotating sphere (128x128) Slide 16 of 23
66 Experimental Results A rotating sphere (128x128) the flow field Slide 16 of 23
67 Experimental Results the flow field after scaling Slide 17 of 23
68 Experimental Results The marble sequence (512x512) Slide 18 of 23
69 Experimental Results The marble sequence (512x512) the flow field Slide 18 of 23
70 Experimental Results the flow field after scaling Slide 19 of 23
71 Experimental Results The Hamburg taxi sequence (256x190) Slide 20 of 23
72 Experimental Results The Hamburg taxi sequence (256x190) the flow field Slide 20 of 23
73 Experimental Results the flow field after scaling Slide 21 of 23
74 Experimental Results α = 5 and u 0 = v 0 = 0 on P4 2.4GHz Sphere Marble Taxi ρ CPU ρ CPU ρ CPU Horn-Schunck VMG V(2,1) Slide 22 of 23
75 Experimental Results α = 5 and u 0 = v 0 = 0 on P4 2.4GHz Sphere Marble Taxi ρ CPU ρ CPU ρ CPU Horn-Schunck VMG V(2,1) Robust (but not yet optimal convergence) Slide 22 of 23
76 Experimental Results α = 5 and u 0 = v 0 = 0 on P4 2.4GHz Sphere Marble Taxi ρ CPU ρ CPU ρ CPU Horn-Schunck VMG V(2,1) Robust (but not yet optimal convergence) * need matrix dependent transfer grid operators Slide 22 of 23
77 Experimental Results α = 5 and u 0 = v 0 = 0 on P4 2.4GHz Sphere Marble Taxi ρ CPU ρ CPU ρ CPU Horn-Schunck VMG V(2,1) Robust (but not yet optimal convergence) * need matrix dependent transfer grid operators Galerkin leads to high memory costs search for better representation of the CGOs Slide 22 of 23
78 Future work Application in medicine imaging Consider other regularization (e.g. by Nagel, Weickert...) Slide 23 of 23
79 Future work Application in medicine imaging Consider other regularization (e.g. by Nagel, Weickert...) PLEASE WAKE UP! Slide 23 of 23
80 Future work Application in medicine imaging Consider other regularization (e.g. by Nagel, Weickert...) PLEASE WAKE UP! THANKS FOR YOUR ATTENTION Slide 23 of 23
Motion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationDense Optical Flow Estimation from the Monogenic Curvature Tensor
Dense Optical Flow Estimation from the Monogenic Curvature Tensor Di Zang 1, Lennart Wietzke 1, Christian Schmaltz 2, and Gerald Sommer 1 1 Department of Computer Science, Christian-Albrechts-University
More informationMotion Estimation (I)
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationLucas-Kanade Optical Flow. Computer Vision Carnegie Mellon University (Kris Kitani)
Lucas-Kanade Optical Flow Computer Vision 16-385 Carnegie Mellon University (Kris Kitani) I x u + I y v + I t =0 I x = @I @x I y = @I u = dx v = dy I @y t = @I dt dt @t spatial derivative optical flow
More informationElaborazione delle Immagini Informazione multimediale - Immagini. Raffaella Lanzarotti
Elaborazione delle Immagini Informazione multimediale - Immagini Raffaella Lanzarotti OPTICAL FLOW Thanks to prof. Mubarak Shah,UCF 2 Video Video: sequence of frames (images) catch in the time Data: function
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationCS4495/6495 Introduction to Computer Vision. 6B-L1 Dense flow: Brightness constraint
CS4495/6495 Introduction to Computer Vision 6B-L1 Dense flow: Brightness constraint Motion estimation techniques Feature-based methods Direct, dense methods Motion estimation techniques Direct, dense methods
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationVideo and Motion Analysis Computer Vision Carnegie Mellon University (Kris Kitani)
Video and Motion Analysis 16-385 Computer Vision Carnegie Mellon University (Kris Kitani) Optical flow used for feature tracking on a drone Interpolated optical flow used for super slow-mo optical flow
More informationIntroduction to motion correspondence
Introduction to motion correspondence 1 IPAM - UCLA July 24, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Why estimate visual motion? 2 Tracking Segmentation Structure
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationT H E S I S. Computer Engineering May Smoothing of Matrix-Valued Data. Thomas Brox
D I P L O M A Computer Engineering May 2002 T H E S I S Smoothing of Matrix-Valued Data Thomas Brox Computer Vision, Graphics, and Pattern Recognition Group Department of Mathematics and Computer Science
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationUniversity of Illinois at Urbana-Champaign. Multigrid (MG) methods are used to approximate solutions to elliptic partial differential
Title: Multigrid Methods Name: Luke Olson 1 Affil./Addr.: Department of Computer Science University of Illinois at Urbana-Champaign Urbana, IL 61801 email: lukeo@illinois.edu url: http://www.cs.uiuc.edu/homes/lukeo/
More informationMethods in Computer Vision: Introduction to Optical Flow
Methods in Computer Vision: Introduction to Optical Flow Oren Freifeld Computer Science, Ben-Gurion University March 22 and March 26, 2017 Mar 22, 2017 1 / 81 A Preliminary Discussion Example and Flow
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationAn Adaptive Confidence Measure for Optical Flows Based on Linear Subspace Projections
An Adaptive Confidence Measure for Optical Flows Based on Linear Subspace Projections Claudia Kondermann, Daniel Kondermann, Bernd Jähne, Christoph Garbe Interdisciplinary Center for Scientific Computing
More informationFast Space Charge Calculations with a Multigrid Poisson Solver & Applications
Fast Space Charge Calculations with a Multigrid Poisson Solver & Applications Gisela Pöplau Ursula van Rienen Rostock University DESY, Hamburg, April 26, 2005 Bas van der Geer Marieke de Loos Pulsar Physics
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationIntroduction to Scientific Computing II Multigrid
Introduction to Scientific Computing II Multigrid Miriam Mehl Slide 5: Relaxation Methods Properties convergence depends on method clear, see exercises and 3), frequency of the error remember eigenvectors
More informationA multigrid method for large scale inverse problems
A multigrid method for large scale inverse problems Eldad Haber Dept. of Computer Science, Dept. of Earth and Ocean Science University of British Columbia haber@cs.ubc.ca July 4, 2003 E.Haber: Multigrid
More informationA Distributed-Parameter Approach to Dynamic Image Motion
A Distributed-Parameter Approach to Dynamic Image Motion Annette Stahl, Paul Ruhnau, Christoph Schnörr Computer Vision, Graphics, and Pattern Recognition Group Department of Mathematics and Computer Science;
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationMultigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme Accelerated by Krylov Subspace
201, TextRoad Publication ISSN: 2090-27 Journal of Applied Environmental and Biological Sciences www.textroad.com Multigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationA Multilevel Method for Image Registration
A Multilevel Method for Image Registration Eldad Haber Jan Modersitzki May 1, 004 Abstract In this paper we introduce a new framework for image registration. Our formulation is based on consistent discretization
More informationSecond Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form
Second Order ODEs Second Order ODEs In general second order ODEs contain terms involving y, dy But here only consider equations of the form A d2 y dx 2 + B dy dx + Cy = 0 dx, d2 y dx 2 and F(x). where
More informationOptic Flow Computation with High Accuracy
Cognitive Computer Vision Colloquium Prague, January 12 13, 2004 Optic Flow Computation with High ccuracy Joachim Weickert Saarland University Saarbrücken, Germany joint work with Thomas Brox ndrés Bruhn
More informationPerformance of Optical Flow Techniques. Abstract. While dierent optical ow techniques continue to appear, there has been a
1 Performance of Optical Flow Techniques J.L. Barron 1, D.J. Fleet 2 and S.S. Beauchemin 1 1 Dept. of Computer Science University of Western Ontario London, Ontario, N6A 5B7 2 Dept. of Computing Science
More informationGross Motion Planning
Gross Motion Planning...given a moving object, A, initially in an unoccupied region of freespace, s, a set of stationary objects, B i, at known locations, and a legal goal position, g, find a sequence
More informationSection 4.3 Vector Fields
Section 4.3 Vector Fields DEFINITION: A vector field in R n is a map F : A R n R n that assigns to each point x in its domain A a vector F(x). If n = 2, F is called a vector field in the plane, and if
More informationAggregation-based algebraic multigrid
Aggregation-based algebraic multigrid from theory to fast solvers Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire CEMRACS, Marseille, July 18, 2012 Supported by the Belgian FNRS
More informationANALYSIS AND COMPARISON OF DIFFERENT APPROXIMATIONS TO NONLOCAL DIFFUSION AND LINEAR PERIDYNAMIC EQUATIONS
SIAM J. NUMER. ANAL. Vol. xx, No. x, pp. to appear c 214 Society for Industrial and Applied Mathematics ANALYSIS AND COMPARISON OF DIFFERENT APPROXIMATIONS TO NONLOCAL DIFFUSION AND LINEAR PERIDYNAMIC
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 68 (2014) 1151 1160 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A GPU
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationRevisiting Horn and Schunck: Interpretation as Gauß-Newton Optimisation
ZIKIC, KAMEN, AND NAVAB: REVISITING HORN AND SCHUNCK 1 Revisiting Horn and Schunck: Interpretation as Gauß-Newton Optimisation Darko Zikic 1 zikic@in.tum.de Ali Kamen ali.kamen@siemens.com Nassir Navab
More informationRegular article Nonlinear multigrid methods for total variation image denoising
Comput Visual Sci 7: 199 06 (004) Digital Object Identifier (DOI) 10.1007/s00791-004-0150-3 Computing and Visualization in Science Regular article Nonlinear multigrid methods for total variation image
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More information6. Multigrid & Krylov Methods. June 1, 2010
June 1, 2010 Scientific Computing II, Tobias Weinzierl page 1 of 27 Outline of This Session A recapitulation of iterative schemes Lots of advertisement Multigrid Ingredients Multigrid Analysis Scientific
More informationAN AGGREGATION MULTILEVEL METHOD USING SMOOTH ERROR VECTORS
AN AGGREGATION MULTILEVEL METHOD USING SMOOTH ERROR VECTORS EDMOND CHOW Abstract. Many algebraic multilevel methods for solving linear systems assume that the slowto-converge, or algebraically smooth error
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationGoal-oriented error control of the iterative solution of finite element equations
J. Numer. Math., Vol. 0, No. 0, pp. 1 31 (2009) DOI 10.1515/ JNUM.2009.000 c de Gruyter 2009 Prepared using jnm.sty [Version: 20.02.2007 v1.3] Goal-oriented error control of the iterative solution of finite
More informationLucas/Kanade Meets Horn/Schunck: Combining Local and Global Optic Flow Methods
International Journal of Computer Vision 61(3), 211 231, 2005 c 2005 Springer Science + Business Media, Inc. Manufactured in The Netherlands. Lucas/Kanade Meets Horn/Schunck: Combining Local and Global
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 141 Adaptive Structure Tensors and their Applications Thomas Brox, Rein van den Boomgaard,
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationRobust and Adaptive Multigrid Methods: comparing structured and algebraic approaches
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 0000; 00:1 25 Published online in Wiley InterScience www.interscience.wiley.com). Robust and Adaptive Multigrid Methods: comparing
More informationA Generative Model Based Approach to Motion Segmentation
A Generative Model Based Approach to Motion Segmentation Daniel Cremers 1 and Alan Yuille 2 1 Department of Computer Science University of California at Los Angeles 2 Department of Statistics and Psychology
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationGlobal parametric image alignment via high-order approximation
Global parametric image alignment via high-order approximation Y. Keller, A. Averbuch 2 Electrical & Computer Engineering Department, Ben-Gurion University of the Negev. 2 School of Computer Science, Tel
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationLecture 8: Interest Point Detection. Saad J Bedros
#1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Last Lecture : Edge Detection Preprocessing of image is desired to eliminate or at least minimize noise effects There is always tradeoff
More informationMULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW
MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW VAN EMDEN HENSON CENTER FOR APPLIED SCIENTIFIC COMPUTING LAWRENCE LIVERMORE NATIONAL LABORATORY Abstract Since their early application to elliptic
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 25, No. 6, pp. 982 23 c 24 Society for Industrial and Applied Mathematics A BOUNDARY CONDITION CAPTURING MULTIGRID APPROAC TO IRREGULAR BOUNDARY PROBLEMS JUSTIN W. L. WAN AND
More informationGoal-oriented error control of the iterative solution of finite element equations
J. Numer. Math., Vol. 17, No. 2, pp. 143 172 (2009) DOI 10.1515/ JNUM.2009.009 c de Gruyter 2009 Goal-oriented error control of the iterative solution of finite element equations D. MEIDNER, R. RANNACHER,
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationIn practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points,
1 Interpolation 11 Introduction In practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points, {x i, y i = f(x i ) i = 1 n, obtained,
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 4, No. 3, pp. 6 9 c 004 Society for Industrial and Applied Mathematics CONVERGENCE ANALYSIS OF A MULTIGRID METHOD FOR A CONVECTION-DOMINATED MODEL PROBLEM MAXIM A. OLSHANSKII
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2013 14 CALCULUS AND MULTIVARIABLE CALCULUS MTHA4005Y Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationFinite Element Methods with Linear B-Splines
Finite Element Methods with Linear B-Splines K. Höllig, J. Hörner and M. Pfeil Universität Stuttgart Institut für Mathematische Methoden in den Ingenieurwissenschaften, Numerik und geometrische Modellierung
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationPlane Smoothers for Multiblock Grids: Computational Aspects
NASA/CR-999-20933 ICASE Report No. 99-7 Plane Smoothers for Multiblock Grids: Computational Aspects Ignacio M. Llorente Universidad Complutense, Madrid, Spain Boris Diskin ICASE, Hampton, Virginia N. Duane
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationA fast method for the solution of the Helmholtz equation
A fast method for the solution of the Helmholtz equation Eldad Haber and Scott MacLachlan September 15, 2010 Abstract In this paper, we consider the numerical solution of the Helmholtz equation, arising
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationMotivation: Sparse matrices and numerical PDE's
Lecture 20: Numerical Linear Algebra #4 Iterative methods and Eigenproblems Outline 1) Motivation: beyond LU for Ax=b A little PDE's and sparse matrices A) Temperature Equation B) Poisson Equation 2) Splitting
More informationClassification of Optical Flow by Constraints
Classification of Optical Flow by Constraints Yusuke Kameda 1 and Atsushi Imiya 2 1 School of Science and Technology, Chiba University, Japan 2 Institute of Media and Information Technology, Chiba University,
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationBootstrap AMG. Kailai Xu. July 12, Stanford University
Bootstrap AMG Kailai Xu Stanford University July 12, 2017 AMG Components A general AMG algorithm consists of the following components. A hierarchy of levels. A smoother. A prolongation. A restriction.
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More information3D Space Charge Routines: The Software Package MOEVE and FFT Compared
3D Space Charge Routines: The Software Package MOEVE and FFT Compared Gisela Pöplau DESY, Hamburg, December 4, 2007 Overview Algorithms for 3D space charge calculations Properties of FFT and iterative
More informationADAPTIVE ALGEBRAIC MULTIGRID
ADAPTIVE ALGEBRAIC MULTIGRID M. BREZINA, R. FALGOUT, S. MACLACHLAN, T. MANTEUFFEL, S. MCCORMICK, AND J. RUGE Abstract. Efficient numerical simulation of physical processes is constrained by our ability
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 06-09 Multigrid for High Dimensional Elliptic Partial Differential Equations on Non-equidistant Grids H. bin Zubair, C.W. Oosterlee, R. Wienands ISSN 389-6520 Reports
More informationLehrstuhl für Informatik 10 (Systemsimulation)
FRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG INSTITUT FÜR INFORMATIK (MATHEMATISCHE MASCHINEN UND DATENVERARBEITUNG) Lehrstuhl für Informatik 10 (Systemsimulation) Comparison of two implementations
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationStochastic uncertainty models for the luminance consistency assumption
JOURNAL OF L A TEX CLASS FILES, VOL., NO., MONTH 1 Stochastic uncertainty models for the luminance consistency assumption Thomas Corpetti 1 and Etienne Mémin 2 1 : CNRS/LIAMA, Institute of Automation,
More informationAlgebraic multigrid and multilevel methods A general introduction. Outline. Algebraic methods: field of application
Algebraic multigrid and multilevel methods A general introduction Yvan Notay ynotay@ulbacbe Université Libre de Bruxelles Service de Métrologie Nucléaire May 2, 25, Leuven Supported by the Fonds National
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More informationCreated by T. Madas VECTOR OPERATORS. Created by T. Madas
VECTOR OPERATORS GRADIENT gradϕ ϕ Question 1 A surface S is given by the Cartesian equation x 2 2 + y = 25. a) Draw a sketch of S, and describe it geometrically. b) Determine an equation of the tangent
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More information