Search

Our underlying motivation is the iterative method for the implementation of the limit representation of the Moore-Penrose inverse lim αI0 (αI + A*A)-1 A* from [Žukovski, Lipcer, On recurent computation of normal solutions of linear algebraic equations, Ž. Vicisl. Mat. i Mat. Fiz. 12 (1972), 843-857] and [Žukovski, Lipcer, On computation pseudoinverse matrices, Ž. Vicisl. Mat. i Mat. Fiz. 15 (1975), 489-492]. The iterative process for the implementation of the general limit formula lim αI0 (αI + R*S)-1R* was defined in [P.S. Stanimirović, Limit representations of gen- eralized inverses and related methods, Appl. Math. Comput. 103 (1999), 51-68]. In this paper we develop an improvement of this iterative process. The iterative method defined in such a way is able to produce the result in a predefined number of iterative steps. Convergence properties of defined iterations are further investigated. [ABSTRACT FROM AUTHOR]

Abstract: The hybrid algorithm that combined particle swarm optimization with simulated annealing behavior (SA-PSO) is proposed in this paper. The SA-PSO algorithm takes both of the advantages of good solution quality in simulated annealing and fast searching ability in particle swarm optimization. As stochastic optimization algorithms are sensitive to their parameters, proper procedure for parameters selection is introduced in this paper to improve solution quality. To verify the usability and effectiveness of the proposed algorithm, simulations are performed using 20 different mathematical optimization functions with different dimensions. The comparative works have also been conducted among different algorithms under the criteria of quality of the solution, the efficiency of searching for the solution and the convergence characteristics. According to the results, the SA-PSO could have higher efficiency, better quality and faster convergence speed than compared algorithms. [Copyright &y& Elsevier]

As a solution to the frequently reported problem of the high frequencies reconstructed sometimes many iterations after the low frequencies in tomography, we proposed in a recent paper [V. Israel-Jost, Ph. Choquet, and A. Constantinesco, Int. J. Biomed. Imaging, 2006, paper 34043.] a general adaptation that can apply to many algorithms: the incomplete backprojection. We stressed the fact that this type of frequency adapted (FA) algorithm works only when the linear problem of tomography to be solved involves deblurring, i.e., when the point spread functions (PSFs) associated with each voxel are large. Thus, depending on the desired level of detail, a parameter is set to achieve the reconstruction within a small number of iterations. The aim of this paper is to give a mathematical analysis of both the acceleration obtained and the convergence of an FA algorithm derived from the simultaneous algebraic reconstruction technique (SART) algorithm of Andersen and Kak, the so-called FA-SART algorithm. [ABSTRACT FROM AUTHOR]

In this paper, a finite filter is used in the structure of the Barzilai–Browein (BB) gradient method in order to propose a new modified BB algorithm for solving large-scale unconstrained optimization problems. Our algorithm is equipped with a relaxed nonmonotone line search technique which allows the algorithm to enjoy the nonmonotonicity properties from scratch. Under some suitable conditions, the global convergence property of the new proposed algorithm is established. Numerical results on some test problems in CUTEr library show the efficiency and effectiveness of the new algorithm in practice too. [ABSTRACT FROM PUBLISHER]

Abstract: In this paper, we discuss the convergence of the double-shift and multi-shift QR algorithms for symmetric tridiagonal matrices. We analyze how to choose multi-shifts by comparing the relationships between the number of iterations, CPU time and the number of multi-shifts. Numerical tests and figures are performed. [Copyright &y& Elsevier]

This paper first proposes a trust region algorithm to obtain a stationary point of unconstrained multiobjective optimisation problem. Under suitable assumptions, the global convergence of the new algorithm is established. We then extend the trust region method to solve the non-smooth multiobjective optimisation problem. [ABSTRACT FROM PUBLISHER]

