Davod Khojasteh Salkuyeh

Mohaddese Kaveh Shaldehi, Davod Khojasteh Salkuyeh

This paper introduces inexact versions of several block-splitting preconditioners for solving the three-by-three block linear systems arising from a special class of indefinite least squares problems. We first establish the convergence conditions for the corresponding stationary iterative methods. Then, it follows that under these conditions, all eigenvalues of the preconditioned matrices are contained within a circle centered at $(1,0)$ with radius $1$. This property implies that these preconditioners are effective in accelerating the convergence of the GMRES method. Furthermore, we analyze the eigenpairs of the preconditioned matrices in detail and derive a theoretical upper bound on the number of GMRES iterations for solving the preconditioned systems. Ultimately, numerical experiments reveal the efficacy of the proposed preconditioners.

Davod Khojasteh Salkuyeh, Mohsen Masoudi, Davod Hezari

In this paper, the generalized shift-splitting preconditioner is implemented for saddle point problems with symmetric positive definite (1,1)-block and symmetric positive semidefinite (2,2)-block. The proposed preconditioner is extracted form a stationary iterative method which is unconditionally convergent. Moreover, a relaxed version of the proposed preconditioner is presented and some properties of the eigenvalues distribution of the corresponding preconditioned matrix are studied. Finally, some numerical experiments on test problems arisen from finite element discretization of the Stokes problem are given to show the effectiveness of the preconditioners.

Davod Khojasteh Salkuyeh, Tahereh Salimi Siahkolaei

We introduce a two-parameter version of the two-step scale-splitting iteration method, called TTSCSP, for solving a broad class of complex symmetric system of linear equations. We present some conditions for the convergence of the method. An upper bound for the spectral radius of the method is presented and optimal parameters which minimize this bound are given. Inexact version of the TTSCSP iteration method (ITTSCSP) is also presented. Some numerical experiments are reported to verify the effectiveness of the TTSCSP iteration method and the numerical results are compared with those of the TSCSP, the SCSP and the PMHSS iteration methods. Numerical comparison of the ITTSCSP method with the inexact version of TSCSP, SCSP and PMHSS are presented. We also compare the numerical results of the BiCGSTAB method in conjunction with the TTSCSP and the ILU preconditioners.

Davod Khojasteh Salkuyeh, Mohsen Masoudi

We present a preconditioner for saddle point problems. The proposed preconditioner is extracted from a stationary iterative method which is convergent under a mild condition. Some properties of the preconditioner as well as the eigenvalues distribution of the preconditioned matrix are presented. The preconditioned system is solved by a Krylov subspace method like restarted GMRES. Finally, some numerical experiments on test problems arisen from finite element discretization of the Stokes problem are given to show the effectiveness of the preconditioner.

Davod Khojasteh Salkuyeh

Based on the Scale-Splitting (SCSP) iteration method presented by Hezari et al. in (A new iterative method for solving a class of complex symmetric system linear of equations, Numerical Algorithms 73 (2016) 927-955), we present a new two-step iteration method, called TSCSP, for solving the complex symmetric system of linear equations $(W+iT)x=b$, where $W$ and $T$ are symmetric positive definite and symmetric positive semidefinite matrices, respectively. It is shown that if the matrices $W$ and $T$ are symmetric positive definite, then the method is unconditionally convergent. The optimal value of the parameter, which minimizes the spectral radius of the iteration matrix is also computed. Numerical {comparisons} of the TSCSP iteration method with the SCSP, the MHSS, the PMHSS and the GSOR methods are given to illustrate the effectiveness of the method.

Davod Khojasteh Salkuyeh, Ali Tavakoli

In order to solve an initial value problem by the variational iteration method, a sequence of functions is produced which converges to the solution under some suitable conditions. In the nonlinear case, after a few iterations the terms of the sequence become complicated, and therefore, computing a highly accurate solution would be difficult or even impossible. In this paper, for one-dimensional initial value problems, we propose a new approach which is based on approximating each term of the sequence by a piecewise linear function. Moreover, the convergence of the method is proved. Three illustrative examples are given to show the superiority of the proposed method over the classical variational iteration method.

