Angles Between Subspaces

Total Page:16

File Type:pdf, Size:1020Kb

Angles Between Subspaces Angles between subspaces Eckhard Hitzer, Department of Applied Physics, University of Fukui, 910-8507 Japan June 10, 2013 Abstract We first review the definition of the angle between subspaces and how it is computed using matrix algebra. Then we introduce the Grassmann and Clifford algebra description of subspaces. The geometric product of two subspaces yields the full relative angular information in an explicit manner. We explain and interpret the result of the geometric product of subspaces gaining thus full access to the relative orientation information. Keywords: Clifford geometric algebra, subspaces, relative angle, principal angles, principal vectors. AMS Subj. Class.: 15A66. 1 Introduction lines, planes, circles ans spheres of three dimen- sional Euclidean geometry [7] and was able to find I first came across Clifford’s geometric algebra in the one general formula fully expressing the relative early 90ies in papers on gauge field theory of gravity orientation of any two of these objects. Yet L. by J.S.R. Chisholm, struck by the seamlessly com- Dorst (Amsterdam) later asked me if this formula pact, elegant, and geometrically well interpretable could be generalized to any dimension, because expressions for elementary particle fields subject to by his experience formulas that work dimension Einstein's gravity. Later I became familiar with D. independent are right. I had no immediate answer, Hestenes' excellent modern formulation of geometric it seemed to complicated to me, having to deal with algebra, which explicitly shows how the geometric too many possible cases. product of two vectors encodes their complete rela- But when I prepared for December 2009 a presen- tive orientation in the scalar inner product part (co- tation on neural computation and Clifford algebra, I sine) and in the bivector outer product part (sine). came across a 1983 paper by Per Ake Wedin on an- Geometric algebra can be viewed as an algebra of gles between subspaces of finite dimensional inner a vector space and all its subspaces, represented by product spaces [2], which taught me the classical socalled blades. I therefore often wondered if the approach. In addition it had a very interesting note geometric product of subspace blades also encodes on solving the problem, essentially using Grassmann their complete relative orientation, and how this is algebra with an additional canonically defined inner done? What is the form of the result, how can it product. After that the various bits and pieces came together and began to show the whole picture, the arXiv:1306.1629v1 [math.MG] 7 Jun 2013 be interpreted and put to further use? I learned more about this problem, when I worked on the picture which I want to explain in this contribution. conformal representation of points, point pairs, Permission to make digital or hard copies of all or part of this work for personal or classroom use is 2 The angle between two lines granted without fee provided that copies are not made or distributed for profit or commercial ad- vantage and that copies bear this notice and the To begin with let us look (see Fig. 1) at two lines A; B in a vector space Rn, which are spanned by two full citation on the first page. To copy otherwise, n or republish, to post on servers or to redistribute (unit) vectors a; b 2 R ; a · a = b · b = 1: to lists, requires prior specific permission and/or a fee. A = spanfag; B = spanfbg: (1) Example: 2 Lines tween the subspaces A, B is characterized by a set of r principal angles θk; 1 ≤ k ≤ r, as indicated in b Fig. 2. A principal angle is the angle between two principal vectors a k 2 A and bk 2 B. The spanning sets of vectors fa 1;:::; a rg, and fb1;:::; brg can A,B a be chosen such that pairs of vectors a k; bk either • agree a k = bk, θk = 0, • or enclose a finite angle 0 < θk ≤ π=2. Figure 1:• Two Angle linesθA;B betweengiven b twoy two lines unitA, B vectors, spanned by unit vectors a, b, respectively. Angles– Line between A: spanned Subs by vectorpaces a A,BIn addition the pairs of vectors fa k; bkg; 1 ≤ k ≤ r span mutually orthogonal lines (for θk = 0) and – Line B: spanned by vector b (principal) planes ik (for 0 < θk ≤ π=2). These mu- tually orthogonal planes ik are indicated in Fig. 2. • Angle 0 ≤ A,B ≤ between lines A,B b1 br Therefore if a k , bk and a l , bl for 1 ≤ k 6= l ≤ r, b2 b3cos A,B = a·b then plane ik is orthogonal to plane il. The cosines B of the socalled principal angles θk may therefore 1 r be cos θk = 1 (for a k = bk), or cos θk = 0 (for 2 3 a k ? bk), or any value 0 < cos θk < 1. The total a1 ar angle between the two subspaces A, B is defined as A a the product i1 a2 3 ir i3 cos θA;B = cos θ1 cos θ2 ::: cos θr: (4) i2 In this definition cos θA;B will automatically be zero A = span(a ,a ,a ,…,a ), B=span(b ,b , b , …if, any b ) pair of principal vectors fa k; bkg; 1 ≤ k ≤ r is Figure 2: Angular1 2 relationship3 r of two subspaces1 2A, 3 perpendicular.r Then the two subspaces are said to B, spanned by two sets of vectors fa ;:::; a g, and cos A,B = cos 1 cos 2 cos1 3 … cosr r be perpendicular A ? B, a familiar notion from three fb1;:::; brg, respectively. dimensions, where two perpendicular planes A, B share a common line spanned by a 1 = b1, and have two mutually orthogonal principal vectors a 2 ? b2, The angle 0 ≤ θA;B ≤ π=2 between lines A and B is which are both in turn orthogonal to the common simply given by line vector a 1. It is further possible to choose the indexes of the vector pairs fa k; bkg; 1 ≤ k ≤ r such cos θA;B = a · b: (2) that the principle angles θk appear ordered by mag- nitude θ1 ≥ θ2 ≥ ::: ≥ θr: (5) 3 Angles between two sub- spaces (described by princi- 4 Matrix algebra computation pal vectors) of angle between subspaces Next let us examine the case of two r-dimensional The conventional method of computing the angle (r ≤ n) subspaces A; B of an n-dimensional Eu- θA;B between two r-dimensional subspaces A; B ⊂ n n 0 0 0 clidean vector space R . The situation is depicted R spanned by two sets of vectors fa 1; a 2;::: a rg 0 0 0 in Fig. 2. Each subspace A, B is spanned by a set and fb1; b2;::: brg is to first arrange these vectors of r linearly independent vectors as column vectors into two n × r matrices n 0 0 0 0 A = spanfa 1;:::; a rg ⊂ R ; MA = [a 1;:::; a r];MB = [b1;:::; br]: (6) n B = spanfb1;:::; brg ⊂ R : (3) Then standard matrix algebra methods of QR de- Using Fig. 2 we introduce the following notation composition and singular value decomposition are for principal vectors. The angular relationship be- applied to obtain • r pairs of singular unit vectorsVectorsa k; bk and and Bivectors b a • r singular values σ = cos θ = a · b . oriented k k k k unit area = This approach is very computation1+a2 e2, b intensive.b1e1+b2e2 a b (a1b2 a2b1) e1e2 = |a||b| sin e1e2 5 Even more subtle b∧a …oriented ways area spanned by a,b 0 ⇔ a || b ⇔ b = a , ∊R Per Ake Wedin in his 1983th direction contribution a : {x [2] ∊R ton | ax∧a= 0} conference on Matrix Pencils entitled On Angles between Subspaces of ar Finites can be Dimensional freely reshaped Inner Product Space first carefully treatsa∧ theb = abovea∧(b men-+ a) a tioned matrix algebra approacha product to of computing orthogonal the vectors angle θA;B in great detail and clarity. Towards the b’ end of his paper he dedicatesa∧b less= a thanb’ , onea pageb’ (a to · b’ = 0 ) mentioning an alternative method starting out with a the words: But there are even more subtle ways to define angle functions. Figure 3: Bivectors a ^ b as oriented area elements There he essentially reviews how r-dimensional can be reshaped (e.g. by b ! b + µa; µ 2 R) with- subspaces A; B ⊂ Rn can be represented by r-vectors out changing their value (area and orientation). The (blades) in Grassmann algebra A; B 2 Λ(Rn): bottom figure shows orthogonal reshaping into the form of an oriented rectangle. n A = fx 2 R jx ^ A = 0g; n B = fx 2 R jx ^ B = 0g: (7) n 6 Clifford (geometric) algebra The angle θA;B between the two subspaces A; B 2 R can then be computed in a single step Clifford (geometric) algebra is based on the geomet- ric product of vectors a; b 2 p;q; p + q = n A · Be R cos θA;B = = cos θ1 cos θ2 ::: cos θr; (8) jAjjBj ab = a · b + a ^ b; (9) where the inner product is canonically defined on and the associative algebra Cl thus generated the Grassmann algebra Λ( n) corresponding to the p;q R with and p;q as subspaces of Cl . a · b is geometry of n. The tilde operation is the reverse R R p;q R the symmetric inner product of vectors and a ^ b operation representing a dimension dependent sign r(r−1) is Grassmann's outer product of vectors represent- change Be = (−1) 2 B, and jAj represents the ing the oriented parallelogram area spanned by a; b, 2 norm of blade A, i.e.
Recommended publications
  • Quaternions and Cli Ord Geometric Algebras
    Quaternions and Cliord Geometric Algebras Robert Benjamin Easter First Draft Edition (v1) (c) copyright 2015, Robert Benjamin Easter, all rights reserved. Preface As a rst rough draft that has been put together very quickly, this book is likely to contain errata and disorganization. The references list and inline citations are very incompete, so the reader should search around for more references. I do not claim to be the inventor of any of the mathematics found here. However, some parts of this book may be considered new in some sense and were in small parts my own original research. Much of the contents was originally written by me as contributions to a web encyclopedia project just for fun, but for various reasons was inappropriate in an encyclopedic volume. I did not originally intend to write this book. This is not a dissertation, nor did its development receive any funding or proper peer review. I oer this free book to the public, such as it is, in the hope it could be helpful to an interested reader. June 19, 2015 - Robert B. Easter. (v1) [email protected] 3 Table of contents Preface . 3 List of gures . 9 1 Quaternion Algebra . 11 1.1 The Quaternion Formula . 11 1.2 The Scalar and Vector Parts . 15 1.3 The Quaternion Product . 16 1.4 The Dot Product . 16 1.5 The Cross Product . 17 1.6 Conjugates . 18 1.7 Tensor or Magnitude . 20 1.8 Versors . 20 1.9 Biradials . 22 1.10 Quaternion Identities . 23 1.11 The Biradial b/a .
    [Show full text]
  • Introduction to Linear Bialgebra
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of New Mexico University of New Mexico UNM Digital Repository Mathematics and Statistics Faculty and Staff Publications Academic Department Resources 2005 INTRODUCTION TO LINEAR BIALGEBRA Florentin Smarandache University of New Mexico, [email protected] W.B. Vasantha Kandasamy K. Ilanthenral Follow this and additional works at: https://digitalrepository.unm.edu/math_fsp Part of the Algebra Commons, Analysis Commons, Discrete Mathematics and Combinatorics Commons, and the Other Mathematics Commons Recommended Citation Smarandache, Florentin; W.B. Vasantha Kandasamy; and K. Ilanthenral. "INTRODUCTION TO LINEAR BIALGEBRA." (2005). https://digitalrepository.unm.edu/math_fsp/232 This Book is brought to you for free and open access by the Academic Department Resources at UNM Digital Repository. It has been accepted for inclusion in Mathematics and Statistics Faculty and Staff Publications by an authorized administrator of UNM Digital Repository. For more information, please contact [email protected], [email protected], [email protected]. INTRODUCTION TO LINEAR BIALGEBRA W. B. Vasantha Kandasamy Department of Mathematics Indian Institute of Technology, Madras Chennai – 600036, India e-mail: [email protected] web: http://mat.iitm.ac.in/~wbv Florentin Smarandache Department of Mathematics University of New Mexico Gallup, NM 87301, USA e-mail: [email protected] K. Ilanthenral Editor, Maths Tiger, Quarterly Journal Flat No.11, Mayura Park, 16, Kazhikundram Main Road, Tharamani, Chennai – 600 113, India e-mail: [email protected] HEXIS Phoenix, Arizona 2005 1 This book can be ordered in a paper bound reprint from: Books on Demand ProQuest Information & Learning (University of Microfilm International) 300 N.
    [Show full text]
  • Multivector Differentiation and Linear Algebra 0.5Cm 17Th Santaló
    Multivector differentiation and Linear Algebra 17th Santalo´ Summer School 2016, Santander Joan Lasenby Signal Processing Group, Engineering Department, Cambridge, UK and Trinity College Cambridge [email protected], www-sigproc.eng.cam.ac.uk/ s jl 23 August 2016 1 / 78 Examples of differentiation wrt multivectors. Linear Algebra: matrices and tensors as linear functions mapping between elements of the algebra. Functional Differentiation: very briefly... Summary Overview The Multivector Derivative. 2 / 78 Linear Algebra: matrices and tensors as linear functions mapping between elements of the algebra. Functional Differentiation: very briefly... Summary Overview The Multivector Derivative. Examples of differentiation wrt multivectors. 3 / 78 Functional Differentiation: very briefly... Summary Overview The Multivector Derivative. Examples of differentiation wrt multivectors. Linear Algebra: matrices and tensors as linear functions mapping between elements of the algebra. 4 / 78 Summary Overview The Multivector Derivative. Examples of differentiation wrt multivectors. Linear Algebra: matrices and tensors as linear functions mapping between elements of the algebra. Functional Differentiation: very briefly... 5 / 78 Overview The Multivector Derivative. Examples of differentiation wrt multivectors. Linear Algebra: matrices and tensors as linear functions mapping between elements of the algebra. Functional Differentiation: very briefly... Summary 6 / 78 We now want to generalise this idea to enable us to find the derivative of F(X), in the A ‘direction’ – where X is a general mixed grade multivector (so F(X) is a general multivector valued function of X). Let us use ∗ to denote taking the scalar part, ie P ∗ Q ≡ hPQi. Then, provided A has same grades as X, it makes sense to define: F(X + tA) − F(X) A ∗ ¶XF(X) = lim t!0 t The Multivector Derivative Recall our definition of the directional derivative in the a direction F(x + ea) − F(x) a·r F(x) = lim e!0 e 7 / 78 Let us use ∗ to denote taking the scalar part, ie P ∗ Q ≡ hPQi.
    [Show full text]
  • Università Degli Studi Di Trieste a Gentle Introduction to Clifford Algebra
    Università degli Studi di Trieste Dipartimento di Fisica Corso di Studi in Fisica Tesi di Laurea Triennale A Gentle Introduction to Clifford Algebra Laureando: Relatore: Daniele Ceravolo prof. Marco Budinich ANNO ACCADEMICO 2015–2016 Contents 1 Introduction 3 1.1 Brief Historical Sketch . 4 2 Heuristic Development of Clifford Algebra 9 2.1 Geometric Product . 9 2.2 Bivectors . 10 2.3 Grading and Blade . 11 2.4 Multivector Algebra . 13 2.4.1 Pseudoscalar and Hodge Duality . 14 2.4.2 Basis and Reciprocal Frames . 14 2.5 Clifford Algebra of the Plane . 15 2.5.1 Relation with Complex Numbers . 16 2.6 Clifford Algebra of Space . 17 2.6.1 Pauli algebra . 18 2.6.2 Relation with Quaternions . 19 2.7 Reflections . 19 2.7.1 Cross Product . 21 2.8 Rotations . 21 2.9 Differentiation . 23 2.9.1 Multivectorial Derivative . 24 2.9.2 Spacetime Derivative . 25 3 Spacetime Algebra 27 3.1 Spacetime Bivectors and Pseudoscalar . 28 3.2 Spacetime Frames . 28 3.3 Relative Vectors . 29 3.4 Even Subalgebra . 29 3.5 Relative Velocity . 30 3.6 Momentum and Wave Vectors . 31 3.7 Lorentz Transformations . 32 3.7.1 Addition of Velocities . 34 1 2 CONTENTS 3.7.2 The Lorentz Group . 34 3.8 Relativistic Visualization . 36 4 Electromagnetism in Clifford Algebra 39 4.1 The Vector Potential . 40 4.2 Electromagnetic Field Strength . 41 4.3 Free Fields . 44 5 Conclusions 47 5.1 Acknowledgements . 48 Chapter 1 Introduction The aim of this thesis is to show how an approach to classical and relativistic physics based on Clifford algebras can shed light on some hidden geometric meanings in our models.
    [Show full text]
  • The Making of a Geometric Algebra Package in Matlab Computer Science Department University of Waterloo Research Report CS-99-27
    The Making of a Geometric Algebra Package in Matlab Computer Science Department University of Waterloo Research Report CS-99-27 Stephen Mann∗, Leo Dorsty, and Tim Boumay [email protected], [email protected], [email protected] Abstract In this paper, we describe our development of GABLE, a Matlab implementation of the Geometric Algebra based on `p;q (where p + q = 3) and intended for tutorial purposes. Of particular note are the C matrix representation of geometric objects, effective algorithms for this geometry (inversion, meet and join), and issues in efficiency and numerics. 1 Introduction Geometric algebra extends Clifford algebra with geometrically meaningful operators, and its purpose is to facilitate geometrical computations. Present textbooks and implementation do not always convey this geometrical flavor or the computational and representational convenience of geometric algebra, so we felt a need for a computer tutorial in which representation, computation and visualization are combined to convey both the intuition and the techniques of geometric algebra. Current software packages are either Clifford algebra only (such as CLICAL [9] and CLIFFORD [1]) or do not include graphics [6], so we decide to build our own. The result is GABLE (Geometric Algebra Learning Environment) a hands-on tutorial on geometric algebra that should be accessible to the second year student in college [3]. The GABLE tutorial explains the basics of Geometric Algebra in an accessible manner. It starts with the outer product (as a constructor of subspaces), then treats the inner product (for perpendilarity), and moves via the geometric product (for invertibility) to the more geometrical operators such as projection, rotors, meet and join, and end with to the homogeneous model of Euclidean space.
    [Show full text]
  • Introduction to Clifford's Geometric Algebra
    解説:特集 コンピューテーショナル・インテリジェンスの新展開 —クリフォード代数表現など高次元表現を中心として— Introduction to Clifford’s Geometric Algebra Eckhard HITZER* *Department of Applied Physics, University of Fukui, Fukui, Japan *E-mail: [email protected] Key Words:hypercomplex algebra, hypercomplex analysis, geometry, science, engineering. JL 0004/12/5104–0338 C 2012 SICE erty 15), 21) that any isometry4 from the vector space V Abstract into an inner-product algebra5 A over the field6 K can Geometric algebra was initiated by W.K. Clifford over be uniquely extended to an isometry7 from the Clifford 130 years ago. It unifies all branches of physics, and algebra Cl(V )intoA. The Clifford algebra Cl(V )is has found rich applications in robotics, signal process- the unique associative and multilinear algebra with this ing, ray tracing, virtual reality, computer vision, vec- property. Thus if we wish to generalize methods from tor field processing, tracking, geographic information algebra, analysis, calculus, differential geometry (etc.) systems and neural computing. This tutorial explains of real numbers, complex numbers, and quaternion al- the basics of geometric algebra, with concrete exam- gebra to vector spaces and multivector spaces (which ples of the plane, of 3D space, of spacetime, and the include additional elements representing 2D up to nD popular conformal model. Geometric algebras are ideal subspaces, i.e. plane elements up to hypervolume ele- to represent geometric transformations in the general ments), the study of Clifford algebras becomes unavoid- framework of Clifford groups (also called versor or Lip- able. Indeed repeatedly and independently a long list of schitz groups). Geometric (algebra based) calculus al- Clifford algebras, their subalgebras and in Clifford al- lows e.g., to optimize learning algorithms of Clifford gebras embedded algebras (like octonions 17))ofmany neurons, etc.
    [Show full text]
  • A Bit About Hilbert Spaces
    A Bit About Hilbert Spaces David Rosenberg New York University October 29, 2016 David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 1 / 9 Inner Product Space (or “Pre-Hilbert” Spaces) An inner product space (over reals) is a vector space V and an inner product, which is a mapping h·,·i : V × V ! R that has the following properties 8x,y,z 2 V and a,b 2 R: Symmetry: hx,yi = hy,xi Linearity: hax + by,zi = ahx,zi + b hy,zi Postive-definiteness: hx,xi > 0 and hx,xi = 0 () x = 0. David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 2 / 9 Norm from Inner Product For an inner product space, we define a norm as p kxk = hx,xi. Example Rd with standard Euclidean inner product is an inner product space: hx,yi := xT y 8x,y 2 Rd . Norm is p kxk = xT y. David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 3 / 9 What norms can we get from an inner product? Theorem (Parallelogram Law) A norm kvk can be generated by an inner product on V iff 8x,y 2 V 2kxk2 + 2kyk2 = kx + yk2 + kx - yk2, and if it can, the inner product is given by the polarization identity jjxjj2 + jjyjj2 - jjx - yjj2 hx,yi = . 2 Example d `1 norm on R is NOT generated by an inner product. [Exercise] d Is `2 norm on R generated by an inner product? David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 4 / 9 Pythagorean Theroem Definition Two vectors are orthogonal if hx,yi = 0.
    [Show full text]
  • True False Questions from 4,5,6
    (c)2015 UM Math Dept licensed under a Creative Commons By-NC-SA 4.0 International License. Math 217: True False Practice Professor Karen Smith 1. For every 2 × 2 matrix A, there exists a 2 × 2 matrix B such that det(A + B) 6= det A + det B. Solution note: False! Take A to be the zero matrix for a counterexample. 2. If the change of basis matrix SA!B = ~e4 ~e3 ~e2 ~e1 , then the elements of A are the same as the element of B, but in a different order. Solution note: True. The matrix tells us that the first element of A is the fourth element of B, the second element of basis A is the third element of B, the third element of basis A is the second element of B, and the the fourth element of basis A is the first element of B. 6 6 3. Every linear transformation R ! R is an isomorphism. Solution note: False. The zero map is a a counter example. 4. If matrix B is obtained from A by swapping two rows and det A < det B, then A is invertible. Solution note: True: we only need to know det A 6= 0. But if det A = 0, then det B = − det A = 0, so det A < det B does not hold. 5. If two n × n matrices have the same determinant, they are similar. Solution note: False: The only matrix similar to the zero matrix is itself. But many 1 0 matrices besides the zero matrix have determinant zero, such as : 0 0 6.
    [Show full text]
  • A Guided Tour to the Plane-Based Geometric Algebra PGA
    A Guided Tour to the Plane-Based Geometric Algebra PGA Leo Dorst University of Amsterdam Version 1.15{ July 6, 2020 Planes are the primitive elements for the constructions of objects and oper- ators in Euclidean geometry. Triangulated meshes are built from them, and reflections in multiple planes are a mathematically pure way to construct Euclidean motions. A geometric algebra based on planes is therefore a natural choice to unify objects and operators for Euclidean geometry. The usual claims of `com- pleteness' of the GA approach leads us to hope that it might contain, in a single framework, all representations ever designed for Euclidean geometry - including normal vectors, directions as points at infinity, Pl¨ucker coordinates for lines, quaternions as 3D rotations around the origin, and dual quaternions for rigid body motions; and even spinors. This text provides a guided tour to this algebra of planes PGA. It indeed shows how all such computationally efficient methods are incorporated and related. We will see how the PGA elements naturally group into blocks of four coordinates in an implementation, and how this more complete under- standing of the embedding suggests some handy choices to avoid extraneous computations. In the unified PGA framework, one never switches between efficient representations for subtasks, and this obviously saves any time spent on data conversions. Relative to other treatments of PGA, this text is rather light on the mathematics. Where you see careful derivations, they involve the aspects of orientation and magnitude. These features have been neglected by authors focussing on the mathematical beauty of the projective nature of the algebra.
    [Show full text]
  • The Determinant Inner Product and the Heisenberg Product of Sym(2)
    INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 14 NO. 1 PAGE 1–12 (2021) DOI: HTTPS://DOI.ORG/10.32323/IEJGEO.602178 The determinant inner product and the Heisenberg product of Sym(2) Mircea Crasmareanu (Dedicated to the memory of Prof. Dr. Aurel BEJANCU (1946 - 2020)Cihan Ozgur) ABSTRACT The aim of this work is to introduce and study the nondegenerate inner product < ·; · >det induced by the determinant map on the space Sym(2) of symmetric 2 × 2 real matrices. This symmetric bilinear form of index 2 defines a rational symmetric function on the pairs of rays in the plane and an associated function on the 2-torus can be expressed with the usual Hopf bundle projection 3 2 1 S ! S ( 2 ). Also, the product < ·; · >det is treated with complex numbers by using the Hopf invariant map of Sym(2) and this complex approach yields a Heisenberg product on Sym(2). Moreover, the quadratic equation of critical points for a rational Morse function of height type generates a cosymplectic structure on Sym(2) with the unitary matrix as associated Reeb vector and with the Reeb 1-form being half of the trace map. Keywords: symmetric matrix; determinant; Hopf bundle; Hopf invariant AMS Subject Classification (2020): Primary: 15A15 ; Secondary: 15A24; 30C10; 22E47. 1. Introduction This short paper, whose nature is mainly expository, is dedicated to the memory of Professor Dr. Aurel Bejancu, who was equally a gifted mathematician and a dedicated professor. As we hereby present an in memoriam paper, we feel it would be better if we choose an informal presentation rather than a classical structured discourse, where theorems and propositions are followed by their proofs.
    [Show full text]
  • Inner Product Spaces
    CHAPTER 6 Woman teaching geometry, from a fourteenth-century edition of Euclid’s geometry book. Inner Product Spaces In making the definition of a vector space, we generalized the linear structure (addition and scalar multiplication) of R2 and R3. We ignored other important features, such as the notions of length and angle. These ideas are embedded in the concept we now investigate, inner products. Our standing assumptions are as follows: 6.1 Notation F, V F denotes R or C. V denotes a vector space over F. LEARNING OBJECTIVES FOR THIS CHAPTER Cauchy–Schwarz Inequality Gram–Schmidt Procedure linear functionals on inner product spaces calculating minimum distance to a subspace Linear Algebra Done Right, third edition, by Sheldon Axler 164 CHAPTER 6 Inner Product Spaces 6.A Inner Products and Norms Inner Products To motivate the concept of inner prod- 2 3 x1 , x 2 uct, think of vectors in R and R as x arrows with initial point at the origin. x R2 R3 H L The length of a vector in or is called the norm of x, denoted x . 2 k k Thus for x .x1; x2/ R , we have The length of this vector x is p D2 2 2 x x1 x2 . p 2 2 x1 x2 . k k D C 3 C Similarly, if x .x1; x2; x3/ R , p 2D 2 2 2 then x x1 x2 x3 . k k D C C Even though we cannot draw pictures in higher dimensions, the gener- n n alization to R is obvious: we define the norm of x .x1; : : : ; xn/ R D 2 by p 2 2 x x1 xn : k k D C C The norm is not linear on Rn.
    [Show full text]
  • Geometric Algebra Techniques for General Relativity
    Geometric Algebra Techniques for General Relativity Matthew R. Francis∗ and Arthur Kosowsky† Dept. of Physics and Astronomy, Rutgers University 136 Frelinghuysen Road, Piscataway, NJ 08854 (Dated: February 4, 2008) Geometric (Clifford) algebra provides an efficient mathematical language for describing physical problems. We formulate general relativity in this language. The resulting formalism combines the efficiency of differential forms with the straightforwardness of coordinate methods. We focus our attention on orthonormal frames and the associated connection bivector, using them to find the Schwarzschild and Kerr solutions, along with a detailed exposition of the Petrov types for the Weyl tensor. PACS numbers: 02.40.-k; 04.20.Cv Keywords: General relativity; Clifford algebras; solution techniques I. INTRODUCTION Geometric (or Clifford) algebra provides a simple and natural language for describing geometric concepts, a point which has been argued persuasively by Hestenes [1] and Lounesto [2] among many others. Geometric algebra (GA) unifies many other mathematical formalisms describing specific aspects of geometry, including complex variables, matrix algebra, projective geometry, and differential geometry. Gravitation, which is usually viewed as a geometric theory, is a natural candidate for translation into the language of geometric algebra. This has been done for some aspects of gravitational theory; notably, Hestenes and Sobczyk have shown how geometric algebra greatly simplifies certain calculations involving the curvature tensor and provides techniques for classifying the Weyl tensor [3, 4]. Lasenby, Doran, and Gull [5] have also discussed gravitation using geometric algebra via a reformulation in terms of a gauge principle. In this paper, we formulate standard general relativity in terms of geometric algebra. A comprehensive overview like the one presented here has not previously appeared in the literature, although unpublished works of Hestenes and of Doran take significant steps in this direction.
    [Show full text]