In this paper, equivalence between the Mann and Ishikawa iterations for a generalized contraction mapping in cone subset of a real Banach space is discussed.
This note explores the relations between two different methods. The first one is the Alternating Least Squares (ALS) method for calculating a rank<em>-k</em> approximation of a real <em>m</em>&...This note explores the relations between two different methods. The first one is the Alternating Least Squares (ALS) method for calculating a rank<em>-k</em> approximation of a real <em>m</em>×<em>n</em> matrix, <em>A</em>. This method has important applications in nonnegative matrix factorizations, in matrix completion problems, and in tensor approximations. The second method is called Orthogonal Iterations. Other names of this method are Subspace Iterations, Simultaneous Iterations, and block-Power method. Given a real symmetric matrix, <em>G</em>, this method computes<em> k</em> dominant eigenvectors of <em>G</em>. To see the relation between these methods we assume that <em>G </em>=<em> A</em><sup>T</sup> <em>A</em>. It is shown that in this case the two methods generate the same sequence of subspaces, and the same sequence of low-rank approximations. This equivalence provides new insight into the convergence properties of both methods.展开更多
A new restoration algorithm based on double loops and alternant iterations is proposed to restore the object image effectively from a few frames of turbulence-degraded images, Based on the double loops, the iterative ...A new restoration algorithm based on double loops and alternant iterations is proposed to restore the object image effectively from a few frames of turbulence-degraded images, Based on the double loops, the iterative relations for estimating the turbulent point spread function PSF and object image alternately are derived. The restoration experiments have been made on computers, showing that the proposed algorithm can obtain the optimal estimations of the object and the point spread function, with the feasibility and practicality of the proposed algorithm being convincing.展开更多
In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduce...In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduced that enables successive decomposition of high Eigen values of the iteration matrix to enable convergence. While extrapolation of a convergent fixed point iteration using a geometric series sum is a known form of Aitken acceleration, it is shown that in this paper, the same formula can be used to estimate the solution of sets of linear equations from diverging Gauss-Seidel iterations. In both convergent and divergent iterations, the ratios of differences among the consecutive values of iteration eventually form a convergent (divergent) series with a factor equal to the largest Eigen value of the iteration matrix. Higher order Aitken extrapolation is shown to eliminate the influence of dominant Eigen values of the iteration matrix in successive order until the iteration is determined by the lowest possible Eigen values. For the convergent part of the Gauss-Seidel iteration, further acceleration is made possible by coupling of the extrapolation technique with the successive over relaxation (SOR) method. Application examples from both convergent and divergent iterations have been provided. Coupling of the extrapolation with the SOR technique is also illustrated for a steady state two dimensional heat flow problem which was solved using MATLAB programming.展开更多
The simultaneous iterations rithms of the ART family. It is used reconstruction technique (SIRT) widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matr...The simultaneous iterations rithms of the ART family. It is used reconstruction technique (SIRT) widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matrices. Its theoretical background and iteration model are discussed at the beginning of this paper. Then, the implementation of the SIRT to reconstruct the three-dimensional distribution of water vapor by simulation is discussed. The results show that the SIRT can function effectively in water vapor tomography, obtain rapid convergence, and be implemented more easily than inversion.展开更多
Under the weak Lipschitz condition about the solution of the equation, convergence theorems for a family of iterations with one parameter are obtained. An estimation of the radius of the attraction ball is shown. At l...Under the weak Lipschitz condition about the solution of the equation, convergence theorems for a family of iterations with one parameter are obtained. An estimation of the radius of the attraction ball is shown. At last two examples are given.展开更多
Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phas...Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.展开更多
In order to improve the performance of time difference of arrival(TDOA)localization,a nonlinear least squares algorithm is proposed in this paper.Firstly,based on the criterion of the minimized sum of square error of ...In order to improve the performance of time difference of arrival(TDOA)localization,a nonlinear least squares algorithm is proposed in this paper.Firstly,based on the criterion of the minimized sum of square error of time difference of arrival,the location estimation is expressed as an optimal problem of a non-linear programming.Then,an initial point is obtained using the semi-definite programming.And finally,the location is extracted from the local optimal solution acquired by Newton iterations.Simulation results show that when the number of anchor nodes is large,the performance of the proposed algorithm will be significantly better than that of semi-definite programming approach with the increase of measurement noise.展开更多
This paper considers Stokes and Newton iterations to solve stationary Navier- Stokes equations based on the finite element discretization. We obtain new sufficient conditions of stability and convergence for the two i...This paper considers Stokes and Newton iterations to solve stationary Navier- Stokes equations based on the finite element discretization. We obtain new sufficient conditions of stability and convergence for the two iterations. Specifically, when 0 〈 σ =N||f||-1/v2≤1/√2+1 , the Stokes iteration is stable and convergent, where N is defined in the paper. When 0 〈 σ ≤5/11, the Newton iteration is stable and convergent. This work gives a more accurate admissible range of data for stability and convergence of the two schemes, which improves the previous results. A numerical test is given to verify the theory.展开更多
We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-st...We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-step iterative method of Khan et al. [1] as a special case. Our results are new in hyperbolic spaces and generalize many known results in Banach spaces and CAT(0) spaces, simultaneously.展开更多
The iteration maps of Euler family for finding zeros of an operatorf in Banach spaces is defined as the partial sum of Taylor expansion of the local inversef z -1 off atz. The unified convergence theorem is establishe...The iteration maps of Euler family for finding zeros of an operatorf in Banach spaces is defined as the partial sum of Taylor expansion of the local inversef z -1 off atz. The unified convergence theorem is established for the iterations of Euler family under the assumption that $\alpha \leqslant 3 - 2\sqrt 2 $ , while the strong condition thatf is analytic in Smale’s criterion α is replaced by the weak condition thatf is of finite order derivative.展开更多
In the sense of the nonlinear multisplitting and based on the principle of suffi-ciently using the delayed information, we propose models of asynchronous parallelaccelerated overrelaxation iteration methods for solvin...In the sense of the nonlinear multisplitting and based on the principle of suffi-ciently using the delayed information, we propose models of asynchronous parallelaccelerated overrelaxation iteration methods for solving large scale system of non-linear equations. Under proper conditions, we set up the local convergence theoriesof these new method models.展开更多
Presents the first estimation conditions for Nourein iterations for simultaneous finding all zeros of a polynomial under which the iteration processes are guaranteed to converge. Computational formulas; Theorems and p...Presents the first estimation conditions for Nourein iterations for simultaneous finding all zeros of a polynomial under which the iteration processes are guaranteed to converge. Computational formulas; Theorems and proofs.展开更多
In this paper, we suggest and analyse a three-step iterative scheme with errors for solving nonlinear strongly accretive operator equation Tx = f without the Lipshitz condition. The results presented in this paper imp...In this paper, we suggest and analyse a three-step iterative scheme with errors for solving nonlinear strongly accretive operator equation Tx = f without the Lipshitz condition. The results presented in this paper improve and extend current results in the more general setting.展开更多
Presents preconditioning matrices having parallel computing function for the coefficient matrix and a class of parallel hybrid algebraic multilevel iteration methods for solving linear equations. Solution to elliptic ...Presents preconditioning matrices having parallel computing function for the coefficient matrix and a class of parallel hybrid algebraic multilevel iteration methods for solving linear equations. Solution to elliptic boundary value problem; Discussion on symmetric positive definite matrix; Computational complexities.展开更多
SINCE there are many difficulties in finding a single root of a polynomial, it becomes moreand more important for parallel iterations to determine all roots simultaneously. Among theproposed methods, the iterative fam...SINCE there are many difficulties in finding a single root of a polynomial, it becomes moreand more important for parallel iterations to determine all roots simultaneously. Among theproposed methods, the iterative fami1y produced by paralleling the iterative family of Halleyusing Bell’s polynomial appears the most systematic and richest. In fact the paper becomesthe main contents of the monograph.展开更多
In this paper, we establish the dynamical systems of iterated entire aigebroid functions. According to dynamics, we give a classification theorem of entire aigebroid functions. Some typical properties on Julia set and...In this paper, we establish the dynamical systems of iterated entire aigebroid functions. According to dynamics, we give a classification theorem of entire aigebroid functions. Some typical properties on Julia set and Fatou set are proved. And we obtain a similar property on the dstribution of J(f) and Vf.展开更多
We establish a general convergence theory of the Shift-Invert Residual Arnoldi(SIRA)method for computing a simple eigenvalue nearest to a given targetσand the associated eigenvector.In SIRA,a subspace expansion vecto...We establish a general convergence theory of the Shift-Invert Residual Arnoldi(SIRA)method for computing a simple eigenvalue nearest to a given targetσand the associated eigenvector.In SIRA,a subspace expansion vector at each step is obtained by solving a certain inner linear system.We prove that the inexact SIRA method mimics the exact SIRA well,i.e.,the former uses almost the same outer iterations to achieve the convergence as the latter does if all the inner linear systems are iteratively solved with low or modest accuracy during outer iterations.Based on the theory,we design practical stopping criteria for inner solves.Our analysis is on one step expansion of subspace and the approach applies to the Jacobi-Davidson(JD)method with the fixed targetσas well,and a similar general convergence theory is obtained for it.Numerical experiments confirm our theory and demonstrate that the inexact SIRA and JD are similarly effective and are considerably superior to the inexact SIA.展开更多
A class of asynchronous matrix multi-splitting multi-parameter relaxation methods, including the asynchronous matrix multisplitting SAOR, SSOR and SGS methods as well. as the known asynchronous matrix multisplitting A...A class of asynchronous matrix multi-splitting multi-parameter relaxation methods, including the asynchronous matrix multisplitting SAOR, SSOR and SGS methods as well. as the known asynchronous matrix multisplitting AOR, SOR and GS methods, etc., is proposed for solving the large sparse systems of linear equations by making use of the principle of sufficiently using the delayed information. These new methods can greatly execute the parallel computational efficiency of the MIMD-systems, and are shown to be convergent when the coefficient matrices are H-matrices. Moreover, necessary and sufficient conditions ensuring the convergence of these methods are concluded for the case that the coefficient matrices are L-matrices.展开更多
Based on the Homotopy Analysis Method, a direct numerical method for strongly nonlinear problems was proposed. The 2-D laminar flow over semi-infinite plate was used. The method can give the accurate enough approximat...Based on the Homotopy Analysis Method, a direct numerical method for strongly nonlinear problems was proposed. The 2-D laminar flow over semi-infinite plate was used. The method can give the accurate enough approximations of a strongly nonlinear problem by means of no iteration and can provide a family of iterative formulas with traditional approaches.展开更多
文摘In this paper, equivalence between the Mann and Ishikawa iterations for a generalized contraction mapping in cone subset of a real Banach space is discussed.
文摘This note explores the relations between two different methods. The first one is the Alternating Least Squares (ALS) method for calculating a rank<em>-k</em> approximation of a real <em>m</em>×<em>n</em> matrix, <em>A</em>. This method has important applications in nonnegative matrix factorizations, in matrix completion problems, and in tensor approximations. The second method is called Orthogonal Iterations. Other names of this method are Subspace Iterations, Simultaneous Iterations, and block-Power method. Given a real symmetric matrix, <em>G</em>, this method computes<em> k</em> dominant eigenvectors of <em>G</em>. To see the relation between these methods we assume that <em>G </em>=<em> A</em><sup>T</sup> <em>A</em>. It is shown that in this case the two methods generate the same sequence of subspaces, and the same sequence of low-rank approximations. This equivalence provides new insight into the convergence properties of both methods.
文摘A new restoration algorithm based on double loops and alternant iterations is proposed to restore the object image effectively from a few frames of turbulence-degraded images, Based on the double loops, the iterative relations for estimating the turbulent point spread function PSF and object image alternately are derived. The restoration experiments have been made on computers, showing that the proposed algorithm can obtain the optimal estimations of the object and the point spread function, with the feasibility and practicality of the proposed algorithm being convincing.
文摘In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduced that enables successive decomposition of high Eigen values of the iteration matrix to enable convergence. While extrapolation of a convergent fixed point iteration using a geometric series sum is a known form of Aitken acceleration, it is shown that in this paper, the same formula can be used to estimate the solution of sets of linear equations from diverging Gauss-Seidel iterations. In both convergent and divergent iterations, the ratios of differences among the consecutive values of iteration eventually form a convergent (divergent) series with a factor equal to the largest Eigen value of the iteration matrix. Higher order Aitken extrapolation is shown to eliminate the influence of dominant Eigen values of the iteration matrix in successive order until the iteration is determined by the lowest possible Eigen values. For the convergent part of the Gauss-Seidel iteration, further acceleration is made possible by coupling of the extrapolation technique with the successive over relaxation (SOR) method. Application examples from both convergent and divergent iterations have been provided. Coupling of the extrapolation with the SOR technique is also illustrated for a steady state two dimensional heat flow problem which was solved using MATLAB programming.
基金supported by the National Natural Science Foundation of China(40974018)Nationa l863 Plan Projects(2009AA12Z307)
文摘The simultaneous iterations rithms of the ART family. It is used reconstruction technique (SIRT) widely in tomography because of is one of several reconstruction algoits convenience in dealing with large sparse matrices. Its theoretical background and iteration model are discussed at the beginning of this paper. Then, the implementation of the SIRT to reconstruct the three-dimensional distribution of water vapor by simulation is discussed. The results show that the SIRT can function effectively in water vapor tomography, obtain rapid convergence, and be implemented more easily than inversion.
基金Supported by Shanghai Municipal Foundation of Selected Academic Research and the National Natural Science Foundation of China(10571059,10571060).
文摘Under the weak Lipschitz condition about the solution of the equation, convergence theorems for a family of iterations with one parameter are obtained. An estimation of the radius of the attraction ball is shown. At last two examples are given.
基金Project supported by the National Natural Science Foundation of China(Nos.11971410 and12071404)the Natural Science Foundation of Hunan Province of China(No.2019JJ40279)+2 种基金the Excellent Youth Program of Scientific Research Project of Hunan Provincial Department of Education(Nos.18B064 and 20B564)the China Postdoctoral Science Foundation(Nos.2018T110073 and 2018M631402)the International Scientific and Technological Innovation Cooperation Base of Hunan Province for Computational Science(No.2018WK4006)。
文摘Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.
基金This study was supported by the“High level research and training project for professional leaders of teachers in Higher Vocational Colleges in Jiangsu Province”.
文摘In order to improve the performance of time difference of arrival(TDOA)localization,a nonlinear least squares algorithm is proposed in this paper.Firstly,based on the criterion of the minimized sum of square error of time difference of arrival,the location estimation is expressed as an optimal problem of a non-linear programming.Then,an initial point is obtained using the semi-definite programming.And finally,the location is extracted from the local optimal solution acquired by Newton iterations.Simulation results show that when the number of anchor nodes is large,the performance of the proposed algorithm will be significantly better than that of semi-definite programming approach with the increase of measurement noise.
基金supported by the National Natural Science Foundation of China(No.11271298)
文摘This paper considers Stokes and Newton iterations to solve stationary Navier- Stokes equations based on the finite element discretization. We obtain new sufficient conditions of stability and convergence for the two iterations. Specifically, when 0 〈 σ =N||f||-1/v2≤1/√2+1 , the Stokes iteration is stable and convergent, where N is defined in the paper. When 0 〈 σ ≤5/11, the Newton iteration is stable and convergent. This work gives a more accurate admissible range of data for stability and convergence of the two schemes, which improves the previous results. A numerical test is given to verify the theory.
文摘We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-step iterative method of Khan et al. [1] as a special case. Our results are new in hyperbolic spaces and generalize many known results in Banach spaces and CAT(0) spaces, simultaneously.
文摘The iteration maps of Euler family for finding zeros of an operatorf in Banach spaces is defined as the partial sum of Taylor expansion of the local inversef z -1 off atz. The unified convergence theorem is established for the iterations of Euler family under the assumption that $\alpha \leqslant 3 - 2\sqrt 2 $ , while the strong condition thatf is analytic in Smale’s criterion α is replaced by the weak condition thatf is of finite order derivative.
文摘In the sense of the nonlinear multisplitting and based on the principle of suffi-ciently using the delayed information, we propose models of asynchronous parallelaccelerated overrelaxation iteration methods for solving large scale system of non-linear equations. Under proper conditions, we set up the local convergence theoriesof these new method models.
基金National Natural Science Foundation of China Natural ScienceFoundation of Zhejiang Province.
文摘Presents the first estimation conditions for Nourein iterations for simultaneous finding all zeros of a polynomial under which the iteration processes are guaranteed to converge. Computational formulas; Theorems and proofs.
文摘In this paper, we suggest and analyse a three-step iterative scheme with errors for solving nonlinear strongly accretive operator equation Tx = f without the Lipshitz condition. The results presented in this paper improve and extend current results in the more general setting.
基金Subsidized by the Special Funds for Major State Basic Research Projects G1999032803 and Suported bythe National Natural Scienc
文摘Presents preconditioning matrices having parallel computing function for the coefficient matrix and a class of parallel hybrid algebraic multilevel iteration methods for solving linear equations. Solution to elliptic boundary value problem; Discussion on symmetric positive definite matrix; Computational complexities.
文摘SINCE there are many difficulties in finding a single root of a polynomial, it becomes moreand more important for parallel iterations to determine all roots simultaneously. Among theproposed methods, the iterative fami1y produced by paralleling the iterative family of Halleyusing Bell’s polynomial appears the most systematic and richest. In fact the paper becomesthe main contents of the monograph.
基金Project supported by the Tianyuan Foundation of China
文摘In this paper, we establish the dynamical systems of iterated entire aigebroid functions. According to dynamics, we give a classification theorem of entire aigebroid functions. Some typical properties on Julia set and Fatou set are proved. And we obtain a similar property on the dstribution of J(f) and Vf.
基金supported by National Basic Research Program of China(Grant No.2011CB302400)National Natural Science Foundation of China(Grant No.11071140)
文摘We establish a general convergence theory of the Shift-Invert Residual Arnoldi(SIRA)method for computing a simple eigenvalue nearest to a given targetσand the associated eigenvector.In SIRA,a subspace expansion vector at each step is obtained by solving a certain inner linear system.We prove that the inexact SIRA method mimics the exact SIRA well,i.e.,the former uses almost the same outer iterations to achieve the convergence as the latter does if all the inner linear systems are iteratively solved with low or modest accuracy during outer iterations.Based on the theory,we design practical stopping criteria for inner solves.Our analysis is on one step expansion of subspace and the approach applies to the Jacobi-Davidson(JD)method with the fixed targetσas well,and a similar general convergence theory is obtained for it.Numerical experiments confirm our theory and demonstrate that the inexact SIRA and JD are similarly effective and are considerably superior to the inexact SIA.
基金Project 19601036 supported by the National Natural Science Foundation of China.
文摘A class of asynchronous matrix multi-splitting multi-parameter relaxation methods, including the asynchronous matrix multisplitting SAOR, SSOR and SGS methods as well. as the known asynchronous matrix multisplitting AOR, SOR and GS methods, etc., is proposed for solving the large sparse systems of linear equations by making use of the principle of sufficiently using the delayed information. These new methods can greatly execute the parallel computational efficiency of the MIMD-systems, and are shown to be convergent when the coefficient matrices are H-matrices. Moreover, necessary and sufficient conditions ensuring the convergence of these methods are concluded for the case that the coefficient matrices are L-matrices.
文摘Based on the Homotopy Analysis Method, a direct numerical method for strongly nonlinear problems was proposed. The 2-D laminar flow over semi-infinite plate was used. The method can give the accurate enough approximations of a strongly nonlinear problem by means of no iteration and can provide a family of iterative formulas with traditional approaches.