22,367 research outputs found
A High-Order Scheme for Image Segmentation via a modified Level-Set method
In this paper we propose a high-order accurate scheme for image segmentation
based on the level-set method. In this approach, the curve evolution is
described as the 0-level set of a representation function but we modify the
velocity that drives the curve to the boundary of the object in order to obtain
a new velocity with additional properties that are extremely useful to develop
a more stable high-order approximation with a small additional cost. The
approximation scheme proposed here is the first 2D version of an adaptive
"filtered" scheme recently introduced and analyzed by the authors in 1D. This
approach is interesting since the implementation of the filtered scheme is
rather efficient and easy. The scheme combines two building blocks (a monotone
scheme and a high-order scheme) via a filter function and smoothness indicators
that allow to detect the regularity of the approximate solution adapting the
scheme in an automatic way. Some numerical tests on synthetic and real images
confirm the accuracy of the proposed method and the advantages given by the new
velocity.Comment: Accepted version for publication in SIAM Journal on Imaging Sciences,
86 figure
Can local single-pass methods solve any stationary Hamilton-Jacobi-Bellman equation?
The use of local single-pass methods (like, e.g., the Fast Marching method)
has become popular in the solution of some Hamilton-Jacobi equations. The
prototype of these equations is the eikonal equation, for which the methods can
be applied saving CPU time and possibly memory allocation. Then, some natural
questions arise: can local single-pass methods solve any Hamilton-Jacobi
equation? If not, where the limit should be set? This paper tries to answer
these questions. In order to give a complete picture, we present an overview of
some fast methods available in literature and we briefly analyze their main
features. We also introduce some numerical tools and provide several numerical
tests which are intended to exhibit the limitations of the methods. We show
that the construction of a local single-pass method for general Hamilton-Jacobi
equations is very hard, if not impossible. Nevertheless, some special classes
of problems can be actually solved, making local single-pass methods very
useful from the practical point of view.Comment: 19 page
A Method to Study Complex Enzyme Kinetics Involving Numerical Analysis of Enzymatic Schemes. The Mannitol Permease of Escherichia coli as an Example
An analysis of complex kinetic mechanisms is proposed that consists of two steps, (i) building of an kinetic scheme from experimental data other than steady-state kinetics and (ii) numerical simulation and analysis of the kinetics of the proposed scheme in relation to the experimental kinetics. Procedures are introduced to deal with large numbers of enzymatic states and rate constants, and numerical tools are defined to support the analysis of the scheme.
The approach is explored by taking the mannitol permease of Escherichia coli as an example. This enzyme catalyzes both the transport of mannitol across the cytoplasmic membrane and the phosphorylation of mannitol. The challenge is to deduce the transport properties of this dimeric enzyme from the phosphorylation kinetics. It is concluded that (i) the steady-state kinetic behavior is largely consistent with the proposed catalytic cycle of the monomeric subunit, (ii) the kinetics provide no direct support but also do not disprove a coupled translocation of the binding sites on the two monomeric subunits. The approach reveals the need for further experimentation where the implementation of experimental results in the scheme conflict with the experimental kinetics and where specific experimental characteristics do not show up in the simulations of the proposed kinetic scheme.
Stochastic methods for solving high-dimensional partial differential equations
We propose algorithms for solving high-dimensional Partial Differential
Equations (PDEs) that combine a probabilistic interpretation of PDEs, through
Feynman-Kac representation, with sparse interpolation. Monte-Carlo methods and
time-integration schemes are used to estimate pointwise evaluations of the
solution of a PDE. We use a sequential control variates algorithm, where
control variates are constructed based on successive approximations of the
solution of the PDE. Two different algorithms are proposed, combining in
different ways the sequential control variates algorithm and adaptive sparse
interpolation. Numerical examples will illustrate the behavior of these
algorithms
A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators
We develop an algorithm for computing the solution of a large system of
linear ordinary differential equations (ODEs) with polynomial inhomogeneity.
This is equivalent to computing the action of a certain matrix function on the
vector representing the initial condition. The matrix function is a linear
combination of the matrix exponential and other functions related to the
exponential (the so-called phi-functions). Such computations are the major
computational burden in the implementation of exponential integrators, which
can solve general ODEs. Our approach is to compute the action of the matrix
function by constructing a Krylov subspace using Arnoldi or Lanczos iteration
and projecting the function on this subspace. This is combined with
time-stepping to prevent the Krylov subspace from growing too large. The
algorithm is fully adaptive: it varies both the size of the time steps and the
dimension of the Krylov subspace to reach the required accuracy. We implement
this algorithm in the Matlab function phipm and we give instructions on how to
obtain and use this function. Various numerical experiments show that the phipm
function is often significantly more efficient than the state-of-the-art.Comment: 20 pages, 3 colour figures, code available from
http://www.maths.leeds.ac.uk/~jitse/software.html . v2: Various changes to
improve presentation as suggested by the refere
A Parallel Decomposition Scheme for Solving Long-Horizon Optimal Control Problems
We present a temporal decomposition scheme for solving long-horizon optimal
control problems. In the proposed scheme, the time domain is decomposed into a
set of subdomains with partially overlapping regions. Subproblems associated
with the subdomains are solved in parallel to obtain local primal-dual
trajectories that are assembled to obtain the global trajectories. We provide a
sufficient condition that guarantees convergence of the proposed scheme. This
condition states that the effect of perturbations on the boundary conditions
(i.e., initial state and terminal dual/adjoint variable) should decay
asymptotically as one moves away from the boundaries. This condition also
reveals that the scheme converges if the size of the overlap is sufficiently
large and that the convergence rate improves with the size of the overlap. We
prove that linear quadratic problems satisfy the asymptotic decay condition,
and we discuss numerical strategies to determine if the condition holds in more
general cases. We draw upon a non-convex optimal control problem to illustrate
the performance of the proposed scheme
- …