AN EFFECTIVE SEQUENTIAL QUADRATIC PROGRAMMING ALGORITHM FOR NONLINEAR OPTIMIZATION PROBLEMS

In this paper,a new globally convergent algorithm for nonlinear optimization prablems with equality and inequality constraints is presented. The new algorithm is of SQP type which determines a search direction by solving a quadratic programming subproblem per itera-tion. Some revisions on the quadratic programming subproblem have been made in such a way that the associated constraint region is nonempty for each point x generated by the algorithm, i. e. , the subproblems always have optimal solutions. The new algorithm has two important properties. The computation of revision parameter for guaranteeing the consistency of quadratic sub-problem and the computation of the second order correction step for superlinear convergence use the same inverse of a matrix per iteration, so the computation amount of the new algorithm will not be increased much more than other SQP type algorithms; Another is that the new algorithm can give automatically a feasible point as a starting point for the quadratic subproblems pe     

Parareal算法的均方稳定性分析

Parareal算法是一种非常有效的实时并行计算方法.与传统的并行计算方法相比,该算法的显著特点是它的时间并行性-先将整个计算时间划分成若干个子区间,然后在每个子区间内同时进行计算.Parareal算法收敛速度快,并行效率高,且易于编程实现,从2001年由Lions,Maday和Turinici等人首次提出至今,在短短...     

A two step algorithm for solving a large scale semi-definite logit model

This paper proposes a two step algorithm for solving a large scale semi-definite logit model, which is appreciated as a powerful model in failure discriminant analysis. This problem has been successfully solved by a cutting plane (outer approximation) algorithm. However, it requires much more computation time than the corresponding linear logit model. A two step algorithm to be proposed in this paper is intended to reduce the amount of computation time by eliminating a certain portion of the data based on the information obtained by solving an associated linear logit model. It will be shown that this algorithm can generate a solution with almost the same quality as the solution obtained by solving the original large scale semi-definite model within a fraction of computation time.     

零误差计算

研究采用有误差的数值计算来获得无误差的准确值具有重要的理论价值和应用价值.这种通过近似的数值方法获得准确结果的计算被称为零误差计算.本文首先指出,只有一致离散集合中的数才能够开展零误差计算,即有非零隔离界的数集,这也是"数"可以进行零误差计算的一个充要条件.以此为基本出发点,本文分析代数数零误差计算的最低理论精度,该精度对应于恢复近似代数数的准确值时必要的误差控制条件,但由于所采用恢复算法的局限性,这一理论精度往往不能保证成功恢复出代数数的准确值.为此,本文给出采用PSLQ (partial-sum-LQ-decomposition)算法进行代数数零误差计算所需的精度控制条件,与基于LLL (Lenstra-Lenstra-Lovász)算法相比,该精度控制条件关于代数数次数的依赖程度由二次降为拟线性,从而可降低相应算法的复杂度.最后探讨零误差计算未来的发展趋势.     

The Algorithm of Sequential KKT Equations by Nonmonotone Search for Arbitrary Initial Point

In this paper, a new algorithm for arbitrary initial point is proposed. The feature of the algorithm is that only KKT equations are solved in each iteration which greatly decreases the amount of computation. Under some proper assumptions, the algorithm is globally and superlinearly convergent.     

Computing generalized inverses of a rational matrix and applications

In this paper we investigate symbolic implementation of two modifications of the Leverrier-Faddeev algorithm, which are applicable in computation of the Moore-Penrose and the Drazin inverse of rational matrices. We introduce an algorithm for computation of the Drazin inverse of rational matrices. This algorithm represents an extension of the papers [11] and [14]. and a continuation of the papers [15, 16]. The symbolic implementation of these algorithms in the package mathEmatica is developed. A few matrix equations are solved by means of the Drazin inverse and the Moore-Penrose inverse of rational matrices.     

An efficient algorithm for the computation of the metric average of two intersecting convex polygons, with application to morphing

Motivated by the method for the reconstruction of 3D objects from a set of parallel cross sections, based on the binary operation between 2D sets termed “metric average”, we developed an algorithm for the computation of the metric average between two intersecting convex polygons in 2D. For two 1D sets there is an algorithm for the computation of the metric average, with linear time in the number of intervals in the two 1D sets. The proposed algorithm has linear computation time in the number of vertices of the two polygons. As an application of this algorithm, a new technique for morphing between two convex polygons is developed. The new algorithm performs morphing in a non-intuitive way.     

The Cascade Algorithm for the Numerical Computation of Refinable Functions

The cascade algorithm is a method which can be used for the computation of refinable functions. We prove here that in general the cascade algorithm will not inherit the accuracy and refinability of the limit refinable function. We provide conditions for which the cascade iteration does indeed preserve accuracy.     

Recursive computational formulas of the least squares criterion functions for scalar system identification

The paper discusses recursive computation problems of the criterion functions of several least squares type parameter estimation methods for linear regression models, including the well-known recursive least squares (RLS) algorithm, the weighted RLS algorithm, the forgetting factor RLS algorithm and the finite-data-window RLS algorithm without or with a forgetting factor. The recursive computation formulas of the criterion functions are derived by using the recursive parameter estimation equations. The proposed recursive computation formulas can be extended to the estimation algorithms of the pseudo-linear regression models for equation error systems and output error systems. Finally, the simulation example is provided.     

Computing the Volume of n-Dimensional Copulas

Abstract

A problem that is very relevant in applications of copula functions to finance is the computation of the survival copula, which is applied to enforce multivariate put–call parity. This may be very complex for large dimensions. The problem is a special case of the more general problem of volume computation in high-dimensional copulas. We provide an algorithm for the exact computation of the volume of copula functions in cases where the copula function is computable in closed form. We apply the algorithm to the problem of computing the survival of a copula function in the pricing problem of a multivariate digital option, and we provide evidence that this is feasible for baskets of up to 20 underlying assets, with acceptable CPU time performance.     

