The restarted FOM method presented by Simoncini[7]according to the natural collinearity of all residuals is an efficient method for solving shifted systems,which generates the same Krylov subspace when the shifts are ...The restarted FOM method presented by Simoncini[7]according to the natural collinearity of all residuals is an efficient method for solving shifted systems,which generates the same Krylov subspace when the shifts are handled simultaneously.However,restarting slows down the convergence.We present a practical method for solving the shifted systems by adding some Ritz vectors into the Krylov subspace to form an augmented Krylov subspace. Numerical experiments illustrate that the augmented FOM approach(restarted version)can converge more quickly than the restarted FOM method.展开更多
The seed method is used for solving multiple linear systems A (i)x (i) =b (i) for 1≤i≤s, where the coefficient matrix A (i) and the right-hand side b (i) are different in general. It is known that the CG meth...The seed method is used for solving multiple linear systems A (i)x (i) =b (i) for 1≤i≤s, where the coefficient matrix A (i) and the right-hand side b (i) are different in general. It is known that the CG method is an effective method for symmetric coefficient matrices A (i). In this paper, the FOM method is employed to solve multiple linear sy stems when coefficient matrices are non-symmetric matrices. One of the systems is selected as the seed system which generates a Krylov subspace, then the resi duals of other systems are projected onto the generated Krylov subspace to get t he approximate solutions for the unsolved ones. The whole process is repeated u ntil all the systems are solved.展开更多
In this paper, we propose another version of the full orthogonalization method (FOM) for the resolution of linear system Ax = b, based on an extended definition of Sturm sequence in the calculation of the determinant ...In this paper, we propose another version of the full orthogonalization method (FOM) for the resolution of linear system Ax = b, based on an extended definition of Sturm sequence in the calculation of the determinant of an upper hessenberg matrix in o(n2). We will also give a new version of Givens method based on using a tensor product and matrix addition. This version can be used in parallel calculation.展开更多
基金Supported-by the National Natural Science Foundation of China(19971057)the Science and Technology Developing Foundation of University in Shanghai of China(02AK41).
文摘The restarted FOM method presented by Simoncini[7]according to the natural collinearity of all residuals is an efficient method for solving shifted systems,which generates the same Krylov subspace when the shifts are handled simultaneously.However,restarting slows down the convergence.We present a practical method for solving the shifted systems by adding some Ritz vectors into the Krylov subspace to form an augmented Krylov subspace. Numerical experiments illustrate that the augmented FOM approach(restarted version)can converge more quickly than the restarted FOM method.
基金Project supported by the National Natural Science Foundation of China (Grant No.10271075)
文摘The seed method is used for solving multiple linear systems A (i)x (i) =b (i) for 1≤i≤s, where the coefficient matrix A (i) and the right-hand side b (i) are different in general. It is known that the CG method is an effective method for symmetric coefficient matrices A (i). In this paper, the FOM method is employed to solve multiple linear sy stems when coefficient matrices are non-symmetric matrices. One of the systems is selected as the seed system which generates a Krylov subspace, then the resi duals of other systems are projected onto the generated Krylov subspace to get t he approximate solutions for the unsolved ones. The whole process is repeated u ntil all the systems are solved.
文摘In this paper, we propose another version of the full orthogonalization method (FOM) for the resolution of linear system Ax = b, based on an extended definition of Sturm sequence in the calculation of the determinant of an upper hessenberg matrix in o(n2). We will also give a new version of Givens method based on using a tensor product and matrix addition. This version can be used in parallel calculation.
