Real Root Isolation of Polynomial Equations Ba. Equations Based on Hybrid Computation
|
|
- Esmond Newman
- 5 years ago
- Views:
Transcription
1 Real Root Isolation of Polynomial Equations Based on Hybrid Computation Fei Shen 1 Wenyuan Wu 2 Bican Xia 1 LMAM & School of Mathematical Sciences, Peking University Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences Beijing, Oct. 27, 2012
2 Outline 1 Introduction The problem Related work Construct initial boxes An empirical estimation 5 6
3 The problem Related work The problem Let F = (f 1,..., f n ) T be a polynomial system defined on R n where f i R[x 1,..., x n ]. Suppose F(x) = 0 has only finite many real roots, say ξ (1),..., ξ (m). The target of real root isolation is to compute a family of regions S 1,..., S m, S j R n (1 j m), such that ξ (j) S j and S i S j = (1 i, j m). Hypothesis on the problem. 1 The system is square. 2 The system has only finite many roots. 3 The Jacobian matrix of F is non-singular at each root of F(x) = 0.
4 The problem Related work The problem Let F = (f 1,..., f n ) T be a polynomial system defined on R n where f i R[x 1,..., x n ]. Suppose F(x) = 0 has only finite many real roots, say ξ (1),..., ξ (m). The target of real root isolation is to compute a family of regions S 1,..., S m, S j R n (1 j m), such that ξ (j) S j and S i S j = (1 i, j m). Hypothesis on the problem. 1 The system is square. 2 The system has only finite many roots. 3 The Jacobian matrix of F is non-singular at each root of F(x) = 0.
5 The problem Related work Related work Univariate case Based on Descartes rule and bisection: Collins-Akritas(1976), Collins-Loos(1983), Collins-Johnson(1989), Johnson-Krandick(1997), von zur Gathen-Gerhard(1997), Johnson(1998), Rouillier-Zimmermann(2004), Johnson-Krandicket. al.(2005, 2006), Sagraloff(2012),... Based on Vincent s Theorem and continued fractions: Akritas(1980), Akritas-Strzeboński-et. al.(2005, 2006, 2008), Ştefănescu(2005), Sharma(2008),...
6 The problem Related work form. Multivariate case Xia-Yang(2002), Zhang-Xiao-Xia(2005), Xia-Zhang(2006), Cheng- Gao-Yap(2007), Boulier-Chen-Lemaire-Moreno Maza(2009), Mourrain-Pavone(2009), Cheng-Gao-Guo(2012),... The system should be in some special shapes, e.g., triangular symbolically. Or, the system has to be transformed into such forms Motivation Can we develop a method which avoids pre-processing the system symbolically?
7 The problem Related work form. Multivariate case Xia-Yang(2002), Zhang-Xiao-Xia(2005), Xia-Zhang(2006), Cheng- Gao-Yap(2007), Boulier-Chen-Lemaire-Moreno Maza(2009), Mourrain-Pavone(2009), Cheng-Gao-Guo(2012),... The system should be in some special shapes, e.g., triangular symbolically. Or, the system has to be transformed into such forms Motivation Can we develop a method which avoids pre-processing the system symbolically?
8 The problem Related work form. Multivariate case Xia-Yang(2002), Zhang-Xiao-Xia(2005), Xia-Zhang(2006), Cheng- Gao-Yap(2007), Boulier-Chen-Lemaire-Moreno Maza(2009), Mourrain-Pavone(2009), Cheng-Gao-Guo(2012),... The system should be in some special shapes, e.g., triangular symbolically. Or, the system has to be transformed into such forms Motivation Can we develop a method which avoids pre-processing the system symbolically?
9 Interval arithmetic Let X = [x, x], Y = [y, y] I (R), - X + Y = [x + y, x + y] - X Y = [x y, x y] - X Y = [min(xy, xy, xy, xy), max(xy, xy, xy, xy)] - X/Y = [x, x] [1/y, 1/y], 0 Y
10 Matlab package Intlab (Rump(1999, 2010)) Interval arithmetic operations, Jacobian matrix and Hessian matrix calculations, etc. Using floating-point arithmetic but the results are rigorous. Theoretically, the width of intervals for some special problems can be very small. Hence, we assume that the accuracy of numerical computation in this paper can be arbitrarily high (For arbitrarily high accuracy, we can call Matlab s vpa). However, it is also important to point out that such case rarely happens and double-precision is usually enough to obtain very small intervals in practice.
11 Matlab package Intlab (Rump(1999, 2010)) Interval arithmetic operations, Jacobian matrix and Hessian matrix calculations, etc. Using floating-point arithmetic but the results are rigorous. Theoretically, the width of intervals for some special problems can be very small. Hence, we assume that the accuracy of numerical computation in this paper can be arbitrarily high (For arbitrarily high accuracy, we can call Matlab s vpa). However, it is also important to point out that such case rarely happens and double-precision is usually enough to obtain very small intervals in practice.
12 Krawczyk operator Definition (Moore(1966)) Suppose f : D R n R n is continuous differentiable on D. Consider the equation f (x) = 0. Let f be the Jacobi matrix of f, F and F be the interval expand of f and f with inclusive monotonicity, respectively. For X I (D) and any y X, define the Krawczyk operator (K-operator) as: K(y, X) = y Y f (y) + (I Y F (X))(X y) (1) where Y is any n n non-singular matrix. K(X) = mid(x) mid(f (X)) 1 f (mid(x)) +(I mid(f (X)) 1 F (X))rad(X)[ 1, 1] (2)
13 Krawczyk operator Definition (Moore(1966)) Suppose f : D R n R n is continuous differentiable on D. Consider the equation f (x) = 0. Let f be the Jacobi matrix of f, F and F be the interval expand of f and f with inclusive monotonicity, respectively. For X I (D) and any y X, define the Krawczyk operator (K-operator) as: K(y, X) = y Y f (y) + (I Y F (X))(X y) (1) where Y is any n n non-singular matrix. K(X) = mid(x) mid(f (X)) 1 f (mid(x)) +(I mid(f (X)) 1 F (X))rad(X)[ 1, 1] (2)
14 Properties of K-operator Proposition (Moore(1966)) Suppose K(y, X) is defined as Formula (1), then 1 If x X is a root of F(x) = 0, then for any y X, we have x K(y, X); 2 y X, if X K(y, X) =, there is no roots in X; 3 y X and any non-singular matrix Y, if K(y, X) X, F(x) = 0 has a root in X; 4 y X and any non-singular matrix Y, if K(y, X) X, F(x) = 0 has only one root in X.
15 For a bounded region, say a box (i.e., interval vector), if we apply K-operator and subdivision strategy, it is possible to isolate all the roots in the box of F(x) = 0. The problem is: if the box is not small, the efficiency is low. Especially, if we want to isolate all the real roots of F(x) = 0, the estimated box is usually too big.
16 For a bounded region, say a box (i.e., interval vector), if we apply K-operator and subdivision strategy, it is possible to isolate all the roots in the box of F(x) = 0. The problem is: if the box is not small, the efficiency is low. Especially, if we want to isolate all the real roots of F(x) = 0, the estimated box is usually too big.
17 Homotopy continuation method is an important numerical computation method, which is used in various fields. We only treat it as an algorithm black box here, where the input is a polynomial system, and the output is its approximate roots. The details about the theory can be found in Li(1997) and Sommese-Wampler(2005). Homotopy software Hom4ps-2.0 (Li(2008)), PHCpack (Verschelde(1999)), HomLab (Wampler(2005)), etc. In our implementation, we use Hom4ps-2.0, which could return all the approximate complex roots of a given polynomial system efficiently, along with residues and condition numbers.
18 Homotopy continuation method is an important numerical computation method, which is used in various fields. We only treat it as an algorithm black box here, where the input is a polynomial system, and the output is its approximate roots. The details about the theory can be found in Li(1997) and Sommese-Wampler(2005). Homotopy software Hom4ps-2.0 (Li(2008)), PHCpack (Verschelde(1999)), HomLab (Wampler(2005)), etc. In our implementation, we use Hom4ps-2.0, which could return all the approximate complex roots of a given polynomial system efficiently, along with residues and condition numbers.
19 Homotopy continuation method is an important numerical computation method, which is used in various fields. We only treat it as an algorithm black box here, where the input is a polynomial system, and the output is its approximate roots. The details about the theory can be found in Li(1997) and Sommese-Wampler(2005). Homotopy software Hom4ps-2.0 (Li(2008)), PHCpack (Verschelde(1999)), HomLab (Wampler(2005)), etc. In our implementation, we use Hom4ps-2.0, which could return all the approximate complex roots of a given polynomial system efficiently, along with residues and condition numbers.
20 Construct initial boxes An empirical estimation Basic idea By Homotopy continuation, we can obtain all the approximate (complex) roots of the given polynomial system efficiently. Then, it is possible to get very small initial boxes. For each small initial box, K-operator can serve as an efficient tool verifying whether or not there is a real root in it. Key problem How to construct those initial boxes?
21 Construct initial boxes An empirical estimation Basic idea By Homotopy continuation, we can obtain all the approximate (complex) roots of the given polynomial system efficiently. Then, it is possible to get very small initial boxes. For each small initial box, K-operator can serve as an efficient tool verifying whether or not there is a real root in it. Key problem How to construct those initial boxes?
22 Construct initial boxes An empirical estimation Suppose z is an approximate root corresponding to the accurate root z of the system. Let I be an initial box to be constructed containing z. We must guarantee that I contains z. We give a rigorous result to estimate the distance between z and z based on the Kantorovich Theorem. For the sake of efficiency, we also give an empirical estimation to detect those roots with big imaginary parts which are highly impossible to be real.
23 Construct initial boxes An empirical estimation Suppose z is an approximate root corresponding to the accurate root z of the system. Let I be an initial box to be constructed containing z. We must guarantee that I contains z. We give a rigorous result to estimate the distance between z and z based on the Kantorovich Theorem. For the sake of efficiency, we also give an empirical estimation to detect those roots with big imaginary parts which are highly impossible to be real.
24 Construct initial boxes An empirical estimation Suppose z is an approximate root corresponding to the accurate root z of the system. Let I be an initial box to be constructed containing z. We must guarantee that I contains z. We give a rigorous result to estimate the distance between z and z based on the Kantorovich Theorem. For the sake of efficiency, we also give an empirical estimation to detect those roots with big imaginary parts which are highly impossible to be real.
25 Construct initial boxes An empirical estimation Construct initial boxes Theorem (Kantorovich) Let X and Y be Banach spaces and F : D X Y be an operator, which is Fréchet differentiable on an open convex set D 0 D. equation F(x) = 0, if the given approximate zero x 0 D 0 meets the following three conditions: F (x 0 ) 1 For exists, and there are real numbers B and η such that F (x 0 ) 1 B, F (x 0 ) 1 F(x 0 ) η, F satisfies the Lipschitz condition on D 0 : F (x) F (y) K x y, x, y D 0, h = BKη 1 2, O(x 0, 1 1 2h h η) D 0,
26 Construct initial boxes An empirical estimation Construct initial boxes Theorem (Kantorovich(continued)) then we claim that: F(x) = 0 has a root x in O(x 0, 1 1 2h h η) D 0, and the sequence {x k : x k+1 = x k F (x k ) 1 F(x k )} of Newton method converges to x ; For the convergence of x, we have: x x k+1 θ2k+1 (1 θ 2 ) θ(1 θ 2k+1 ) η where θ = 1 1 2h h ; The root x is unique in D 0 O(x 0, h h η). The proof can be found in Gragg-Tapia(1974).
27 Construct initial boxes An empirical estimation Construct initial boxes Proposition Let F = (f 1,..., f n) T be a polynomial system, where f i R[x 1,..., x n]. Denote by J the Jacobian matrix of F. For an approximation x 0 C n, if the following conditions hold: 1 J 1 (x 0) exists, and there are real numbers B and η such that J 1 (x 0) B, J 1 (x 0)F(x 0) η, 2 There exists a neighborhood O(x 0, ω) such that J(x) satisfies the Lipschitz condition on it: J(x) J(y) K x y 3 Let h = BKη, h 1 2, and ω 1 1 2h h η, then F(x) = 0 has only one root x in O(x 0, ω) O(x 0, h η). h
28 Construct initial boxes An empirical estimation Construct initial boxes Proof. We consider F as an operator on C n C n, obviously it is Fréchet differentiable, and from F(x +h) = F(x)+J(x)h +o(h) we can get F(x + h) F(x) J(x)h lim = 0. h 0 h Thus the first order Fréchet derivative of F is just the Jacobian matrix J, i.e. F (x) = J(x). So by the Kantorovich Theorem, the proof is completed immediately after checking the situation of x x 0.
29 Construct initial boxes An empirical estimation Construct initial boxes It is easy to know that 1 1 2h h 2. So we can just assign ω = 2η. Then we need to check whether BKη 1 2 in the neighborhood O(x 0, 2η). Even though the initial x 0 does not satisfy the conditions, we can still find a proper x k after several Newton iterations, since B and K are bounded and η will approach zero. And we only need to find an upper bound for the Lipschitz constant K.
30 Construct initial boxes An empirical estimation Construct initial boxes Bound for the Lipschitz constant K Let J ij = f i/ x j, apply mean value inequality (Mujica(1986)) to each element of J on O(x 0, ω) to get J ij(y) J ij(x) sup J ij(κ ij) y x, x, y O(x 0, ω) (3) κ ij line(x,y) where = ( / x 1, / x 2,..., / x n) is the gradient operator and line(x, y) refers to the line connecting x with y. By estimating the inequality, we finally get that K = max 1 i n n j=1 h (i) j (X 0) max, (4) where h (i) j are the column vectors of the Hessian matrix of f i.
31 Construct initial boxes An empirical estimation Construct initial boxes Bound for the Lipschitz constant K Let J ij = f i/ x j, apply mean value inequality (Mujica(1986)) to each element of J on O(x 0, ω) to get J ij(y) J ij(x) sup J ij(κ ij) y x, x, y O(x 0, ω) (3) κ ij line(x,y) where = ( / x 1, / x 2,..., / x n) is the gradient operator and line(x, y) refers to the line connecting x with y. By estimating the inequality, we finally get that K = max 1 i n n j=1 h (i) j (X 0) max, (4) where h (i) j are the column vectors of the Hessian matrix of f i.
32 Construct initial boxes An empirical estimation Algorithm: init width Input: Output: Equation F; Approximation x 0 ; Number of variables n Initial interval s radius r 1: repeat 2: x 0 = x 0 J 1 (x 0 )F(x 0 ); 3: B = J 1 (x 0 ) ; η = J 1 (x 0 )F(x 0 ) ; 4: ω = 2η; 5: X 0 = midrad(x 0, ω); 6: K = 0; 7: for i = 1 to n do 8: Compute the Hessian matrix H i = (h (i) 1, h(i) 2,..., h(i) n ) of F on X 0; 9: n if j=1 h(i) (X j 0 ) max > K then n j=1 h(i) (X j 0 ) max; 10: K = 11: end if 12: end for 13: h = BKη; 14: until h 1/2 15: return r = 1 1 2h η h
33 Construct initial boxes An empirical estimation An empirical estimation For F = 0, let x be an accurate root and x 0 be its approximation. we obtain the empirical estimation x x 0 J 1 (ξ)f(x 0) J 1 (x 0) F(x 0) 1 λ J 1 (x 0) 2 F(x 0), n where λ = max 1 i n j=1 h(i) j (x 0) max. Notice that the inequality is only a non-rigorous estimation. All the computation in it are carried out in a point-wise way, so it is faster than the rigorous result. In the numerical experiments later we will see that this empirically estimated radius performs very well. So we can use it to detect those non-real roots rather than the interval arithmetic.
34 Construct initial boxes An empirical estimation An empirical estimation For F = 0, let x be an accurate root and x 0 be its approximation. we obtain the empirical estimation x x 0 J 1 (ξ)f(x 0) J 1 (x 0) F(x 0) 1 λ J 1 (x 0) 2 F(x 0), n where λ = max 1 i n j=1 h(i) j (x 0) max. Notice that the inequality is only a non-rigorous estimation. All the computation in it are carried out in a point-wise way, so it is faster than the rigorous result. In the numerical experiments later we will see that this empirically estimated radius performs very well. So we can use it to detect those non-real roots rather than the interval arithmetic.
35 Our algorithms have been implemented as a Matlab program. We use Intlab of Version 6 for interval calculation and Hom4ps-2.0 as our homotopy continuation tool. All the experiments are undertaken in Matlab2008b on a PC (OS: Windows Vista, CPU: Inter R Core 2 Duo T GHz, Memory: 2G.). For arbitrarily high accuracy, we can call Matlab s vpa (variable precision arithmetic), but in fact all the real roots of the examples below are isolated by using Matlab s default double-precision floating point.
36 Our algorithms have been implemented as a Matlab program. We use Intlab of Version 6 for interval calculation and Hom4ps-2.0 as our homotopy continuation tool. All the experiments are undertaken in Matlab2008b on a PC (OS: Windows Vista, CPU: Inter R Core 2 Duo T GHz, Memory: 2G.). For arbitrarily high accuracy, we can call Matlab s vpa (variable precision arithmetic), but in fact all the real roots of the examples below are isolated by using Matlab s default double-precision floating point.
37 Equation # roots # real roots DISCOVERER complex roots detected barry cyclic cyclic des eco eco geneig kinema reimer reimer virasoro Table: Real root isolation results comparison
38 Equations Total time Homotopy time Interval time DISCOVERER barry cyclic cyclic > 1000 des > 1000 eco eco > 1000 geneig > 1000 kinema > 1000 reimer reimer > 1000 virasoro > 1000 Table: Execution time comparison, unit:s
39 Equation avg.rad.of intv avg.rad.of init. avg.time of iter max time of iter barry e e cyclic e e cyclic e e des e e eco e e eco e e geneig e e kinema e e reimer e e reimer e e virasoro e e Table: Detail data for each iteration, unit:s
40 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
41 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
42 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
43 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
44 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
45 For non-singular and zero-dimensional polynomial systems we propose a new algorithm for real root isolation based on hybrid computation. Our experiments on many benchmarks show that the new hybrid method is more efficient, compared with the traditional symbolic approaches. Future work: Systems with multiple roots or positive dimensional solutions (ongoing); Better construction of initial intervals; Parallel computation; Subdivision strategy and other techniques; Better implementation...
46 Thanks!
Real Root Isolation of Regular Chains.
Real Root Isolation of Regular Chains. François Boulier 1, Changbo Chen 2, François Lemaire 1, Marc Moreno Maza 2 1 University of Lille I (France) 2 University of London, Ontario (Canada) ASCM 2009 (Boulier,
More informationA New Method for Real Root Isolation of Univariate
A New Method for Real Root Isolation of Univariate Polynomials Ting Zhang and Bican Xia Abstract. A new algorithm for real root isolation of univariate polynomials is proposed, which is mainly based on
More informationQuaternion Data Fusion
Quaternion Data Fusion Yang Cheng Mississippi State University, Mississippi State, MS 39762-5501, USA William D. Banas and John L. Crassidis University at Buffalo, State University of New York, Buffalo,
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More information17 Solution of Nonlinear Systems
17 Solution of Nonlinear Systems We now discuss the solution of systems of nonlinear equations. An important ingredient will be the multivariate Taylor theorem. Theorem 17.1 Let D = {x 1, x 2,..., x m
More informationA Blackbox Polynomial System Solver on Parallel Shared Memory Computers
A Blackbox Polynomial System Solver on Parallel Shared Memory Computers Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science The 20th Workshop on
More informationA derivative-free nonmonotone line search and its application to the spectral residual method
IMA Journal of Numerical Analysis (2009) 29, 814 825 doi:10.1093/imanum/drn019 Advance Access publication on November 14, 2008 A derivative-free nonmonotone line search and its application to the spectral
More informationTutorial: Numerical Algebraic Geometry
Tutorial: Numerical Algebraic Geometry Back to classical algebraic geometry... with more computational power and hybrid symbolic-numerical algorithms Anton Leykin Georgia Tech Waterloo, November 2011 Outline
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationNumerical Analysis: Solving Nonlinear Equations
Numerical Analysis: Solving Nonlinear Equations Mirko Navara http://cmp.felk.cvut.cz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office 104a
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationSTOP, a i+ 1 is the desired root. )f(a i) > 0. Else If f(a i+ 1. Set a i+1 = a i+ 1 and b i+1 = b Else Set a i+1 = a i and b i+1 = a i+ 1
53 17. Lecture 17 Nonlinear Equations Essentially, the only way that one can solve nonlinear equations is by iteration. The quadratic formula enables one to compute the roots of p(x) = 0 when p P. Formulas
More informationFactoring Solution Sets of Polynomial Systems in Parallel
Factoring Solution Sets of Polynomial Systems in Parallel Jan Verschelde Department of Math, Stat & CS University of Illinois at Chicago Chicago, IL 60607-7045, USA Email: jan@math.uic.edu URL: http://www.math.uic.edu/~jan
More informationAn Alternative Three-Term Conjugate Gradient Algorithm for Systems of Nonlinear Equations
International Journal of Mathematical Modelling & Computations Vol. 07, No. 02, Spring 2017, 145-157 An Alternative Three-Term Conjugate Gradient Algorithm for Systems of Nonlinear Equations L. Muhammad
More informationExactness in numerical algebraic computations
Exactness in numerical algebraic computations Dan Bates Jon Hauenstein Tim McCoy Chris Peterson Andrew Sommese Wednesday, December 17, 2008 MSRI Workshop on Algebraic Statstics Main goals for this talk
More informationarxiv: v1 [math.na] 25 Sep 2012
Kantorovich s Theorem on Newton s Method arxiv:1209.5704v1 [math.na] 25 Sep 2012 O. P. Ferreira B. F. Svaiter March 09, 2007 Abstract In this work we present a simplifyed proof of Kantorovich s Theorem
More informationAn Accelerated Block-Parallel Newton Method via Overlapped Partitioning
An Accelerated Block-Parallel Newton Method via Overlapped Partitioning Yurong Chen Lab. of Parallel Computing, Institute of Software, CAS (http://www.rdcps.ac.cn/~ychen/english.htm) Summary. This paper
More informationarxiv: v1 [math.na] 28 Jun 2013
A New Operator and Method for Solving Interval Linear Equations Milan Hladík July 1, 2013 arxiv:13066739v1 [mathna] 28 Jun 2013 Abstract We deal with interval linear systems of equations We present a new
More informationA Blackbox Polynomial System Solver on Parallel Shared Memory Computers
A Blackbox Polynomial System Solver on Parallel Shared Memory Computers Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science 851 S. Morgan Street
More informationSolving Parametric Polynomial Systems by RealComprehensiveTriangularize
Solving Parametric Polynomial Systems by RealComprehensiveTriangularize Changbo Chen 1 and Marc Moreno Maza 2 1 Chongqing Key Laboratory of Automated Reasoning and Cognition, Chongqing Institute of Green
More informationVerified computations for hyperbolic 3-manifolds
Verified computations for hyperbolic 3-manifolds Neil Hoffman Kazuhiro Ichihara Masahide Kashiwagi Hidetoshi Masai Shin ichi Oishi Akitoshi Takayasu (U. Melbourne) (Nihon U.) (Waseda U.) (U. Tokyo) (Waseda
More informationComplete Numerical Isolation of Real Zeros in General Triangular Systems
Complete Numerical Isolation of Real Zeros in General Triangular Systems Jin-San Cheng 1, Xiao-Shan Gao 1 and Chee-Keng Yap 2,3 1 Key Lab of Mathematics Mechanization Institute of Systems Science, AMSS,
More informationNewton-type Methods for Solving the Nonsmooth Equations with Finitely Many Maximum Functions
260 Journal of Advances in Applied Mathematics, Vol. 1, No. 4, October 2016 https://dx.doi.org/10.22606/jaam.2016.14006 Newton-type Methods for Solving the Nonsmooth Equations with Finitely Many Maximum
More informationRoot Clustering for Analytic and their Complexity
Root Clustering for Analytic and their Complexity Chee Yap Courant Institute, NYU August 5, 2015 Mini-Symposium: Algorithms and Complexity in Polynomial System Solving SIAM AG15, Daejeon, Korea Joint work
More informationN. Hoffman (The University of Melbourne) ( CREST, JST) 1/15
3 N. Hoffman (The University of Melbourne) ( ) ( ) ( ) ( CREST, JST) ( ) 2014 @ 3 15 1/15 Gluing equations [W. Thurston] M: triangulated 3-manifold, possibly with torus boundary. equation s.t. whose solution
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationCache Complexity and Multicore Implementation for Univariate Real Root Isolation
Cache Complexity and Multicore Implementation for Univariate Real Root Isolation Changbo Chen, Marc Moreno Maza and Yuzhen Xie University of Western Ontario, London, Ontario, Canada E-mail: cchen5,moreno,yxie@csd.uwo.ca
More informationSOLVING VIA MODULAR METHODS
SOLVING VIA MODULAR METHODS DEEBA AFZAL, FAIRA KANWAL, GERHARD PFISTER, AND STEFAN STEIDEL Abstract. In this article we present a parallel modular algorithm to compute all solutions with multiplicities
More informationComputing Real Roots of Real Polynomials
Computing Real Roots of Real Polynomials and now For Real! Alexander Kobel Max-Planck-Institute for Informatics, Saarbrücken, Germany Fabrice Rouillier INRIA & Université Pierre et Marie Curie, Paris,
More informationComplete Numerical Isolation of Real Zeros in General Triangular Systems 1)
MM Research Preprints, 1 31 KLMM, AMSS, Academia Sinica Vol. 25, December 2006 1 Complete Numerical Isolation of Real Zeros in General Triangular Systems 1) Jin-San Cheng 1, Xiao-Shan Gao 1 and Chee-Keng
More informationInterval Arithmetic An Elementary Introduction and Successful Applications
Interval Arithmetic An Elementary Introduction and Successful Applications by R. Baker Kearfott Department of Mathematics University of Southwestern Louisiana U.S.L. Box 4-1010 Lafayette, LA 70504-1010
More informationResearch Article Strong Convergence of a Projected Gradient Method
Applied Mathematics Volume 2012, Article ID 410137, 10 pages doi:10.1155/2012/410137 Research Article Strong Convergence of a Projected Gradient Method Shunhou Fan and Yonghong Yao Department of Mathematics,
More informationOn the complexity of the hybrid proximal extragradient method for the iterates and the ergodic mean
On the complexity of the hybrid proximal extragradient method for the iterates and the ergodic mean Renato D.C. Monteiro B. F. Svaiter March 17, 2009 Abstract In this paper we analyze the iteration-complexity
More informationResearch Article A Novel Differential Evolution Invasive Weed Optimization Algorithm for Solving Nonlinear Equations Systems
Journal of Applied Mathematics Volume 2013, Article ID 757391, 18 pages http://dx.doi.org/10.1155/2013/757391 Research Article A Novel Differential Evolution Invasive Weed Optimization for Solving Nonlinear
More informationConvex Optimization on Large-Scale Domains Given by Linear Minimization Oracles
Convex Optimization on Large-Scale Domains Given by Linear Minimization Oracles Arkadi Nemirovski H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Joint research
More informationSecond Order Optimization Algorithms I
Second Order Optimization Algorithms I Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye Chapters 7, 8, 9 and 10 1 The
More informationNumerical Computations in Algebraic Geometry. Jan Verschelde
Numerical Computations in Algebraic Geometry Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science http://www.math.uic.edu/ jan jan@math.uic.edu AMS
More informationConvexity in R n. The following lemma will be needed in a while. Lemma 1 Let x E, u R n. If τ I(x, u), τ 0, define. f(x + τu) f(x). τ.
Convexity in R n Let E be a convex subset of R n. A function f : E (, ] is convex iff f(tx + (1 t)y) (1 t)f(x) + tf(y) x, y E, t [0, 1]. A similar definition holds in any vector space. A topology is needed
More informationCertified predictor-corrector tracking for Newton homotopies
Certified predictor-corrector tracking for Newton homotopies Jonathan D. Hauenstein Alan C. Liddell, Jr. March 11, 2014 Abstract We develop certified tracking procedures for Newton homotopies, which are
More information1 Newton s Method. Suppose we want to solve: x R. At x = x, f (x) can be approximated by:
Newton s Method Suppose we want to solve: (P:) min f (x) At x = x, f (x) can be approximated by: n x R. f (x) h(x) := f ( x)+ f ( x) T (x x)+ (x x) t H ( x)(x x), 2 which is the quadratic Taylor expansion
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationINERTIAL ACCELERATED ALGORITHMS FOR SOLVING SPLIT FEASIBILITY PROBLEMS. Yazheng Dang. Jie Sun. Honglei Xu
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX INERTIAL ACCELERATED ALGORITHMS FOR SOLVING SPLIT FEASIBILITY PROBLEMS Yazheng Dang School of Management
More informationCycling in Newton s Method
Universal Journal of Computational Mathematics 1(3: 73-77, 2013 DOI: 10.13189/ujcmj.2013.010302 http://www.hrpub.org Cycling in Newton s Method Mikheev Serge E. Faculty of Applied Mathematics & Control
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationNonlinear Equations and Continuous Optimization
Nonlinear Equations and Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems
More informationA Smoothing Newton Method for Solving Absolute Value Equations
A Smoothing Newton Method for Solving Absolute Value Equations Xiaoqin Jiang Department of public basic, Wuhan Yangtze Business University, Wuhan 430065, P.R. China 392875220@qq.com Abstract: In this paper,
More information0, otherwise. Find each of the following limits, or explain that the limit does not exist.
Midterm Solutions 1, y x 4 1. Let f(x, y) = 1, y 0 0, otherwise. Find each of the following limits, or explain that the limit does not exist. (a) (b) (c) lim f(x, y) (x,y) (0,1) lim f(x, y) (x,y) (2,3)
More informationA Posteriori Existence in Numerical Computations
Oxford University Computing Laboratory OXMOS: New Frontiers in the Mathematics of Solids OXPDE: Oxford Centre for Nonlinear PDE January, 2007 Introduction Non-linear problems are handled via the implicit
More informationAttitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation
Attitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation Yang Cheng Mississippi State University, Mississippi State, MS 3976 John L. Crassidis University at Buffalo, State
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationOn Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials. University of Bucharest CADE 2007.
On Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials Doru Ştefănescu University of Bucharest CADE 2007 21 February 2007 Contents Approximation of the Real Roots Bounds
More informationComputing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm
Computing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm Poteaux Adrien XLIM-DMI UMR CNRS 6172 Université de Limoges, France SNC'07 University of
More informationx 2 x n r n J(x + t(x x ))(x x )dt. For warming-up we start with methods for solving a single equation of one variable.
Maria Cameron 1. Fixed point methods for solving nonlinear equations We address the problem of solving an equation of the form (1) r(x) = 0, where F (x) : R n R n is a vector-function. Eq. (1) can be written
More informationEXISTENCE VERIFICATION FOR SINGULAR ZEROS OF REAL NONLINEAR SYSTEMS
EXISTENCE VERIFICATION FOR SINGULAR ZEROS OF REAL NONLINEAR SYSTEMS JIANWEI DIAN AND R BAKER KEARFOTT Abstract Traditional computational fixed point theorems, such as the Kantorovich theorem (made rigorous
More informationExplicit Criterion to Determine the Number of Positive Roots of a Polynomial 1)
MM Research Preprints, 134 145 No. 15, April 1997. Beijing Explicit Criterion to Determine the Number of Positive Roots of a Polynomial 1) Lu Yang and Bican Xia 2) Abstract. In a recent article, a complete
More informationarxiv:cs/ v1 [cs.cc] 16 Aug 2006
On Polynomial Time Computable Numbers arxiv:cs/0608067v [cs.cc] 6 Aug 2006 Matsui, Tetsushi Abstract It will be shown that the polynomial time computable numbers form a field, and especially an algebraically
More informationMaria Cameron. f(x) = 1 n
Maria Cameron 1. Local algorithms for solving nonlinear equations Here we discuss local methods for nonlinear equations r(x) =. These methods are Newton, inexact Newton and quasi-newton. We will show that
More informationHeterogenous Parallel Computing with Ada Tasking
Heterogenous Parallel Computing with Ada Tasking Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science http://www.math.uic.edu/ jan jan@math.uic.edu
More informationHOM4PS-2.0: A software package for solving polynomial systems by the polyhedral homotopy continuation method
HOM4PS-2.0: A software package for solving polynomial systems by the polyhedral homotopy continuation method Tsung-Lin Lee, Tien-Yien Li and Chih-Hsiung Tsai January 15, 2008 Abstract HOM4PS-2.0 is a software
More informationLecture 7 Monotonicity. September 21, 2008
Lecture 7 Monotonicity September 21, 2008 Outline Introduce several monotonicity properties of vector functions Are satisfied immediately by gradient maps of convex functions In a sense, role of monotonicity
More informationAn Accelerated Hybrid Proximal Extragradient Method for Convex Optimization and its Implications to Second-Order Methods
An Accelerated Hybrid Proximal Extragradient Method for Convex Optimization and its Implications to Second-Order Methods Renato D.C. Monteiro B. F. Svaiter May 10, 011 Revised: May 4, 01) Abstract This
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationRemarks On Second Order Generalized Derivatives For Differentiable Functions With Lipschitzian Jacobian
Applied Mathematics E-Notes, 3(2003), 130-137 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Remarks On Second Order Generalized Derivatives For Differentiable Functions
More informationNew Inexact Line Search Method for Unconstrained Optimization 1,2
journal of optimization theory and applications: Vol. 127, No. 2, pp. 425 446, November 2005 ( 2005) DOI: 10.1007/s10957-005-6553-6 New Inexact Line Search Method for Unconstrained Optimization 1,2 Z.
More informationNumerical Irreducible Decomposition
Numerical Irreducible Decomposition Jan Verschelde Department of Math, Stat & CS University of Illinois at Chicago Chicago, IL 60607-7045, USA e-mail: jan@math.uic.edu web: www.math.uic.edu/ jan CIMPA
More informationRecent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables
Recent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables Department of Systems Engineering and Engineering Management The Chinese University of Hong Kong 2014 Workshop
More informationIntroduction to Numerical Algebraic Geometry
Introduction to Numerical Algebraic Geometry What is possible in Macaulay2? Anton Leykin Georgia Tech Applied Macaulay2 tutorials, Georgia Tech, July 2017 Polynomial homotopy continuation Target system:
More informationINTRODUCTION TO COMPUTATIONAL MATHEMATICS
INTRODUCTION TO COMPUTATIONAL MATHEMATICS Course Notes for CM 271 / AMATH 341 / CS 371 Fall 2007 Instructor: Prof. Justin Wan School of Computer Science University of Waterloo Course notes by Prof. Hans
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 18, No. 1, pp. 106 13 c 007 Society for Industrial and Applied Mathematics APPROXIMATE GAUSS NEWTON METHODS FOR NONLINEAR LEAST SQUARES PROBLEMS S. GRATTON, A. S. LAWLESS, AND N. K.
More informationSolving Polynomial Systems in the Cloud with Polynomial Homotopy Continuation
Solving Polynomial Systems in the Cloud with Polynomial Homotopy Continuation Jan Verschelde joint with Nathan Bliss, Jeff Sommars, and Xiangcheng Yu University of Illinois at Chicago Department of Mathematics,
More informationContinuous Optimization
Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2009 Outline 1 Introduction 2 Golden Section Search 3 Multivariate Functions Steepest Descent Method
More informationj=1 r 1 x 1 x n. r m r j (x) r j r j (x) r j (x). r j x k
Maria Cameron Nonlinear Least Squares Problem The nonlinear least squares problem arises when one needs to find optimal set of parameters for a nonlinear model given a large set of data The variables x,,
More informationIterative Reweighted Minimization Methods for l p Regularized Unconstrained Nonlinear Programming
Iterative Reweighted Minimization Methods for l p Regularized Unconstrained Nonlinear Programming Zhaosong Lu October 5, 2012 (Revised: June 3, 2013; September 17, 2013) Abstract In this paper we study
More informationOPER 627: Nonlinear Optimization Lecture 9: Trust-region methods
OPER 627: Nonlinear Optimization Lecture 9: Trust-region methods Department of Statistical Sciences and Operations Research Virginia Commonwealth University Sept 25, 2013 (Lecture 9) Nonlinear Optimization
More informationConvex Optimization Notes
Convex Optimization Notes Jonathan Siegel January 2017 1 Convex Analysis This section is devoted to the study of convex functions f : B R {+ } and convex sets U B, for B a Banach space. The case of B =
More informationNumerical Integration over unit sphere by using spherical t-designs
Numerical Integration over unit sphere by using spherical t-designs Congpei An 1 1,Institute of Computational Sciences, Department of Mathematics, Jinan University Spherical Design and Numerical Analysis
More informationToolboxes and Blackboxes for Solving Polynomial Systems. Jan Verschelde
Toolboxes and Blackboxes for Solving Polynomial Systems Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science http://www.math.uic.edu/ jan jan@math.uic.edu
More informationarxiv: v2 [math.ag] 24 Jun 2015
TRIANGULATIONS OF MONOTONE FAMILIES I: TWO-DIMENSIONAL FAMILIES arxiv:1402.0460v2 [math.ag] 24 Jun 2015 SAUGATA BASU, ANDREI GABRIELOV, AND NICOLAI VOROBJOV Abstract. Let K R n be a compact definable set
More informationMotivation: We have already seen an example of a system of nonlinear equations when we studied Gaussian integration (p.8 of integration notes)
AMSC/CMSC 460 Computational Methods, Fall 2007 UNIT 5: Nonlinear Equations Dianne P. O Leary c 2001, 2002, 2007 Solving Nonlinear Equations and Optimization Problems Read Chapter 8. Skip Section 8.1.1.
More informationThe Inverse Function Theorem via Newton s Method. Michael Taylor
The Inverse Function Theorem via Newton s Method Michael Taylor We aim to prove the following result, known as the inverse function theorem. Theorem 1. Let F be a C k map (k 1) from a neighborhood of p
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationA globally and R-linearly convergent hybrid HS and PRP method and its inexact version with applications
A globally and R-linearly convergent hybrid HS and PRP method and its inexact version with applications Weijun Zhou 28 October 20 Abstract A hybrid HS and PRP type conjugate gradient method for smooth
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationAN ITERATIVE METHOD FOR COMPUTING ZEROS OF OPERATORS SATISFYING AUTONOMOUS DIFFERENTIAL EQUATIONS
Electronic Journal: Southwest Journal of Pure Applied Mathematics Internet: http://rattler.cameron.edu/swjpam.html ISBN 1083-0464 Issue 1 July 2004, pp. 48 53 Submitted: August 21, 2003. Published: July
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationRoot Finding (and Optimisation)
Root Finding (and Optimisation) M.Sc. in Mathematical Modelling & Scientific Computing, Practical Numerical Analysis Michaelmas Term 2018, Lecture 4 Root Finding The idea of root finding is simple we want
More informationSolving non-linear systems of equations
Solving non-linear systems of equations Felix Kubler 1 1 DBF, University of Zurich and Swiss Finance Institute October 7, 2017 Felix Kubler Comp.Econ. Gerzensee, Ch2 October 7, 2017 1 / 38 The problem
More informationTopology of implicit curves and surfaces
Topology of implicit curves and surfaces B. Mourrain, INRIA, BP 93, 06902 Sophia Antipolis mourrain@sophia.inria.fr 24th October 2002 The problem Given a polynomial f(x, y) (resp. f(x, y, z)), Compute
More informationA Regularized Directional Derivative-Based Newton Method for Inverse Singular Value Problems
A Regularized Directional Derivative-Based Newton Method for Inverse Singular Value Problems Wei Ma Zheng-Jian Bai September 18, 2012 Abstract In this paper, we give a regularized directional derivative-based
More informationMethods that avoid calculating the Hessian. Nonlinear Optimization; Steepest Descent, Quasi-Newton. Steepest Descent
Nonlinear Optimization Steepest Descent and Niclas Börlin Department of Computing Science Umeå University niclas.borlin@cs.umu.se A disadvantage with the Newton method is that the Hessian has to be derived
More informationSolving Polynomial Systems Symbolically and in Parallel
Solving Polynomial Systems Symbolically and in Parallel Marc Moreno Maza & Yuzhen Xie Ontario Research Center for Computer Algebra University of Western Ontario, London, Canada MITACS - CAIMS, June 18,
More informationConvergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces
Applied Mathematical Sciences, Vol. 6, 212, no. 63, 319-3117 Convergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces Nguyen Buong Vietnamese
More informationAcceleration Method for Convex Optimization over the Fixed Point Set of a Nonexpansive Mapping
Noname manuscript No. will be inserted by the editor) Acceleration Method for Convex Optimization over the Fixed Point Set of a Nonexpansive Mapping Hideaki Iiduka Received: date / Accepted: date Abstract
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 14, 017 One of the most frequently occurring problems in scientific work is to find the roots of equations of the form f(x) = 0. (1)
More informationParallel and Distributed Stochastic Learning -Towards Scalable Learning for Big Data Intelligence
Parallel and Distributed Stochastic Learning -Towards Scalable Learning for Big Data Intelligence oé LAMDA Group H ŒÆOŽÅ Æ EâX ^ #EâI[ : liwujun@nju.edu.cn Dec 10, 2016 Wu-Jun Li (http://cs.nju.edu.cn/lwj)
More informationComparison between Interval and Fuzzy Load Flow Methods Considering Uncertainty
Comparison between Interval and Fuzzy Load Flow Methods Considering Uncertainty T.Srinivasarao, 2 P.Mallikarajunarao Department of Electrical Engineering, College of Engineering (A), Andhra University,
More informationMotivation Subgradient Method Stochastic Subgradient Method. Convex Optimization. Lecture 15 - Gradient Descent in Machine Learning
Convex Optimization Lecture 15 - Gradient Descent in Machine Learning Instructor: Yuanzhang Xiao University of Hawaii at Manoa Fall 2017 1 / 21 Today s Lecture 1 Motivation 2 Subgradient Method 3 Stochastic
More informationEconomics 204 Fall 2013 Problem Set 5 Suggested Solutions
Economics 204 Fall 2013 Problem Set 5 Suggested Solutions 1. Let A and B be n n matrices such that A 2 = A and B 2 = B. Suppose that A and B have the same rank. Prove that A and B are similar. Solution.
More informationIs the Hénon map chaotic
Is the Hénon map chaotic Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology, Poland, galias@agh.edu.pl International Workshop on Complex Networks and Applications
More information