Abstract: In this paper, we consider numerical simulation to a class of constrained dynamic problems where the overall dynamics are determined by the interactions between two sub-systems. We present an iterative algorithm that naturally decouples the computation of the two sub-systems and that ensures an accurate and efficient solution procedure. We conduct rigorous error analysis for the convergence of the iterative algorithm, and verify the analytical results through careful numerical tests. [Copyright &y& Elsevier]

Abstract: In this paper, we propose a new penalty-free-type method for nonlinear equality constrained problems. The new algorithm uses trust region framework and feasibility safeguarding technique. Moreover, it has no choice of penalty parameter and penalty function as a merit function, and it does not use the filter technique to avoid the penalty function either. We analyze the global convergence of the main algorithm under the standard assumptions. The preliminary numerical tests are reported. [Copyright &y& Elsevier]

Recently, Xiao et al. proposed a nonsmooth equations-based method to solve the l₁-norm minimization problem (2011). The advantage of this method is its simplicity and lower storage. In this paper, based on new nonsmooth equations reformulation, we investigate new nonsmooth equations-based algorithms for solving l₁-norm minimization problems. Under mild conditions, we show that the proposed algorithms are globally convergent. The preliminary numerical results demonstrate the effectiveness of the proposed algorithms. [ABSTRACT FROM AUTHOR]

Abstract: In this paper it is presented a compensated de Casteljau algorithm to accurately evaluate a bivariate polynomial in Bernstein–Bézier form. The principle is to apply error-free transformations to improve the traditional de Casteljau algorithm. A forward error and a running error analysis are performed. Finally, some numerical experiments illustrate the accuracy of the proposed algorithm in ill-conditioned problems. [Copyright &y& Elsevier]

Abstract: This paper proposes a novel hybrid immune algorithm (HIA) that can overcome the typical drawback of the artificial immune algorithm (AIA), which runs slowly and experiences slow convergence. The HIA combines the adaptive AIA based on the steepest descent algorithm. The HIA fully displays global search ability and the global convergence of the immune algorithm. At the same time, it inserts a quasi-descent operator to strengthen its local search ability. A good convergence of the HIA with the quasi-descent idea is shown as well. Numerical experiment results show that the HIA successfully improves running speed and convergence performance. [Copyright &y& Elsevier]

Abstract: In this paper, an active set limited BFGS algorithm is proposed for bound constrained optimization. The global convergence will be established under some suitable conditions. Numerical results show that the given method is effective. [Copyright &y& Elsevier]

Abstract: In this paper, the issue of multi-degree reduction of Bézier curves with C 1 and G 2-continuity at the end points of the curve is considered. An iterative method, which is the first of this type, is derived. It is shown that this algorithm converges and can be applied iteratively to get the required accuracy. Some examples and figures are given to demonstrate the efficiency of this method. [Copyright &y& Elsevier]

Abstract: In this paper, we propose an interval version of the generalized accelerated overrelaxation methods, which we refer to as IGAOR, for solving the linear complementarity problems, LCP (M, q), and develop a class of multisplitting IGAOR methods which can be easily implemented in parallel. In addition, in regards to the H-matrix with positive diagonal elements, we prove the convergence of these algorithms and illustrate their efficiency through our numerical results. [Copyright &y& Elsevier]

Abstract: The aim of this paper is to incorporate the preconditioned gradient path in a nonmonotone stabilization algorithm for unconstrained optimization. The global convergence and locally superlinear convergence are established for this class of algorithms. Finally, we report in details the numerical results which show the effectiveness of the proposed algorithm. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we use the auxiliary principle technique in conjunction with the Bregman function to suggest and analyze a predictor–corrector method for solving invex equilibrium problems. We also study the convergence criteria of this new method under some mild conditions. We obtain, as special cases, new and known methods for solving variational-like inequalities and equilibrium problems. [Copyright &y& Elsevier]

