Related papers: Lamperti Semi-Discrete method
We study numerical methods for solving a system of quasilinear stochastic partial differential equations known as the stochastic Landau-Lifshitz-Bloch (LLB) equation on a bounded domain in $\mathbb R^d$ for $d=1,2$. Our main results are…
This article is concerned with the multilevel Monte Carlo (MLMC) methods for approximating expectations of some functions of the solution to the Heston 3/2-model from mathematical finance, which takes values in $(0, \infty)$ and possesses…
Discrete gradient methods are a powerful tool for the time discretization of dynamical systems, since they are structure-preserving regardless of the form of the total energy. In this work, we discuss the application of discrete gradient…
We review our recent work leading to steady-state solutions of the semiclassical (Maxwell-Bloch) equations of a laser. These are coupled non-linear partial differential equations in space and time which have previously been solved either by…
This paper provides insight into the estimation and asymptotic behavior of parameters in interest rate models, focusing primarily on the Cox-Ingersoll-Ross (CIR) process and its extension -- the more general Chan-Karolyi-Longstaff-Sanders…
We introduce a new class of numerical methods for solving McKean-Vlasov stochastic differential equations, which are relevant in the context of distribution-dependent or mean-field models, under super-linear growth conditions for both the…
Semi-supervised regression (SSR), which aims to predict continuous scores for samples while reducing the reliance on large-scale labeled data, has recently attracted considerable attention across various applications, including computer…
The nonlinear Schr\"{o}dinger (NLS) equation possesses an infinite hierarchy of conserved densities and the numerical preservation of some of these quantities is critical for accurate long-time simulations, particularly for multi-soliton…
We study the asymptotic behavior of solution of semi-linear PDEs. Neither periodicity nor ergodicity will be assumed. In return, we assume that the coefficients admit a limit in \`{C}esaro sense. In such a case, the averaged coefficients…
In this paper, we study the semiclassical Schr\"odinger equation with random parameters and develop several robust multi-fidelity methods. We employ the time-splitting Fourier pseudospectral (TSFP) method for the high-fidelity solver, and…
Anomalous diffusions are ubiquitous in nature, whose functional distributions are governed by the backward Feynman-Kac equation. In this paper, the local discontinuous Galerkin (LDG) method is used to solve the 2D backward Feynman-Kac…
Density-adaptive domain discretization is essential for high-utility privacy-preserving analytics but remains challenging under Local Differential Privacy (LDP) due to the privacy-budget costs associated with iterative refinement. We…
We present the Latent Sequence Decompositions (LSD) framework. LSD decomposes sequences with variable lengthed output units as a function of both the input sequence and the output sequence. We present a training algorithm which samples…
In this paper, we propose an RADI-type method for large-scale stochastic continuous-time algebraic Riccati equations with sparse and low-rank matrices. This new variant of RADI-type methods is developed by integrating the core concept of…
In this paper, we present a sparse grid-based Monte Carlo method for solving high-dimensional semi-linear nonlocal diffusion equations with volume constraints. The nonlocal model is governed by a class of semi-linear partial…
This paper deals with the error processing problem of sparse identification of nonlinear dynamical systems(SINDy) through introducing the $L_\infty$ approximation to take place of the former $L_2$ approximation. The motivation is that the…
Many stochastic differential equations (SDEs) in the literature have a superlinearly growing nonlinearity in their drift or diffusion coefficient. Unfortunately, moments of the computationally efficient Euler-Maruyama approximation method…
Priors with non-smooth log-densities, such as the l1-prior, are widely used in Bayesian inverse problems for their sparsity-inducing properties. Existing Langevin-based sampling methods typically rely on proximal mappings or smooth…
We consider the one dimensional focusing (cubic) Nonlinear Schr\"odinger equation (NLS) in the semiclassical limit with exponentially decaying complex-valued initial data, whose phase is multiplied by a real parameter. We prove smooth…
A system of semi-discrete coupled nonlinear Schr\"{o}dinger equations is studied. To show the complete integrability of the model with multiple components, we extend the discrete version of the inverse scattering method for the…