Davod Khojasteh Salkuyeh, Maryam Rahimian

A modification of the generalized shift-splitting (GSS) method is presented for solving singular saddle point problems. In this kind of modification, the diagonal shift matrix is replaced by a block diagonal matrix which is symmetric positive definite. Semi-convergence of the proposed method is investigated. The induced preconditioner is applied to the saddle point problem and the preconditioned system is solved by the restarted generalized minimal residual method. Eigenvalue distribution of the preconditioned matrix is also discussed. Finally some numerical experiments are given to show the effectiveness and robustness of the new preconditioner. Numerical results show that the modified GSS method is superior to the classical GSS method.

Davod Khojasteh Salkuyeh, Mohsen Hasani, Fatemeh Panjeh Ali Beik

Several preconditioned AOR methods have been proposed to solve system of linear equations $Ax=b$, where $A \in \mathbb{R}^{n \times n}$ is a unit Z-matrix. The aim of this paper is to give a comparison result for a class of preconditioners $P$, where $P\in \mathbb{R}^{n\times n}$ is nonsingular, nonnegative and has unit diagonal entries. Numerical results for corresponding preconditioned GMRES methods are given to illustrate the theoretical results.

Davod Khojasteh Salkuyeh, Amin Rafiei, Hadi Roohani

A technique for computing an ILU preconditioner based on the FAPINV algorithm is presented. We show that this algorithm is well-defined for H-matrices. Moreover, when used in conjunction with Krylov-subspace-based iterative solvers such as the GMRES algorithm, results in reliable solvers. Numerical experiments on some test matrices are given to show the efficiency of the new ILU preconditioner.

Davod Khojasteh Salkuyeh, Faezeh Toutounian

In this paper, a method via sparse-sparse iteration for computing a sparse incomplete factorization of the inverse of a symmetric positive definite matrix is proposed. The resulting factorized sparse approximate inverse is used as a preconditioner for solving symmetric positive definite linear systems of equations by using the preconditioned conjugate gradient algorithm. Some numerical experiments on test matrices from the Harwell-Boeing collection for comparing the numerical performance of the presented method with one available well-known algorithm are also given.

Davod Khojasteh Salkuyeh, Tahereh Salimi Siahkalaei

We present a block lower triangular (BLT) preconditioner to accelerate the convergence of nthe Krylov subspace iterative methods, such as generalized minimal residual (GMRES), for solving a broad class of complex symmetric system of linear equations. We analyze the eigenvalues distribution of preconditioned coefficient matrix. Numerical experiments are given to demonstrate the effectiveness of the BLT preconditioner.

Davod Khojasteh Salkuyeh, Mohsen Masoudi, Davod Hezari

In this paper, we propose a preconditioner based on the shift-splitting method for generalized saddle point problems with nonsymmetric positive definite (1,1)-block and symmetric positive semidefinite $(2,2)$-block. The proposed preconditioner is obtained from an basic iterative method which is unconditionally convergent. We also present a relaxed version of the proposed method. Some numerical experiments are presented to show the effectiveness of the method.

Davod Hezari, Vahid Edalatpour, Davod Khojasteh Salkuyeh

In this paper, we present a preconditioned variant of the generalized successive overrelaxation (GSOR) iterative method for solving a broad class of complex symmetric linear systems. We study conditions under which the spectral radius of the iteration matrix of the preconditioned GSOR method is smaller than that of the GSOR method and determine the optimal values of iteration parameters. Numerical experiments are given to verify the validity of the presented theoretical results and the effectiveness of the preconditioned GSOR method.

Davod Khojasteh Salkuyeh

In [7], a new iterative method for solving linear system of equations was presented which can be considered as a modification of the Gauss-Seidel method. Then in [4] a different approach, say 2D-DSPM, and more effective one was introduced. In this paper, we improve this method and give a generalization of it. Convergence properties of this kind of generalization are also discussed. We finally give some numerical experiments to show the efficiency of the method and compare with 2D-DSPM.