Geometric Algebra (GA)
Total Page:16
File Type:pdf, Size:1020Kb
Geometric Algebra (GA) Werner Benger, 2007 CCT@LSU SciViz 1 Abstract • Geometric Algebra (GA) denotes the re-discovery and geometrical interpretation of the Clifford algebra applied to real fields. Hereby the so-called „geometrical product“ allows to expand linear algebra (as used in vector calculus in 3D) by an invertible operation to multiply and divide vectors. In two dimenions, the geometric algebra can be interpreted as the algebra of complex numbers. In extends in a natural way into three dimensions and corresponds to the well-known quaternions there, which are widely used to describe rotations in 3D as an alternative superior to matrix calculus. However, in contrast to quaternions, GA comes with a direct geometrical interpretation of the respective operations and allows a much finer differentation among the involved objects than is achieveable via quaternions. Moreover, the formalism of GA is independent from the dimension of space. For instance, rotations and reflections of objects of arbitrary dimensions can be easily described intuitively and generic in spaces of arbitrary higher dimensions. • Due to the elegance of the GA and its wide applicabililty it is sometimes denoted as a new „fundamental language of mathematics“. Its unified formalism covers domains such as differential geometry (relativity theory), quantum mechanics, robotics and last but not least computer graphics in a natural way. • This talk will present the basics of Geometric Algebra and specifically emphasizes on the visualization of its elementary operations. Furthermore, the potential of GA will be demonstrated via usage in various application domains. 2 Motivation of GA • Unification of many domains: quantum mechanics, computer graphics, general relativity, robotics… • Completing algebraic operations on vectors • Unified concept for geometry and algebra • Superior formalism for rotations in arbitrary dimensions • Explicit geometrical interpretation of the involved objects and operations on them 3 Definition: “Algebra” • Vector space V over field K with multiplication “ ” • Null-element, One-element, Inverse • Commutative? a b = b a • Associative? (a b) c = a (b c) • Division algebra? • a≠0 a-1 such that a a-1 = 1 = a-1 a • Alternatively: a b=0 a=0 or b=0 4 Historical Roots • Complex Plane (Gauss ~1800) • Real/Imaginary part: a+ib where i2= -1 • Associative, commutative division-algebra • Polar representation: r ei = r ( cos + i sin ) • Multiplication corresponds to rotation in the plane i sin cos 5 Historical Roots, II William Rowan Hamilton (1805-65) invents Quaternions (1844): – Generalization of complex numbers: • 4 components, non commutative: ab ba (in general) • Basic idea: ii=jj=kk= ijk = -1 • Alternative to younger vector- and matrix algebra (Josiah Willard Gibbs, 1839-1903) • p=(p,p), q=(q,q), p q=(pq - p q , pq + pq + p q) • rotation in R3 are around axis of the vector component: v’ = q v q-1 6 Historical Roots, III • Construction by Cayley-Dickson (a,b)(c,d) = (ac-d *b, *a d+cb) – hypercomplex numbers: • octaves/octonions (8 components) • sedenions/hexadekanions (16 components) • … – incremental loss of • commutativity (quaternions,…) • associativity (octonions,…) • division algebra (sedenions,…) 7 Renaissance of the GA 1878: Clifford introduces “geometric algebra”, but dies at age 34 superseded by Gibb’s vector calculus 1920er: Renaissance in quantum mechanics (Pauli, Dirac) algebra on complex fields no geometrical interpretation 1966-2005 David Hestenes (Arizona State University) revives the geometrical interpretation 1997: Gravitation theory using GA (Lasenby, Doran, Gull; Cambridge) 2001: Geometric Algebra at SIGGRAPH (L. Dorst, S. Mann) 8 9 Geometry and Vectors • Geometric interpretation of a vector – Directed line segment or tangent • Vector-algebra in Euclidean Geometry or Tp(M) • Addition / subtraction of vectors a+b • Multiplication / division by scalars a • Multiplication / Division of vectors?? Multiplication of vectors 10 Complete Vector-algebra? • Invertible product of vectors? • What means vector-division “a/b” ? • ab=C b=a-1C • Note: C not necessarily a vector! • Inner product (not associative): a b Skalar – Not invertible e.g. a b =0 with a≠0, b≠0 but orthogonal • Outer product (associative): a b Bivektor – Generalized cross-product from 3D: a b – Not invertible e.g. a b =0 with a≠0, b≠0 but parallel Multiplikation von Vectoren 11 Bivector a b Describes the plane spun by a and b, sign is orientation a b b a = -a b Defined in arbitrary dimensions, anti-symmetric ( not commutative), associative, distributive, spans a vector space, does not require additional structures Multiplikation von Vektoren 12 Constructing Bivectors No unique determination of the generating vectors possible = = a b = (a+λb) b b b b =0 Basis-element |a| |b| sin a+λb Multiplikation von Vektoren 13 Bivectors in R3 • 3 Basis-elements ex ey, ey ez, ez ex • Generalization: ex ey ez is a volume Multiplikation von Vektoren 14 Vectorspace of Bivectors Linear combinations possible e.g.: ex ey, ez ex Multiplikation von Vektoren 15 Coordinate representation of “ ”-product in R3 • Generic Bivector: A = Axy ex ey + Ayz ey ez + Azx ez ex • (axex + ayey + azez) (bxex + byey + bzez)= axex bxex + axex byey + axex bzez + ayey bxex + ayey byey + ayey bzez + azez bxex + azez byey + azez bzez = (axby - aybx)exy+(aybz-azby)eyz+(axbz- azbx)exz Multiplikation von Vektoren 16 Inner product a b • Describes projections a b = |a| |b| cos = b a Symmetric (commutative), requires quadratic form (Metric) as additional structure, not associative (a b) c a (b c) Multiplikation von Vektoren 17 Comparing the products • Inner product • Outer product – Not associative – Associative • (a b) c ≠ a (b c) • (a b) c= a (b c) – Commutative – Not commutative • a b = b a • a b ≠ b a – Not invertible – Not invertible – Yields a scalar – Yields a bivector 18 Geometric Product 1. Requirements and definition 2. Structure of the operands 3. Calculus using GP 4. Rotations using GP Das Geometrische Produkt 19 Requirements to GP • For elements A,B,C of a vector space with quadratic form Q(v) [i.e. a metric g(u,v) = Q(u+v) - Q(u) – Q(v)] we require: 1. Associative: (AB)C = A(BC) 2. Left-distributive: A(B+C) = AB+AC 3. Right-distributive: (B+C)A= BA+CA 4. Scalar product: A2 = Q(A) = |A|2 Das Geometrische Produkt 20 Properties of the GP • Right-angled triangle |a+b|2 = |a|2+|b|2 (A+B)(A+B) = AA+BA+AB+BB = A2 + B2 AB = -BA for A B = 0 anti-symm if orthogonal • However: not purely anti-symmetric |AB|2 =|A|2 |B|2 for A B = 0 (i.e. A,B parallel: B= A) Das Geometrische Produkt 21 Geometric Product • William Kingdon Clifford (1845-79): • Combine inner and outer product to defined the geometric product AB (1878): AB := A B A B • Result is not a vector, but the sum of a scalar + bivector! • Operates on “multivectors” • Subset of the tensoralgebra • Geometric Product is invertible! Das Geometrische Produkt 22 Multi-vector components 2 • R : A = A0 + A1 e0 + A2 e1 + A3 e0 e1 2.7819… + • R3: A = A0 + + + + A1 e0 + A2 e1 + A3 e2 + + + A4 e0 e1+A5 e1 e2+A6 e0 e2 + A e e e 7 0 1 2 + Struktur von Multivektoren 23 Structure of Multi-vectors Linear combination of anti-symmetric basis elements 2n components 0D 1 Scalar 1D 1 Scalar, 1 Vector 2D 1 Scalar, 2 Vectors, 1 Bivector 3D 1 Scalar, 3 Vectors, 3 Bivectors, 1 Volume 4D 1 Scalar, 4 Vectors, 6 Bivectors, 4 Volume, 1 Hyper-volume 5D … Struktur von Multivektoren 24 Inversion • Given vectors a,b: a b = ½ (ab + ba) symmetric part a b = ½ (ab - ba) anti-symmetric part a b = -(a b) (ex ey ez) Dual in 3D Rechnen mit Multivektoren 25 Reflection at a Vector • Unit vector n, arbitrary vector v Vector v projected to n: v║=(v n) n Reflected vector w = v┴ – v║ = v – 2v║ thus w = v – 2(v n) n with GP w = v – 2[½(vn+nv) ] n = v – vnn – nvn w = -nvn Rechnen mit Multivektoren 26 Rotations 1. Identification with Quaternions 2. Rotation in 2D 3. Rotation in nD 4. Rotation of arbitrary Multivectors in nD Rotation 27 Geometrical Quadrate Consider (AB)2 of Bivector-basis element where |A|=1, |B|=1, A B = 0 AB=A B=-BA Basiselement (AB)2 = (AB) (AB) = -(AB) (BA)=-A(BB) A= -1 Rotation 28 Quaternion Algebra • 2D: complex numbers 2 • i:= exey, i = -1 • 3D: quaternions • i:= ex ey= exey, j:= ey ez = eyez, k:=ex ez=exez • i2 = -1, j2 = -1 , k2 = -1 • ijk = (exey)(eyez)(exez) = -1 • 4D: Biquaternions (complex quaternions, spacetime algebra) Rotation 29 Rotation and GA Right-multiplication of Vectors by Bivectors ex i = ex (exey) = (exex ) ey= ey = ey i = ey(exey)=-ey(eyex)= -ex = Rotation 30 Generic Rotation in 2D • Multiple Rotation ex i i = (ex i) i = ey i = -ex = -1 ex • Arbitrary vector A = Ax ex + Ay ey A i = Ax ex i + Ay ey i = Ax ey - Ay ex • Rotation by arbitrary angle: A cos + A i sin ≡ “A e i ” rotates vector A by angle in plane i Inverse rotation: Ai = -iA : - A ei = e-i A Rotation 31 Rotor in 2D • Rotor R := e i = cos + i sin mit i² = -1 A ei = e-i A = e-i /2 A e i/2 = R A R-1 With R=e-i /2 “Rotor” R-1=ei /2 “inverse Rotor” A R-2 = R2 A = R A R-1 • Product of rotors is multiple rotation R=ABCD, R-1=DCBA is “reverse” R Rotation 32 Rotor in nD • Rotor in plane U, Vektor v: R = cos + sin U U² = -1 Expect: Rv or vR-1 or R v R-1 • Problem: With arbitrary vector v there would be a tri-vector component: Rv = v cos + sin (U v + U v ) iff U v ≠ 0 ( v not coplanar with U) Rotation 33 Rotation in nD -1 Consider: R v R mit v =v┴ + v║ : – We have: U v┴ = 0 d.h.