The Flexible, Extensible and Efficient Toolbox of Level Set Methods

Abstract Level set methods are a popular and powerful class of numerical algorithms for dynamic implicit surfaces and solution of Hamilton-Jacobi PDEs. While the advanced level set schemes combine both efficiency and accuracy, their implementation complexity makes it difficult for the community to reproduce new results and make quantitative comparisons between methods. This paper describes the Toolbox of Level Set Methods, a collection of Matlab routines implementing the basic level set algorithms on fixed Cartesian grids for rectangular domains in arbitrary dimension. The Toolbox’s code and interface are designed to permit flexible combinations of different schemes and PDE forms, allow easy extension through the addition of new algorithms, and achieve efficient execution despite the fact that the code is entirely written as m-files. The current contents of the Toolbox and some coding patterns important to achieving its flexibility, extensibility and efficiency are briefly explained, as is the process of adding two new algorithms. Code for both the Toolbox and the new algorithms is available from the Web.

[1]  Alexandre M. Bayen,et al.  A time-dependent Hamilton-Jacobi formulation of reachable sets for continuous dynamic games , 2005, IEEE Transactions on Automatic Control.

[2]  D. Talay.  Numerical solution of stochastic differential equations , 1994 .

[3]  G. Barles,et al.  Convergence of approximation schemes for fully nonlinear second order equations , 1991 .

[4]  Adam M. Oberman.  A convergent monotone difference scheme for motion of level sets by mean curvature , 2004, Numerische Mathematik.

[5]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[6]  Bernt Øksendal,et al.  Stochastic differential equations (3rd ed.): an introduction with applications , 1992 .

[7]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[8]  Ian M. Mitchell,et al.  A Toolbox of Hamilton-Jacobi Solvers for Analysis of Nondeterministic Continuous and Hybrid Systems , 2005, HSCC.

[9]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[10]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[11]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[12]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[13]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[14]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[15]  P. Saint-Pierre,et al.  Set-Valued Numerical Analysis for Optimal Control and Differential Games , 1999 .

[16]  Ryo Takei.  MODERN THEORY OF NUMERICAL METHODS FOR MOTION BY MEAN CURVATUR,E , 2007 .

[17]  Wang Hai-bing,et al.  High-order essentially non-oscillatory schemes for Hamilton-Jacobi equations , 2006 .

[18]  J. Sethian.  Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[19]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[20]  P. Lions,et al.  User’s guide to viscosity solutions of second order partial differential equations , 1992, math/9207212.

[21]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[22]  S. Osher,et al.  Geometric Level Set Methods in Imaging, Vision, and Graphics , 2011, Springer New York.

[23]  Alexander B. Kurzhanski,et al.  Control Synthesis for State Constrained Systems and Obstacle Problems , 2004 .

[24]  P. Smereka,et al.  A Remark on Computing Distance Functions , 2000 .

[25]  B. Øksendal.  Stochastic differential equations : an introduction with applications , 1987 .

[26]  P. Lions,et al.  Viscosity solutions of Hamilton-Jacobi equations , 1983 .

[27]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[28]  S. Osher.  A level set formulation for the solution of the Dirichlet problem for Hamilton-Jacobi equations , 1993 .

[29]  Ian M. Mitchell,et al.  A Toolbox of Level Set Methods , 2005 .

[30]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[31]  Adam M. Oberman,et al.  Convergent Difference Schemes for Degenerate Elliptic and Parabolic Equations: Hamilton-Jacobi Equations and Free Boundary Problems , 2006, SIAM J. Numer. Anal..

[32]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[33]  P. Smereka.  Spiral crystal growth , 2000 .

[34]  S. Osher,et al.  Lax-Friedrichs sweeping scheme for static Hamilton-Jacobi equations , 2004 .

[35]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[36]  Ian M. Mitchell,et al.  Level Set Methods for Computation in Hybrid Systems , 2000, HSCC.

[37]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[38]  P. Lions,et al.  Two approximations of solutions of Hamilton-Jacobi equations , 1984 .

[39]  S. Osher,et al.  An improved level set method for incompressible two-phase flows , 1998 .

[40]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .