In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order...In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.展开更多
In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-m...In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.展开更多
In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of th...In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 (1997), 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67(2007), 1758-1763].展开更多
In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is...In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.展开更多
Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterativ...Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.展开更多
Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed ...Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform(FFT)algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.展开更多
In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validi...In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validity and performance of these iterative methods, we have applied to solve some nonlinear problems.展开更多
Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial ...Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial guess or the subsequent updated solution.In the second type,the NBE is rearranged into a quadratic form of the absolute vorticity with the positive root of this quadratic form used in the form of a Poisson equation to solve NBE iteratively.The two methods are rederived by expanding the solution asymptotically upon a small Rossby number,and a criterion for optimally truncating the asymptotic expansion is proposed to obtain the super-asymptotic approximation of the solution.For each rederived method,two iterative procedures are designed using the integral-form Poisson solver versus the over-relaxation scheme to solve the boundary value problem in each iteration.Upon testing with analytically formulated wavering jet flows on the synoptic,sub-synoptic and meso-αscales,the iterative procedure designed for the first method with the Poisson solver,named M1a,is found to be the most accurate and efficient.For the synoptic wavering jet flow in which the NBE is entirely elliptic,M1a is extremely accurate.For the sub-synoptic wavering jet flow in which the NBE is mostly elliptic,M1a is sufficiently accurate.For the meso-αwavering jet flow in which the NBE is partially hyperbolic so its boundary value problem becomes seriously ill-posed,M1a can effectively reduce the solution error for the cyclonically curved part of the wavering jet flow,but not for the anti-cyclonically curved part.展开更多
This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods...This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.展开更多
There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,...There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,namely the basins of attraction of the method.The purpose of this study is to compare several iterative schemes for nonlinear equations.All the selected schemes are of the third-order of convergence and most of them have the same efficiency index.The comparison depends on the basins of attraction of the iterative techniques when applied on several polynomials of different degrees.As a comparison,we determine the CPU time(in seconds)needed by each scheme to obtain the basins of attraction,besides,we illustrate the area of convergence of these schemes by finding the number of convergent and divergent points in a selected range for all methods.Comparisons confirm the fact that basins of attraction differ for iterative methods of different orders,furthermore,they vary for iterative methods of the same order even if they have the same efficiency index.Consequently,this leads to the need for a new index that reflects the real efficiency of the iterative scheme instead of the commonly used efficiency index.展开更多
Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced...Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced to linear systems.Due to the typical ill-posedness of inverse problems,the reduced linear systems are often illposed,especially when their scales are large.This brings great computational difficulty.Particularly,a small perturbation in the right side of an ill-posed linear system may cause a dramatical change in the solution.Therefore,regularization methods should be adopted for stable solutions.In this paper,a new class of accelerated iterative regularization methods is applied to solve this kind of large-scale ill-posed linear systems.An iterative scheme becomes a regularization method only when the iteration is early terminated.And a Morozov’s discrepancy principle is applied for the stop criterion.Compared with the conventional Landweber iteration,the new methods have acceleration effect,and can be compared to the well-known acceleratedν-method and Nesterov method.From the numerical results,it is observed that using appropriate discretization schemes,the proposed methods even have better behavior when comparing withν-method and Nesterov method.展开更多
In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Conver...In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.展开更多
In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and di...In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.展开更多
Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods....Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods. The remarkable feature of iterative methods based on the new class of matrices is that the methods are easily implemented for parallel computation. In particular, a generalization of the accelerated overrelaxation method (GAOR) is introduced. Some theories of the AOR method are extended to the generalized method to include a wide class of matrices. The convergence of the new method is derived for Hermitian positive definite matrices. Finally, some examples are given in order to show the superiority of the new method.展开更多
In this paper, we present and analyze modified families of predictor-corrector iterative methods for finding simple zeros of univariate nonlinear equations, permitting near the root. The main advantage of our methods ...In this paper, we present and analyze modified families of predictor-corrector iterative methods for finding simple zeros of univariate nonlinear equations, permitting near the root. The main advantage of our methods is that they perform better and moreover, have the same efficiency indices as that of existing multipoint iterative methods. Furthermore, the convergence analysis of the new methods is discussed and several examples are given to illustrate their efficiency.展开更多
A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-or...A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.展开更多
Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is ...Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is discussed. The second order iterative methods behave quite similar to first order methods and the development of efficient preconditioners for solving the original linear system is a decisive factor for making the second order iterative methods superior to the first order iterative methods. Adaptive preconditioned Conjugate Gradient methods using explicit approximate preconditioners for solving efficiently large sparse systems of algebraic equations are also presented. The generalized Approximate Inverse Matrix techniques can be efficiently used in conjunction with explicit iterative schemes leading to effective composite semi-direct solution methods for solving large linear systems of algebraic equations.展开更多
Iterative methods for solving discrete optimal control problems are constructed and investigated. These discrete problems arise when approximating by finite difference method or by finite element method the optimal co...Iterative methods for solving discrete optimal control problems are constructed and investigated. These discrete problems arise when approximating by finite difference method or by finite element method the optimal control problems which contain a linear elliptic boundary value problem as a state equation, control in the righthand side of the equation or in the boundary conditions, and point-wise constraints for both state and control functions. The convergence of the constructed iterative methods is proved, the implementation problems are discussed, and the numerical comparison of the methods is executed.展开更多
Large catalogues of classified galaxy images have been useful in many studies of the universe in astronomy. There are too many objects to classify manually in the Sloan Digital Sky Survey, one of the premier data sour...Large catalogues of classified galaxy images have been useful in many studies of the universe in astronomy. There are too many objects to classify manually in the Sloan Digital Sky Survey, one of the premier data sources in astronomy. Therefore, efficient machine learning and classification algorithms are required to automate the classifying process. We propose to apply the Support Vector Machine (SVM) algorithm to classify galaxy morphologies and Krylov iterative methods to improve runtime of the classification. The accuracy of the classification is measured on various categories of galaxies from the survey. A three-class algorithm is presented that makes use of multiple SVMs. This algorithm is used to assign the categories of spiral, elliptical, and irregular galaxies. A selection of Krylov iterative solvers are compared based on their efficiency and accuracy of the resulting classification. The experimental results demonstrate that runtime can be significantly improved by utilizing Krylov iterative methods without impacting classification accuracy. The generalized minimal residual method (GMRES) is shown to be the most efficient solver to classify galaxy morphologies.展开更多
This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of init...This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of initially known and unknown multiplicity. Efficient methods are presented in this note for the evaluation of the multiplicity index of the root being sought. Also reviewed here are super-linear and super-cubic methods that converge contrarily or alternatingly, enabling us, not only to approach the root briskly and confidently but also to actually bound and bracket it as we progress.展开更多
文摘In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.
文摘In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.
基金the Thailand Research Fund for financial support under Grant BRG5280016
文摘In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 (1997), 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67(2007), 1758-1763].
文摘In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.
基金Supported by the National Natural Science Foundation of China(61272300)
文摘Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.
基金supported by the National Natural Science Foundation of China(41304022,41174026,41104047)the National 973 Foundation(61322201,2013CB733303)+1 种基金the Key laboratory Foundation of Geo-space Environment and Geodesy of the Ministry of Education(13-01-08)the Youth Innovation Foundation of High Resolution Earth Observation(GFZX04060103-5-12)
文摘Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform(FFT)algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.
文摘In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validity and performance of these iterative methods, we have applied to solve some nonlinear problems.
基金the NSF of China Grants 91937301 and 41675060,the National Key Scientific and Technological Infrastructure Project"EarthLab",and the ONR Grants N000141712375 and N000142012449 to the University of Oklahoma(OU)The numerical experiments were performed at the OU supercomputer SchoonerCIMMS by NOAA/Office of Oceanic and Atmospheric Research under NOAA-OU Cooperative Agreement#NA110AR4320072,U.S.Department of Commerce.
文摘Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial guess or the subsequent updated solution.In the second type,the NBE is rearranged into a quadratic form of the absolute vorticity with the positive root of this quadratic form used in the form of a Poisson equation to solve NBE iteratively.The two methods are rederived by expanding the solution asymptotically upon a small Rossby number,and a criterion for optimally truncating the asymptotic expansion is proposed to obtain the super-asymptotic approximation of the solution.For each rederived method,two iterative procedures are designed using the integral-form Poisson solver versus the over-relaxation scheme to solve the boundary value problem in each iteration.Upon testing with analytically formulated wavering jet flows on the synoptic,sub-synoptic and meso-αscales,the iterative procedure designed for the first method with the Poisson solver,named M1a,is found to be the most accurate and efficient.For the synoptic wavering jet flow in which the NBE is entirely elliptic,M1a is extremely accurate.For the sub-synoptic wavering jet flow in which the NBE is mostly elliptic,M1a is sufficiently accurate.For the meso-αwavering jet flow in which the NBE is partially hyperbolic so its boundary value problem becomes seriously ill-posed,M1a can effectively reduce the solution error for the cyclonically curved part of the wavering jet flow,but not for the anti-cyclonically curved part.
基金This work was supported by the National Natural Science Foundation of China
文摘This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.
基金We are grateful for the financial support from UKM’s research Grant GUP-2019-033。
文摘There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,namely the basins of attraction of the method.The purpose of this study is to compare several iterative schemes for nonlinear equations.All the selected schemes are of the third-order of convergence and most of them have the same efficiency index.The comparison depends on the basins of attraction of the iterative techniques when applied on several polynomials of different degrees.As a comparison,we determine the CPU time(in seconds)needed by each scheme to obtain the basins of attraction,besides,we illustrate the area of convergence of these schemes by finding the number of convergent and divergent points in a selected range for all methods.Comparisons confirm the fact that basins of attraction differ for iterative methods of different orders,furthermore,they vary for iterative methods of the same order even if they have the same efficiency index.Consequently,this leads to the need for a new index that reflects the real efficiency of the iterative scheme instead of the commonly used efficiency index.
基金supported by the Natural Science Foundation of China (Nos. 11971230, 12071215)the Fundamental Research Funds for the Central Universities(No. NS2018047)the 2019 Graduate Innovation Base(Laboratory)Open Fund of Jiangsu Province(No. Kfjj20190804)
文摘Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced to linear systems.Due to the typical ill-posedness of inverse problems,the reduced linear systems are often illposed,especially when their scales are large.This brings great computational difficulty.Particularly,a small perturbation in the right side of an ill-posed linear system may cause a dramatical change in the solution.Therefore,regularization methods should be adopted for stable solutions.In this paper,a new class of accelerated iterative regularization methods is applied to solve this kind of large-scale ill-posed linear systems.An iterative scheme becomes a regularization method only when the iteration is early terminated.And a Morozov’s discrepancy principle is applied for the stop criterion.Compared with the conventional Landweber iteration,the new methods have acceleration effect,and can be compared to the well-known acceleratedν-method and Nesterov method.From the numerical results,it is observed that using appropriate discretization schemes,the proposed methods even have better behavior when comparing withν-method and Nesterov method.
基金the Natural Science Foundation of China(Grant Nos.61673169,11301127,11701176,11626101,and 11601485)The Natural Science Foundation of Huzhou City(Grant No.2018YZ07).
文摘In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.
基金national natural science foundation natural science foundation of Gansu province.
文摘In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.
基金Project supported by the Natural Science Foundation of Liaoning Province of China (No.20022021)
文摘Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods. The remarkable feature of iterative methods based on the new class of matrices is that the methods are easily implemented for parallel computation. In particular, a generalization of the accelerated overrelaxation method (GAOR) is introduced. Some theories of the AOR method are extended to the generalized method to include a wide class of matrices. The convergence of the new method is derived for Hermitian positive definite matrices. Finally, some examples are given in order to show the superiority of the new method.
文摘In this paper, we present and analyze modified families of predictor-corrector iterative methods for finding simple zeros of univariate nonlinear equations, permitting near the root. The main advantage of our methods is that they perform better and moreover, have the same efficiency indices as that of existing multipoint iterative methods. Furthermore, the convergence analysis of the new methods is discussed and several examples are given to illustrate their efficiency.
基金Supported by National Science Foundation of China(Grant 10871179)the National Basic Research Programme of China(Grant 2008CB717806)the Department of Education of Zhejiang Province(GrantY200803559).
文摘A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.
文摘Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is discussed. The second order iterative methods behave quite similar to first order methods and the development of efficient preconditioners for solving the original linear system is a decisive factor for making the second order iterative methods superior to the first order iterative methods. Adaptive preconditioned Conjugate Gradient methods using explicit approximate preconditioners for solving efficiently large sparse systems of algebraic equations are also presented. The generalized Approximate Inverse Matrix techniques can be efficiently used in conjunction with explicit iterative schemes leading to effective composite semi-direct solution methods for solving large linear systems of algebraic equations.
文摘Iterative methods for solving discrete optimal control problems are constructed and investigated. These discrete problems arise when approximating by finite difference method or by finite element method the optimal control problems which contain a linear elliptic boundary value problem as a state equation, control in the righthand side of the equation or in the boundary conditions, and point-wise constraints for both state and control functions. The convergence of the constructed iterative methods is proved, the implementation problems are discussed, and the numerical comparison of the methods is executed.
文摘Large catalogues of classified galaxy images have been useful in many studies of the universe in astronomy. There are too many objects to classify manually in the Sloan Digital Sky Survey, one of the premier data sources in astronomy. Therefore, efficient machine learning and classification algorithms are required to automate the classifying process. We propose to apply the Support Vector Machine (SVM) algorithm to classify galaxy morphologies and Krylov iterative methods to improve runtime of the classification. The accuracy of the classification is measured on various categories of galaxies from the survey. A three-class algorithm is presented that makes use of multiple SVMs. This algorithm is used to assign the categories of spiral, elliptical, and irregular galaxies. A selection of Krylov iterative solvers are compared based on their efficiency and accuracy of the resulting classification. The experimental results demonstrate that runtime can be significantly improved by utilizing Krylov iterative methods without impacting classification accuracy. The generalized minimal residual method (GMRES) is shown to be the most efficient solver to classify galaxy morphologies.
文摘This paper considers practical, high-order methods for the iterative location of the roots of nonlinear equations, one at a time. Special attention is being paid to algorithms also applicable to multiple roots of initially known and unknown multiplicity. Efficient methods are presented in this note for the evaluation of the multiplicity index of the root being sought. Also reviewed here are super-linear and super-cubic methods that converge contrarily or alternatingly, enabling us, not only to approach the root briskly and confidently but also to actually bound and bracket it as we progress.