<?xml version="1.0" encoding="UTF-8" standalone="yes"?><oembed><version><![CDATA[1.0]]></version><provider_name><![CDATA[Geometry and the imagination]]></provider_name><provider_url><![CDATA[https://lamington.wordpress.com]]></provider_url><author_name><![CDATA[Danny Calegari]]></author_name><author_url><![CDATA[https://lamington.wordpress.com/author/dannycaltech/]]></author_url><title><![CDATA[Ellipsoids and KAK]]></title><type><![CDATA[link]]></type><html><![CDATA[<p>As many readers are no doubt aware, the title of this blog comes from the famous book <a href="http://www.amazon.com/Geometry-Imagination-CHEL-Chelsea-Publishing/dp/0821819984/ref=sr_1_1?ie=UTF8&amp;s=books&amp;qid=1244736666&amp;sr=1-1"><em>Geometry and the Imagination</em></a> by Hilbert and Cohn-Vossen (based on lectures given by Hilbert). One of the first things discussed in that book is the geometry of conics, especially in two and three dimensions. An <em>ellipsoid</em> is a certain kind of (real) quadric surface, i.e. a surface in <img src="https://s0.wp.com/latex.php?latex=%5Cmathbb%7BR%7D%5En&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;mathbb{R}^n" title="&#92;mathbb{R}^n" class="latex" /> defined by a single quadratic equation of the co-ordinates. It may also be defined as the image of the unit <img src="https://s0.wp.com/latex.php?latex=%28n-1%29&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="(n-1)" title="(n-1)" class="latex" />-dimensional sphere under an affine self-map of <img src="https://s0.wp.com/latex.php?latex=%5Cmathbb%7BR%7D%5En&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;mathbb{R}^n" title="&#92;mathbb{R}^n" class="latex" />. After composing with a translation, one may imagine an ellipsoid centered at the origin, and think of it as the image of the unit sphere under a <em>linear</em> automorphism of <img src="https://s0.wp.com/latex.php?latex=%5Cmathbb%7BR%7D%5En&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;mathbb{R}^n" title="&#92;mathbb{R}^n" class="latex" /> &#8212; i.e. transformation by a nonsingular matrix <img src="https://s0.wp.com/latex.php?latex=M&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="M" title="M" class="latex" />.</p>
<p>A (generic) ellipsoid has <img src="https://s0.wp.com/latex.php?latex=n&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="n" title="n" class="latex" /> axes; in dimension three, these are the &#8220;major axis&#8221;, the &#8220;minor axis&#8221; and the &#8220;mean axis&#8221;. Distance to the origin is a Morse function on a generic ellipsoid; the symmetry of an ellipsoid under the antipodal map means that critical points occur in antipodal pairs. There are a pair of critical points of each index between <img src="https://s0.wp.com/latex.php?latex=0&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="0" title="0" class="latex" /> and <img src="https://s0.wp.com/latex.php?latex=n&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="n" title="n" class="latex" />. There is a gradient flow line of this Morse function between each pair of critical points whose index differs by <img src="https://s0.wp.com/latex.php?latex=1&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="1" title="1" class="latex" />, and the union of these flowlines are the (<img src="https://s0.wp.com/latex.php?latex=2&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="2" title="2" class="latex" />-dimensional) ellipse obtained by intersecting the ellipsoid with the plane spanned by the pair of axes in question. This shows that these axes are mutually perpendicular.</p>
<p>One may use this geometric picture to &#8220;see&#8221; the <img src="https://s0.wp.com/latex.php?latex=KAK&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="KAK" title="KAK" class="latex" /> decomposition of <img src="https://s0.wp.com/latex.php?latex=%5Ctext%7BGL%7D%28n%2C%5Cmathbb%7BR%7D%29&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;text{GL}(n,&#92;mathbb{R})" title="&#92;text{GL}(n,&#92;mathbb{R})" class="latex" /> as follows, where <img src="https://s0.wp.com/latex.php?latex=K&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="K" title="K" class="latex" /> denotes the orthogonal subgroup <img src="https://s0.wp.com/latex.php?latex=%5Ctext%7BO%7D%28n%2C%5Cmathbb%7BR%7D%29&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;text{O}(n,&#92;mathbb{R})" title="&#92;text{O}(n,&#92;mathbb{R})" class="latex" />, and <img src="https://s0.wp.com/latex.php?latex=A&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="A" title="A" class="latex" /> denotes the subgroup of diagonal matrices with positive entries. Let <img src="https://s0.wp.com/latex.php?latex=M&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="M" title="M" class="latex" /> be a linear map of <img src="https://s0.wp.com/latex.php?latex=%5Cmathbb%7BR%7D%5En&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;mathbb{R}^n" title="&#92;mathbb{R}^n" class="latex" />, and let <img src="https://s0.wp.com/latex.php?latex=E&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="E" title="E" class="latex" /> be the ellipsoid which is the image of the unit sphere under <img src="https://s0.wp.com/latex.php?latex=M&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="M" title="M" class="latex" />. Let <img src="https://s0.wp.com/latex.php?latex=%5Cxi_i&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;xi_i" title="&#92;xi_i" class="latex" /> be the axes of <img src="https://s0.wp.com/latex.php?latex=E&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="E" title="E" class="latex" /> of index <img src="https://s0.wp.com/latex.php?latex=i&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="i" title="i" class="latex" />. There is a unique orthogonal matrix <img src="https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="O" title="O" class="latex" /> taking the <img src="https://s0.wp.com/latex.php?latex=%5Cxi_i&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;xi_i" title="&#92;xi_i" class="latex" /> to the co-ordinate axes. There is a unique diagonal matrix <img src="https://s0.wp.com/latex.php?latex=D&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="D" title="D" class="latex" /> taking <img src="https://s0.wp.com/latex.php?latex=O%28E%29&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="O(E)" title="O(E)" class="latex" /> to the round sphere. Hence the composition <img src="https://s0.wp.com/latex.php?latex=ODM&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="ODM" title="ODM" class="latex" /> is orthogonal, and we can express <img src="https://s0.wp.com/latex.php?latex=M&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="M" title="M" class="latex" /> as a product of an orthogonal matrix, a diagonal matrix, and another orthogonal matrix.</p>
<p>One can use ellipsoids to visualize another less standard matrix decomposition as follows. For simplicity we concentrate on the case of dimension <img src="https://s0.wp.com/latex.php?latex=3&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="3" title="3" class="latex" />. The minor and mean axis span a plane <img src="https://s0.wp.com/latex.php?latex=%5Cpi&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;pi" title="&#92;pi" class="latex" /> which intersects the ellipsoid in the &#8220;smallest&#8221; possible ellipse. Rotate this plane by keeping the mean axis fixed, and tilting the minor axis towards the major axis. At some unique point one obtains a plane <img src="https://s0.wp.com/latex.php?latex=%5Cpi%27&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="&#92;pi&#039;" title="&#92;pi&#039;" class="latex" /> that intersects the ellipsoid in a round circle. One may shear the ellipsoid, keeping this plane fixed, into an ellipsoid of rotation. This describes a way to factorize <img src="https://s0.wp.com/latex.php?latex=M&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="M" title="M" class="latex" /> as a product of a shear, a diagonal matrix with two equal eigenvalues, and a rotation.</p>
<p><strong>Question:</strong> What is the generalization of the &#8220;shear, dilate, rotate&#8221; factorization in higher dimensions?</p>
<p><strong>Question:</strong> Is there a way to see the Iwasawa (<img src="https://s0.wp.com/latex.php?latex=KAN&#038;bg=ffffff&#038;fg=000&#038;s=0" alt="KAN" title="KAN" class="latex" />) decomposition geometrically, by using ellipsoids?</p>
]]></html></oembed>