A numerical implementation of a predictor-corrector algorithm for sufcient linear complementarity problem
|
|
- Bathsheba Nichols
- 5 years ago
- Views:
Transcription
1 A numerical imlementation of a redictor-corrector algorithm for sufcient linear comlementarity roblem BENTERKI DJAMEL University Ferhat Abbas of Setif-1 Faculty of science Laboratory of fundamental and numerical Mathematics ALGERIA dj benterki@yahoo.fr BOULOUDENINE NADJIBA University Ferhat Abbas of Setif-1 Faculty of science Deartment of Mathematics ALGERIA nadjiba010@hotmail.fr Abstract: In this aer, we imlement a version of a Mizuno-Todd-Ye redictor-corrector interior oint algorithm for the P*k)-matrix linear comlementarity roblem LCP) resented by T. Illes and M. Nagy. We make some numerical develoments to calculate the directions and the ste length in the redictor case. Our results shows that the roosed algorithm achieve high, accuracy when it alies to monotone linear comlementarity roblem. Key Words: Linear comlementarity roblem, P*k)-matrix, Predictor-corrector methods., 1 Introduction We consider the following linear comlementarity roblem LCP ) : nd vectors x; y R n which veries : 8 < LCP ) : Mx + y x; y 0 xy = 0: = q Where M R nn and q R n. Comlementarity roblems are known in 1940, they have many alications in mathematical rogramming and equilibrium roblems. Indeed, it is known that by exloiting the rst-order otimality conditions of the otimization roblem, any differentiable convex quadratic rogram can be formulated into a monotone LCP ), and vice versa [16]. Variational inequality roblems are widely used in the study of equilibriums in economics, transortation lanning, and game theory. They also have a close connection to the LCP )s. The reader can refer to [] for the basic theory. Although it's earlier aearance, it was not until 1960 that researchers interested in these issues. This certainly comes to the lack of numerical tools for solving such roblem. Since then, several algorithms are roosed in the literature, esecially after the aearance of the famous method of Karmarkar in The interior oint methods are more suitable for solving the LCP ) roblem, in articular, the central ath methods have become the methods of choice thanks to their good roerties rimal-dual methods, olynomial convergence). Many alternatives are available in the literature, the work of Kojima, Mizuno and Yoshise [9], Achache [1] are cited to solve a linear comlementarity roblem LCP ) monotone M is ositive semidenite ). Desite the remarkable success of these methods, the central ath methods suffer esecially from the initial oint. Much more, in general, they don't have a feasible initial oint. Worse, even if this is known, the convergence of the central ath methods is not guaranteed unless this oint is close to the central ath. To address this issue, different rocedures are roosed. This is the method of weight with central ath [3, 6, 14], which serves to return the iterated neighboring to the central ath. Moreover, these methods are not achievable the algorithm starts from any oint x > 0 ) and rove theoretically converge. Unfortunately, these statements are not erformed numerically. Moreover, such methods did not matter thereafter. Currently, the redictor-corrector methods are imortant to answer this question. The rst version of the redictor-corrector algorithm was initiated by Sonnevend, Stor and Zhao [14] for solving a linear rogramming roblems. This algorithm needs more corrector stes after each redictor ste in order to return to the aroriate neighborhood of the central ath. Mizuno, Todd and Ye [11] have roosed a redictor-corrector algorithm for solving a linear rograming in a olynomial time which each redictor ste is followed by a single corrector ste. This algo- ISBN:
2 rithm has oened the way for other more comlicated alication to solve linear comlementarity roblems. Indeed Anstricher and Ye [17] extended this result to the linear comlementarity roblem with a ositive semidenite matrix with the same iteration comlexity. Potra have roosed a Mizuno-Todd-Ye algorithm in a larger neighborhood of the central ath [1]. In 1991, Kojima and al have ublished their book on interior oint method for LCP )s [8]. Since then the quality of a variant of an interior oint algorithm is measured by the fact whether it can be generalized to the P k)-matrix linear comlementarity roblem or not. Several variants of Muzino -Todd-Ye tye redictor-corrector interior oint algorithm are known in the literature. In 007, Illes and Nagy [5] have roosed a redictor-corrector algorithm tye Mizuno - Todd -Ye for sufcient linear comlementarity roblems. Our work focuses on the comrehensive theoretical and numerical study of this algorithm [5], where we have made theoretical and numerical develoments in order to make the algorithm The work is divided into three sections organized as follows : In the rst section, we mention some matrix classes denitions, and some results. The second section contains a summary of Illes and Nagy's [5] algorithm, followed by the descrition of the algorithm itself, and the last section resent our numerical exeriments of the algorithm and a conclusion. We denote by: e = 1; 1; :::; 1) R n : xy = x 1 y 1 ; x y ; :::; x n y n ) R n Hadamard roduct of vectors x and y. x = x 1 ; x ; :::; x n ) R n. x y = x 1 y 1 ; x y ; :::; xn y n ) R n : X = diagx 1 ; x ; ::; x n ): F = x; y) R n + : y = Mx + q ; the feasible set of the LCP ) roblem. F + = x; y) R n ++ : y = Mx + q ; the strict feasible set of the LCP ) roblem. F = x; y) R n + : y = Mx + q; xy = 0, the solution set of the LCP ) roblem. T = fx; y) F + : xy = e; > 0g ; the central ath. The P k)-matrices and their roerties Denition 1 [8] Let k R ++. A matrix M R nn is called P k)-matrix if for all x R n we have: 1 + 4k) X x i Mx) i + X x i Mx) i 0; ii + x) ii where I + x) = f1 i n : x i Mx) i > 0 g and I x) = f1 i n : x i Mx) i < 0 g : - If k = 0; then M is called a semidenite matrix M is a P 0)-matrix). - Lets k 1 and k tow ositif real numbers. If k 1 k then P k 1 ) P k ): Denition A matrix M R nn is called a P - matrix if it is a P k)-matrix for some k > 0, i.e., P = [P k) k0 Denition 3 A matrix M R nn is called a column sufcient matrix if: 8x R n ; x i Mx) i 0 =) x i Mx) i = 0; 8i = 1; :::; n: Denition 4 A matrix M R nn is called row sufcient if M T is a column sufcient matrix. Denition 5 A matrix M R nn is called sufcient if it is both row and column sufcient. P matrices are sufcient. Kojima and al [8] roved that a P k) matrix is column sufcient and Guu and Cottle [4] roved that it is row sufcient, too. Valiaho [15] roved that each sufcient matrix is a P matrix, so the class of P -matrices is equal to the class of sufcient matrices. Proosition 6 [8] If M R nn is a P k)-matrix then M 0 M I = Y X is a nonsingular matrix for any ositive diagonal matrices X; Y R nn. Corollary 7 Let M R nn be a P k) matrix, x; y R n + then for all a R n the system x) M4x + 4y = 0 Y 4x + X4y = e xy has a unique solution 4x; 4y): ISBN:
3 Theorem 8 [8] Let a linear comlementarity roblem with a P k)-matrix M be given. Then the following statements are equivalent: 1 F + 6=?. 8w R n +; 9!x; y) F + : xy = w: 3 8 > 0; 9!x; y) F + : xy = e: Throughout the aer we make the following assumtions: 1. F + 6=?:. An initial oint x 0 ; y 0 ) F + exists. 3. M is a P k)-matrix. 3 Descrition of the method After the relaxation of the LCP ) roblem, we obtain the following central ath roblem: 8 < Mx + y = q LCP ) x; y > 0 : xy = e When! 0; the sequence x); y)) solutions of the central ath LCP ) roblem aroach the solution x; y) of the LCP ) roblem. By using the Newton method to solve the LCP ) roblem, we obtain the following system: M4x + 4y = 0 1) Y 4x + X4y = e xy Let the roximity measure be x; y; ) = v) = v 1 v ; q where v = xy. For the oint x; y) of the central ath, the roximity must be zero, i.e., xy = e for some > 0. The roximity measure aroaches the innity if some coordinates of x or y tends to 0 or innity. Let 0 < < 0 be the suitable roximity arameters dening the neighborhood of the central ath. In redictor ste, from a given initial oint x 0 ; y 0 ; 0 ) F + and satises x 0 ; y 0 ; 0 ) <, we solve the Newton system with = 0 and we denote the solution by 4x; 4y). The new oint will be: x = x 0 +4x; y = y 0 +4y; = 1 ) 0 : Where 0; 1] is the maximal real number for which x ; y ; ) F + and x ; y ; ) 0. In the corrector ste, we return from the 0 neighborhood to the neighborhood of the central ath. Starting from the oint x ; y ; ); we solve the Newton system and we get 4ex; 4ey). Let the new oint be x c = x + 4ex; y c = y + 4ey; c = : The oint x c ; y c ; c ) will be in the neighborhood again, so we can continue the iteration from this oint as an initial oint x 0 ; y 0 ) while the duality ga is not small enough. Algorithm Inut An accuracy arameter " > 0; Proximity arameters ; 0 ; 0 < < 0 ; An initial oint x 0 ; y 0 ; 0 ); such that x 0 ; y 0 ; 0 ) ; Begin x = x 0 ; y = y 0 ; = 0 ; While x T y " do Predictor ste Determine 4x; 4y) by solving 1) with = 0 and let x = x 0 +4x; y = y 0 +4y; = 1 ); Determine : max f > 0 : x ; y ; ) F + and x ; y ; ) 0 g Corrector ste Determine 4ex; 4ey) by solving 1) with x = x ; y = y ; = and let x c = x +4ex; y c = y +4ey; c = ; x = x c ; y = y c ; = c ; End while End. 4 The algorithm iteration In this section, we describe the comutation of the redictor and corrector directions, and the ste-length = us it was q resented in q [5]: Let v = xy ; d = x y ; d 1 = e d ; d x = d 1 4x = v4x x ; d y = d 1 4y = v4y y. After rescaling the Newton-system we have Mdx + d y = 0 d x + d y = v 1 v; where M = DMD; and D = diagd i ): ISBN:
4 4.1 The redictor ste Starting from a oint x; y) F + where x; y; ) <, we solve the Newton system with arameter = 0; that means to solve the following system: Mdx + d y = 0 d x + d y = v; which has the following solution: d x = I + M) 1 v; and d y = MI + M) 1 v: We choose the ste-length so that the following two conditions are satised: 1. x ; y ) F +, i.e., x y > 0,. x ; y ; ) < 0, where x = x + 4x; y = y + 4y; and = 1 ); i.e., the ositivity condition is satised at the new oint and we do not deviate from the central ath too much. We have: and x y = 1 )v + d x d y ; v ) = x y 1 ) = v + 1 d xd y : For the rst condition, x y > 0 imlies v + 1 d xd y > 0: Here we have tow cases: Case 1: If d x d y ) i 0; then obviously v ) i > 0 is true for all 0:1): Case : If d x d y ) i < 0; then the condition 1 < vi d xd y) i must hold for the ste-length. Let v i b' = min : d x d y ) i < 0 d x d y ) i and ' = 1 ; then ' [0; b') is the necessary condition so that x ; y ) F + is satised. To examine the second condition about choosing the ste-length, we comute the roximity measure of the vector v so we have: v ) = v ) 1 v = v) + e T e v + 'd x d y e v + 'd xd y ): Let f') : [0; b')! R, therefore f') = v ) v) = e T e v + 'd x d y e v + 'd xd y ): The function f is strictly convex on the interval [0; b') and veries, f0) = 0 and lim '!b' f') = +1: From the roerties of function f it follows that the equation f') = 0 v) has a unique solution ' on interval [0; b'): We can ut ' = ) 1 ; ) then the ste length can be comuted as = ' + ' ) + 4' : Remark 9 To comute the ste length ), we can use a numerical method such as for examle Newton method. Proosition 10 [5] is the maximal feasible ste length such that v ) 0 and x ; y ) F + : 4. The corrector ste In this ste, our aim is to return from the 0 - neighborhood of the central ath. To obtain the corrector directions g 4x; f 4y); we solve the following Newton-system: M f d x + e d y = 0 fd x + e d y = v ) 1 v where fd x = v g 4x x ; e dy = v f 4y y ; e d = Let r x y ; f M = e DM e D: x c = x + g 4x; y c = y + f 4y; c = = 1 ): The new oint x c ; y c ; c ) must be feasible witch means the following tow conditions are satised: 1. x c ; y c ) F +,. x c ; y c ; c ) <. ISBN:
5 For that, we will choose a suitable value of and 0 where the full Newton ste is feasible. r q Proosition 11 [5]If < , 1+4k) 1+4k 4 + then there exists a 0 that satises 0 < and the full Newton ste is feasible at corrector ste. 4.3 Iteration comlexity Theorem 1 [5] Let the linear comlementarity roblem for any P k)-matrix M be given, where k 0 and let 0 = 1. Then the Mizuno-Todd-Ye algorithm generates an x; y; ) oint satisfying x T y < " in at most [ Where: n 1 log 4n+ 0 4" ] iterations. 1 = ) 0 = m)1 + 4k) ) = min 0 3 ; h) 8m) ; r m) = h) = 8 m ) m ) ) Corollary 13 [5] Let the linear comlementarity roblem for any P k)-matrix M be given, where k 0 and let 0 = 1; = 1 1+4k and 0 = 1+4k. Then the Mizuno-Todd-Ye algorithm generates an x; y; ) oint satisfying x T y < " in at most O1 + k) 3 n log n " ) iterations. 5 Numerical exeriments In this section, we resent some numerical exeriments for a sufcient linear comlementarity roblems. Unfortunately, because the comutation of the arameter k is difcult those exeriments are limited on the monotone linear comlementarity roblems where the matrix M is a ositive semidenite. Our imlementation was tested on a Pentium 3 and was rogrammed in Dev-Pascal on a number of roblems taken from the literature and reorted in [1, 7]. In all the following roblems, we take " = 10 6, = 1:, 0 = 1:76, and we denote by exn) the taille of different examles and Iter reresents the number of iterations necessary to obtain the otimal solution. Examle 14 with xed size n) [1, 7] examle exn) iter ex) 4 ex3) 4 ex4) 5 ex5) 4 ex10) 7 Examle 15 with variable size n) [7] 0 1 ::: 0 1 ::: ::: M = : : : : : B : : : ::: : : : : : ::: 1 q = 1; :::; 1 ) T ; x 0 = 0:05; :::; 0:05; 1:05 ) T ; Conclusion x = 0 ; :::; 0 ; 1 ) T : Size n Iter ; C A The theoretical study on redictor-corrector algorithms for sufcient linear comlementarity roblems has allowed us to imrove the behavior of this tye of algorithm acting on the best values of and 0 rediction and correction neighborhood), numerical imlementation conrmed that. Our imlementation suose treated sufcient linear comlementarity roblem, because the calculation of the arameter k is difcult. The imlementation was limited on monotone linear comlementarity roblem. Our testes roved that the redictor-corrector algorithm worked very satisfactorily. References [1] M. Achache, Comlexity analysis and numerical imlementation of a short-ste rimal-dual algorithm for linear comlementarity roblems, Alied M athematics and Comutation 16, 010, ISBN:
6 [] R.W. Cottle, J.S. Pang, R.E. Stone, The linear comlementarity roblem, Academic P ress Inc, San Diego, CA, 199. [3] J. Ding and T. Y. Li, An algorithm based on weighted logarithmic barrier function for linear comlementarity roblem, T he Arabian Journal f or Science and Engineering 15, N4B [4] S. M. Guu, R. W. Cottle, On a subclass of P 0 linear algebra and its alications 3/4; [5] T. Illés, M. Nagy, A Mizuno-Todd-Ye tye redictor-corrector algorithm for sufcient linear comlementarity roblems, Euroean Journal of Oerational Research 181, 007, [6] Z. Kebbiche, Etude et extension d'algorithmes de oints intérieurs our la rogrammation non linéaire, T hese de doctorat; U niversite F erhat Abbas de Setif; 007. [13] E. M. Simantiraki et D. F. Shanno, An infeasible-interior oint method for linear comlementarity roblems, SIAM: J: otim:.7, 1997, : [14] Gy. Sonnevend, J. Stoer, G. Zhao, On the comlexity of following the central ath of linear rograms by linear extraolation, M ethods of Oerations Research; 63, 1989, [15] H. Valiaho, P k) matrices are just sufcient, Linear Algebra and its Alications; 39, 1996, [16] S. J. Wright, Primal-dual interior oint methods, SIAM:P hiladelhia [17] Y. Ye, K. Anstreicher, On quadratic and O nl) convergence of a redictor-corrector algorithm for LCP ), M athematical P rogramming 6, 1993, [7] R. Khebchache, Etude théorique et numérique de quelques méthodes de oints intérieurs de trajectoire avec oids our le comlémentarité linéaire monotone, T hese de magisterse; Universite F erhat Abbas de Setif 010. [8] M. Kojima, N. Megiddo, T. Noma and A. Yoshise, A unied aroach to interior oint algorithms for linear comlementarity roblems, Lecture N otes in Comuter Science; vol: 538; Sringer V erlag; Berlin; Germany; [9] M. Kojima, S. Mizuno and A. Yoshise, A olynomial time algorithm for class of linear comlementarity roblems, M athematical P rogramming; 44, 1989,. 1-6: [10] J. Miao, A quadratically convergent Ok + 1) nl) -iteration algorithm for the P k) matrix linear comlementarity roblem, Mathematical Programming ) [11] S. Mizuno, M. J. Todd, Y. Ye, On adative-ste rimal-dual interior-oint algorithms for linear rogramming, M athematics of Oerations Research; 18,4, 1993, [1] F. A. Potra, The Mizuno-Todd-Ye algorithm in a larger neighborhood of the central ath, Euroean Journal of Oerational Research; , ISBN:
Tibor Illes, Jiming Peng 2, Kees Roos, Tamas Terlaky Faculty of Information Technology and Systems Subfaculty of Technical Mathematics and Informatics
DELFT UNIVERSITY OF TECHNOLOGY REPORT 98{5 A Strongly Polynomial Rounding Procedure Yielding A Maximally Comlementary Solution for P () Linear Comlementarity Problems T. Illes, J. Peng, C. Roos, T. Terlaky
More informationAn Infeasible Interior Point Method for the Monotone Linear Complementarity Problem
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 17, 841-849 An Infeasible Interior Point Method for the Monotone Linear Complementarity Problem Z. Kebbiche 1 and A. Keraghel Department of Mathematics,
More informationSome results of convex programming complexity
2012c12 $ Ê Æ Æ 116ò 14Ï Dec., 2012 Oerations Research Transactions Vol.16 No.4 Some results of convex rogramming comlexity LOU Ye 1,2 GAO Yuetian 1 Abstract Recently a number of aers were written that
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
More informationA full-newton step feasible interior-point algorithm for P (κ)-lcp based on a new search direction
Croatian Operational Research Review 77 CRORR 706), 77 90 A full-newton step feasible interior-point algorithm for P κ)-lcp based on a new search direction Behrouz Kheirfam, and Masoumeh Haghighi Department
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationA Generalized Homogeneous and Self-Dual Algorithm. for Linear Programming. February 1994 (revised December 1994)
A Generalized Homogeneous and Self-Dual Algorithm for Linear Programming Xiaojie Xu Yinyu Ye y February 994 (revised December 994) Abstract: A generalized homogeneous and self-dual (HSD) infeasible-interior-point
More informationAn EP theorem for dual linear complementarity problems
An EP theorem for dual linear complementarity problems Tibor Illés, Marianna Nagy and Tamás Terlaky Abstract The linear complementarity problem (LCP ) belongs to the class of NP-complete problems. Therefore
More informationEfficient algorithms for the smallest enclosing ball problem
Efficient algorithms for the smallest enclosing ball roblem Guanglu Zhou, Kim-Chuan Toh, Jie Sun November 27, 2002; Revised August 4, 2003 Abstract. Consider the roblem of comuting the smallest enclosing
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More informationOn Isoperimetric Functions of Probability Measures Having Log-Concave Densities with Respect to the Standard Normal Law
On Isoerimetric Functions of Probability Measures Having Log-Concave Densities with Resect to the Standard Normal Law Sergey G. Bobkov Abstract Isoerimetric inequalities are discussed for one-dimensional
More informationA PREDICTOR-CORRECTOR PATH-FOLLOWING ALGORITHM FOR SYMMETRIC OPTIMIZATION BASED ON DARVAY'S TECHNIQUE
Yugoslav Journal of Operations Research 24 (2014) Number 1, 35-51 DOI: 10.2298/YJOR120904016K A PREDICTOR-CORRECTOR PATH-FOLLOWING ALGORITHM FOR SYMMETRIC OPTIMIZATION BASED ON DARVAY'S TECHNIQUE BEHROUZ
More informationNONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS. The Goldstein-Levitin-Polyak algorithm
- (23) NLP - NONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS The Goldstein-Levitin-Polya algorithm We consider an algorithm for solving the otimization roblem under convex constraints. Although the convexity
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More information1 1 c (a) 1 (b) 1 Figure 1: (a) First ath followed by salesman in the stris method. (b) Alternative ath. 4. D = distance travelled closing the loo. Th
18.415/6.854 Advanced Algorithms ovember 7, 1996 Euclidean TSP (art I) Lecturer: Michel X. Goemans MIT These notes are based on scribe notes by Marios Paaefthymiou and Mike Klugerman. 1 Euclidean TSP Consider
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationLecture 10. Primal-Dual Interior Point Method for LP
IE 8534 1 Lecture 10. Primal-Dual Interior Point Method for LP IE 8534 2 Consider a linear program (P ) minimize c T x subject to Ax = b x 0 and its dual (D) maximize b T y subject to A T y + s = c s 0.
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationOn the Chvatál-Complexity of Knapsack Problems
R u t c o r Research R e o r t On the Chvatál-Comlexity of Knasack Problems Gergely Kovács a Béla Vizvári b RRR 5-08, October 008 RUTCOR Rutgers Center for Oerations Research Rutgers University 640 Bartholomew
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationEnlarging neighborhoods of interior-point algorithms for linear programming via least values of proximity measure functions
Enlarging neighborhoods of interior-point algorithms for linear programming via least values of proximity measure functions Y B Zhao Abstract It is well known that a wide-neighborhood interior-point algorithm
More informationApplicable Analysis and Discrete Mathematics available online at HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS
Alicable Analysis and Discrete Mathematics available online at htt://efmath.etf.rs Al. Anal. Discrete Math. 4 (010), 3 44. doi:10.98/aadm1000009m HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS Zerzaihi
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationA full-newton step infeasible interior-point algorithm for linear complementarity problems based on a kernel function
Algorithmic Operations Research Vol7 03) 03 0 A full-newton step infeasible interior-point algorithm for linear complementarity problems based on a kernel function B Kheirfam a a Department of Mathematics,
More informationFinding Shortest Hamiltonian Path is in P. Abstract
Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationarxiv: v2 [math.na] 6 Apr 2016
Existence and otimality of strong stability reserving linear multiste methods: a duality-based aroach arxiv:504.03930v [math.na] 6 Ar 06 Adrián Németh January 9, 08 Abstract David I. Ketcheson We rove
More informationAn Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization
An Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization H. Mansouri M. Zangiabadi Y. Bai C. Roos Department of Mathematical Science, Shahrekord University, P.O. Box 115, Shahrekord,
More informationLEIBNIZ SEMINORMS IN PROBABILITY SPACES
LEIBNIZ SEMINORMS IN PROBABILITY SPACES ÁDÁM BESENYEI AND ZOLTÁN LÉKA Abstract. In this aer we study the (strong) Leibniz roerty of centered moments of bounded random variables. We shall answer a question
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationA FULL-NEWTON STEP INFEASIBLE-INTERIOR-POINT ALGORITHM COMPLEMENTARITY PROBLEMS
Yugoslav Journal of Operations Research 25 (205), Number, 57 72 DOI: 0.2298/YJOR3055034A A FULL-NEWTON STEP INFEASIBLE-INTERIOR-POINT ALGORITHM FOR P (κ)-horizontal LINEAR COMPLEMENTARITY PROBLEMS Soodabeh
More information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationNew weighing matrices and orthogonal designs constructed using two sequences with zero autocorrelation function - a review
University of Wollongong Research Online Faculty of Informatics - Paers (Archive) Faculty of Engineering and Information Sciences 1999 New weighing matrices and orthogonal designs constructed using two
More informationE-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula
e-comanion to Han Du and Zuluaga: Etension of Scarf s ma-min order formula ec E-comanion to A risk- and ambiguity-averse etension of the ma-min newsvendor order formula Qiaoming Han School of Mathematics
More informationMetrics Performance Evaluation: Application to Face Recognition
Metrics Performance Evaluation: Alication to Face Recognition Naser Zaeri, Abeer AlSadeq, and Abdallah Cherri Electrical Engineering Det., Kuwait University, P.O. Box 5969, Safat 6, Kuwait {zaery, abeer,
More informationRobust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System
Vol.7, No.7 (4),.37-38 htt://dx.doi.org/.457/ica.4.7.7.3 Robust Predictive Control of Inut Constraints and Interference Suression for Semi-Trailer System Zhao, Yang Electronic and Information Technology
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More informationAn efficient method for generalized linear multiplicative programming problem with multiplicative constraints
DOI 0.86/s40064-06-2984-9 RESEARCH An efficient method for generalized linear multilicative rogramming roblem with multilicative constraints Yingfeng Zhao,2* and Sanyang Liu Oen Access *Corresondence:
More informationApproximation of the Euclidean Distance by Chamfer Distances
Acta Cybernetica 0 (0 399 47. Aroximation of the Euclidean Distance by Chamfer Distances András Hajdu, Lajos Hajdu, and Robert Tijdeman Abstract Chamfer distances lay an imortant role in the theory of
More informationLineáris komplementaritási feladatok: elmélet, algoritmusok, alkalmazások
General LCPs Lineáris komplementaritási feladatok: elmélet, algoritmusok, alkalmazások Illés Tibor BME DET BME Analízis Szeminárium 2015. november 11. Outline Linear complementarity problem: M u + v =
More informationConvexification of Generalized Network Flow Problem with Application to Power Systems
1 Convexification of Generalized Network Flow Problem with Alication to Power Systems Somayeh Sojoudi and Javad Lavaei + Deartment of Comuting and Mathematical Sciences, California Institute of Technology
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationCCO Commun. Comb. Optim.
Communications in Combinatorics and Optimization Vol. 3 No., 08 pp.5-70 DOI: 0.049/CCO.08.580.038 CCO Commun. Comb. Optim. An infeasible interior-point method for the P -matrix linear complementarity problem
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More informationAn Overview of Witt Vectors
An Overview of Witt Vectors Daniel Finkel December 7, 2007 Abstract This aer offers a brief overview of the basics of Witt vectors. As an alication, we summarize work of Bartolo and Falcone to rove that
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. III Stability Theory - Peter C. Müller
STABILITY THEORY Peter C. Müller University of Wuertal, Germany Keywords: Asymtotic stability, Eonential stability, Linearization, Linear systems, Lyaunov equation, Lyaunov function, Lyaunov stability,
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationSecond Order Symmetric and Maxmin Symmetric Duality with Cone Constraints
International Journal of Oerations Research International Journal of Oerations Research Vol. 4, No. 4, 99 5 7) Second Order Smmetric Mamin Smmetric Dualit with Cone Constraints I. Husain,, Abha Goel, M.
More informationMultivariable Generalized Predictive Scheme for Gas Turbine Control in Combined Cycle Power Plant
Multivariable Generalized Predictive Scheme for Gas urbine Control in Combined Cycle Power Plant L.X.Niu and X.J.Liu Deartment of Automation North China Electric Power University Beiing, China, 006 e-mail
More informationProximal methods for the latent group lasso penalty
Proximal methods for the latent grou lasso enalty The MIT Faculty has made this article oenly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher Villa,
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationRotations in Curved Trajectories for Unconstrained Minimization
Rotations in Curved rajectories for Unconstrained Minimization Alberto J Jimenez Mathematics Deartment, California Polytechnic University, San Luis Obiso, CA, USA 9407 Abstract Curved rajectories Algorithm
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationarxiv: v1 [quant-ph] 20 Jun 2017
A Direct Couling Coherent Quantum Observer for an Oscillatory Quantum Plant Ian R Petersen arxiv:76648v quant-h Jun 7 Abstract A direct couling coherent observer is constructed for a linear quantum lant
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationAsymptotic Properties of the Markov Chain Model method of finding Markov chains Generators of..
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, -ISSN: 319-765X. Volume 1, Issue 4 Ver. III (Jul. - Aug.016), PP 53-60 www.iosrournals.org Asymtotic Proerties of the Markov Chain Model method of
More informationA Note on the Positive Nonoscillatory Solutions of the Difference Equation
Int. Journal of Math. Analysis, Vol. 4, 1, no. 36, 1787-1798 A Note on the Positive Nonoscillatory Solutions of the Difference Equation x n+1 = α c ix n i + x n k c ix n i ) Vu Van Khuong 1 and Mai Nam
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationAN INTERIOR POINT METHOD, BASED ON RANK-ONE UPDATES, Jos F. Sturm 1 and Shuzhong Zhang 2. Erasmus University Rotterdam ABSTRACT
October 13, 1995. Revised November 1996. AN INTERIOR POINT METHOD, BASED ON RANK-ONE UPDATES, FOR LINEAR PROGRAMMING Jos F. Sturm 1 Shuzhong Zhang Report 9546/A, Econometric Institute Erasmus University
More informationSOME TRACE INEQUALITIES FOR OPERATORS IN HILBERT SPACES
Kragujevac Journal of Mathematics Volume 411) 017), Pages 33 55. SOME TRACE INEQUALITIES FOR OPERATORS IN HILBERT SPACES SILVESTRU SEVER DRAGOMIR 1, Abstract. Some new trace ineualities for oerators in
More informationParallel Quantum-inspired Genetic Algorithm for Combinatorial Optimization Problem
Parallel Quantum-insired Genetic Algorithm for Combinatorial Otimization Problem Kuk-Hyun Han Kui-Hong Park Chi-Ho Lee Jong-Hwan Kim Det. of Electrical Engineering and Comuter Science, Korea Advanced Institute
More informationApplied Mathematics and Computation
Alied Mathematics and Comutation 217 (2010) 1887 1895 Contents lists available at ScienceDirect Alied Mathematics and Comutation journal homeage: www.elsevier.com/locate/amc Derivative free two-oint methods
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationNetwork Configuration Control Via Connectivity Graph Processes
Network Configuration Control Via Connectivity Grah Processes Abubakr Muhammad Deartment of Electrical and Systems Engineering University of Pennsylvania Philadelhia, PA 90 abubakr@seas.uenn.edu Magnus
More information1 Extremum Estimators
FINC 9311-21 Financial Econometrics Handout Jialin Yu 1 Extremum Estimators Let θ 0 be a vector of k 1 unknown arameters. Extremum estimators: estimators obtained by maximizing or minimizing some objective
More informationChapter 3. GMM: Selected Topics
Chater 3. GMM: Selected oics Contents Otimal Instruments. he issue of interest..............................2 Otimal Instruments under the i:i:d: assumtion..............2. he basic result............................2.2
More informationPositive Definite Uncertain Homogeneous Matrix Polynomials: Analysis and Application
BULGARIA ACADEMY OF SCIECES CYBEREICS AD IFORMAIO ECHOLOGIES Volume 9 o 3 Sofia 009 Positive Definite Uncertain Homogeneous Matrix Polynomials: Analysis and Alication Svetoslav Savov Institute of Information
More informationImproved Full-Newton-Step Infeasible Interior- Point Method for Linear Complementarity Problems
Georgia Southern University Digital Commons@Georgia Southern Mathematical Sciences Faculty Publications Mathematical Sciences, Department of 4-2016 Improved Full-Newton-Step Infeasible Interior- Point
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationThe Algebraic Structure of the p-order Cone
The Algebraic Structure of the -Order Cone Baha Alzalg Abstract We study and analyze the algebraic structure of the -order cones We show that, unlike oularly erceived, the -order cone is self-dual for
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationElements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley
Elements of Asymtotic Theory James L. Powell Deartment of Economics University of California, Berkeley Objectives of Asymtotic Theory While exact results are available for, say, the distribution of the
More informationConvex Analysis and Economic Theory Winter 2018
Division of the Humanities and Social Sciences Ec 181 KC Border Conve Analysis and Economic Theory Winter 2018 Toic 16: Fenchel conjugates 16.1 Conjugate functions Recall from Proosition 14.1.1 that is
More informationOn Doob s Maximal Inequality for Brownian Motion
Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t
More informationBest approximation by linear combinations of characteristic functions of half-spaces
Best aroximation by linear combinations of characteristic functions of half-saces Paul C. Kainen Deartment of Mathematics Georgetown University Washington, D.C. 20057-1233, USA Věra Kůrková Institute of
More informationModel checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]
Chater 5 Model checking, verification of CTL One must verify or exel... doubts, and convert them into the certainty of YES or NO. [Thomas Carlyle] 5. The verification setting Page 66 We introduce linear
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationMonopolist s mark-up and the elasticity of substitution
Croatian Oerational Research Review 377 CRORR 8(7), 377 39 Monoolist s mark-u and the elasticity of substitution Ilko Vrankić, Mira Kran, and Tomislav Herceg Deartment of Economic Theory, Faculty of Economics
More informationfrom the primal-dual interior-point algorithm (Megiddo [16], Kojima, Mizuno, and Yoshise
1. Introduction The primal-dual infeasible-interior-point algorithm which we will discuss has stemmed from the primal-dual interior-point algorithm (Megiddo [16], Kojima, Mizuno, and Yoshise [7], Monteiro
More informationOptimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment
Neutrosohic Sets and Systems Vol 14 016 93 University of New Mexico Otimal Design of Truss Structures Using a Neutrosohic Number Otimization Model under an Indeterminate Environment Wenzhong Jiang & Jun
More informationx 2 a mod m. has a solution. Theorem 13.2 (Euler s Criterion). Let p be an odd prime. The congruence x 2 1 mod p,
13. Quadratic Residues We now turn to the question of when a quadratic equation has a solution modulo m. The general quadratic equation looks like ax + bx + c 0 mod m. Assuming that m is odd or that b
More informationFigure : An 8 bridge design grid. (a) Run this model using LOQO. What is the otimal comliance? What is the running time?
5.094/SMA53 Systems Otimization: Models and Comutation Assignment 5 (00 o i n ts) Due Aril 7, 004 Some Convex Analysis (0 o i n ts) (a) Given ositive scalars L and E, consider the following set in three-dimensional
More informationOn Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.**
On Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.** * echnology and Innovation Centre, Level 4, Deartment of Electronic and Electrical Engineering, University of Strathclyde,
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationOn Superlinear Convergence of Infeasible Interior-Point Algorithms for Linearly Constrained Convex Programs *
Computational Optimization and Applications, 8, 245 262 (1997) c 1997 Kluwer Academic Publishers. Manufactured in The Netherlands. On Superlinear Convergence of Infeasible Interior-Point Algorithms for
More informationCombinatorics of topmost discs of multi-peg Tower of Hanoi problem
Combinatorics of tomost discs of multi-eg Tower of Hanoi roblem Sandi Klavžar Deartment of Mathematics, PEF, Unversity of Maribor Koroška cesta 160, 000 Maribor, Slovenia Uroš Milutinović Deartment of
More informationA WIDE NEIGHBORHOOD PRIMAL-DUAL INTERIOR-POINT ALGORITHM WITH ARC-SEARCH FOR LINEAR COMPLEMENTARITY PROBLEMS 1. INTRODUCTION
J Nonlinear Funct Anal 08 (08), Article ID 3 https://doiorg/0395/jnfa083 A WIDE NEIGHBORHOOD PRIMAL-DUAL INTERIOR-POINT ALGORITHM WITH ARC-SEARCH FOR LINEAR COMPLEMENTARITY PROBLEMS BEIBEI YUAN, MINGWANG
More information1. Introduction A number of recent papers have attempted to analyze the probabilistic behavior of interior point algorithms for linear programming. Ye
Probabilistic Analysis of an Infeasible-Interior-Point Algorithm for Linear Programming Kurt M. Anstreicher 1, Jun Ji 2, Florian A. Potra 3, and Yinyu Ye 4 Final Revision June, 1998 Abstract We consider
More informationOtimal exercise boundary for an American ut otion Rachel A. Kuske Deartment of Mathematics University of Minnesota-Minneaolis Minneaolis, MN 55455 e-mail: rachel@math.umn.edu Joseh B. Keller Deartments
More information