The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with carefully chosen implicit-explicit (IMEX) Runge-Kutta time discretization up to third order accuracy, for solving one-dimensional linear advection-diffusion equations. In the time discretization the advection term is treated explicitly and the diffusion term implicitly. There are three highlights of this work. The first is that we establish an important relationship between the gradient and interface jump of the numerical solution with the independent numerical solution of the gradient in the LDG methods. The second is that, by aid of the aforementioned relationship and the energy method, we show that the IMEX LDG schemes are unconditionally stable for the linear problems, in the sense that the time-step $\dt$ is only required to be upper-bounded by a constant which depends on the ratio of the diffusion and the square of the advection coefficients and is independent of the mesh-size $h$, even though the advection term is treated explicitly. The last is that under this time step condition, we obtain optimal error estimates in both space and time for the third order IMEX Runge-Kutta time-marching coupled with LDG spatial discretization. Numerical experiments are also given to verify the main results.

A new Tikhonov regularization method of Fuhry and Reichel [A new Tikhonov regularization method, Numerical Algorithms, 59:433-445, 2011] exhibits the excellent properties for ill-posed problems, but it can only deal with small or moderate size problems because of the expensive computation of singular value decomposition (SVD). In this paper, we extend the above new Tikhonov regularization method to solve large-scale problems, e.g., image restoration problem with periodic boundary conditions, and realize this extending by applying Fast Fourier Transformation (FFT) algorithm to the spectral decomposition of the block circulant with circulant blocks (BCCB) matrices. Experimental results confirm the superiority of our new method.

We propose a computational efficient yet simple numerical algorithm to solve the surface eikonal equation on general implicit surfaces. The method is developed based on the embedding idea and the fast sweeping methods. We first approximate the solution to the surface eikonal equation by the Euclidean weighted distance function defined in a tubular neighbourhood of the implicit surface, and then apply the efficient fast sweeping method to numerically compute the corresponding viscosity solution. Unlike some other embedding methods which require the radius of the computational tube satisfies h = O ( x ) for some h = O ( x ) , our approach allows h = O ( x ) . This implies that the total number of grid points in the computational tube is optimal and is given by h = O ( x ) for a co-dimensional one surface in h = O ( x ) . The method can be easily extended to general static HamiltonJacobi equation defined on implicit

No abstract uploaded!

In this article, we propose two numerical methods, the Gaussian Process (GP) method and the Fourier Features (FF) algorithm, to solve mean field games (MFGs). The GP algorithm approximates the solution of a MFG with maximum a posteriori probability estimators of GPs conditioned on the partial differential equation (PDE) system of the MFG at a finite number of sample points. The main bottleneck of the GP method is to compute the inverse of a square gram matrix, whose size is proportional to the number of sample points. To improve the performance, we introduce the FF method, whose insight comes from the recent trend of approximating positive definite kernels with random Fourier features. The FF algorithm seeks approximated solutions in the space generated by sampled Fourier features. In the FF method, the size of the matrix to be inverted depends only on the number of Fourier features selected, which is much less than the size of sample points. Hence, the FF method reduces the precomputation time, saves the memory, and achieves comparable accuracy to the GP method. We give the existence and the convergence proofs for both algorithms. The convergence argument of the GP method does not depend on any monotonicity condition, which suggests the potential applications of the GP method to solve MFGs with non-monotone couplings in future work. We show the efficacy of our algorithms through experiments on a stationary MFG with a non-local coupling and on a time-dependent planning problem. We believe that the FF method can also serve as an alternative algorithm to solve general PDEs.