Matching parentheses in parallel

Parallel algorithms for evaluating arithmetic expressions generally assume the computation tree form to be at hand. The computation tree form can be generated within the same resource bounds as the parenthesis matching problem can be solved. We provide a new cost optimal parallel algorithm for the latter problem, which runs in time O(log n) using O(n/log n) processors on an . We also prove that the algorithm is the fastest possible independently of the number of processors available.     

RETRACTED ARTICLE: Discrete optimization method based on grassmannian parameterization in multidimensional dichotomic data structuring

A solution of the large computational time problem arising in multidimensional data structuring is addressed by employing algebraic properties of finite geometries. A vector parameterization of the Grassmannian Gr₂(k, n) is proposed which makes it possible to minimize the amount of memory and reduce the number of operations required to solve the problem. An algorithm based on this parameterization and the Gray codes is constructed; the algorithm is suitable for parallel computation, which further reduces computation time.     

Modified screening-based Kriging method with cross validation and application to engineering design

In this paper, a basis screening Kriging method using cross validation error is proposed to alleviate computational burden of the dynamic Kriging while maintaining its accuracy. Metamodeling is widely used for design optimization of complex engineering applications where considerable computation time is required. The Kriging method is one of popular metamodeling methods due to its accuracy and efficiency. There have been many attempts to improve accuracy of the Kriging method, and the dynamic Kriging method using cross-validation error, which selects adequate basis functions to best describe the mean structure of a response using a genetic algorithm, achieves outstanding performance in terms of accuracy. However, despite its accuracy, the dynamic Kriging requires very large amounts of computation because of the genetic algorithm and no limitation for order of basis functions. In the proposed method, a basis function set is determined by screening each basis function instead of using the genetic algorithm, which has advantages in computation for high dimensional metamodels or repeated metamodel generation. Numerical studies with four mathematical examples and two engineering applications verify that the proposed basis screening Kriging significantly reduces computation time with similar accuracy as the dynamic Kriging.     

连续代数扩域上多项式因式分解的Trager算法

多项式的因式分解是符号计算中最基本的算法,二十世纪六十年代开始出现的关于多项式因式分解的工作被认为是符号计算领域的起源．目前多项式的因式分解已经成熟,并已在Maple等符号计算软件中实现,但代数扩域上的因式分解算法还有待进一步改进．代数扩域上的基本算法是Trager算法．Weinberger等提出了基于Hensel提升的算法．这些算法是在单个扩域上做因式分解．而在吴零点分解定理中,多个代数扩域上的因式分解是非常基本的一步,主要用于不可约升列的计算．为了解决这一问题,吴文俊,胡森、王东明分别提出了基于方程求解的多个扩域上的因式分解算法．王东明、林东岱提出了另外一个算法Trager算法相似,将问题化为有理数域上的分解．他们应用了吴的三角化算法,因此算法的终止性依赖于吴方法的计算．支丽红则将提升技巧用于多个扩域上的因式分解算法．本文将Trager的算法直接推广为连续扩域上的因式分解,只涉及结式计算与有理数域上的因式分解,给出了多个代数扩域上的因式分解一个直接的算法．     

The arithmetic fourier transform and the visual pathway

The Arithmetic Fourier Transform is an algorithm for the computation of Fourier coefficients. In this paper it is extended to the computation of double Fourier coefficients by a repeated sum. We discuss the relevance of the repeated sum algorithm to signal processing by neurons in the visual pathway.     

A bisection method for measuring the distance of a quadratic matrix polynomial to the quadratic matrix polynomials that are singular on the unit circle

The computation of the distance of a quadratic matrix polynomial to the quadratic matrix polynomials that are singular on the unit circle is investigated. The emphasis is placed on backward stable methods that transform the computation of the distance to a palindromic eigenvalue problem for which structure-preserving eigensolvers can be utilized in conjunction with a bisection algorithm. Reliability of the suggested methods is guaranteed by a novel error analysis.     

Clique analysis of a tolerance relation

The Harary‐Ross clique detection algorithm has often been neglected in analyses of symmetric graphs and tolerance relations because of presumed storage and computation problems. This paper presents a modification of the original algorithm which overcomes these difficulties. The modified algorithm is shown to perform as well as other recent structurally different algorithms. Indeed, for arc densities which might be expected in typical sociometric analyses, the modified algorithm has certain conceptual and operational advantages. The paper also indicates possible applications of the algorithm in the analysis of non‐symmetric functional relations and points out the conceptual similarity with simplicial connectivity analysis.     

A truncated Newton method in an augmented Lagrangian framework for nonlinear programming

In this paper we propose a primal-dual algorithm for the solution of general nonlinear programming problems. The core of the method is a local algorithm which relies on a truncated procedure for the computation of a search direction, and is thus suitable for large scale problems. The truncated direction produces a sequence of points which locally converges to a KKT pair with superlinear convergence rate.     

Complexity of linear programming

The complexity of linear programming is discussed in the “integer” and “real number” models of computation. Even though the integer model is widely used in theoretical computer science, the real number model is more useful for estimating an algorithm's running time in actual computation.Although the ellipsoid algorithm is a polynomial-time algorithm in the integer model, we prove that it has unbounded complexity in the real number model. We conjecture that there exists no polynomial-time algorithm for the linear inequalities problem in the real number model. We also conjecture that linear inequalities are strictly harder than linear equalities in all “reasonable” models of computation.