Related papers: A high-order exponential integrator for nonlinear …
For a linear equality constrained convex optimization problem involving two objective functions with a ``nonsmooth" + ``nonsmooth" composite structure, we study two algorithms derived from a mixed-order dynamical system which incorporates…
We present a derivative-based algorithm for nonlinearly constrained optimization problems that is tolerant of inaccuracies in the data. The algorithm solves a semi-smooth set of nonlinear equations that are equivalent to the first-order…
In this paper, we develop a numerical scheme for the space-time fractional parabolic equation, i.e., an equation involving a fractional time derivative and a fractional spatial operator. Both the initial value problem and the…
Successive quadratic approximations, or second-order proximal methods, are useful for minimizing functions that are a sum of a smooth part and a convex, possibly nonsmooth part that promotes regularization. Most analyses of iteration…
In this paper, we present new types of exponential integrators for Stochastic Differential Equations (SDEs) that take the advantage of the exact solution of (generalised) geometric Brownian motion. We examine both Euler and Milstein…
In this paper a variant of nonlinear exponential Euler scheme is proposed for solving nonlinear heat conduction problems. The method is based on nonlinear iterations where at each iteration a linear initial-value problem has to be solved.…
A sequential quadratic optimization algorithm is proposed for solving smooth nonlinear equality constrained optimization problems in which the objective function is defined by an expectation of a stochastic function. The algorithmic…
In the framework of a Particle-In-Cell scheme for some 1D Vlasov-Poisson system depending on a small parameter, we propose a time-stepping method which is numerically uniformly accurate when the parameter goes to zero. Based on an…
We will discuss an extension of the pseudospectral method developed by Wineberg, McGrath, Gabl, and Scott for the numerical integration of the KdV initial value problem. Our generalization of their algorithm can be used to solve initial…
It is classical that, when the small deformation is assumed, the incremental analysis problem of an elastoplastic structure with a piecewise-linear yield condition and a linear strain hardening model can be formulated as a convex quadratic…
In this study, we propose high-order implicit and semi-implicit schemes for solving ordinary differential equations (ODEs) based on Taylor series expansion. These methods are designed to handle stiff and non-stiff components within a…
The strong numerical approximation of semilinear stochastic partial differential equations (SPDEs) driven by infinite dimensional Wiener processes is investigated. There are a number of results in the literature that show that Euler-type…
By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for…
In this work, we propose an adaptive spectral element algorithm for solving nonlinear optimal control problems. The method employs orthogonal collocation at the shifted Gegenbauer-Gauss points combined with very accurate and stable…
We present a new and relatively elementary method for studying the solution of the initial-value problem for dispersive linear and integrable equations in the large-$t$ limit, based on a generalization of steepest descent techniques for…
We introduce low regularity exponential-type integrators for nonlinear Schr\"odinger equations for which first-order convergence only requires the boundedness of one additional derivative of the solution. More precisely, we will prove…
Exponential integrators are a well-known class of time integration methods that have been the subject of many studies and developments in the past two decades. Surprisingly, there have been limited efforts to analyze their stability and…
A cutting-plane model for a nonsmooth function is the maximum of several first-order expansions centered at different points. Using such a model in a bundle method leads to linear convergence (of serious steps) to a minimum. In smooth…
In this work, we present a first-order unfiltered exponential integrator for the one-dimensional derivative nonlinear Schr\"odinger equation with low regularity. Our analysis shows that for any $s>\frac12$, the method converges with…
In this paper, we propose a $\mu$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a $d$-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix…