Related papers: Symplectic finite-difference methods for solving p…
The Lyapunov exponents of a chaotic system quantify the exponential divergence of initially nearby trajectories. For Hamiltonian systems the exponents are related to the eigenvalues of a symplectic matrix. We make use of this fact to…
The rapid advancements in high-dimensional statistics and machine learning have increased the use of first-order methods. Many of these methods can be regarded as instances of the proximal point algorithm. Given the importance of the…
We develop a hybrid spatial discretization for the wave equation in second order form, based on high-order accurate finite difference methods and discontinuous Galerkin methods. The hybridization combines computational efficiency of finite…
In this paper we propose a new method to stabilise non-symmetric indefinite problems. The idea is to solve a forward and an adjoint problem simultaneously using a suitable stabilised finite element method. Both stabilisation of the element…
The computation time required by standard finite difference methods with fixed timesteps for solving fractional diffusion equations is usually very large because the number of operations required to find the solution scales as the square of…
We present a finite-difference integration algorithm for solution of a system of differential equations containing a diffusion equation with nonlinear terms. The approach is based on Crank-Nicolson method with predictor-corrector algorithm…
Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a…
In this paper, we present numerical procedures to compute solutions of partial differential equations posed on fractals. In particular, we consider the strong form of the equation using standard graph Laplacian matrices and also weak forms…
Symplectic integrators for Hamiltonian systems have been quite successful for studying few-body dynamical systems. These integrators are frequently derived using a formalism built on symplectic maps. There have been recent efforts to extend…
A kind of spatial fractional diffusion equations in this paper are studied. Firstly, an L1 formula is employed for the spatial discretization of the equations. Then, a second order scheme is derived based on the resulting semi-discrete…
This paper develops and analyzes a fully discrete finite element method for a class of semilinear stochastic partial differential equations (SPDEs) with multiplicative noise. The nonlinearity in the diffusion term of the SPDEs is assumed to…
On this paper, we have proposed an approach to observe the time-centered difference scheme for dissipative mechanical systems from a Hamiltonian perspective and to introduce the idea of symplectic algorithm to dissipative systems. The…
Fractional differential equations (FDEs) are an extension of the theory of fractional calculus. However, due to the difficulty in finding analytical solutions, there have not been extensive applications of FDEs until recent decades. With…
Accelerated gradient methods have had significant impact in machine learning -- in particular the theoretical side of machine learning -- due to their ability to achieve oracle lower bounds. But their heuristic construction has hindered…
We present a discrete total variation calculus in Hamiltonian formalism in this paper. Using this discrete variation calculus and generating functions for the flows of Hamiltonian systems, we derive two-step symplectic-energy integrators of…
Symplectic integrators can be excellent for Hamiltonian initial value problems. Reasons for this include their preservation of invariant sets like tori, good energy behaviour, nonexistence of attractors, and good behaviour of statistical…
A novel second order family of explicit stabilized Runge-Kutta-Chebyshev methods for advection-diffusion-reaction equations is introduced. The new methods outperform existing schemes for relatively high Peclet number due to their favorable…
We build convergent discretizations and semi-implicit solvers for the Infinity Laplacian and the game theoretical $p$-Laplacian. The discretizations simplify and generalize earlier ones. We prove convergence of the solution of the Wide…
We develop a fully discrete scheme for time-fractional diffusion equations by using a finite difference method in time and a finite element method in space. The fractional derivatives are used in Caputo sense. Stability and error estimates…
The evolution of any factorized time-reversible symplectic integrators, when applied to the harmonic oscillator, can be exactly solved in a closed form. The resulting modified Hamiltonians demonstrate the convergence of the Lie series…