Abstract: This paper presents an affine scaling optimal path approach in association with nonmonotonic interior backtracking line search technique for nonlinear optimization subject to linear constraints. We shall employ the eigensystem decomposition and affine scaling mapping to form affine scaling optimal curvilinear path very easily. By using interior backtracking line search technique, each iterate switches to trial step of strict interior feasibility. The nonmonotone criterion is used to speed up the convergence progress in the contours of objective function with large curvature. Theoretical analysis is given which prove that the proposed algorithm is globally convergent and has a local superlinear convergence rate under some reasonable conditions. The results of numerical experiments are reported to show the effectiveness of the proposed algorithm. [Copyright &y& Elsevier]

In this paper we classify several algorithmic problems in group theory in the complexity classes PZK and SZK (problems with perfect/statistical zero-knowledge proofs respectively). Prior to this, these problems were known to be in $\mbox {\rm AM}\cap \mbox {\rm coAM}$ . As $\mbox {\rm PZK}\subseteq \mbox {\rm SZK}\subseteq \mbox {\rm AM}\cap \mbox {\rm coAM}$ , we have a tighter upper bound for these problems. Specifically: [ABSTRACT FROM AUTHOR]

The rates of convergence of plug-in kernel, partitioning, and nearest neighbors classification rules are analyzed. A margin condition, which measures how quickly the a posteriori probabilities cross the decision boundary, smoothness conditions on the a posteriori probabilities, and boundedness of the feature vector are imposed. The rates of convergence of the plug-in classifiers shown in this paper are faster than previously known. [ABSTRACT FROM AUTHOR]

Abstract: In this paper extended algorithm of Karmarkar’s interior point algorithm is considered. Theoretical aspects of the new algorithm are discussed. We show that the new iterative algorithm converges faster than former Karmarkar’s algorithms. Successful convergence for problems of different sizes is obtained. Numerical results show that when we use a new parameter in the classical Karmarkar’s algorithm, number of iterations is less than the number of iterations for two mentioned methods. [Copyright &y& Elsevier]

In this paper, we reinvestigate the trust-region method by reformulating its subproblem: the trust-region radius is guided by gradient information at the current iteration and is self-adaptively adjusted. A trust-region algorithm based on the proposed subproblem is proved to be globally convergent. Moreover, the superlinear convergence of the new algorithm is shown without the condition that the Hessian of the objective function at the solution be positive definite. Preliminary numerical results indicate that the performance of the new method is notable. [ABSTRACT FROM AUTHOR]

In this paper we propose an improved quasi-Monte Carlo method for solving Linear Algebra problems. We show that by using low-discrepancy sequences both the convergence and the CPU time of the algorithm are improved. Two parallelization schemes using the Message Passing Interface with static and dynamic load balancing are proposed. The dynamic scheme is useful for computing in the GRID environment. [ABSTRACT FROM AUTHOR]

For the problem of minimizing the sum of Euclidean norms (MSN), most existing quadratically convergent algorithms require a strict complementarity assumption. However, this assumption is not satisfied for a number of MSN problems. In this paper, we present a globally and quadratically convergent algorithm for the MSN problem. In particular, the quadratic convergence result is obtained without assuming strict complementarity. Examples without strictly complementary solutions are given to show that our algorithm can indeed achieve quadratic convergence. Preliminary numerical results are reported. [ABSTRACT FROM AUTHOR]

This paper presents a new trust-region algorithm for n-dimension nonlinear optimization subject to m nonlinear inequality constraints. Equivalent KKT conditions are derived, which is the basis for constructing the new algorithm. Global convergence of the algorithm to a first-order KKT point is established under mild conditions on the trial steps, local quadratic convergence theorem is proved for nondegenerate minimizer point. Numerical experiment is presented to show the effectiveness of our approach. [ABSTRACT FROM AUTHOR]

In this paper we analyse the Kaczmarz projection algorithm with remotest set control of projection indices. According to this procedure, at each iteration the projection index is one which gives the maximal absolute value of the corresponding residual. We prove that for under-determined full row rank systems and under some assumptions valid for problems arising in algebraic reconstruction of images in computerized tomography, this selection procedure has the property that each row index is selected at least once during the Kaczmarz algorithm iterations. [ABSTRACT FROM AUTHOR]