The Semi Implicit Gradient Augmented Level Set Method

Here a semi-implicit formulation of the gradient augmented level set method is presented. By tracking both the level set and it's gradient accurate subgrid information is provided,leading to highly accurate descriptions of a moving interface. The result is a hybrid Lagrangian-Eulerian method that may be easily applied in two or three dimensions. The new approach allows for the investigation of interfaces evolving by mean curvature and by the intrinsic Laplacian of the curvature. In this work the algorithm, convergence and accuracy results are presented. Several numerical experiments in both two and three dimensions demonstrate the stability of the scheme.Comment: 19 Pages, 14 Figure

The level set method for the two-sided eigenproblem

We consider the max-plus analogue of the eigenproblem for matrix pencils Ax=lambda Bx. We show that the spectrum of (A,B) (i.e., the set of possible values of lambda), which is a finite union of intervals, can be computed in pseudo-polynomial number of operations, by a (pseudo-polynomial) number of calls to an oracle that computes the value of a mean payoff game. The proof relies on the introduction of a spectral function, which we interpret in terms of the least Chebyshev distance between Ax and lambda Bx. The spectrum is obtained as the zero level set of this function.Comment: 34 pages, 4 figures. Changes with respect to the previous version: we explain relation to mean-payoff games and discrete event systems, and show that the reconstruction of spectrum is pseudopolynomia

A deep level set method for image segmentation

This paper proposes a novel image segmentation approachthat integrates fully convolutional networks (FCNs) with a level setmodel. Compared with a FCN, the integrated method can incorporatesmoothing and prior information to achieve an accurate segmentation.Furthermore, different than using the level set model as a post-processingtool, we integrate it into the training phase to fine-tune the FCN. Thisallows the use of unlabeled data during training in a semi-supervisedsetting. Using two types of medical imaging data (liver CT and left ven-tricle MRI data), we show that the integrated method achieves goodperformance even when little training data is available, outperformingthe FCN or the level set model alone

Stochastic level-set method for shape optimisation

We present a new method for stochastic shape optimisation of engineering structures. The method generalises an existing deterministic scheme, in which the structure is represented and evolved by a level-set method coupled with mathematical programming. The stochastic element of the algorithm is built on the methods of statistical mechanics and is designed so that the system explores a Boltzmann-Gibbs distribution of structures. In non-convex optimisation problems, the deterministic algorithm can get trapped in local optima: the stochastic generalisation enables sampling of multiple local optima, which aids the search for the globally-optimal structure. The method is demonstrated for several simple geometrical problems, and a proof-of-principle calculation is shown for a simple engineering structure.Comment: 17 pages, 10 fig

A level set based method for fixing overhangs in 3D printing

3D printers based on the Fused Decomposition Modeling create objects layer-by-layer dropping fused material. As a consequence, strong overhangs cannot be printed because the new-come material does not find a suitable support over the last deposed layer. In these cases, one can add some support structures (scaffolds) which make the object printable, to be removed at the end. In this paper we propose a level set method to create object-dependent support structures, specifically conceived to reduce both the amount of additional material and the printing time. We also review some open problems about 3D printing which can be of interests for the mathematical community

Level set method for motion by mean curvature

Modeling of a wide class of physical phenomena, such as crystal growth and flame propagation, leads to tracking fronts moving with curvature-dependent speed. When the speed is the curvature this leads to one of the classical degenerate nonlinear second order differential equations on Euclidean space. One naturally wonders "what is the regularity of solutions?" A priori solutions are only defined in a weak sense, but it turns out that they are always twice differentiable classical solutions. This result is optimal; their second derivative is continuous only in very rigid situations that have a simple geometric interpretation. The proof weaves together analysis and geometry. Without deeply understanding the underlying geometry, it is impossible to prove fine analytical properties.Comment: To appear in Notices of the AM