Source author record

Xian-Ming Gu

Xian-Ming Gu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

7works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Fast second-order implicit difference schemes for time distributed-order and Riesz space fractional diffusion-wave equations

In this paper, fast numerical methods are established for solving a class of time distributed-order and Riesz space fractional diffusion-wave equations. We derive new difference schemes by the weighted and shifted Gr$\ddot{\rm{u}}$nwald formula in time and the fractional centered difference formula in space. The unconditional stability and second-order convergence in time, space and distributed-order of the difference schemes are analyzed. In the one-dimensional case, the Gohberg-Semencul formula utilizing the preconditioned Krylov subspace method is developed to solve the symmetric positive definite Toeplitz linear systems derived from the proposed difference scheme. In the two-dimensional case, we also design a global preconditioned conjugate gradient method with a truncated preconditioner to solve the discretized Sylvester matrix equations. We prove that the spectrums of the preconditioned matrices in both cases are clustered around one, such that the proposed numerical methods with preconditioners converge very quickly. Some numerical experiments are carried out to demonstrate the effectiveness of the proposed difference schemes and show that the performances of the proposed fast solution algorithms are better than other numerical methods.

preprint2020arXiv

Two CSCS-based iteration methods for solving absolute value equations

Recently, two families of HSS-based iteration methods are constructed for solving the system of absolute value equations (AVEs), which is a class of non-differentiable NP-hard problems. In this study, we establish the Picard-CSCS iteration method and the nonlinear CSCS-like iteration method for AVEs involving the Toeplitz matrix. Then, we analyze the convergence of the Picard-CSCS iteration method for solving AVEs. By using the theory about nonsmooth analysis, we particularly prove the convergence of the nonlinear CSCS-like iterationsolver for AVEs. The advantage of these methods is that they do not require the storage of coefficient matrices at all, and the sub-system of linear equations can be solved efficiently via the fast Fourier transforms (FFTs). Therefore, computational cost and storage can be saved in practical implementations. Numerical examples including numerical solutions of nonlinear fractional diffusion equations are reported to show the effectiveness of the proposed methods in comparison with some existing methods.

preprint2019arXiv

An implicit integration factor method for a kind of spatial fractional diffusion equations

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 ordinary differential system by using the implicit integration factor method, which is a class of efficient semi-implicit temporal scheme. Numerical results show that the proposed scheme is accurate even for the discontinuous coefficients.

preprint2018arXiv

A second-order accurate implicit difference scheme for time fractional reaction-diffusion equation with variable coefficients and time drift term

An implicit finite difference scheme based on the $L2$-$1_σ$ formula is presented for a class of one-dimensional time fractional reaction-diffusion equations with variable coefficients and time drift term. The unconditional stability and convergence of this scheme are proved rigorously by the discrete energy method, and the optimal convergence order in the $L_2$-norm is $\mathcal{O}(τ^2 + h^2)$ with time step $τ$ and mesh size $h$. Then, the same measure is exploited to solve the two-dimensional case of this problem and a rigorous theoretical analysis of the stability and convergence is carried out. Several numerical simulations are provided to show the efficiency and accuracy of our proposed schemes and in the last numerical experiment of this work, three preconditioned iterative methods are employed for solving the linear system of the two-dimensional case.

preprint2016arXiv

Block variants of the COCG and COCR methods for solving complex symmetric linear systems with multiple right-hand sides

In the present study, we establish two new block variants of the Conjugate Orthogonal Conjugate Gradient (COCG) and the Conjugate A-Orthogonal Conjugate Residual (COCR) Krylov subspace methods for solving complex symmetric linear systems with multiple right hand sides. The proposed Block iterative solvers can fully exploit the complex symmetry property of coefficient matrix of the linear system. We report on extensive numerical experiments to show the favourable convergence properties of our newly developed Block algorithms for solving realistic electromagnetic simulations.

preprint2015arXiv

Numerical Gradient Schemes for Heat Equations Based on the Collocation Polynomial and Hermite Interpolation

As is well-known, the advantage of the high-order compact difference scheme (H-OCD) is unconditionally stable and convergent with the order $O(τ^2+h^4)$ under the maximum norm. In this article, a new numerical gradient scheme based on the collocation polynomial and Hermite interpolation is presented. Moreover, the convergence order of this kind of method is also $O(τ^2+h^4)$ under the discrete maximum norm when the space step size is just twice the one of H-OCD method, which accelerates the computational process and makes the result much smoother to some extent. In addition, some corresponding analyses are made and the Richardson extrapolation technique is also considered in time direction. The results of numerical experiments are also consistent with these theoretical analysis.

preprint2013arXiv

Strang-type preconditioners for solving fractional diffusion equations by boundary value methods

The finite difference scheme with the shifted Grünwarld formula is employed to semi-discrete the fractional diffusion equations. This spatial discretization can reduce to the large system of ordinary differential equations (ODEs) with initial values. Recently, boundary value method (BVM) was developed as a popular algorithm for solving large systems of ODEs. This method requires the solutions of one or more nonsymmetric, large and sparse linear systems. In this paper, the GMRES method with the block circulant preconditioner is proposed for solving these linear systems. One of the main results is that if an $A_{ν_1,ν_2}$-stable boundary value method is used for an m-by-m system of ODEs, then the preconditioner is invertible and the preconditioned matrix can be decomposed as I+L, where I is the identity matrix and the rank of L is at most $2m(ν_1+ν_2)$. It means that when the GMRES method is applied to solve the preconditioned linear systems, the method will converge in at most $2m(ν_1+ν_2)+1$ iterations.Finally, extensive numerical experiments are reported to illustrate the effectiveness of our methods for solving the fractional diffusion equations.