Related papers: Positivity-preserving methods for population model…
Recently a majorization method for optimizing partition functions of log-linear models was proposed alongside a novel quadratic variational upper-bound. In the batch setting, it outperformed state-of-the-art first- and second-order…
In this paper, we develop parametrized positivity satisfying flux limiters for the high order finite difference Runge-Kutta weighted essentially non-oscillatory (WENO) scheme solving compressible Euler equations to maintain positive density…
We study the second-order boundary value problem \begin{equation*} \begin{cases} \, -u''=a_{\lambda,\mu}(t) \, u^{2}(1-u), & t\in(0,1), \\ \, u'(0)=0, \quad u'(1)=0, \end{cases} \end{equation*} where $a_{\lambda,\mu}$ is a step-wise…
In this paper, we propose new structured second-order methods and structured adaptive-gradient methods obtained by performing natural-gradient descent on structured parameter spaces. Natural-gradient descent is an attractive approach to…
We are interested in developing a numerical method for capturing stationary sheaths, that a plasma forms in contact with a metallic wall. This work is based on a bi-species (ion/electron) Vlasov-Amp{\`e}re model proposed in [3]. The main…
In the present paper, a class of stochastic Runge-Kutta methods containing the second order stochastic Runge-Kutta scheme due to E. Platen for the weak approximation of It\^o stochastic differential equation systems with a multi-dimensional…
Fitting a matrix of a given rank to data in a least squares sense can be done very effectively using 2nd order methods such as Levenberg-Marquardt by explicitly optimizing over a bilinear parameterization of the matrix. In contrast, when…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
In this paper, a class of high-order compact finite difference Hermite scheme is presented for the simulation of double-diffusive convection. To maintain linear stability, the convective fluxes are split into positive and negative parts,…
We propose a monotone, and consistent numerical scheme for the approximation of the Dirichlet problem for the normalized Infinity Laplacian, which could be related to the family of so--called two--scale methods. We show that this method is…
First-order methods have been popularly used for solving large-scale problems. However, many existing works only consider unconstrained problems or those with simple constraint. In this paper, we develop two first-order methods for…
We study the existence of infinitely many positive homoclinic solutions to a second-order difference equation on integers with $p_k$-Laplacian. To achieve our goal we use the critical point theory and the general variational principle of…
Multiphysics systems are driven by multiple processes acting simultaneously, and their simulation leads to partitioned systems of differential equations. This paper studies the solution of partitioned systems of differential equations using…
The exponential fitting technique uses information on the expected behaviour of the solution of a differential problem to define accurate and efficient numerical methods. In particular, exponentially fitted methods are very effective when…
In this paper, a family of arbitrarily high-order structure-preserving exponential Runge-Kutta methods are developed for the nonlinear Schr\"odinger equation by combining the scalar auxiliary variable approach with the exponential…
Runge-Kutta methods have an irreplaceable position among numerical methods designed to solve ordinary differential equations. Especially, implicit ones are suitable for approximating solutions of stiff initial value problems. We propose a…
Concern about how to aggregate sensitive user data without compromising individual privacy is a major barrier to greater availability of data. The model of differential privacy has emerged as an accepted model to release sensitive…
This work proposes a model-reduction methodology that preserves Lagrangian structure (equivalently Hamiltonian structure) and achieves computational efficiency in the presence of high-order nonlinearities and arbitrary parameter dependence.…
This paper presents two kinds of strategies to construct structure-preserving algorithms with homogeneous Neumann boundary conditions for the sine-Gordon equation, while most existing structure-preserving algorithms are only valid for zero…
In this paper, we develop a higher order symmetric partitioned Runge-Kutta method for a coupled system of differential equations on Lie groups. We start with a discussion on partitioned Runge-Kutta methods on Lie groups of arbitrary order.…