Point (geometry) - Revision history

en>ClueBot NG: Reverting possible vandalism by 123.231.21.238 to version by D.Lazard. False positive? Report it. Thanks, ClueBot NG. (2049609) (Bot)

2014-12-03T08:28:43Z

Reverting possible vandalism by 123.231.21.238 to version by D.Lazard. False positive? Report it. Thanks, ClueBot NG. (2049609) (Bot)

← Older revision		Revision as of 10:28, 3 December 2014
Line 1:		Line 1:
	~~Common history~~ of the ~~writer was Adrianne Quesada. Her husband~~ [http://~~Photo~~.~~net/gallery/tag-search~~/search?~~query_string~~=~~doesn%27t doesn't~~] ~~like~~ the ~~idea~~ the ~~way she truly but what she relatively likes doing is~~ to ~~bake but she's pondering on starting something fresh~~, ~~new. After being out~~ of the ~~boy~~'s ~~job~~ for ~~years~~ the ~~guy became an order maid~~ of ~~honor. Vermont~~ will ~~where her~~ house ~~has always been~~. ~~Her husband~~ and ~~her care~~ for a ~~website~~. ~~You might wish~~ to ~~check it out:~~ http://circuspartypanama.com<br><br>~~Here~~ is ~~my homepage ::~~ [http://~~circuspartypanama~~.com ~~clash of clans cheat gems~~]		Regardless of the Clash of Clans hack tool; there might be also [http://search.un.org/search?ie=utf8&site=un_org&output=xml_no_dtd&client=UN_Website_en&num=10&lr=lang_en&proxystylesheet=UN_Website_en&oe=utf8&q=hack+tools&Submit=Go hack tools] by other games. Americans can check out all of the hacks and obtain dozens of which they need. It is sure that will have lost of fun once they have the hack tool saved.<br><br>Its upsides of video flash games can include fun, recreational and even education. The downsides range by means of addictive game play and younger individuals seeing and so hearing things they continue to be not old enough in order for. With luck, all of the ideas presented within the following article can help customers manage video games amazingly well within your home to gain everyone's benefit.<br><br>Circumvent purchasing big title events near their launch dates. Waiting means that you're prone to acquire clash of clans cheats after developing a patch or two is carrying emerge to mend obtrusive holes and bugs which will impact your pleasure also game play. Near keep an eye through for titles from studios which are understood our bodies and wellbeing patching and support.<br><br>Normally clash of clans get into tool no survey builds believe in among people. Society may be definitely powered by professional pressure, one of the most powerful forces during the planet. Whereas long as peer power utilizes its power for many good, clash of clans hack tool no feedback survey will have its house in community.<br><br>Your site is offering Clash along with Clans hack tool trainer to users who would like to. The website offering this tool is safe and it guarantees finist quality software. There likewise other sites which secure the tool. But would be the either incomplete or connected with bad quality. When users download these just a few hack tools, instead to complete well they end along in trouble. So, players are advised pick the tool from an internet site that offers complete computer software programs.Users who are finding it tough to cross-stitching the hurdles can look for a site that allows folks to download the cheats. Most of the world-wide-web websites allow free download a few websites charge fees. If you have any kind of concerns relating to where and how you can utilize [http://circuspartypanama.com clash of clans cheat gems], you can call us at the web page. Users can locate an online site from where they can obtain good quality software.<br><br>Purchase are playing a dressing in game, and you have no any experience with it, set the difficulty skill level to rookie. This is considered help you pick in on the unique presents of the game learn your way in the vicinity of the field. If you set it higher than that, you commonly tend to get frustrated and thus not have any pleasure.<br><br>Now you have read this composition, you need to a great easier time locating as well as a loving video games in your daily life. Notwithstanding your favored platform, from your cellphone in your own own computer, playing as well as a enjoying video gaming can allow you to take the benefit of the worries of the actual busy week get [http://Browse.deviantart.com/?qh=&section=&global=1&q=specifics specifics].

en>Wavelength: inserting 2 hyphens: —> "two-dimensional" and "three-dimensional"—wikt:two-dimensional—wikt:three-dimensional

2014-02-27T00:08:25Z

inserting 2 hyphens: —> "two-dimensional" and "three-dimensional"—wikt:two-dimensional—wikt:three-dimensional

← Older revision		Revision as of 02:08, 27 February 2014
Line 1:		Line 1:
	~~In [[numerical linear algebra]], the '''QR algorithm''' is an [[eigenvalue algorithm]]: that is, a procedure to calculate~~ the ~~[[eigenvalues and eigenvectors]] of a [[Matrix (mathematics)\|matrix]]. The QR transformation~~ was ~~developed in the late 1950s by [[John G.F. Francis]] (England) and by [[Vera N. Kublanovskaya]] (USSR), working independently.<ref>~~		Common history of the writer was Adrianne Quesada. Her husband [http://Photo.net/gallery/tag-search/search?query_string=doesn%27t doesn't] like the idea the way she truly but what she relatively likes doing is to bake but she's pondering on starting something fresh, new. After being out of the boy's job for years the guy became an order maid of honor. Vermont will where her house has always been. Her husband and her care for a website. You might wish to check it out: http://circuspartypanama.com<br><br>Here is my homepage :: [http://circuspartypanama.com clash of clans cheat gems]
	~~J.G.F. Francis, "The QR Transformation, I", ''The Computer Journal'', vol. 4, no~~. ~~3, pages 265-271 (1961, received Oct 1959)~~ [http://~~comjnl~~.~~oxfordjournals.org~~/~~cgi/content/abstract/4/3/265 online at oxfordjournals.org];<br~~ />
	~~J.G.F. Francis, "The QR Transformation, II" ''The Computer Journal'', vol. 4, no. 4, pages 332~~-~~345 (1962) [http:~~/~~/comjnl.oxfordjournals.org/cgi/content/abstract/4/4/332 online at oxfordjournals.org].<br />~~
	Vera N. Kublanovskaya, "On some algorithms for the solution of the complete eigenvalue problem," ''USSR Computational Mathematics and Mathematical Physics'', vol. 1, no. 3, pages 637–657 (1963, received Feb 1961). Also published in: ''Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki'', vol.1, no. 4, pages 555–570 (1961).</ref> The basic idea is to perform a [[QR decomposition]], writing the matrix as a product of an [[orthogonal matrix]] and an upper [[triangular matrix]], multiply the factors in the reverse order, and iterate.

	~~==The practical QR algorithm==~~
	Formally, let ''A'' be a real matrix of which we want to compute the eigenvalues, and let ''A''<sub>0</sub>:=''A''. At the ''k''-th step (starting with ''k'' = 0), we compute the QR decomposition ''A''<sub>''k''</sub>='~~'Q''<sub>''k''</sub>''R''<sub>''k''</sub> where ''Q''<sub>''k''</sub> is an [[orthogonal matrix~~]~~] and ''R''<sub>''k''</sub> is an upper triangular matrix. We then form ''A''<sub>''k''+1</sub> = ''R''<sub>''k''</sub>''Q''<sub>''k''</sub>. Note that~~
	~~:<math> A_{k+1} = R_k Q_k = Q_k^T Q_k R_k Q_k = Q_k^T A_k Q_k = Q_k^{-1} A_k Q_k, </math>~~
	~~so all~~ the ~~''A''<sub>''k''</sub> are [[Similar matrix\|similar]] and hence they have~~ the ~~same eigenvalues. The algorithm~~ is ~~[[numerical stability\|numerically stable]] because it proceeds by ''orthogonal'' similarity transforms.~~

	Under certain conditions,<ref>Golub, G. H. and Van Loan, C. F.: Matrix Computations, 3rd ed., Johns Hopkins University Press, Baltimore, 1996, ISBN 0-8018-5414-8.</ref> the matrices ''A''<sub>''k''</sub> converge to a triangular matrix, the [[Schur form]] of ''A''. The eigenvalues of a triangular matrix are listed on the diagonal, and the eigenvalue problem is solved. In testing for convergence it is impractical to require exact zeros, but ~~the [[Gershgorin circle theorem]] provides a bound on the error.~~

	~~In this crude form the iterations are relatively expensive. This can be mitigated by first bringing the matrix~~ ''A'' to upper [[Hessenberg form]] (which costs <math>\begin{matrix}\frac{10}{3}\end{matrix} n^3 + O(n^2)</math> arithmetic operations using a technique based on [[Householder transformation\|Householder reduction]]), with a finite sequence of orthogonal similarity transforms, somewhat like a two-sided QR decomposition.<ref name=Demmel>[[James W. Demmel]], ''Applied Numerical Linear Algebra'' (SIAM, 1997).</ref><ref name=Trefethen>[[Lloyd N. Trefethen]] and David Bau, ''Numerical Linear Algebra'' (SIAM, 1997).</ref> (For QR decomposition, the Householder reflectors are multiplied only on the left, but for the Hessenberg case they are multiplied on both left and right.) Determining the QR decomposition of an upper Hessenberg matrix costs <math>6 n^2 + O(n)</math> arithmetic operations. Moreover, because the Hessenberg form is already nearly upper-triangular (it has just one nonzero entry below each diagonal), using it as a starting ~~point reduces the number of steps required for convergence of the QR algorithm.~~

	~~If the original matrix is [[symmetric matrix\|symmetric]], then the upper Hessenberg matrix is also symmetric and thus [[tridiagonal matrix\|tridiagonal]]~~, ~~and so are all the ''A''<sub>''k''</sub>~~. This procedure costs <math>\begin{matrix}\frac{4}{3}\end{matrix} n^3 + O(n^2)</math> arithmetic operations using a technique based on Householder reduction.<ref name=Demmel/><ref name=Trefethen/> Determining the QR decomposition of a symmetric tridiagonal matrix costs <math>O(n)</math> operations.<ref>James M. Ortega and Henry F. Kaiser, "The ''LL<sup>T</sup>'' and ''QR'' methods for symmetric tridiagonal matrices," ''The Computer Journal'' '''6''' (1), 99–101 (1963).</ref>

	The rate of convergence depends on the separation between eigenvalues, so a practical algorithm will use shifts, either explicit or implicit, to increase separation and accelerate convergence. A typical symmetric QR algorithm isolates each eigenvalue (then reduces the size of the matrix) with only one or two iterations, making it efficient as well as robust.{{clarify\|date=June 2012}}

	~~== The implicit QR algorithm ==~~
	In modern computational practice,<ref>Golub and Van Loan, chapter 7</ref> the QR algorithm is performed in an implicit version which makes the use of multiple shifts easier to introduce. The matrix is first brought to upper Hessenberg form <math>A_0=QAQ^T</math> as in the explicit version; then, at each step, the first column of <math>A_k</math> is transformed via a small-size Householder similarity transformation to the first column of <math>p(A_k)</math> (or <math>p(A_k)e_1</math>), where <math>p(A_k)</math>, of degree <math>r</math>, is the polynomial that defines the shifting strategy (often <math>p(x)=(x-\lambda)(x-\bar\lambda)</math>, where <math>\lambda</math> and <math>\bar\lambda</math> are the two eigenvalues of the trailing <math>2 \times 2</math> principal submatrix of ~~<math>A_k</math>,~~ the ~~so-called ''implicit double-shift~~''). Then successive Householder transformation of size <math>r+1</math> are performed in order to return the working matrix <math>A_k</math> to upper Hessenberg form. This operation is known as ''bulge chasing'', due to the peculiar shape of the non-zero entries of the matrix along the steps of the algorithm. As in the first version, deflation is performed as soon as one of the sub-diagonal entries of <math>A_k</math> is sufficiently small.

	~~===Renaming proposal===~~
	~~Since in the modern implicit version of the procedure no [[QR decomposition]]~~s ~~are explicitly performed, some authors,~~ for instance Watkins,<ref>Watkins, David S.: The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods, SIAM, Philadelphia, PA, USA, 2007, ISBN 0-89871-641-1, ISBN 978-0-89871-641-2</ref> suggested changing its name to ''Francis algorithm''. [[Gene H. Golub\|Golub]] and [[Charles F. Van Loan\|Van Loan]] use the term ''Francis QR step''.

	~~== Interpretation and convergence ==~~
	~~The QR algorithm can be seen as a more sophisticated variation of~~ the ~~basic "power" [[eigenvalue algorithm]]. Recall that the power algorithm repeatedly multiplies ''A'' times a single vector, normalizing after each iteration. The vector converges to~~ an ~~eigenvector~~ of ~~the largest eigenvalue~~. ~~Instead, the QR algorithm works with a complete basis of vectors, using QR decomposition to renormalize (and orthogonalize). For a symmetric matrix ''A'', upon convergence, ''AQ'' = ''QΛ'',~~ where ''Λ'' is the diagonal matrix of eigenvalues to which ''A'' converged, and where ''Q'' is a composite of all the orthogonal similarity transforms required to get there. Thus the columns of ''Q'' are the eigenvectors.

	~~== History ==~~
	~~The QR algorithm was preceded by the ''LR algorithm'', which uses the [[LU decomposition]] instead of the QR decomposition. The QR algorithm is more stable, so the LR algorithm is rarely used nowadays~~. ~~However, it represents an important step in the development of the QR algorithm.~~

	The LR algorithm was developed in the early 1950s by [[Heinz Rutishauer]], who worked at that time as a research assistant of [[Eduard Stiefel]] at [[ETH Zurich]]. Stiefel suggested that Rutishauer use the sequence of moments ''y''<sub>0</sub><sup>T</sup> ''A''<sup>''k''</sup> ''x''<sub>0</sub>, ''k'' = 0, 1, … (where ''x''<sub>0</sub> and ~~''y''<sub>0</sub> are arbitrary vectors) to find the eigenvalues of ''A''. Rutishauer took an algorithm of [[Alexander Aitken]]~~ for ~~this task and developed it into the ''quotient–difference algorithm'' or ''qd algorithm''. After arranging the computation in~~ a suitable shape, he discovered that the qd algorithm is in fact the iteration ''A''<sub>''k''</sub> = ''L''<sub>''k''</sub>''U''<sub>''k''</sub> (LU decomposition), ''A''<sub>''k''+1</sub> = ''U''<sub>''k''</sub>''L''<sub>''k''</sub>, applied on a tridiagonal matrix, from which the LR algorithm follows.~~<ref>{{Citation \| last1=Parlett \| first1=Beresford N. \| last2=Gutknecht \| first2=Martin H. \| title=From qd~~ to ~~LR, or, how were the qd and LR algorithms discovered? \| doi=10.1093/imanum/drq003 \| year=2011 \| journal=IMA Journal of Numerical Analysis \| issn=0272-4979 \| volume=31 \| pages=741–754}}</ref>~~

	~~== Other variants ==~~
	One variant of the ''QR algorithm'', ''the Golub-Kahan-Reinsch'' algorithm starts with reducing a general matrix into a bidiagonal one.<ref>Bochkanov Sergey Anatolyevich. ALGLIB User Guide - General Matrix operations - Singular value decomposition . ALGLIB Project. 2010-12-11. URL:http://~~www~~.~~alglib.net/matrixops/general/svd.php. Accessed: 2010-12-11. (Archived by WebCite at http://www.webcitation.org/5utO4iSnR)~~<~~/ref~~> ~~This variant of the ''QR algorithm'' for the computation of eigenvalues, which was first described by {{harvtxt\|Golub\|Kahan\|1965}}. The [[LAPACK]] subroutine [http~~://www.netlib.org/lapack/double/dbdsqr.f DBDSQR] implements this iterative method, with some modifications to cover the case where the singular values are very small {{harv\|Demmel\|Kahan\|1990}}. Together with a first step using Householder reflections and, if appropriate, [[QR decomposition]], this forms the [http:~~//www.netlib.org/lapack/double/dgesvd.f DGESVD] routine for the computation of the [[singular value decomposition]].~~

	~~== References ==~~
	~~<references/>~~

	~~== External links ==~~
	* {{planetmath reference\|id=1474\|title=Eigenvalue problem}}
	* [http://~~www.math.umn.edu/~olver/aims_/qr~~.~~pdf Prof. Peter Olver's notes on orthogonal bases and the workings~~ of ~~the QR algorithm]~~
	* [http://math.fullerton.edu/mathews/n2003/QRMethodMod.html Module for the QR Method]

	~~{{Numerical linear algebra}}~~

	~~{{DEFAULTSORT:Qr Algorithm}}~~
	~~[[Category:Numerical linear algebra]~~]

en>Twsx: unneeded / unhelpful wikilink (HG)

2014-01-16T10:18:10Z

unneeded / unhelpful wikilink (HG)

← Older revision		Revision as of 12:18, 16 January 2014
Line 1:		Line 1:
	~~My name~~: ~~Henry Dawkins~~<br>~~My age~~: 33<br>~~Country~~: ~~Sweden~~<br>~~Home town~~: ~~Gullholmen~~ <br>~~ZIP~~: ~~440 84~~<br>~~Street~~: ~~Villagatan 48~~<br><br>~~my web site~~ :: [~~http:~~//~~ebook-pdfree~~.~~com~~/~~free~~-~~mark~~-~~athena~~-~~pdf~~-~~heroes~~-~~olympus~~-~~book~~-3-~~download~~/ The ~~Mark~~ of ~~Athena~~ pdf]		In [[numerical linear algebra]], the '''QR algorithm''' is an [[eigenvalue algorithm]]: that is, a procedure to calculate the [[eigenvalues and eigenvectors]] of a [[Matrix (mathematics)\|matrix]]. The QR transformation was developed in the late 1950s by [[John G.F. Francis]] (England) and by [[Vera N. Kublanovskaya]] (USSR), working independently.<ref>
			J.G.F. Francis, "The QR Transformation, I", ''The Computer Journal'', vol. 4, no. 3, pages 265-271 (1961, received Oct 1959) [http://comjnl.oxfordjournals.org/cgi/content/abstract/4/3/265 online at oxfordjournals.org];<br />
			J.G.F. Francis, "The QR Transformation, II" ''The Computer Journal'', vol. 4, no. 4, pages 332-345 (1962) [http://comjnl.oxfordjournals.org/cgi/content/abstract/4/4/332 online at oxfordjournals.org].<br />
			Vera N. Kublanovskaya, "On some algorithms for the solution of the complete eigenvalue problem," ''USSR Computational Mathematics and Mathematical Physics'', vol. 1, no. 3, pages 637–657 (1963, received Feb 1961). Also published in: ''Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki'', vol.1, no. 4, pages 555–570 (1961).</ref> The basic idea is to perform a [[QR decomposition]], writing the matrix as a product of an [[orthogonal matrix]] and an upper [[triangular matrix]], multiply the factors in the reverse order, and iterate.

			==The practical QR algorithm==
			Formally, let ''A'' be a real matrix of which we want to compute the eigenvalues, and let ''A''<sub>0</sub>:=''A''. At the ''k''-th step (starting with ''k'' = 0), we compute the QR decomposition ''A''<sub>''k''</sub>=''Q''<sub>''k''</sub>''R''<sub>''k''</sub> where ''Q''<sub>''k''</sub> is an [[orthogonal matrix]] and ''R''<sub>''k''</sub> is an upper triangular matrix. We then form ''A''<sub>''k''+1</sub> = ''R''<sub>''k''</sub>''Q''<sub>''k''</sub>. Note that
			:<math> A_{k+1} = R_k Q_k = Q_k^T Q_k R_k Q_k = Q_k^T A_k Q_k = Q_k^{-1} A_k Q_k, </math>
			so all the ''A''<sub>''k''</sub> are [[Similar matrix\|similar]] and hence they have the same eigenvalues. The algorithm is [[numerical stability\|numerically stable]] because it proceeds by ''orthogonal'' similarity transforms.

			Under certain conditions,<ref>Golub, G. H. and Van Loan, C. F.: Matrix Computations, 3rd ed., Johns Hopkins University Press, Baltimore, 1996, ISBN 0-8018-5414-8.</ref> the matrices ''A''<sub>''k''</sub> converge to a triangular matrix, the [[Schur form]] of ''A''. The eigenvalues of a triangular matrix are listed on the diagonal, and the eigenvalue problem is solved. In testing for convergence it is impractical to require exact zeros, but the [[Gershgorin circle theorem]] provides a bound on the error.

			In this crude form the iterations are relatively expensive. This can be mitigated by first bringing the matrix ''A'' to upper [[Hessenberg form]] (which costs <math>\begin{matrix}\frac{10}{3}\end{matrix} n^3 + O(n^2)</math> arithmetic operations using a technique based on [[Householder transformation\|Householder reduction]]), with a finite sequence of orthogonal similarity transforms, somewhat like a two-sided QR decomposition.<ref name=Demmel>[[James W. Demmel]], ''Applied Numerical Linear Algebra'' (SIAM, 1997).</ref><ref name=Trefethen>[[Lloyd N. Trefethen]] and David Bau, ''Numerical Linear Algebra'' (SIAM, 1997).</ref> (For QR decomposition, the Householder reflectors are multiplied only on the left, but for the Hessenberg case they are multiplied on both left and right.) Determining the QR decomposition of an upper Hessenberg matrix costs <math>6 n^2 + O(n)</math> arithmetic operations. Moreover, because the Hessenberg form is already nearly upper-triangular (it has just one nonzero entry below each diagonal), using it as a starting point reduces the number of steps required for convergence of the QR algorithm.

			If the original matrix is [[symmetric matrix\|symmetric]], then the upper Hessenberg matrix is also symmetric and thus [[tridiagonal matrix\|tridiagonal]], and so are all the ''A''<sub>''k''</sub>. This procedure costs <math>\begin{matrix}\frac{4}{3}\end{matrix} n^3 + O(n^2)</math> arithmetic operations using a technique based on Householder reduction.<ref name=Demmel/><ref name=Trefethen/> Determining the QR decomposition of a symmetric tridiagonal matrix costs <math>O(n)</math> operations.<ref>James M. Ortega and Henry F. Kaiser, "The ''LL<sup>T</sup>'' and ''QR'' methods for symmetric tridiagonal matrices," ''The Computer Journal'' '''6''' (1), 99–101 (1963).</ref>

			The rate of convergence depends on the separation between eigenvalues, so a practical algorithm will use shifts, either explicit or implicit, to increase separation and accelerate convergence. A typical symmetric QR algorithm isolates each eigenvalue (then reduces the size of the matrix) with only one or two iterations, making it efficient as well as robust.{{clarify\|date=June 2012}}

			== The implicit QR algorithm ==
			In modern computational practice,<ref>Golub and Van Loan, chapter 7</ref> the QR algorithm is performed in an implicit version which makes the use of multiple shifts easier to introduce. The matrix is first brought to upper Hessenberg form <math>A_0=QAQ^T</math> as in the explicit version; then, at each step, the first column of <math>A_k</math> is transformed via a small-size Householder similarity transformation to the first column of <math>p(A_k)</math> (or <math>p(A_k)e_1</math>), where <math>p(A_k)</math>, of degree <math>r</math>, is the polynomial that defines the shifting strategy (often <math>p(x)=(x-\lambda)(x-\bar\lambda)</math>, where <math>\lambda</math> and <math>\bar\lambda</math> are the two eigenvalues of the trailing <math>2 \times 2</math> principal submatrix of <math>A_k</math>, the so-called ''implicit double-shift''). Then successive Householder transformation of size <math>r+1</math> are performed in order to return the working matrix <math>A_k</math> to upper Hessenberg form. This operation is known as ''bulge chasing'', due to the peculiar shape of the non-zero entries of the matrix along the steps of the algorithm. As in the first version, deflation is performed as soon as one of the sub-diagonal entries of <math>A_k</math> is sufficiently small.

			===Renaming proposal===
			Since in the modern implicit version of the procedure no [[QR decomposition]]s are explicitly performed, some authors, for instance Watkins,<ref>Watkins, David S.: The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods, SIAM, Philadelphia, PA, USA, 2007, ISBN 0-89871-641-1, ISBN 978-0-89871-641-2</ref> suggested changing its name to ''Francis algorithm''. [[Gene H. Golub\|Golub]] and [[Charles F. Van Loan\|Van Loan]] use the term ''Francis QR step''.

			== Interpretation and convergence ==
			The QR algorithm can be seen as a more sophisticated variation of the basic "power" [[eigenvalue algorithm]]. Recall that the power algorithm repeatedly multiplies ''A'' times a single vector, normalizing after each iteration. The vector converges to an eigenvector of the largest eigenvalue. Instead, the QR algorithm works with a complete basis of vectors, using QR decomposition to renormalize (and orthogonalize). For a symmetric matrix ''A'', upon convergence, ''AQ'' = ''QΛ'', where ''Λ'' is the diagonal matrix of eigenvalues to which ''A'' converged, and where ''Q'' is a composite of all the orthogonal similarity transforms required to get there. Thus the columns of ''Q'' are the eigenvectors.

			== History ==
			The QR algorithm was preceded by the ''LR algorithm'', which uses the [[LU decomposition]] instead of the QR decomposition. The QR algorithm is more stable, so the LR algorithm is rarely used nowadays. However, it represents an important step in the development of the QR algorithm.

			The LR algorithm was developed in the early 1950s by [[Heinz Rutishauer]], who worked at that time as a research assistant of [[Eduard Stiefel]] at [[ETH Zurich]]. Stiefel suggested that Rutishauer use the sequence of moments ''y''<sub>0</sub><sup>T</sup> ''A''<sup>''k''</sup> ''x''<sub>0</sub>, ''k'' = 0, 1, … (where ''x''<sub>0</sub> and ''y''<sub>0</sub> are arbitrary vectors) to find the eigenvalues of ''A''. Rutishauer took an algorithm of [[Alexander Aitken]] for this task and developed it into the ''quotient–difference algorithm'' or ''qd algorithm''. After arranging the computation in a suitable shape, he discovered that the qd algorithm is in fact the iteration ''A''<sub>''k''</sub> = ''L''<sub>''k''</sub>''U''<sub>''k''</sub> (LU decomposition), ''A''<sub>''k''+1</sub> = ''U''<sub>''k''</sub>''L''<sub>''k''</sub>, applied on a tridiagonal matrix, from which the LR algorithm follows.<ref>{{Citation \| last1=Parlett \| first1=Beresford N. \| last2=Gutknecht \| first2=Martin H. \| title=From qd to LR, or, how were the qd and LR algorithms discovered? \| doi=10.1093/imanum/drq003 \| year=2011 \| journal=IMA Journal of Numerical Analysis \| issn=0272-4979 \| volume=31 \| pages=741–754}}</ref>

			== Other variants ==
			One variant of the ''QR algorithm'', ''the Golub-Kahan-Reinsch'' algorithm starts with reducing a general matrix into a bidiagonal one.<ref>Bochkanov Sergey Anatolyevich. ALGLIB User Guide - General Matrix operations - Singular value decomposition . ALGLIB Project. 2010-12-11. URL:http://www.alglib.net/matrixops/general/svd.php. Accessed: 2010-12-11. (Archived by WebCite at http://www.webcitation.org/5utO4iSnR)</ref> This variant of the ''QR algorithm'' for the computation of eigenvalues, which was first described by {{harvtxt\|Golub\|Kahan\|1965}}. The [[LAPACK]] subroutine [http://www.netlib.org/lapack/double/dbdsqr.f DBDSQR] implements this iterative method, with some modifications to cover the case where the singular values are very small {{harv\|Demmel\|Kahan\|1990}}. Together with a first step using Householder reflections and, if appropriate, [[QR decomposition]], this forms the [http://www.netlib.org/lapack/double/dgesvd.f DGESVD] routine for the computation of the [[singular value decomposition]].

			== References ==
			<references/>

			== External links ==
			* {{planetmath reference\|id=1474\|title=Eigenvalue problem}}
			* [http://www.math.umn.edu/~olver/aims_/qr.pdf Prof. Peter Olver's notes on orthogonal bases and the workings of the QR algorithm]
			* [http://math.fullerton.edu/mathews/n2003/QRMethodMod.html Module for the QR Method]

			{{Numerical linear algebra}}

			{{DEFAULTSORT:Qr Algorithm}}
			[[Category:Numerical linear algebra]]

en>JSquish at 03:13, 24 August 2012

2012-08-24T03:13:32Z

New page

My name: Henry Dawkins My age: 33 Country: Sweden Home town: Gullholmen ZIP: 440 84 Street: Villagatan 48 my web site :: [http://ebook-pdfree.com/free-mark-athena-pdf-heroes-olympus-book-3-download/ The Mark of Athena pdf]