This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with...This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.展开更多
Inexact Newton methods are constructed by combining Newton's method with another iterative method that is used to solve the Newton equations inexactly. In this paper, we establish two semilocal convergence theorems f...Inexact Newton methods are constructed by combining Newton's method with another iterative method that is used to solve the Newton equations inexactly. In this paper, we establish two semilocal convergence theorems for the inexact Newton methods. When these two theorems are specified to Newton's method, we obtain a different Newton-Kantorovich theorem about Newton's method. When the iterative method for solving the Newton equations is specified to be the splitting method, we get two estimates about the iteration steps for the special inexact Newton methods.展开更多
We propose an inexact Newton method with a filter line search algorithm for nonconvex equality constrained optimization. Inexact Newton's methods are needed for large-scale applications which the iteration matrix can...We propose an inexact Newton method with a filter line search algorithm for nonconvex equality constrained optimization. Inexact Newton's methods are needed for large-scale applications which the iteration matrix cannot be explicitly formed or factored. We incorporate inexact Newton strategies in filter line search, yielding algorithm that can ensure global convergence. An analysis of the global behavior of the algorithm and numerical results on a collection of test problems are presented.展开更多
In this paper,we introduce a physics-based nonlinear preconditioned Inexact Newton Method(INB)for the multiphysical simulation of fractured reservoirs.Instead of solving the partial differential equations(PDE)exactly,...In this paper,we introduce a physics-based nonlinear preconditioned Inexact Newton Method(INB)for the multiphysical simulation of fractured reservoirs.Instead of solving the partial differential equations(PDE)exactly,Inexact Newton method finds a direction for the iteration and solves the equations inexactly with fewer iterations.However,when the equations are not smooth enough,especially when lo-cal discontinuities exits,and when proper preconditioning operations are not adopted,the Inexact Newton method may be slow or even stagnant.As pointed out by Keyes et al.[1],multi-physical numerical simulation faces several challenges,one of which is the local-scale nonlinearity and discontinuity.In this work,we have proposed and studied a nonlinear preconditioner to improve the performance of Inexact Newton Method.The nonlinear preconditioner is essentially a physics-based strategy to adaptively identify and eliminate the highly nonlinear zones.The proposed algorithm has been implemented into our fully coupled,fully implicit THM reservoir simulator(Wang et al.[2,3])to study the effects of cold water injection on fractured petroleum reservoirs.The results of this work show that after the implementation of this nonlinear preconditioner,the iterative solver has become significantly more robust and efficient.展开更多
We first propose a new class of smoothing functions for the non- linear complementarity function which contains the well-known Chen-Harker- Kanzow-Smale smoothing function and Huang-Han-Chen smoothing function as spec...We first propose a new class of smoothing functions for the non- linear complementarity function which contains the well-known Chen-Harker- Kanzow-Smale smoothing function and Huang-Han-Chen smoothing function as special cases, and then present a smoothing inexact Newton algorithm for the P0 nonlinear complementarity problem. The global convergence and local superlinear convergence are established. Preliminary numerical results indicate the feasibility and efficiency of the algorithm.展开更多
A conic Newton method is attractive because it converges to a local minimizzer rapidly from any sufficiently good initial guess. However, it may be expensive to solve the conic Newton equation at each iterate. In this...A conic Newton method is attractive because it converges to a local minimizzer rapidly from any sufficiently good initial guess. However, it may be expensive to solve the conic Newton equation at each iterate. In this paper we consider an inexact conic Newton method, which solves the couic Newton equation oldy approximately and in sonm unspecified manner. Furthermore, we show that such method is locally convergent and characterizes the order of convergence in terms of the rate of convergence of the relative residuals.展开更多
In this paper, the non-quasi-Newton's family with inexact line search applied to unconstrained optimization problems is studied. A new update formula for non-quasi-Newton's family is proposed. It is proved that the ...In this paper, the non-quasi-Newton's family with inexact line search applied to unconstrained optimization problems is studied. A new update formula for non-quasi-Newton's family is proposed. It is proved that the constituted algorithm with either Wolfe-type or Armijotype line search converges globally and Q-superlinearly if the function to be minimized has Lipschitz continuous gradient.展开更多
In full waveform inversion (FWI), Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the He...In full waveform inversion (FWI), Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the Hessian matrix and its inverse. Although the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) or Hessian-free inexact Newton (HFN) methods are able to use approximate Hessian information, the information they collect is limited. The two methods can be interlaced because they are able to provide Hessian information for each other; however, the performance of the hybrid iterative method is dependent on the effective switch between the two methods. We have designed a new scheme to realize the dynamic switch between the two methods based on the decrease ratio (DR) of the misfit function (objective function), and we propose a modified hybrid iterative optimization method. In the new scheme, we compare the DR of the two methods for a given computational cost, and choose the method with a faster DR. Using these steps, the modified method always implements the most efficient method. The results of Marmousi and overthrust model testings indicate that the convergence with our modified method is significantly faster than that in the L-BFGS method with no loss of inversion quality. Moreover, our modified outperforms the enriched method by a little speedup of the convergence. It also exhibits better efficiency than the HFN method.展开更多
基金Project supported by the National Natural Science Foundation of China (No. 10871130)the Ph. D.Programs Foundation of Ministry of Education of China (No. 20093127110005)the Shanghai Leading Academic Discipline Project (No. T0401)
文摘This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.
基金Supported by State Key Laboratory of Scientific/Engineering Computing,Chinese Academy of Sciencesthe National Natural Science Foundation of China (10571059,10571060).
文摘Inexact Newton methods are constructed by combining Newton's method with another iterative method that is used to solve the Newton equations inexactly. In this paper, we establish two semilocal convergence theorems for the inexact Newton methods. When these two theorems are specified to Newton's method, we obtain a different Newton-Kantorovich theorem about Newton's method. When the iterative method for solving the Newton equations is specified to be the splitting method, we get two estimates about the iteration steps for the special inexact Newton methods.
基金Supported in part by the National Natural Science Foundation of China under Grant No.11371253Natural Science Foundation of Hunan Province under Grant No.2016JJ2038the project of Scientific Research Fund of Hunan Provincial Education Department under Grant No.14B044
文摘We propose an inexact Newton method with a filter line search algorithm for nonconvex equality constrained optimization. Inexact Newton's methods are needed for large-scale applications which the iteration matrix cannot be explicitly formed or factored. We incorporate inexact Newton strategies in filter line search, yielding algorithm that can ensure global convergence. An analysis of the global behavior of the algorithm and numerical results on a collection of test problems are presented.
文摘In this paper,we introduce a physics-based nonlinear preconditioned Inexact Newton Method(INB)for the multiphysical simulation of fractured reservoirs.Instead of solving the partial differential equations(PDE)exactly,Inexact Newton method finds a direction for the iteration and solves the equations inexactly with fewer iterations.However,when the equations are not smooth enough,especially when lo-cal discontinuities exits,and when proper preconditioning operations are not adopted,the Inexact Newton method may be slow or even stagnant.As pointed out by Keyes et al.[1],multi-physical numerical simulation faces several challenges,one of which is the local-scale nonlinearity and discontinuity.In this work,we have proposed and studied a nonlinear preconditioner to improve the performance of Inexact Newton Method.The nonlinear preconditioner is essentially a physics-based strategy to adaptively identify and eliminate the highly nonlinear zones.The proposed algorithm has been implemented into our fully coupled,fully implicit THM reservoir simulator(Wang et al.[2,3])to study the effects of cold water injection on fractured petroleum reservoirs.The results of this work show that after the implementation of this nonlinear preconditioner,the iterative solver has become significantly more robust and efficient.
文摘We first propose a new class of smoothing functions for the non- linear complementarity function which contains the well-known Chen-Harker- Kanzow-Smale smoothing function and Huang-Han-Chen smoothing function as special cases, and then present a smoothing inexact Newton algorithm for the P0 nonlinear complementarity problem. The global convergence and local superlinear convergence are established. Preliminary numerical results indicate the feasibility and efficiency of the algorithm.
文摘A conic Newton method is attractive because it converges to a local minimizzer rapidly from any sufficiently good initial guess. However, it may be expensive to solve the conic Newton equation at each iterate. In this paper we consider an inexact conic Newton method, which solves the couic Newton equation oldy approximately and in sonm unspecified manner. Furthermore, we show that such method is locally convergent and characterizes the order of convergence in terms of the rate of convergence of the relative residuals.
文摘In this paper, the non-quasi-Newton's family with inexact line search applied to unconstrained optimization problems is studied. A new update formula for non-quasi-Newton's family is proposed. It is proved that the constituted algorithm with either Wolfe-type or Armijotype line search converges globally and Q-superlinearly if the function to be minimized has Lipschitz continuous gradient.
基金financially supported by the National Important and Special Project on Science and Technology(2011ZX05005-005-007HZ)the National Natural Science Foundation of China(No.41274116)
文摘In full waveform inversion (FWI), Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the Hessian matrix and its inverse. Although the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) or Hessian-free inexact Newton (HFN) methods are able to use approximate Hessian information, the information they collect is limited. The two methods can be interlaced because they are able to provide Hessian information for each other; however, the performance of the hybrid iterative method is dependent on the effective switch between the two methods. We have designed a new scheme to realize the dynamic switch between the two methods based on the decrease ratio (DR) of the misfit function (objective function), and we propose a modified hybrid iterative optimization method. In the new scheme, we compare the DR of the two methods for a given computational cost, and choose the method with a faster DR. Using these steps, the modified method always implements the most efficient method. The results of Marmousi and overthrust model testings indicate that the convergence with our modified method is significantly faster than that in the L-BFGS method with no loss of inversion quality. Moreover, our modified outperforms the enriched method by a little speedup of the convergence. It also exhibits better efficiency than the HFN method.
基金Project supported by Key Industrial Projects of Major Science and Technology Projects of Zhejiang(No.2009C11023)Foundation of Zhejiang Educational Committee(No.Y200907886)Major High-Tech Industrialization Project of Jiaxing(No.2009BY10004)