Clifford Algebras
Clifford Algebras
Clifford Algebras
Daniel Klawitter
Clifford Algebras
Geometric Modelling
and Chain Geometries
with Application in Kinematics
Foreword by Prof. Dr. Gunter Weiss
ISBN 978-3-658-07617-7
DOI 10.1007/978-3-658-07618-4
To my parents,
Manuela and Burghard
Foreword
What is the right mathematical model to a real phenomenon of our
world? Do there exist criteria whether a model can be called elegant as well as practically ecient? These questions will surely be
answered dierently by e.g. a pure mathematician on the one hand
and an engineer on the other. But both will have to start with abstracting real world phenomena to objects of a more or less platonic
world. In this idealized world one constructs the structured geometric image of real world processes and objects in consideration. Thus,
on the way to constructing explicit mathematical models, geometric
representations, and geometric reasoning will always play an essential
part and provides insight into perhaps not that obvious interrelations
and structures.
The idea for this book took its origins from space kinematics and
robotics and the standardized way. Its aim is to be understood and
handled practically. There are two essentially dierent ways to model
space kinematics, although both ways are based on an Euclidean motion as the core element of the manifold of motions.
Loosely spoken, a seemingly clear model will interpret such a basic
element as a point and the manifold of motions as a point manifold.
(It turns out, that, via a clever transfer mapping, this point manifold becomes a hyper-quadric in a seven-dimensional real projective
space. This hyper-quadric is referred to as Studys quadric.) The
second model interprets a motion as a number rather than a point on
a line: In order to describe the position of it on that line, one needs a
pair of Hamilton quaternions as coordinates, which in addition fulll
conditions. Both models relate to well understood families of models.
The rst one can be seen e.g. as a higher dimensional analogue to the
geometry of lines and screws. Indeed, screws are closely connected
to motions. The second model relates to so-called circle geometries.
VIII
Foreword
Foreword
IX
group structures to model not only Euclidean kinematics and line geometry, but also circle geometries and their generalizations, so-called
chain geometries and, to take it even further, also non-Euclidean
kinematics and line geometry. Projective and chain geometric properties alternately relate to algebraic properties and concepts. Again,
the concept cross ratio becomes such an interface between Geometry and Algebra, being on the one hand visualized as quadruplets
of points and on the other being a number resp. a ring element as
well as responsible for structuring e.g. specic point sets of a line
to chains. The presented method allows to connect disciplines of
(Applied) Mathematics, that, up to now, are treated independently.
Furthermore, it provides an eective tool for solving practical problems of Mechanical Engineering and Robotics, too.
Is this book a Maths book on theoretical and practical applications
of Cliord Algebras? Yes, but in its rich details and examples it aims
at visualizing the more or less hidden Geometries of the Cliord Algebra calculus, and therefore, it could rightly be called a Geometry
book. By its wide range of examples it is protable to a broad span
of readers, be it mechanical engineers or mathematicians. Like me,
they surely consider its publication as an absolute stroke of luck.
Wien
Gunter Weiss
Preface
Usually, dual unit quaternions are used to describe Euclidean displacements in three-dimensional kinematics. Engineers widely use
this elegant calculus next to a real projective geometric model called
Studys quadric. Moreover, there exists a sphere-model in a fourdimensional module space over the ring of dual numbers. However,
there are no investigations connecting these three models up to now.
After recalling Studys quadric, the sphere-model , and other models, we present a survey-like introduction to Cliord algebras and their
Spin and Pin groups. The advantage of Cliord algebras is that geometric objects and transformation may be represented as algebra elements. Using the so-called sandwich operator, it is possible to apply
transformations that are represented by algebra elements directly to
geometric objects, also described by algebra elements. We introduce
the homogeneous and conformal Cliord algebra model. Furthermore, we show how to model special Cayley-Klein geometries and
their isometry groups in a homogeneous Cliord algebra model. As
an example, we focus on the homogeneous Cliord algebra model corresponding to line geometry and derive the correspondence between
projective transformations of three-dimensional projective space and
the Pin group Pin(3,3,0) . Therefore, we introduce the Cliord algebra
C(3,3,0) , constructed over the quadratic space R(3,3) , and describe how
points on Kleins quadric are embedded as null vectors. We discuss
how geometric entities that are known from Kleins model, i.e., linear line manifolds can be transferred to this homogeneous Cliord
algebra model. All entities known from line geometry occur naturally in this model and can be transformed projectively by the application of the sandwich operator. The action of grade-1 elements
corresponds to the action of null polarities on P3 (R), i.e., correlations
that are involutions as the basic elements building up the group of
regular projective transformations. It is proven that every regular
projective transformation of P3 (R) can be expressed as the product
XII
Preface
Preface
XIII
In the last chapter we recall the concept of kinematic mappings. Different kinematic mappings have been found, for example the kinematic mapping of Study that maps the group of Euclidean displacements to Studys quadric, or the kinematic mapping of Blaschke and
Gr
unwald. Furthermore, we use the Cliord algebra calculus to unify
dierent kinematic mappings. With this construction it is possible
to unify the concept of kinematic mappings for isometry groups of
Cayley-Klein geometries and for orthogonal groups SO(p, q). Collineations in any kinematic image and the corresponding Cayley-Klein
space can be derived from the homogeneous Cliord algebra model.
The kinematic images of Pin- and Spin groups are projective varieties. Due to the fact that projective varieties correspond to ideals,
methods of Grobner basis calculus can be applied to kinematic image spaces. We present kinematic mappings for Euclidean spaces of
dimensions two, three and four in detail.
Dresden
Daniel Klawitter
Contents
Foreword
Preface
Introduction
1 Models and Representations
1.1 Description of Displacements . . . . . . . . . . . . . . .
1.1.1 Homogeneous Matrices . . . . . . . . . . . . . .
1.1.2 Dual Quaternions . . . . . . . . . . . . . . . . .
1.1.3 Dual Orthogonal Matrices . . . . . . . . . . . .
1.2 Point Models for Lines and Displacements . . . . . . .
1.2.1 Kleins Quadric . . . . . . . . . . . . . . . . . . .
1.2.2 Studys Quadric . . . . . . . . . . . . . . . . . .
1.2.3 Studys Sphere . . . . . . . . . . . . . . . . . . .
1.3 Geometric Algebras, Cliord Algebras . . . . . . . . . .
1.3.1 Denition of a Geometric Algebra . . . . . . . .
1.3.2 Properties of Cliord Algebras . . . . . . . . . .
1.3.3 Pin and Spin Groups . . . . . . . . . . . . . . .
1.3.4 Matrix Representation of Cliord Algebras . .
1.3.5 Linear Transformation of the Vector Space . .
1.4 The Homogeneous Model . . . . . . . . . . . . . . . . .
1.4.1 The Construction . . . . . . . . . . . . . . . . .
1.4.2 Exterior Algebra . . . . . . . . . . . . . . . . . .
1.4.3 Homogeneous Model via projective Grassmann
Algebra . . . . . . . . . . . . . . . . . . . . . . .
1.5 The Conformal Model . . . . . . . . . . . . . . . . . . .
1.5.1 Construction of Conformal Geometric Algebra
1.5.2 Blades in CGA . . . . . . . . . . . . . . . . . . .
1.5.3 Conformal Transformations . . . . . . . . . . . .
VII
XI
1
5
5
5
6
13
15
15
21
25
26
26
28
33
36
37
38
38
39
41
42
43
45
48
XVI
Contents
1.6
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
50
51
53
55
61
67
72
72
73
74
76
77
78
78
80
84
87
89
94
101
101
101
102
104
105
106
109
114
116
118
121
123
124
126
131
131
Contents
XVII
133
135
136
139
139
141
145
150
151
156
161
161
165
168
169
176
177
181
181
182
184
185
185
185
187
187
Conclusion
201
190
192
192
193
195
197
XVIII
Contents
Index
203
List of Symbols
205
Acknowledgment
209
Bibliography
211
Introduction
Euclidean displacements in three-dimensional kinematics are usually
described by dual unit quaternions. Engineers widely use this elegant
calculus aside a real projective geometric model called Studys quadric. Moreover, there exists a sphere-model in a four-dimensional
module space over the ring of dual numbers. Seemingly there are no
investigations to connect these three models up to now and this was
the rst motivation for the present research, which turned out to lead
to a far more general view. One theme of interest was thereby to nd
interpretations of dierent projective closures of the two ring based
models and the real projective model.
The cross ratio is a fundamental invariant of projective geometry. It
is well known that the cross ratio of four complex numbers is real
if, and only if, they all lie on a M
obius circle, i.e., a circle or a line
augmented by a point at innity. A generalization of the so-called
Mobius geometry is obtained by using dierent algebras instead of
the complex numbers. A cross ratio for dual quaternions was dened
in [58] by
cr(a, b, c, d) := (a b) (b c)1 (c d) (d a)1 .
With this denition a lot of questions arise. What is the cross ratio
of dual unit quaternions and how can it be interpreted geometrically
and kinematically? Under which preliminaries is the cross ratio of
four dual unit quaternions a real or a dual number? Are there points
on a conic section on Studys quadric corresponding to dual unit quaternions with a real or dual cross ratio?
Kinematic and line geometric models that are used today date back
to E. Study [61] and F. Klein [42]. These models nd application in
modern treatises, see for example [13], [38]. Especially line geometry
Introduction
Introduction
geometry as example. Furthermore, a new geometric algebra which allows the description of inversions with respect to quadrics in principal
position as Pin group is presented. This new model is constructed for
dimension two and three in detail. Furthermore, the generalization
to arbitrary dimension is presented.
In chapter 2 we give the answer to the question for a cross ratio
of four dual unit quaternions. This leads to a branch of geometry
called chain geometry. Chains are subsets of the projective line over
a ring that can be parametrised with the cross ratio. The origins
of chain geometry can be found in Benz [5]. Benz dealt with twodimensional algebras over the real numbers, i.e., the double numbers
or split-complex numbers A, the complex numbers C, and the dual
numbers D. He classied the chains in these three algebras. More
recent results can be found in [11] or [33]. Chain geometry is an
essential tool for this work. The connection between quadric models
and chain geometries over kinematic algebras was found by H. Hotje,
see [37]. A. Blunck and M. Stroppel generalized this approach
to so-called Klingenberg chain spaces, see [9] and [12]. Therefore, it
is natural to apply this theory to dual quaternions and to examine
the kinematic and geometric interpretation. The cross ratio of four
dual unit quaternions corresponding to line-symmetric motions with
respect to a hyperboloid of revolution
x2
y2
z2
+
= 1, with a = b,
a2
b2
c2
Introduction
p1
p0 ,
p2
p0 ,
p3
p0
T
1
t
0T
A
p,
A SO(3), t R3 , p = (1, x1 , x2 , x3 )T .
(1.1)
a, b, c, d R,
(1.2)
(1.3)
The multiplication rules for the quaternion units (1.3) are summarized in the following scheme:
i
j
k
i
1
k
j
j
k
k j
.
1
i
i 1
This scheme allows to extend the concept of quaternion multiplication to general quaternions. The real numbers form the center of the
quaternions, i.e., real numbers commute with all other quaternions.
Later, in section 1.3 we study algebras and we see that the quaternions are an example of an R-algebra.
Denition 1.2. The anti-involution
: H H,
a + bi + cj + dk = q q = a bi cj dk
qq =
a 2 + b2 + c 2 + d 2 .
1
q
q
=
.
q2
qq
In the same way we see the action of this operator on the vectorial
quaternions j and k
qjq = 2(bc ad)i + (a2 b2 + c2 d2 )j + 2(cd + ab)k,
qkq = 2(bd + ac)i + 2(cd ab)j + (a2 b2 c2 + d2 )k.
The sandwich operator is linear and the image of a vectorial quaternion is a vectorial quaternion again. Furthermore, the scalar product of two vectors x, y R3 is invariant under the action of the sandwich operator and it is orientation preserving. When collecting the
images of the basis vectors e1 , e2 , e3 in a matrix A, we get
a 2 + b2 c 2 d 2
A=
2(bc + ad)
2(bd ac)
2(bc ad)
a2 b 2 + c 2 d 2
2(cd + ab)
2(bd + ac)
2(cd ab) .
2
a b2 c 2 + d 2
(1.4)
Matrix (1.4) is the well-known form of a rotation matrix. The components of a unit quaternion are the Euler parameters of a rotation.
Dual Numbers Like complex numbers dual numbers are an extension of the real numbers. A dual number has the form z = a + b,
where a and b are real numbers and is the dual unit that squares
to zero, 2 = 0. Addition is dened component-wise. For two dual
numbers the product is dened by
(a1 + b1 )(a2 + b2 ) = a1 a2 + (a1 b2 + a2 b1 ) + 2 (b1 b2 )
= a1 a2 + (a1 b2 + a2 b1 ).
z z =
(a + b)(a b) = |a|.
Any dual number without vanishing real part has an inverse dual
number
z1 = (a + b)1 :=
1
z .
a2
b
1
1
1
2 = 2 (a b) = 2 z .
a
a
a
a
Dual Vectors Later we will use dual vectors. Therefore, we introduce the n-dimensional module
Dn := {v | v = v +
v , 2 = 0, v, v Rn }.
+
v , w
),
where ,
denotes the standard scalar product of Rn . Especially for
the dimension n = 3 we are able to dene a cross product by
v w = v w + (
v w + v w).
(1.5)
10
Multiplication is dened with the relations for quaternion. Furthermore, the dual unit commutes with the quaternion units i = i,
j = j, k = k. Dual quaternions form an eight-dimensional vector
space over the real numbers. The basis elements are 1, i, j, k, , i, j, k.
Displacements Euclidean displacements can be described by dual
unit quaternions. A dual quaternion q = a0 + a1 i + a2 j + a3 k + (c0 +
c1 i + c2 j + c3 k) is normed or a dual unit quaternion, if the norm is
equal to one
N (q ) := q q = a20 + a21 + a22 + a23 + 2(a0 c0 + a1 c1 + a2 c2 + a3 c3 ) = 1,
(1.6)
(1.7)
11
0
a20 +a21 a22 a23
2a0 a3 +2a1 a2
2a1 a3 2a0 a2
0
2a1 a2 2a0 a3
a20 a21 +a22 a23
2a0 a1 +2a2 a3
0
1
x
2a0 a2 +2a1 a3
,
2a2 a3 2a0 a1 y
a20 a21 a22 +a23
z
where
l = 2c1 a0 2c0 a1 + 2c3 a2 2c2 a3 ,
m = 2c2 a0 2c3 a1 2c0 a2 + 2c1 a3 ,
n = 2c3 a0 + 2c2 a1 2c1 a2 2c0 a3 .
v
v
sin sin l3 + cos l0 i
2 2 2
2 2
v
v
sin l4 + cos l1 j sin l5 + cos l2 k,
2
2
2
2
2
2
(1.8)
12
Matrix Representation
A = a0 + a1 i + a2 j + a3 k + (c0 + c1 i + c2 j + c3 k),
X = x0 + x1 i + x2 j + x3 k + (x4 + x5 i + x6 j + x7 k),
a0 a1 a2
a1 a0 a3
A=
a 2 a 3 a 0
a3 a2 a1
A O
x, with
C A
c0 c1 c2 c3
a3
a2
, C = c1 c0 c3 c2 ,
a1
c2 c3 c0 c1
a0
c3 c2 c1 c0
A
C
O
x, with
A
a0 a1 a2 a3
c0 c1 c2 c3
a1 a0 a3 a2
c1 c0 c3 c2
A=
a2 a3 a0 a1 , C = c2 c3 c0 c1 ,
a3 a2 a1 a0
c3 c2 c1 c0
13
(1.9)
Hence, the dual scalar product is invariant with respect to the action of orthogonal dual matrices. Furthermore, the group of dual
orthogonal matrices is a Lie group. Its Lie algebra is the set of skewsymmetric dual matrices. In [57] the authors used this for motion
interpolation purposes. Ruled surfaces are interpreted as curves on
14
the dual unit sphere that we will introduce later on. Therefore, a generalization of the Rodriguez formula for dual skew symmetric 3 3
matrices is used:
O = exp g = I + sin g + (cos 1)g 2 ,
(1.10)
0
g3
1 2 3 T
3
= (g , g , g ) g
0
g2 g1
(1.11)
g2
g1
=g
see right before remark 1.1. With the ad mapping Eq. (1.11) and the
dual Rodriguez formula (Eq. (1.10)) we write
O(t) = exp( (t))ad(n ) = I + sin( (t))ad(n ) + (cos( (t)) 1)(ad(n ))2 ,
where
n =
v w
v w
15
16
x1
y1
x1
y1
x2
y2
x2
y2
x3
y3 x0 x1 x2 x3 x0 x2 x3 x1 x0 x3 x1 x2
+
+
= 0,
=
x3 y0 y1 y2 y3 y0 y2 y3 y1 y0 y3 y1 y2
y3
which yields
M24 : p01 p23 + p02 p31 + p03 p12 = 0.
(1.12)
x Qx = 0, with Q =
O
1
2I
1
2I
, x = (x0 , x1 , x2 , x3 , x4 , x5 )T ,
(1.13)
1
= 0.
64
17
18
Two-dimensional subspaces of M24 There are two types of twospaces that are completely contained by M24 . The rst one is given
by the image of a eld of lines, i.e., all lines contained by the same
plane in P3 (R). The second type is given by all lines concurrent to the
same point, i.e., a bundle of lines. Now we ask for the intersection of
two two-spaces contained in Kleins quadric. Two-spaces of the same
type always intersect in one point. This means two dierent elds
of lines or two dierent bundles of lines contain one common line.
Two two-spaces of dierent type may have either empty intersection
or they intersect in a line. Thus, a eld of lines and a bundle of
lines may have empty intersection or they intersect in a common
pencil of lines. The type of the two-space, i.e., if it corresponds to a
bundle of lines or a eld of lines can be determined by its intersection
with the Klein image of the eld of ideal lines denoted by P . Note
that P corresponds to a eld of lines, and therefore, to a two-space
entirely contained in Kleins quadric. If a two-space corresponds to a
bundle of lines it has either no point or a whole line in common with
P . A two-space corresponding to a eld of lines has one point in
common with P . Fig. 1.1 shows the correspondence between lines
and two-spaces on Kleins quadric and their pre-images under the
Klein mapping. This gure is inspired by a gure from [56, p. 142].
Subspace Intersections of M24 With Kleins mapping sets of lines
can be studied as sets of points in ve-dimensional projective space.
Clearly, all points on the quadric are self-conjugate with respect to
the polarity . It is natural to ask for the intersections of Kleins
quadric with subspaces. In Kleins model k-space intersections, k 4
with M24 P5 (R) dene special sets of lines in P3 (R). These sets of
lines correspond to so-called linear line manifolds, cf. [56].
Conics on Kleins quadric Two-spaces which are neither tangent to
nor belong to Kleins quadric intersect in a conic section. Tangent
two-space intersections with M24 result in degenerate conics, i.e., two
intersecting lines on Kleins quadric that correspond to two pencils
of lines with one line in common in projective space P3 (R). In this
paragraph we discus the non-degenerate case. Therefore, we choose
three skew lines respectively their Pl
ucker coordinates L1 , L2 and L3
M24
19
F
P5 (R)
Figure 1.1: Subspaces on Kleins quadric and their geometric interpretation in P3 (R). A pencil of lines P is mapped to the line
P on M24 . A bundle of lines B and a eld of lines F are
mapped to a two-spaces B and F that are contained entirely in M24 .
20
c c
.
c2
(1.14)
21
m00
m10
M=
m20
m30
0
m11
m21
m31
0
m12
m22
m32
0
m13
m23
m33
22
be a homogeneous transformation matrix of an Euclidean displacement, see Eq. (1.1). The coordinates of a point on S26 \V 3 corresponding to this displacement can now be determined by the relations:
(1.15)
a 0 : a1 : a2 : a3
= m00 +m11 +m22 +m33 : m23 m32 : m31 m13 : m12 m21 ,
= m23 m32 : m00 +m11 m22 m33 : m12 +m21 : m31 +m13 ,
= m31 m13 : m12 +m21 : m00 m11 +m22 m33 : m23 +m31 ,
= m12 m21 : m31 +m13 : m23 +m32 : m00 m11 m22 +m33 ,
and
2m00 c0 = a1 m10 + a2 m20 + a3 m30 ,
(1.16)
These relations can be found in [23] or [61]. Note that each equation
of Eqs. (1.15) is important because the ratio could get (0 : 0 : 0 : 0)
in one equation, but not in all equations, what is not allowed since
the point would lie in the exceptional space V 3 . The entries of the
matrix representing the displacement can be calculated by
m00 = a20 + a21 + a22 + a22 ,
m10 = 2(a2 c3 a3 c2 a0 c1 + a1 c0 ),
m12 = 2(a1 a2 + a0 a3 ),
m13 = 2(a3 a1 a0 a2 ),
m20 = 2(a3 c1 a1 c3 a0 c2 + a2 c0 ),
m21 = 2(a1 a2 a0 a3 ),
m23 = 2(a2 a3 + a0 a1 ),
m30 = 2(a1 c2 a2 c1 a0 c3 + a3 c0 ),
m31 = 2(a3 a1 + a0 a2 ),
m32 = 2(a2 a3 a0 a1 ),
m33 =
a20
a21
a22
(1.17)
a23 .
23
Subspaces contained in S26 A discussion of linear subspaces contained entirely in Studys quadric can be found in [4]. The maximal
possible dimension of subspaces contained in a six-dimensional quadric equals three. In fact S26 carries two three-parameter families of
three-dimensional subspaces. To classify the subspaces that are entirely contained in S26 we follow [24]. First we observe that a point
P S26 corresponds to the identity, a rotation, or a translation if,
and only if, c0 = 0, see [25]. This follows directly with Eq. (1.8).
Thus, all rotations and translations are contained in the hyperplane
x4 = 0, that is the tangent hyperplane at the point PI = (1, 0, . . . , 0)T R,
i.e., the point on S26 corresponding to the image (e) of the identity
element e SE(3) under the Study mapping . All points that are
conjugate to PI are contained in a quadratic cone with equation
C0 : x1 x5 + x2 x6 + x3 x7 = 0.
24
25
(1.18)
Denition 1.8. The mapping that maps an oriented line L of Euclidean three-space to the dual vector v = v + v , where (v, v) are the
normalized Pl
ucker coordinates of L, is called the Study mapping. Its
image space is a model of the set of spears of the three-dimensional
Euclidean space and is called Studys sphere.
Remark 1.4. For dual unit quaternions the same construction can
be applied. In this case we get the set of all dual unit vectors of D4 .
Note, that the direction of the Euler coordinates denes the direction
of the displacement. We have to dier between two displacements
corresponding to dierent oriented Euler coordinates, and therefore,
we consider oriented displacements.
26
1 i, j n.
(1.19)
27
i = j,
(1.20)
The basis {e0 , e1 , e2 , e12 } for C(0,2,0) is called a standard basis for this
geometric algebra. As we already know from the denition e21 = e22 =
1. The square of e12 evaluates to
e212 = e12 e12 = e12 e21 = e11 = 1.
28
0kn
n
i
V of all exterior
i=0
C+
(p,q,r)
C
(p,q,r)
:=
n
i
i=0
i even
n
i
V.
(1.22)
i=0
i odd
(1.23)
For details see [55]. To make this isomorphism more clear we give
another example.
29
Example 1.2. As we already know from Ex. 1.1 the Cliord algebra
C(0,2,0) is isomorphic to H. The even part of this algebra is generated
by {e0 , e12 }, and therefore, isomorphic to C. Furthermore, the Cliord
algebra C(0,1,0) is the algebra with one generator squaring to 1.
Therefore, C(0,1,0)
= C. All in all we have C+
(0,2,0) = C(0,1,0) .
Denition 1.11. The center of a ring R is the set of all elements
that commute with all other elements
C(R) := {c R | cx = xc for all a R} .
(1.24)
30
these observations to the grade-2 basis element e12 . We use the antiinvolution property and the anti-commutativity
e12 = (e1 e2 ) = e2 e1 = (1)2 e2 e1 = e21 = e12 .
n = p + q + r.
(1.25)
(1.26)
The main involution has no eect on the even subalgebra and it commutes with the conjugation. This means for an arbitrary algebra
element x, the equation (x ) = (x) holds.
Cliord Algebra Products We denote grade-1 elements, i.e., vectors, by small gothic letters. Elements of higher grade are denoted
by big gothic letters.
The inner Product The scalar product of two vectors, i.e., grade-1
elements can be written in terms of the geometric product
a b :=
1
(ab + ba).
2
(1.27)
31
A generalization of the
inner product
to blades, see Def. 1.12, can be
found in [36]. For A k V, B l V the generalized inner product
is dened by
A B := [AB]|kl| ,
The geometric Product With the inner (1.27) and the exterior
product (1.28) the geometric product of vectors can be written in
the following form
ab = a b + a b.
(1.29)
is the k-fold exterior product of vectors
Denition 1.12. A k-blade
or grade-1 elements v 1 V . Therefore, a k-blade can be written as
A = a1 a2 . . . ak =:
k
ai .
i=1
32
NI(A) :=
1
V :vA=0 .
k
1
v
V :vA=0 .
The IPNS and the OPNS of k-blades can be used to describe subspaces of the vector space V , see [53]. Furthermore, we have the
property of the outer and inner product null space
NI(a b) = NI(a) NI(b),
(1.30)
p 0
1
2
3
4
5
6
7
q
0
R
R
R(2)
C(2)
H(2)
2
H(2)
H(4)
C(8)
C
R(2)
2
R(2)
R(4)
C(4)
H(4)
2
H(4)
H(8)
H
H
C(2) H(2)
R(4) C(4)
2
R(4) R(8)
R(8) 2 R(8)
C(8) R(16)
H(8) C(16)
2
H(8) H(16)
2
H(2)
H(2)
H(4)
C(8)
R(16)
2
R(16)
R(32)
C(32)
C(4)
H(4)
2
H(4)
H(8)
C(16)
R(32)
2
R(32)
R(64)
R(8)
R(8)
C(8)
R(16)
H(8)
C(16)
2
H(8) H(16)
H(16) 2 H(16)
C(32) H(32)
R(64) C(64)
2
R(64) R(128)
2
Table 1.1 shows the matrix representation of the non-degenerate Clifford algebras C(p,q,0) . In this table R(n) denotes the n n-matrices
over R and 2 R(n) denotes R(n) R(n).
33
0
,
1
e1
1
0
0
,
1
e2
0
1
1
,
0
e12
0
1
1
.
0
(C(p,q,r) ) := g C
|(g)vg
V
v
V
.
(p,q,r)
(1.31)
34
(1.32)
a b =
1
(g)g [g (g)(ab + ba)] .
2
Thus, we get
1
1
((g)g )(g (g)) (ab + ba) = (ab + ba) = a b.
2
2
Linear transformations of 1 V that preserve distances and angles, i.e.,
the scalar product, are elements of the orthogonal group O(p, q, r). In
fact, Pin(p, q, r) is a double cover of the orthogonal group O(p, q, r),
a b =
see [21].
35
g
= g .
gg
In the homogeneous Cliord Algebra models we use g because multiplication with a homogeneous factors does not change the geometric
meaning of the sandwich operator.
Spin(p,q,r) : = g C+
(p,q,r) |N (g) = 1, gvg V v V
(1.33)
= Pin(p,q,r) C+
(p,q,r) .
Note that the main involution has no eect on the even subalgebra,
and therefore, can be neglected.
Remark 1.8. The Spin group is a subgroup of the Pin group. Thus,
the scalar product of vectors is preserved under the action of the Spin
group. Spin group elements are generated by pairs of reections, so
they are rotations.
The Spin group for degenerate Cliord algebras as in the case of
Euclidean geometry are semi-direct products of Spin groups for the
non-degenerate part and an additive matrix group, see [15]. Later we
will construct the Cliord algebra for the three-dimensional Euclidean
space.
Remark 1.9. Due to their denitions the Pin and the Spin group
consist of two connected components. The connected components are
36
,
Pin+
(p,q,r) := g C(p,q,r) |N (g) = +1, (g)vg V v V
Pin
.
(p,q,r) := g C(p,q,r) |N (g) = 1, (g)vg V v V
For the connected components of the Spin group we use the same
notation.
C = B A,
A+ , B ,
(1.34)
Note that these columns are ordered from right to left to match the
chosen ordering of the vector representation for the Cliord algebra.
To obtain the columns of the matrix [B ] we multiply b from the left
side with each basis element
e0 b, e1 b, . . . , en b, e12 b, . . . e(n1)n b, . . . , e1...n b
37
38
39
X Y + Y X = 0.
40
i1 <...<ik
41
1 ,...,ik }
, 0 i1 < . . . < ik n,
42
43
1
(e +e+ ),
2
= (e e+ ).
o e1
0 0
0 1
0 0
0 0
1 0
e2
0
0
1
0
0
e3
0
0
0
1
0
1
0
.
0
0
0
There are also other possibilities to dene the basis of R(4,1) . This
algebra can be thought of as the homogeneous model with one additional restriction in the -direction. The linear subspaces of Rn+1
that model ane subspaces in Rn are intersected with an (n + 1)dimensional paraboloid, called the horosphere, see [20, 35]. The intersection is then projected on the n-dimensional hyperplane spanned
by [e1 , . . . , en ]. The scalar product in conformal geometric algebra
(abbreviated by CGA) helps to give an analogue to the Euclidean distance. Therefore, points are modelled by null vectors, because they
have to satisfy p p = pp = 0, since p p = 0. For a point P Rn the
corresponding algebra element has the form
1
p = o+p+ p2 ,
2
44
(1.35)
With Eq. (1.35) we see that points are described by null vectors, since
pp = p2 = 0 and pp = 0. Furthermore, it now becomes clear that the
Euclidean distance is integrated in the algebra by the scalar product.
Note that the o-component can be understood as the homogeneous
factor. This factor can be extracted from a point by
1
p = (o+p+ p2 ) = .
2
1
R(p,q) , with r + s p + q
Note that we mean the Poincare duality when we talk about duality
in geometric algebra context.
45
1
P o+p+ p2 .
2
(1.36)
46
47
48
1
V v
1
V
Spin(n+1,1,0) := g
C+
(n+1,1,0) | N (g) = 1,
gvg
1
1
V v
V
is a double cover of the orientation preserving conformal transformations. We construct this group with its subgroups by studying the
grade-1 elements, that generate this group. First we identify the
group SE(3) of Euclidean displacements as a subgroup of the conformal group. Therefore, we have to nd all non-null vectors that
leave the point at innity xed. This results in the condition
(g)g1 = .
and see that this transformation leaves the point at innity and the
origin xed. Furthermore, this transformation is a reection in the
plane dened by n. Two reections in parallel planes result in a
49
+ sin (o ) = e 2 o .
2
2
50
Type of Operation
re. in origin plane
re. in real unit sphere
re. in origin
rotation over in R-plane
translation over t
scaling by e
transversion over t
Explicit Form
n
c 12
o
Exp. Form
none
none
none
cos( 12 )sin( 12 )R
e 2 R
1
1
1 12 t
e 2 t
cosh( 21 )sinh( 12 )o
e 2 o
1 + 12 ot
e 2 ot
1
1
Q=
O
I
I
,
O
51
quadric, since multiplication with real scalars has no eect, see Eq.
(1.13). As underlying vector space for the Cliord algebra we take R6
as vector space model for P5 (R). The corresponding Cliord algebra
has signature (p, q, r) = (3, 3, 0) (cf. [47]) and is of dimension 26 =
ucker coordinates, see Def. 1.4,
64. Lines of P3 (R) represented by Pl
correspond to null vectors in this algebra, i.e., vectors that square to
zero. A vector is given by
v = x1 e1 + x2 e2 + x3 e3 + x4 e4 + x5 e5 + x6 e6
(1.37)
52
p01
p23
=
=
x0
y0
x2
y2
x
x1
, p02 = 0
y1
y0
x3
x
,p = 3
y3 31 y3
x
x2
, p03 = 0
y2
y0
x1
x
,p = 1
y1 12 y1
x3
,
y3
x2
.
y2
x x
= i j = xi yj xj yi
yi yj
=
cik xk
cjl yl
cjl xl
cik yk
k
l
l
k
cik cjl (xk yl xl yk ),
=
k,l
where (i, j) is one of (0, 1), (0, 2), (0, 3), (2, 3), (3, 1) or (1, 2), see [56, p.
139]. If we write the action of this transformation on the space of
lines as matrix vector product we get a 6 6 matrix L containing the
coecients from the equations above
(p01 , p02 , p03 , p23 , p31 , p12 )T = L (p01 , p02 , p03 , p23 , p31 , p12 )T .
53
Now we examine geometric entities that are described in this geometric algebra by inner product and outer product null spaces, see
[53].
54
6
xi xj
L = v1 v 2 =
yi yj eij .
i,j=1
i<j
are sucient to describe the inner product null space of L. All null
vectors a = a1 e1 + a2 e2 + a3 e3 + a4 e4 + a5 e5 + a6 e6 satisfying these
conditions correspond to the set of lines contained in a linear line
55
congruence. In a dual way, a linear line congruence can also be represented as the outer product null space of a four-blade that can be
constructed as exterior product of four linearly independent vectors.
1
V | x1 a 4 + x 2 a 5 + x3 a 6 + x 4 a 1 + x 5 a 2 + x6 a 3 = 0 ,
1
V | x1 a 4 + x 2 a 5 + x 3 a 6 + x 4 a 1 + x 5 a 2 + x 6 a 3 = 0 ,
1.6.3 Transformations
To describe the action of the versor group we recall the denition of
a versor, see Def. 1.10. A versor a is an algebra element that can be
expressed as the k-foldgeometric product of non-null vectors. Thus,
1
a = v1 . . . vk with vi
V, i = 1, . . . , k . Transformations are applied
with the so-called sandwich operator.
56
1
v
.
vv
Thus, the inverse element diers from the conjugate element by a real
factor. This is also true for a k-fold product of vectors.
Since we are working in a homogeneous Cliord algebra model multiplication with a real factor does not change the geometric meaning of
an algebra element. Therefore, we use the conjugate element instead
of the inverse element. Hence, the sandwich operator that we use is
given by:
(a)va = aa ((a)va1 ).
This operator does not involve an inverse, and therefore, it can also
be applied if the element a is not invertible.
57
The Cliord group of this Cliord algebra, i.e., the group of projective automorphisms of Kleins quadric is generated by non-null
vectors. Thus, we are interested in the action of non-null vectors on
null vectors. Let
a = a1 e 1 + a 2 e 2 + a 3 e 3 + a 4 e 4 + a 5 e 5 + a 6 e 6 ,
v = x1 e1 + x2 e2 + x3 e3 + x4 e4 + x5 e5 + x6 e6
vv = x1 x4 + x2 x5 + x3 x6 = 0.
The action
of the sandwich operator (a)va is linear on the vector
1
V . The matrix acting on P5 (R) can be represented by
space
k1
a2 a4
a3 a4
M=
a4 a4
a5 a4
a6 a4
a1 a5
k2
a3 a5
a4 a5
a5 a5
a6 a5
a1 a6
a2 a6
k3
a4 a6
a5 a6
a6 a6
a1 a1
a2 a1
a3 a1
k4
a5 a1
a6 a1
a1 a2
a2 a2
a3 a2
a4 a2
k5
a6 a2
a1 a3
a2 a3
a3 a3
,
a4 a3
a5 a3
(1.39)
k6
with
k1 = a5 a2 a6 a3 ,
k2 = a6 a3 a4 a1 ,
k3 = a4 a1 a5 a2 ,
k4 = a5 a2 a6 a3 ,
k5 = a6 a3 a4 a1 ,
k6 = a4 a1 a5 a2 .
Naturally, we now ask for the corresponding projective mapping acting on P3 (R). Therefore, we examine the action of the collineation M
on the space of lines. The action on the base lines corresponding to
the Pl
ucker coordinate vectors
b1 = (1 : 0 : 0 : 0 : 0 : 0),
b2 = (0 : 1 : 0 : 0 : 0 : 0),
b3 = (0 : 0 : 1 : 0 : 0 : 0),
b4 = (0 : 0 : 0 : 1 : 0 : 0),
b5 = (0 : 0 : 0 : 0 : 1 : 0),
b6 = (0 : 0 : 0 : 0 : 0 : 1)
h2 = Mb2 ,
h3 = Mb3 ,
h4 = Mb4 ,
h5 = Mb5 ,
h6 = Mb6 .
If we look at the images of the three ideal lines b4 , b5 , b6 we see, that the
corresponding lines posses linear dependent direction vectors. This
58
means that the three image lines h4 , h5 , h6 intersect in the same ideal
point, and therefore, they belong to a bundle of lines. Since the lines
b4 , b5 , b6 dene a eld of lines and the image is a bundle of lines, the
mapping must be a correlation. With this knowledge we determine
the action of M on the bundle of lines concurrent to the origin spanned
by b1 , b2 , b3 . The image of this bundle of lines is a eld of lines. To
obtain the coordinates of the plane that carries the eld of lines we
intersect the image lines and get three points dening the plane in
P3 (R)
h1 h2 = s1 = (a3 , a5 , a4 , 0)T ,
h1 h3 = s2 = (a2 , a6 , 0, a4 )T ,
h2 h3 = s3 = (a1 , 0, a6 , a5 )T .
h1 h6 = (a3 , a5 , a4 , 0)T ,
h5 h6 = (0, a1 , a2 , a3 )T ,
p2 = R(a4 , 0, a3 , a2 )T .
h2 h6 = (a3 , a5 , a4 , 0)T ,
h4 h6 = (0, a1 , a2 , a3 )T ,
p3 = R(a5 , a3 , 0, a1 )T .
h3 h5 = (a2 , a6 , 0, a4 )T ,
h4 h5 = (0, a1 , a2 , a3 )T ,
p4 = R(a6 , a2 , a1 , 0)T .
59
0 a4 a5 a6
a4 0 a3 a2
.
K=
a5 a3
0 a1
a6 a2 a1
0
1
(aa)2 .
4
(1.41)
60
l2 = e2 + 2e4 ,
l3 = e3 + 2e5 .
1
1
1
e2 +e4 +
e3 +e5 +
e1 +e6 | , , R .
2
2
2
Indeed, the outer product null space of P spans the same two-space
as [l1 , l2 , l3 ]. To get the regulus dened by the three lines Li or their
corresponding null vectors li we have to nd all null vectors in the
two-space p = l1 + l2 + l3 . Therefore, we compute
!
p2 = 4 + 4 + 4 = 0.
(1.42)
The parameters ( : : ) can be interpreted as homogeneous coordinates in a projective plane P12 . Eq. (1.42) denes a conic contained in
this plane. Now we apply an arbitrary null polarity given by
a = a1 e 1 + a 2 e 2 + a 3 e 3 + a 4 e 4 + a 5 e 5 + a 6 e 6 .
61
This equation denes a conic in the image of the projective plane P12
under the null polarity induced by a. Every point on this conic delivers
the parameters for a null vector that corresponds to a line in P3 (R).
62
63
= y03 e123 + y3 y02 e235 y2 y02 e236 + y3 y02 e134 + y32 y0 e345 y2 y3 y0 e346
y02 y1 e136 + y3 y1 y0 e356 + y2 y02 e124 + y2 y3 y0 e245 y22 y0 e246
y0 y1 y2 e146 y02 y1 e125 + y0 y1 y2 e256 + y3 y1 y0 e145 + y0 y12 e156 .
where
v1 = y0 (g1 g32 g20 + g4 g29 + g13 + g9 g24 )
+ 2y1 (g17 + g7 ) + 2y2 (g18 2g3 ) + 2y3 (g2 + 2g19 ) e1
2y0 (g31 + g14 ) + 2y1 (g27 g11 ) 2y2 (g12 y2 + g23 )
+ y3 (g9 g1 + g32 g29 + g4 g13 + g20 g24 ) e5
+ 2y0 (g15 g30 ) + 2y1 (g8 + g28 ) + 2y3 (g10 + g21 )
+ y2 (g9 g24 g4 + g1 g32 + g20 + g29 g13 ) e6 ,
v2 = y0 (g1 g29 + g13 + g9 g24 g32 g20 + g4 )
+ 2y1 (g17 + g7 ) + 2y2 (g18 g3 ) + 2y3 (g19 + g2 ) e2
+ 2y0 (g14 + g31 ) + 2y1 (g27 g11 ) 2y2 (g12 y2 + g23 )
+ y3 (g9 g1 + g32 g29 + g4 g13 + g20 g24 ) e4
2y0 (g16 + g26 ) 2y2 (g5 + g25 ) + 2y3 (g22 g6 )
+ y1 (g13 + g32 g4 + g29 + g9 g24 g20 g1 ) e6 ,
v3 = y0 (g1 g29 + g13 + g9 g24 g32 g20 + g4 )
+ 2y1 (g17 + g7 ) + 2y2 (g18 g3 ) + 2y3 (g2 + g19 ) e3
+ 2y0 (g15 g30 ) + 2y1 (g8 + g28 ) + 2y3 (g10 + g21 )
+ y2 (g1 g24 g4 g32 + g9 + g20 + g29 g13 ) e4
+ 2y0 (g16 + g26 ) + 2y3 (g22 g6 ) 2y2 (g5 + g25 )
+ y1 (g13 + g29 + g32 g4 + g9 g24 g20 g1 ) e5 .
64
(1.44)
L2 X = 0.
65
k1
2(g26 + g16 )
M=
2(g15 g30 )
2(g31 + g14 )
2(g7 + g17 )
k2
2(g8 + g28 )
2(g11 g27 )
2(g18 g3 )
2(g5 + g25 )
k3
2(g23 + g12 )
2(g19 + g2 )
2(g6 g22 )
,
2(g21 + 2g10 )
k4
where
k1 = g1 g20 g24 g32 g29 + g9 + g4 + g13 ,
k2 = g24 g9 + g20 g13 g32 + g1 + g4 g29 ,
k3 = g1 g13 g32 g4 + g29 + g9 g24 + g20 ,
k4 = g24 + g13 + g29 + g1 g4 g9 g20 g32 .
a02 = 2(g18 g3 ),
a03 = 2(g19 + g2 ),
(1.46)
66
1
1
K=
1
1
0
1
2
1
3
0
1
2
0
1
.
0
1
2(g18 g3 ) = 3,
2(g19 + g2 ) = 0,
2(g26 + g16 ) = 1,
2(g5 + g25 ) = 0,
2(g6 g22 ) = 1,
2(g15 g30 ) = 1,
2(g8 + g28 ) = 2,
2(g21 + 2g10 ) = 0,
2(g11 g27 ) = 1,
2(g23 + g12 ) = 2,
2(g31 + g14 ) = 1,
(1.47)
67
68
1
V for all v
1
V.
(1.48)
where
v1 = 2y0 h7 + y1 (h29 h1 h9 h13 ) + y2 (h8 h2 h19 h28 )
+ y3 (h11 + h27 + h18 h3 ) e1
y0 (h3 + h11 h27 + h18 ) + y1 (h25 + h14 + h31 h5 )
+ y2 (h4 + h32 h24 + h20 ) + 2y3 h23 e5
+ y0 (h2 + h8 h19 + h28 ) + y1 (h6 + h22 h15 + h30 ) 2y2 h21
+ y3 (h20 h4 h32 h24 ) e6 ,
69
v2 = 2y0 h7 + y1 (h29 h1 h9 h13 ) + y2 (h8 h2 h28 h19 )
+ y3 (h27 + h18 h3 + h11 ) e2
+ y0 (h3 + h11 + h18 h27 ) + y1 (h25 + h14 + h31 h5 )
+ y2 (h4 + h32 h24 + h20 ) + 2y3 h23 e4
y0 (h1 h29 h13 h9 ) + 2h16 y1 + y2 (h22 h6 h15 h30 )
+ y3 (h25 + h5 h31 + h14 ) e6 ,
v3 = 2y0 h7 + y1 (h29 h13 h1 h9 ) + y2 (h8 h2 h28 h19 )
+ y3 (h11 h3 + h27 + h18 ) e3
y0 (h2 + h8 h19 + h28 ) + y1 (h6 + h22 h15 + h30 ) 2y2 h21
+ y3 (h20 h4 h32 h24 ) e4
+ y0 (h1 h29 h13 h9 ) + 2h16 y1 + y2 (h22 h6 h15 h30 )
+ y3 (h5 + h14 + h25 h31 ) e5 .
70
c01
c11
c21
c31
c02
c12
c22
c32
c03
y0
y1
c13
, with
c23 y2
c33
(1.49)
y3
c30 = h27 h3 h11 h18 , c31 = h5 h31 h25 h14 , c32 = h24 h4 h20 h32 ,
and
c00 = 2h7 ,
c11 = 2h16 ,
c22 = 2h21 ,
c33 = 2h23 .
The matrix C describes the correspondence between regular correlations and elements of C
(3,3,0) . For a given versor that corresponds
to a correlation, we can compute the matrix representation by Eq.
(1.49). Furthermore, we can determine the versor corresponding to
a regular correlation represented by a 4 4 matrix A if we solve the
system of 16 linear equations given by
a01 = h1 +h9 +h13 h29 , a02 = h2 +h19 +h28 h8 , a03 = h3 h18 h27 h11 ,
a10 = h13 h1 +h29 +h9 , a12 = h6 +h30 +h15 h22 , a13 = h31 h25 h14 h5 ,
a20 = h19 h2 h8 h28 , a21 = h15 h30 h6 h22 , a23 = h4 h20 +h32 +h24 ,
a30 = h27 h3 h11 h18 , a31 = h5 h31 h25 h14 , a32 = h24 h4 h20 h32 ,
a00 = 2h7 ,
a11 = 2h16 ,
a22 = 2h21 ,
a33 = 2h23
(1.50)
1
1
K=
1
1
0
1
2
1
3
0
1
2
0
1
,
0
1
71
h2 +h19 +h28 h8 = 3,
h19 h2 h8 h28 = 1,
2h16 = 1,
2h21 = 1,
1
e1 +e2 +e3 +2e4 4e6 2e123 3e124 +e125 +4e126 e134 +e136
8
2e145 2e146 2e156 2e234 e235 +5e236 +2e246 4e245 2e345
6e256 +4e456 2e356 +3e12345 +3e12346 +3e12356 6e23456 .
For hh = 1 we get
h =
1
3e1 +3e2 3e3 6e4 2e123 5e124 +e125 4e126 e134 +e136
8
2e145 +6e146 2e156 +2e234 e235 +3e236 +2e246 2e345 +2e256
+4e346 4e456 2e356 e12345 +e12346 e12356 4e12456 +2e23456 .
72
1
V for all v
1
V,
g C+
(3,3)
V,
h C
(3,3)
1
V for all v
1
for a singular correlation. The matrix representations for general collineations (1.45) and correlations (1.49) are also valid for the singular
case. The system of linear equations
(1.46)
with the constraint equa
tions implied by (g)vg 1 V for all v 1 V with g C+
(3,3) can
not be solved for a singular collineation. For a singular correlation
the system of linear equations
(1.50)
with the constraint equations
implied by (h)vh 1 V for all v 1 V with h C
(3,3) has also
no solution. Nevertheless, we are able to write singular correlations
or collineations as null versors if we know the null polarities that
generate them. These null polarities can be transferred to vectors
in the homogeneous Cliord algebra model, and therefore, the whole
correlation or collineation can be transferred to the algebra model.
73
the power of this calculus we examine Lie sphere geometry in a Clifford algebra context. Lie sphere geometry is the geometry of oriented
spheres. Especially, for the three-dimensional case the set of oriented
spheres can be mapped to a hyperquadric L41 in ve-dimensional projective space P5 (R). The construction goes back to S. Lie and was
treated again by W. Blaschke, cf. [8]. A modern treatment of this
topic can be found in [14]. Moreover, the Lie construction can be
achieved for arbitrary dimension.
Euclidean
points: u Rn
Lie
T
, 1uu
R
2 , u1 , . . . , u n , 0
T
(1, 1, 0, . . . , 0, 0) R
T
1+ppr 2 1pp+r 2
sphere: center p Rn , signed
,
,
p
,
.
.
.
,
p
,
r
R
1
n
2
2
radius r
planes: uN = h, unit normal (h, h, N1 , . . . , Nn , 1)T R
N Rn
1+uu
2
74
75
real vector space R6 as a model for the projective image space together
with the quadratic form of Lies quadric
1
0
0
Q=
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0 0
0 0
0 0
.
0 0
1 0
0 1
76
O
I
I
,
O
77
1
0
2
T
3
SD : x Mx = 0, with X = xD P (D) and M =
0
0
0
1
0
0
0
0
1
0
0
0
.
0
1
The
model is obtained by using the module
ane
1 0 0
M = D3 and Q = 0 1 0. In this case oriented lines are represented
0 0 1
by grade-1
elements
that square to 1 respectively Pin group elements
of the 1 V subspace.
1
0
Q=
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
1
78
as the matrix of the quadratic form for the projective model. Oriented
displacements correspond to null vectors in this model. For the ane
model we need M = D4 and
1
0
Q=
0
0
0
1
0
0
0
0
1
0
0
0
.
0
1
0
0
1
Q=
0
0
0
0
1
0
0
0
0
1
0
0
0
0
0
79
0
0
0
0
0
1
0
0
0
0
1
0
0
0
0
1
0
0
The signature of the resulting algebra is (p, q, r) = (4, 2, 0). For every
is performed,
axis, i.e., the x- and the y -axis a conformal embedding
see [19]. Therefore, we have the embedding : R2 1 V ,
R2 P p
1
V,
1
1
P = (x, y)T e1 + xe2 + x2 e3 + e4 + ye5 + y 2 e6 = p.
2
2
(1.51)
(P )2y = 0.
(1.53)
80
e6 p = 1.
(1.54)
1
u2 = e3 + e4 + y0 e5 + y02 e6 .
2
To calculate the pre-image 1 (p) of p 1 V representing an embedded point, i.e., an algebra element fullling (1.52) and (1.53), we
recall the denition of the GIPNS and the GOPNS, see [53].
The geometric inner product null space (GIPNS) and dual the geo
metric outer product null space (GOPNS) of a k-blade A k V in
Q2GA is dened by
NIG (A) := (x, y)T R2 : (x, y) A = 0 ,
NOG (A) := (x, y)T R2 : (x, y) A = 0 .
81
Remark 1.20. When dealing with an algebra element and the corresponding geometric entity, we mention explicitly what null space is
meant. For example we talk about inner product conics. This means
the inner product null space denes a conic in R2 .
Before we start the examination of geometric objects occurring in this
model we dene special 5-blades that are necessary to change from
inner product to outer product null spaces and vice versa.
Denition 1.19. On the one hand the 5-blade
I = e2 e5 e1 e4 (e3 + e6 )
Note that dualization is realized with the inner product. Now we take
a look at the inner product null space of grade-1 elements that are not
embedded points. Therefore, at least one of the conditions (1.52) or
(1.53) is not satised. Let c = 2a1 e1 +2a2 e2 a3 e3 2a4 e4 +2a5 e5 a6 e6
be a general grade-1 element. The GIPNS results
82
a11
y a12
a13
a12
a22
a23
a13
1
a23 x = 0.
a33
y
a0
a2
a5
a2
a1
0
a5
0 2a1 e1 2a2 e2 +a3 e3 +2a4 e4 2a5 e5 +a6 e6 .
a4
a e6 = 2a4 = 0.
83
Clearly, this entity describes an inner product line and every line
passes through e3 , e6 , and e3 + e6 . An algebra element that contains
just e3 or e6 is a parabola whose axis are parallel to the x-axis or the y axis. An element that contains e3 +e6 , but neither e3 nor e6 , is given by
the condition a (e3 + e6 ) = 2a1 + 2a4 = 0. This means a1 = a4 and the
corresponding conic is an equilateral hyperbola, i.e., the asymptotes
enclose an angle of 90 . All other conics in principal axes position
can be obtained by the wedge product of four embedded points or by
using the bijection between conics and the algebra elements, cf. Eq.
(1.55).
Remark 1.21. Note, that this description of conics also contains
conics without real points.
In the most general case two-blades correspond to inner product point
quadruples. This can be seen from
NIG (a b) = NIG (a) NIG (b),
84
NIG (c) = (x, y)T R2 | x2 + 1 y 2 = 0 .
With the bijection (1.55) we see easily that c is the image of the
conic given by the diagonal matrix diag(1, 1, 1), i.e., the unit circle
centered at the origin.
1.10.3 Transformations
In this section we discuss transformations in this algebra. The Clifford algebra C(4,2,0) corresponds to the quadratic space R(4,2) , and
therefore, the sandwich action of vectors represents reections in hyperplanes in this space. It is clear that the transformations induce
transformations that are not linear in the base space R2 because the
embedding is quadratic. From the last section we know that the
geometric entity corresponding to vectors are conics in principal position. Furthermore, a transformation acts via the sandwich operator
and results again in a k-blade when applied to a k-blade k = 1, . . . , 4.
We begin with an example:
Example 1.12. Let c be the circle from Ex. 1.11
c = 4e1 e3 + 4e4 e6
1
2 e3
and let p = e1 + e2 +
+ e4 + 2e5 + 2e6 be (1, 2)T . Applying the
sandwich operator to p results in
1
p = (c)pc1 = 5e1 + e2 e3 + 5e4 + 2e5 + e6 .
2
p y = 6.
2
85
5
5
1
(x, y)T R2 | x2 + x y 2 + 2y = 0 .
2
2
2
to itself, and nd
(a)aa1 = (a) = a.
86
This two-blade represents the pair of common points of the conic and
l(x, y). The application of the transformation induced by a to l(x, y)a
results in
(a)(l(x, y) a)a1 =4a1 ye12 4a4 ye24 4a1 xe15 + 4(a5 y + a2 x)e25
2a3 ye23 2a3 xe35 2a6 ye26 4a4 xe45 + 2a6 xe56 .
This shows, that all pairs of common points of the pencil of lines with
the conic are xed, and therefore, the whole conic is xed pointwise.
To see that the transformation is an involution we have to apply it
twice to an arbitrary k-blade B
(a)(a)Ba1 a1 = (a2 )B(a2 )
= B.
87
(1.55)
Note that this element represents a pair of points, since lines intersect
also in the point . We parametrize ane lines as the sets of all lines
passing through two points p1 = (x1 , y1 )T and p2 = (x2 , y2 )T . The inner
product line is derived by
l(p1 , p2 ) = I ((p1 ) (p2 ) e3 e6 )
(1.56)
88
just show this theorem for non degenerate conics. A conic with no
terms in x or y is given by
1
1
a = 2a1 e1 + a0 e3 + 2a4 e4 + a0 e6 .
2
2
1 0
0 a1
M = a0
0 0
0
.
a4
a0
e4 +2(x1 x2 )e5 .
a0
(a)l(p1 , p2 )a1 =
a0 (x1 x2 )
0
a0 (y2 y1 )
1
.
N(p1 , p2 ) =
0
a0 (y2 y1 ) 2c1 (x1 y2 y1 x2 )
a0
0
2a4 (x1 y2 y1 x2 )
a0 (x1 x2 )
1 0 0
a (y1y2)
p 1 0 p, with = 2a (x0 y y x ) and =
1 1 2
1 2
0 1
!
a0 (x1x2)
2a4 (x1 y2 y1 x2 ) .
=:T
89
2
2
a2
0 (a4 (y1 y2 ) +a1 (x1 x2 ) )
2a4 a1 (x1 y2 y1 x2 ))
2a1 (x1 y2 y1 x2 )
a0
2a4 (x1 y2 y1 x2 )
a0
If we look at the ane part of the conic, we see that this results in
1
N (p1 , p2 ) = 0
0
0
ka1
0
0
2
4a2
1 (x1 y2 y1 x2 ) a4
.
0 , with k = 2
2
a0 (a4 (y1 y2 ) +a1 (x1 x2 )2 )
ka4
25 2 16 2
x + y = 1,
16
25
3 2
2
a4 : x y = 1.
4
a2 :
1.10.5 Subgroups
In this section we examine some subgroups that are embedded naturally in the Pin and the Spin group of Q2GA.
Rotation First, we concentrate on the group that is generated by
inversions with respect to lines passing through the origin. Therefore,
we study the action of these mappings applied to points embedded
via (1.55). Two lines through the origin may be represented by
l1 = I ((0, 0) (cos , sin ) e3 e6 ) ,
l2 = I ((0, 0) (cos , sin ) e3 e6 ) .
90
l2
l2
c2
c1
e3
e1
c2
e2
o
a1
l1
l3
l1
a2
l3
h3
l2
l2
p2
h1
l1
h2
p1
a3
a4
p3
l3
l1
l3
Furthermore, we are interested in orientation preserving transformations, i.e., elements from the Spin group. The composition of two
reections in l1 and l2 is given by their geometric product
l1 l2 = (sin sin + cos cos ) + (cos sin sin cos )e25
and with the addition theorems for sine and cosine we conclude
l1 l2 = cos( ) + sin( )e25 .
(1.57)
91
The GIPNS of this entity can be computed or we can simply read the
coordinates of the image point.
x = 2 cos()2 x0 + 2 cos sin y0 x0 = cos(2)x0 + sin(2)y0 ,
y = 2 cos()2 y0 2 cos sin x0 y0 = cos(2)y0 sin(2)x0 .
92
T(, t1 , t2 ) = l1 (, t1 )l2 (, t2 )
93
1
1
e25
e23 + e26
e35 e56
e25
e25
-1
(e35 e56 )
e23 + e26
e23 + e26
e23 + e26
e35 e56
e35 e56
e35 e56
(e23 + e26 )
0
0
0
0
1
c0 (e3 + e6 ) + 2c1 e1 + 2c2 e4 .
2
2y0 c1
2x0 c1
2y0 c2
2x0 c2
e12 +
e15 2e25 +
e24 +
e45 .
c0
c0
c0
c0
This is not of the form (1.55), and therefore, it is not the representation of the intersection of two lines. The GIPNS is calculated
NIG (p ) = (x, y)T R2 | 2c1 (xy0 x0 y)e1 (y0 c1 x2 +2yc0 +y0 c2 y 2 )e2
2c2 (xy0 x0 y)e4 + (x0 c1 x2 +x0 c2 y 2 +2xc0 )e5 = 0
2c0 x0
,
c1 x20 + c2 y02
y=
2c0 y0
.
c1 x20 + c2 y02
94
Q=
..
.
D
0 0 1
D = 0 1 0 .
1 0 0
n-times
1
(1.59)
V,
1
1
(x1 , . . . , xn ) e1 +x1 e2 + x21 e3 +. . .+e3n2 +xn e3n1 + x2n e3n .
2
2
(P )2x1 = 0,
(P )2x2 = 0, . . .
(P )2xn = 0.
I=
i
i=1
mod 3=2
n
ei
j
j=1
mod 3=1
n
ej
k
k=1
mod 3=0
ek .
95
Inner product null spaces can be mapped to outer product null spaces
with the blade
n
I :=
i=1
mod 3=2
n
ei
j
n
ej
j=1
mod 3=0
ek .
k=1
mod 3=1
Grade-1 elements correspond to inner product hyperquadrics in principal position. As the dimension is growing, the number of objects
that can be represented grows, too. Blades of grade k, (k n) correspond to the intersection of k hyperquadrics.
Quadrics in three-dimensional Space To construct the quadric geometric algebra for the three-dimensional space we use the quadratic
space R(6,3) given by the nine-dimensional real vector space R9 together with the quadratic form
Q=
D
D
0 0 1
D = 0 1 0 .
1 0 0
For three dimensions the embedding , see Eq (1.59), has the following form
: R3
1
V,
1
1
1
(x, y, z)T e1 + xe2 + x2 e3 + e4 + ye5 + y 2 e6 + e7 + ze8 + z 2 e9 .
2
2
2
(P )2x = 0,
(P )2y = 0,
(P )2z = 0.
96
determined by six values. This can be seen from the symmetric matrix of the equation of the quadric that has, in general, ten free entries.
The fact that we are treating quadrics in principal position reduces
the number of free entries to seven. Furthermore, this matrix representation is homogeneous, and therefore, we have six degrees of
freedom. In analogy to Eq. (1.55) we obtain a bijection that is
dened as
1
1
1
C 2c1 e1 2c2 e2 + c0 e3 +2c4 e4 2c5 e5 + c0 e6 +2c7 e7 2c8 e8 + c0 e9 ,
3
3
3
with
c0
c 2
C=
c5
c8
c2
c1
0
0
c5
0
c4
0
c8
0
.
0
c7
As we did for the planar case, we can now pay our attention on the
intersection of three planes in order to get a pair of points containing one ane and one ideal point. We choose these planes to be
parallel to the coordinate planes and passing through a given point
P = (x0 , y0 , z0 )T . Expressed in terms of the quadric geometric algebra
C(6,3,0) we get
p = ((x0 , y0 , z0 )T (0, 0, z0 )T (x0 , 0, z0 )T e3 e6 e9 ) I
((x0 , y0 , z0 )T (0, y0 , 0)T (x0 , y0 , 0)T e3 e6 e9 ) I
((x0 , y0 , z0 )T (x0 , 0, 0)T (x0 , y0 , 0)T e3 e6 e9 ) I
= 3e258 + (e358 e568 + e589 )x0 + (e238 + e268 e289 )y0
+ (e235 + e256 + e259 )z0 .
Remark 1.26. A representation of the group of Euclidean displacements SE(3) can be obtained by studying the composition of reections in planes. Planes correspond to vectors that are obtained by
(P1 P2 P3 e3 e6 e9 ) I.
Now we dene an inner product inversion quadric with
9
, 0, 0)T ,
10
3
P4 = (0, , 0)T ,
4
P1 = (
9
, 0, 0)T ,
10
5
P5 = (0, 0, )T ,
4
P2 = (
3
P3 = (0, , 0)T
4
5
P6 = (0, 0, )T
4
97
a
a
Figure 1.3: Inversion of the unit cube with respect to an ellipsoid (left),
inversion with respect to a hyperboloid (right).
a = (P1 P2 P3 P4 P5 P6 ) I
3
3
36
3
25
e1 e3 + 4e4 e6 + e7 e9 .
=
9
8
8
25
8
100
16
16
x2 + y 2 + z 2 1 = 0 .
(x, y, z)T R3
81
9
25
In Fig. 1.3 (left) we can see the inversion of the unit cube [1, 1]3
with respect to the ellipsoid dened by a. The image of every face
of the cube, i.e., of every plane is an ellipsoid passing through the
origin. The action of the inversion given by a on pairs of points can
be written as
f (x, y, z) =
452
T
(x, y, z) .
22 (302 y 2 + 252 x2 + 182 z 2 )
Note that the point is mapped to the origin and that the origin is
mapped to . Therefore, we have a map from R3 \ {0} R3 . One
main advantage of this method is that we can calculate the image
ellipsoid of one face of the cube (one plane) directly by applying the
98
sandwich operator to the plane that is expressed as quadric. For example the plane passing through P1 = (1, 1, 1)T , P2 = (1, 1, 1)T , P3 =
(1, 1, 1)T can be expressed as vector by
p = (P1 P2 P3 e3 e6 e9 ) I = 3e2 + e3 + e6 + e9 .
200
32
96
e1 3e2 e4 e7
27
3
25
with GIPNS
NIG (c) =
102
42
42
(x, y, z)T R3 | 2 x2 x + 2 y 2 + 2 z 2 = 0 .
9
3
5
This is one of the ellipsoids displayed in Fig. 1.3 (left). Furthermore, we can intersect two inner product planes p1 , p2 to get an inner
product line that is an edge of the cube. After that we can apply
the sandwich operator to the line and get the intersection curve (an
ellipse) of the two ellipsoids that are the images of p1 , p2 .
Remark 1.27. It is more convenient to compute the sandwich operator by a conjugation instead
inversion. This means we use (a)Xa
k
with a C(p,q,r) and X V . In general, the modied sandwich operator is easier to handle, because computing the inverse of a Cliord
algebra element is extremely expensive. Moreover, we are working in
a projective setting, and therefore, multiplication with a homogeneous
factor does not change the occurring geometric inner product and
outer product null spaces.
The second example in three-dimensional space is a hyperboloid of
two sheets in principal position that is generated by
9
9
a = (P1 P2 P3 P4 P5 P6 ) I = 6e1 + e3 + e4 + e6 + e7 + e9 .
8
8
99
a
a
Here, we have P1 = (1, 0, 0)T, P2 = (1, 0, 0)T, P3 = (2, 0, 4)T, P4 = (2, 0,4)T,
P5 = (2, 4, 0)T, P6 = (2,4, 0)T . We calculate the GIPNS
NIG (a) =
3
3
(x, y, z)T R3 | x2 1 y 2 z 2 = 0 .
16
16
16x2
16
(x, y, z)T .
3y 2 3z 2
The image of the cube [1, 3] [1, 1] [1, 1] under this mapping
is shown in Fig. 1.3 (right). Fig. 1.4 (left) shows the image of an
inversion with respect to a cylinder given by a = 3e12e3+3e42e62e9
applied to a cube. The equation of the cylinder is derived as x2+y 2 = 4.
Planes that are not parallel to the axis of the cylinder are mapped to
paraboloids. Another example is presented in Fig. 1.4 (right). The
inversion quadric is an elliptic paraboloid given by a = 3e12e3+12e4
2e6 +6e8 2e9 respectively by x2 +y 2 4z +4 = 0.
102
103
a
R(a, b) | c, d :
c
b
d
GL(R, 2) .
Such pairs are called admissible, see [32]. Note, that we dene points
of the projective line as left-homogeneous equivalence classes.
Furthermore, we dene equivalence classes of admissible pairs, respectively points.
Denition 2.7. Two points (a, b), (c, d) P1 (R) are equivalent if there
is a unit r R with ra = c and rb = d. This relation is an equivalence
relation and denoted by . For an equivalence class we write R(a, b).
The homogeneous component of the pair is the second coordinate.
This means the algebra element 0 has the projective coordinates U =
R(0, 1), the neutral element 1 has the coordinates V = R(1, 1), and
the projective point W = R(1, 0) corresponds to an ideal point. In the
following we denote points of the projective line P1 (R) with capital
latin letters.
Remark 2.2. If R is a eld, Def. 2.6 describes the classical projective
line.
From Def. 2.6 we see, R(1, 0) is a point, since the 2 2 identity
matrix is in GL(R, 2). If we take a general matrix M GL(R, 2) all
points P P1 (R) can be obtained by R(1, 0) M. Hence, the point set
P1 (R) can also be described as the orbit of R(1, 0) under the action of
GL(R, 2). Two points X = R(x1 , x0 ), Y = R(y1 , y0 ) are called distant
(XY ) if they are complementary, this means if R2 = X Y . Hence,
(P1 (R), ) is a distance space called the projective line over R.
Remark 2.3. We consider only algebras where the left-inverse elements are also right-inverse elements. Thus, every point may be
represented by an admissible pair, see [32].
104
q
0
0
| q C(R) R .
q
Note, that the kernel is equal to the center of the general linear group
ker GL(R, 2) = C(GL(R, 2)).
105
addition
x x + t
1 0
R(x, 1)
t 1
x xq
q 0
R(x, 1)
0 1
x qx
x x1
1 0
0 1
R(x, 1)
R(x, 1)
0 q
1 0
106
1
1
P
(L) of P (R) is called a chain in P (R), if is induced by a matrix
a b
GL(R, 2). The set of all chains in P1 (R) is denoted by C(L, R)
c d
and dened as
C(L, R) := P1 (L) PGL(R, 2) .
We call the incidence structure (L, R) := (P1 (R), C(L, R)) chain geometry over the L-algebra R.
107
108
Remark 2.7. All chain geometries over algebras are also chain spaces,
see [11] or [33]. When we talk about chain spaces we mention chain
geometries over algebras. A chain geometry over an L-algebra R is
denoted by (L, R), where the set of points is given by P = P1 (R) and
the set of chains by
C(L, R) := (P1 (L)) PGL(R, 2) .
109
Furthermore, the cross ratio is invariant under the action of PGL(R, 2):
cr(A, B, C, D) = cr(A , B , C , D ), for all PGL(R, 2).
110
Th. 2.6 allows a parametrisation of chains with the help of the cross
ratio. Therefore, we take a closer look at the cross ratio. The map
PGL(R, 2) that takes three arbitrary mutually distant proper
points A, B, C P1 (R) to the points U, V, W can be constructed in the
following way:
M = M 1 M2 M3 M4 =
tv
1
v 0
1 0
0 1
,
=
(1 ct)v c
0 1
t 1
1 0
1 0
c 1
tv
1
R(a, 1)M = R(a, 1)
(1 ct)v c
= R(atv + (1 ct)v), a c)
P1 (R):
tv
(1 ct)v
1
c
111
R(x, 1)
We repeat this calculation for an ideal point R(1, n), n N (R), where
N (R) denotes the set zero divisors of R.
tv
R(1, n)
(1 ct)v
1
c
= R (tv + n(1 ct)v, 1 nc)
112
(2.2)
With the upper calculation we are able to compute the cross ratio as
f (X) = cr(A, B, C, X) = R(x1 , x0 )M = R(, ).
R(, ) R (x a)(b a)1 (b c)(x c)1 , 1 ,
R(, ) R (1 na)(b a)1 (b c)(1 nc)1 , 1 .
(2.3)
(2.4)
R(, ) R 1, (x c)(b c)1 (b a)(x a)1 ,
R(, ) R 1, (1 nc)(b c)1 (b a)(1 na)1 .
(2.5)
(2.6)
Remark 2.9. For a proper cross ratio R(, 1) we can calculate the
cross ratio directly from the ring elements by
cr(a, b, c, x) = (x a)(b a)1 (b c)(x c)1 .
(2.7)
This is exactly the cross ratio formula that is used for complex numbers
or in classical projective geometry. If C(R) the cross ratio is an
element and not a conjugacy class.
With the help of Th. 2.6 and the formulas (2.3)-(2.6) we can now
parametrise chains in the chain geometry (L, R)
(d a)(b a)1 (b c)(d c)1 =
(d a)(b a)1 (b c) = (d c)
d(b a)1 (b c) a(b a)1 (b c) = d c
d(b a)1 (b c) d = a(b a)1 (b c) c
d (b a)1 (b c) = a(b a)1 (b c) c
1
.
d = a(b a)1 (b c) c (b a)1 (b c)
(2.8)
113
Thus, we can use chains parametrised by the cross ratio for interpolation purposes, since the elements a, b, c are interpolated for the cross
ratio values 0, 1, .
Theorem 2.7. Let a, b, c, d R be four mutually distant elements
incident with one chain. The exchange of two pairs of elements does
not change the value of the cross ratio:
cr(a, b, c, d) = cr(b, a, d, c) = cr(c, d, a, b) = cr(d, c, b, a).
cr(a, b, c, d) = ,
cr(a, c, d, b) = (1 )
cr(a, c, b, d) = ( 1)1
(2.9)
cr(b, a, c, d) = 1
cr(a, d, b, c) = ( 1)1 .
114
115
1
((x + y) (x) (y)).
2
116
and becomes the chain set. Hence, we examine the incidence structure
(Q) := (Q , C(Q)) .
117
118
2.5 Jordan-Systems
Now we take a look at sub structures of chain geometries.
Denition 2.18. Let = (P, C) be a chain space and let S P with
(1) S contains three mutually distant points.
(2) If p, q, r S are mutually distant, then the connecting chain c =
(pqr) is entirely contained in S . The set of all connecting chains
is denoted by C(S).
(3) The incidence structure (S, C(S)) is a chain space.
Under this conditions (S, C(S)) is a subspace of . If (S, C(S)) satises
just condition (2), (S, C(S)) is called a weak subspace of .
Trivial subspaces of = (P, C) are itself and subspaces containing
just a single chain c C. More examples can be generated with
subalgebras.
Denition 2.19. Let R be an algebra over the eld K. A subspace S
of R is called a subalgebra of R, if
(1) 1 S .
(2) For a, b S the product ab is in S .
(3) If a S R , then a1 S , i.e, S R = S .
Since chain geometries over algebras are chain spaces, the chain geometries over subalgebras are also chain spaces, and therefore, subspaces of the chain space. It is not common to demand condition
(3). For the description of subspaces of chain spaces this condition
is important. Just with this condition subalgebras correspond to
Jordan-systems. The following theorem states that subalgebras also
dene chain geometries on their own.
Theorem 2.10. Let R be a K-algebra and let S be a subalgebra of R.
The mapping : P1 (S) P1 (R) with S(a, b) R(a, b) is an injective
morphism of chain spaces (K, S) (K, R). The image ((K, S))
is a subspace of (K, R), that is embedded via and isomorphic to
(K, S).
2.5 Jordan-Systems
119
a, b, c, d R.
A subalgebra is dened by the set {a + cj | a, c R}. In fact this subalgebra is isomorphic to the complex numbers C and it denes a subspace
(R, C) of the chain geometry (R, H).
There are also subspaces of chain spaces that do not correspond to
subalgebras. To describe some of these subspaces we need the denition of a Jordan-system.
Denition 2.20. Let R be a K-algebra and let J be a subspace of
the K-vector space R with 1 J . The subspace J is called
J1 Jordan-system in R, if b J R b1 J for all b J , we
dene J = J R .
J2 Jordan-closed in R, if a, b J aba J for all a, b J .
J3 strong in R, if for all b J the condition |e(b)| > |K\e(b)| is
satised, where
e(x) := f K | x + f R .
120
1 0
Now we generalize the setting and assume that the set S does not
contain the three points R(1, 0), R(0, 1), and R(1, 1). Under these
assumptions S cannot be the projective line over a Jordan-system,
see Def. 2.20. Nevertheless, we can formulate the following theorem,
see [11]:
121
122
Benz [5] dened a contact relation with the help of the residual spaces
without mentioning them explicitly. We also give his denition since
it is more intuitive.
Denition 2.24. Let = (P, C) be a chain geometry, let W =
R(1, 0) P1 (R) and let a, b C. The chains a and b are in contact at the point P if there is a PGL(R, 2), such that P = W and
the lines a \ {W }, b \ {W } are parallel.
We give a theorem from [5]:
Theorem 2.14. Let aP b, a = b. It follows a b = P . The contact
relation p is an equivalence relation, i.e., it is reexive, symmetric,
and transitive.
/ c be points
Contact theorem: Let c C be a chain and let p c, q
with p q . Then there is exactly one chain c incident with p and q
that is in contact with c.
Lemma 2.4. If aP b and PGL(R, 2), then a P b .
The inverse implication is not true in general, but, if there is a socalled beetle-gur (german: Kafergur), then it is true. If = (P, C)
is a chain space, then there is necessarily a beetle-gure. Now we give
two theorems that connect the contact relation with the cross ratio.
Both theorems can be found in [5].
Theorem 2.15. Let R be a K-algebra over a eld and let a, b C
be chains incident with P P1 (R). The following expressions are
equivalent:
(1) aP b,
(2) For two dierent points A, A a\ {P } and for two points B, B
b\ {P }
cr(P, A, A , B) cr(P, A, A , B ) K,
123
A proof can be found in [5]. The proof given there is valid for commutative rings, but every statement can be transferred directly to the
non-commutative case. The proof uses Lemma 2.4 and the invariance
of the cross ration under the action of the group PGL(R, 2). With this
theorem we are able to calculate an explicit parametrisation of the
chain that is in contact to a given one.
a
R(a, b) | c, d :
c
b
d
GL(R, 2) ,
124
where R denotes the Cliord algebra. The chains are dened as Rchains if not stated otherwise:
C(R, R) :=
P1 (R) | PGL(R, 2) .
125
(2.10)
126
e1 0
0 1
=
PGL(C(p,q,r) , 2) that corresponds to a right-multiplication with e1 together with Lemma 2.1 and apply to I . This results in a new subset
that denes a chain geometry and is given by J = I =span(e1 , . . . , en ).
Therefore, I is isomorphic to a Jordan-system J and the isomorphism
is given by .
127
J3 We show |e(x)| > |R\e(x)|. The set R\e(x) is the set of all elements
r R with a r N , a P in(p,q,r) . An element g of a Cliord
algebra is a zero divisor if gg = 0. Hence, we look at the equation
(a r)(a r) = 0.
C
(p,q,r) (g, 1), g Pin(p,q,r) is completely contained in the embedding of
1
Pin
(p,q,r) in the projective line P (C(p,q,r) ).
Proof. We dene a mapping PGL(R, 2) with R(1, 0) R(a, 1),
R(1, 1) R(b, 1), and R(0, 1) R(c, 1). Since the group PGL(R, 2)
acts ag-transitive on the set F, the standard chain R(s, r), s, r
R is mapped to another chain. Thus, we have to show that every
element of the chain P1 (R) corresponds to an element R(g, 1) with
v 1 c v 1
1 tc t
128
0 1
g 0
1 0
,
, and
with g Pin
(p,q,r) . This group acts trans1 0
0 1
0 g
itive on P1 (Pin
(p,q,r) ), and therefore, we can assume that a = 1. To
show that
With Th. 2.18 each connected component of the Pin group denes a
weak subspace and with Th. 2.13 a subspace of (R, C(p,q,r) ).
Remark 2.22. Th. 2.18 also holds for null vectors. Thus, the null
vectors
N 1 :=
V | vv = 0
129
b = e1 ,
c = e2
1
.
2(1)e
+2e
+2(1)e
+
2(1)e
0
1
2
12
22 2+2
We can check, that the chain is entirely contained in the Pin group
cabc ()cabc () = 1,
for all R.
2e0 + 2e12 ,
b = 2e1 ,
c = 2e2
are the elements of Ex. 2.2 multiplied by 2. Note, that the multiplication by 2 can also be expressed as mapping PGL(2, C(0,2,0) ). Now
the elements a, b, c satisfy aa = bb = cc = 4. Again we compute the
chain with Eq. (2.8)
cabc () =
1
.
2(1)e
+2e
+2(1)e
+
2(1)e
0
1
2
12
2 +1
This is exactly the chain of Ex. 2.2 multiplied by 2. Again all points
of cabc () satisfy
cabc ()cabc () = 4,
for all R.
130
q02
1
(q0 e0 q1 e1 q2 e2 ) J .
+ q12 + q22
v 1 c v 1
, with t = (c a)1 , v = (t + (b c)1 )1 ,
1 tc t
131
1
(2.11)
(1q0 +q1 )e0
1+q1 q0 +q22 +q12 +q02
2
(q1 q0 +q22 +q12 +q02 + 2q2 )e2
+ (q0 +q1 )e1 +
2
2
(q1 q0 +q22 +q12 +q02 2q2 )e12 .
+
2
c+(qv 1 t)1 =
c+(qv 1 t)1
c+(qv 1 t)1 = 1, q J .
132
Q=
(2.12)
(2.13)
b = e2 + e6 ,
c = e3 + e 5 .
cabc () =
2 2+2
133
(2.14)
Equating the coecients of Eq. (2.13) and Eq. (2.14) shows that
both
representations describe the same set of points if we identify 1 V
with P5 (R). Furthermore, the chain interpolates the null vectors a, b, c
for the cross ratio values = , 0, 1. A visualization of the resulting
regulus together with its striction curve (blue) is shown in Fig. 2.1.
b
a
Figure 2.1: Regulus corresponding to Eq. (2.13) and Eq. (2.14)
134
b
c
a
0
0
0
,
0
0
1
The subspace 1 V corresponds to the
1
0
0
Q=
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
52
1
2( + 11)e1 + (52 6 20)e2 + 12e3
4 + 2
135
+ 10(1 )e4 + (52 6 + 4)e6 .
52
1
T
(12, 10(1 ), 0)
4 + 2
52 6 + 4
.
52 4 + 2
Q=
O
I
I
,
O
where O is the 44 zero matrix and I is the 44 identity matrix. Similar to the examples for Kleins and Lies quadric we give an example
of a conic on Studys quadric as chain in the homogeneous Cliord
algebra model.
Example 2.7. Let A = (1, 0, 0, 0, 0, 1, 1, 1)T R, B = (0, 1, 0, 0, 0, 0, 1, 1)T R
and C = (0, 0, 1, 0, 1, 0, 0, 0)T R be three points on S26 P7 (R). In analogy
to Eq. (2.12) the conic that interpolates A, B , and C for the parameter
values t = , 0, and 1 can be parametrised by
c(t) = t(A, B)C + (1 t)(A, C)B + t(t 1)(B, C)A,
where denotes the bilinear form belonging to Studys quadric. Especially for the chosen points we derive:
136
0
0
1
t(t 1)
0
1
0 2(1 t)
1
0
0
t
0
0
0
0
.
1
0
0
t
0
0
1 t(t 1)
0
1
1 t2 3t + 2
0
t2 3t + 2
b = e2 +e7 +e8 ,
c = e3 +e5 .
The chain dened by these three null vectors is computed with Eq.
(2.8)
cabc () =
1
2 2+2
+e5 +(1)e6 +(2 3+2)e7 +(2 3+2)e8 .
137
1
Q = 0
0
0
1
0
0
0 .
1
Since we describe the dual unit sphere, the real part has to be equal to
1 and dual part has to vanish. We already showed that the connected
components of the Pin group dene subspaces of (C(C(p,q,r) ), C(p,q,r) ),
see section 2.7.2. Note that we use an ane Cliord algebra model
here. It is also possible to construct a homogeneous. We perform this
construction for the case of dual unit quaternions.
Figure 2.3: Line congruence with constant dual angle to an axis a (red)
138
0
b = 1 ,
c = 0
1
b = e2 + e3 ,
c = e1 + e3 .
Remember that the quadric is constructed over the ring of dual numbers
D. Therefore, the cross ratio may be a dual number = + . The
chain, computed with Eq. (2.8), has the form
1 + (2 1)e1
( 2 1 + 2 2 + + )
+ (3 1)e2
+
( 2 1 + 2 2 + + )
2 + + (2 2 + 1)e3
+
.
2 1 + 2 2 + + )
cabc () = cabc (, ) =
We split the dual numbers into real and dual parts and get
1
(3 22 + + 2)
+
e1
cabc (, ) =
(2.15)
2 + 1
(2 + 1)2
(1 )(22 2 + + 1 + )
+
e2
+
2 + 1
(2 + 1)2
( 1)
3 + 1 + 2
+
e3 .
+
2 + 1
(2 + 1)2
Study [61] and Klein [43] already knew that the intersection of
the dual sphere SD2 with two-dimensional ane submodules results in
congruences of lines, i.e., two-parameter families of oriented lines that
enclose a constant dual angle with an unique oriented line. Therefore,
Eq. (2.15) is a parametrisation for this congruence of lines, see Fig.
2.3 for a schematic view. This congruence can be interpreted as the
set of all tangents of a cylinder of revolution that enclose a constant
angle with its axis.
139
C+
(3,0,1) = C(2,0,1) . Furthermore, the algebra is of the form L + sL,
2
with s = 0, since C(2,0,1)
= C(2,0,0) D. For our investigations it is
important that the norm of a general element A C+
(3,0,1) is contained
+
(2.16)
Hence, we have the two conditions a20 +a21 +a22 +a23 = 1 and c0 a0 a1 c1 +
c2 a2c3 a3 = 0. In analogy to Ex. 2.8 we can construct an ane Cliord
algebra model that describes the group Spin(3,0,1) as the intersection
of the grade-1 subspace with its Pin group.
140
1
0
Q=
0
0
0
1
0
0
0
0
1
0
0
0
.
0
1
1
0
Q=
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
.
0
1
141
C(4,0,0) :
If this is an element of the Spin group, the real part has to equal
1 while the pseudoscalar part has to vanish. The equation in the
pseudoscalar part can, again, be understood as the equation of Studys
quadric. Therefore, the same set of points in projective space can
serve as a point model for dierent geometric entities, see chapter 3.
Ane Cliord Algebra Model Since the center of C(4,0,0) is isomorphic to A we construct the ane model by the use of M = A4
and the quadratic form
1
0
Q=
0
0
0
1
0
0
0
0
1
0
0
0
.
0
1
142
Homogeneous Cliord Algebra Model To construct the homogeneous Cliord algebra model we use M = A5 as model for the projective
space P4 (A) and the quadratic form
1
0
Q=
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
.
0
1
Like in the case of dual quaternions the unit double quaternions are
in this model represented as null vectors.
C(3,1,0) :
1
0
Q=
0
0
0
1
0
0
0
0
1
0
0
0
0
1
1
143
Homogeneous Cliord Algebra Model For the homogeneous quadric model of Spin(3,1,0) we use V = C5 as vector space model for P4 (C).
The corresponding quadratic form is given by
1
0
Q=
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
.
0
1
Remark 2.28. The Cliord algebra C(3,1,0) results in the same model.
We have to use the quadratic form Q = diag(1, 1, 1, 1) for C(3,1,0) .
If we use the standard form Q = diag(1, 1, 1, 1), the equality is not
obvious.
C(2,2,0) :
We proceed analogue to the cases above and determine the ane and
the homogeneous Cliord algebra models for this Spin group.
Ane Cliord Algebra Model As module we use M = A4 . The
corresponding quadratic form is computed as
1
0
Q=
0
0
A general vector v
0 0 0
1 0 0
.
0 1 0
0 0 1
1
144
1
0
Q=
0
0
0
C(2,1,1) :
0 0
1 0
0 1
0 0
0 0
0 0
0 0
0 0
.
1 0
0 1
Ane Cliord Algebra Model The ane Cliord algebra is constructed with M = D4 together with the quadratic form
1
0
Q=
0
0
0
0 0
1 0 0
.
0 1 0
0
0 1
A vector has the form v = (a0 + c0 )e1 + (a1 + c1 )e2 + (a2 c2 )e3 +
(a3 c3 )e4 and its square reads
vv = a20 a21 a22 + a23 + 2(c0 a0 a1 c1 + c2 a2 c3 a3 ).
1 0 0
0 1 0
Q=
0 0 1
0 0 0
0 0 0
0 0
0 0
0 0
.
1 0
0 1
145
a
rk 0
c0
a1
c1
a2
c2
= 3.
146
e12
0
0
1
x1
x2
x3
e0 + 2
e23 2
e13
2
2
2
2
2
x1 + x2 + x3
x1 + x2 + x3
x1 + x22 + x23
x21 x4 2x3 x1 x6 + x4 x22 x23 x4 + 2x3 x5 x2
e24
(x21 + x22 + x23 )2
2x2 x6 x1 + 2x2 x3 x4 x22 x5 + x5 x21 + x23 x5
e14
(x21 + x22 + x22 )2
2x3 x1 x4 + x6 x21 x23 x6 2x2 x5 x1 x22 x6 )e1234
I.
(x21 + x22 + x22 )2
We compute
ABA = (x21 y1 2x3 x1 y4 2x1 x2 y2 x22 y1 x23 y1 )e0
+ (x22 y4 + x21 y4 + 2x3 x1 y1 x23 y4 2x2 x3 y2 )e13
(2x2 x3 y4 + 2x1 x2 y1 x22 y2 + x23 y2 + x21 y2 )e23
+ (2x1 x2 y6 + x23 y5 + x21 y5 2x5 x2 y2 2x4 x3 y2 2x5 x3 y4
2x2 x6 y1 + 2x1 x5 y1 2x1 x6 y2 + 2x2 x3 y3 + 2x4 x2 y4 x22 y5 )e14
(2x1 x4 y1 + 2x3 x1 y6 + 2x3 x6 y1 + x22 y3 + 2x2 x3 y5 2x4 x3 y4
x23 y3 + x21 y3 2x5 x2 y4 + 2x5 x3 y2 + 2x1 x6 y4 2x4 x2 y2 )e24
(2x1 x2 y5 + x21 y6 x23 y6 2x2 x6 y2 2x4 x3 y1 x22 y6 + 2x5 x2 y1
2x1 x4 y4 + 2x1 x5 y2 2x3 x6 y4 2x3 x1 y3 + 2x1 x6 y1 )e1234 I.
147
Remark 2.32. With the identity C(2,0,1) = C(2,0,0) D, see Selig [59,
sect. 9.3], we see that the dual unit quaternions dene a sphere SD3 in
four-dimensional dual space.
PGL(C(3,0,1) , 2) as
148
3
i=0
a2i = 1 and
i=0
149
the dual angle and se(3) denotes the Lie-algebra of SE(3). Then the
chain can be expressed with the sandwich operator
cabc (, d) = ea aea .
Real Cross Ratio Now we ask for a real cross ratio of four dual
unit quaternions. First we restrict ourselves again to the case of linesymmetric displacements. For three given line-symmetric displacements a, b, c respectively lines the corresponding R-chain is a closed
ruled surface, since cr(a, b, c, a) = 0, cr(a, b, c, b) = 1 and cr(a, b, c, c) = .
Every normalized dual vector d of the ruled surface satises l , d
=
, where D is a constant dual angle and l is the axis of the
cylinder that denes the congruence, see Th. 2.20. This means that
the ruled surface is a closed ruled surface of constant slope. In general
its intersection with the ideal plane P is a conic with multiplicity two
that has two complex conjugate contact points with multiplicity four
with the absolute circle. This can be derived with the help of Pl
ucker
coordinates, see section 2.10.2. Hence, the order of the striction curve
is reduced by six, see [51, p. 152]. The parametrisation of the
congruence of lines in terms of Pl
ucker coordinates has degree four.
Therefore, the degree of the striction curve is two,
i.e., it is a planar
curve. In fact, the R-chain dened by a, b, c 2 V corresponds to
the set of line-symmetric displacements with respect to a closed ruled
surface of constant slope whose striction curve is an ellipse, see Fig.
2.4 for example. Note that the striction curve is dened by the planar
intersection of the cylinder that is the envelope of the congruence of
lines. The plane is dened by the contact points of the three given
lines with the cylinder.
150
2(1 + 2)
4(2 + 1)2
e0 +
e14 ,
2+
( + 2)2
151
2(1 + 2)
e14 .
2+
b = e13 + e34 ,
c = e12 + e14 .
Note, that the minus in front of e13 is due to historic reasons, see Ex.
1.6. We compute the chain cabc (), = e0 + e1234 C(C(3,0,1) )
= D,
with Eq. (2.8)
(1 )e23
e13
( 1)e12
+
1 + 2
1 + 2
1 + 2
3
2
( 2 + + 2)e14
(1 2 + 3 + )e34
+
+
(1 + 2 )2
(1 + 2 )2
3
2
2
(2 4 + 3 + 1)e24
.
(1 + 2 )2
cabc () =
In Pl
ucker coordinates the line congruence can be written as
(1 + 2 )(1 )
2
(1 + )
(1 + )( 1)
.
l(, ) =
3
2
( 2 + + 2)
23 + 42 + 2 3 + 1
1 2 + 3 +
152
with the absolute circle determined by x21 + x22 + x23 = 0. This leads to
the equation
1 4 + 102 163 + 194 165 + 106 47 + 8 = 0
with solutions
1 = 2 = 3 = 4 =
i
i
1
1
3, 5 = 6 = 7 = 8 = +
3.
2 2
2 2
Hence, this congruence of lines has two points in common with the
absolute circle. At these points it has contact of order three.
Axis of the Congruence and constant dual Angle We are interested
in the axis of this congruence, i.e., in the dual vector a D3 with
a , cabc ()
= k1 + k2 , for all D
+ c4 l1 (, ) + c5 l2 (, ) + c6 l3 (, ) = 0.
(2.17)
153
1 0 0 0 0 1
0
2 2 1 1 1 2
0
c
1
1 1 0 0 0 0 c2 0
1 4 0 2 1 2 c3 0
MA=
2 0 2 0 0 0 c4 = 0 .
0 3 1 2 1 1 c5 0
0 1 1 0 0 0
0
c6
0 1 1 1 0 0
0
c c
4
= ,
c2
3
2 1 1
(a, a
) = (c, c pc) = (1, 1, 1)T , ( , , )T
3 3 3
=((3,
$
3, 3)T , (2, 1, 1)T ).
154
To compute the dual angle we switch the model to C(3,0,0) (D) and
rewrite the line congruence (see Eq. (2.15)) as:
1
(3 22 + + 2)
+
e1
2 + 1
(2 + 1)2
(1 )(22 2 + + 1 + )
+
+
e2
2 + 1
(2 + 1)2
( 1)
3 + 1 + 2
+
+
e3 .
2 + 1
(2 + 1)2
cabc (, ) =
1
3((3 2)e1 + (3 + )e2 + (3 + )e3 ).
9
1
4
3+
3,
3
9
(1 + 2 )(1 )
2
(1 + )
(1 + 2 )( 2)
.
l() =
3
2
( 2 + )
23 + 42 3 + 1
1 + 3
155
l
l
.
ll
s()
c
a
b
Figure 2.4: Ruled surface of constant slope with striction curve
+22 2+1
23 32 +2
1
.
4 23 +2 +1 +t
f (, t) = 2
3 2 +
( +1)2
3
2
4
3
2
+3 3+1
2 +2
1
,
+1
2 1,
2 2 + 1
T
156
i
c
Figure 2.5: Intersection of the ideal curves with the absolute circle.
B
l
157
B
l
B
l=c
A
C
158
(3) For every point D c(t1 : t0 ) the cross ratio cr(a, b, c, d) is a dual
number.
To show that each case occurs we give examples for every case.
Example 2.10. We take the chain from Ex. 2.8 given by Eq. (2.15).
This chain is generated by
a = e1 + e2 + e3 ,
b = e2 + e3 ,
c = e1 + e3 ,
with corresponding Pl
ucker coordinates given by A = (1, 0, 0, 0, 1, 1)T R,
T
B = (0, 1, 0, 0, 0, 1) R and C = (0, 0, 1, 1, 0, 0)T R. Now we parametrise
the conic on M24 dened by these three points with a homogeneous
parameter (t1 : t0 ).
c(t1 : t0 ) = t1 t0 (A, B)C + (t20 t1 t0 )(A, C)B + (t21 t1 t0 )(C, B)A
= (t21 t0 t1 , 2t20 2t1 t0 , t1 t0 , t1 t0 , t21 t1 t0 , t21 t1 t0 + t20 )T .
1
(3 22 + + 2)
cabc (, ) =
+
e1
2 + 1
(2 + 1)2
(1 )(22 2 + + 1 + )
+
+
e2
2 + 1
(2 + 1)2
( 1)
3 + 1 + 2
+
+
e3 .
2 + 1
(2 + 1)2
If we compare this with the results of section 2.10.2, we see that the
resulting chain is exactly the same. The intersection of c(t1 : t0 ) with
P is given by
c1 (t1 : t0 ) = (t21 t0 t1 , 2t20 2t1 t0 , t1 t0 )T .
159
1
0
0
x0
x0 + x1
y0
0
= y1 .
x1 =
x1
0
x2
x1 + x2
y2
1
1
1
1
With the implicit form of the standard conic we can now determine
the implicit form of c2 (t1 : t0 ) as
y0 y2 y12 = (x0 + x1 )(x1 + x2 ) x21 = x0 x1 + x0 x2 + x1 x2 = 0.
Now we are able to intersect both conics, i.e., we insert the parameter
form of c1 in the implicit form of c2 . This leads to the equation
t20 t21 t0 t31 = 0.
c = cabc (2, 0) =
The corresponding Pl
ucker coordinate vector is given by
T
C = (3, 6, 6, 4, 5, 7) R.
Of course, this results in the same chain, and therefore, in the same
ideal conic of the congruence of lines. The regulus dened by the
points A, B, C can be parametrised by
c(t1 : t0 ) = (4t0 t1 + t21 , 10t1 t0 + 16t20 , 6t1 t0 ,
160
This yields
(t1 : t0 )1 = (0 : 1), (t1 : t0 )2 = (1 : 0), (t1 : t0 )3 = (8 : 1), (t1 : t0 )4 = (1 : 1).
We see that there is a solution that does not correspond to the start
points A, B , or C . Therefore, we have
D = c(8 : 1) = 16(6, 6, 3, 2, 7, 10)T R.
1
((6 2)e1 + (6 + 7)e2 + (3 + 10)e3 )
16
1
((6 2)e1 + (6 + 7)e2 + (3 + 10)e3 ) .
9
2
1
+ .
3
9
161
Remark 2.37. For general conics on S26 there are two possibilities to
perform this calculation. The rst one is based on the transformation
presented in section 2.10. After the transformation is performed we
can proceed with line-symmetric displacements. The second method is
to take the conic that is dened by the rst four Study coordinates and
intersect it with the conic dened by the rst four Study coordinates
corresponding to the chain that is dened by the three points.
(0, 1, 0)T (0, 0, 1)T (0, 0, 1)T I
162
512
64
512
c = (4, 4, 0)T (6, 4, 0)T (5, 5, 0)T
(5, 3, 0)T (5, 4, 1)T (5, 4, 1)T I
= 6e1 30e2 40e3 6e4 24e5 40e6 6e7 40e9 .
The GIPNS of a is the unit sphere centered at the origin, the quadric
corresponding to b is the sphere centered at M = (5, 0, 0)T with radius
r = 23 and c is an inner product sphere with radius r = 1 and center
M = (5, 4, 0)T , see Fig. 2.7 (left). We compute the chain cabc dened
by the three elements:
5 (99 + 65)e2
1 (63 + 2 + 992 )e1
2
3 28 + 33 + 21
3 28 + 332 + 21
1 (63 + 2 + 992 )e4
1 (63 166 + 39602 )e3
18
28 + 332 + 21
3 (28 + 332 + 21)
cabc () =
163
b
a
2
28 + 33 + 21 18
(28 + 332 + 21)
1 (63 + 2 + 992 )e7
1 (63 166 + 39602 )e9
.
2
3 28 + 33 + 21
18
28 + 332 + 21
From Eq. (1.55) we see, that the geometric inner product null space
of a vector determines a sphere if the coecients of e1 , e4 and e7 are
equal. Thus, every entity contained in cabc corresponds to a sphere,
see Fig. 2.7 (right) for the cyclide that is the envelope of this set of
spheres. The GIPNS of the chain cabc () can be determined as
NIG (cabc ()) = (x, y, z)T R3 | 9902 x650x+63x2 +2x2 +992 x2
63166+39602 +63y 2 +2y 2 +992 y 2 7922 y
+792y + 63z 2 +2z 2 +992 z 2 = 0, R .
Furthermore, the chain is dened by three grade-1 elements corresponding to inner product spheres that can also be interpreted as Pin
group elements or as inversions in these spheres. Since the group of
conformal transformations is a subgroup of the group of reections in
quadrics in principal position it also denes a sub chain geometry, see
remark 2.23.
164
Furthermore, we did not normalize the elements since this does not
change the GIPNS and a calculation without normalization is much
faster. The chain cabc results in
15(1997 + 2147)e2
3(19972 817 + 2964)e1
cabc () =
165
2223( 1)e8
3(665 + 1482 + 19972 )e7
+
2
1726 + 1997 + 247
1726 + 19972 + 247
2
2(741 + 4340 + 19970 )e9
.
166
P0
P0
j
L0
j
L0
c0
c1
c2
L1 P1
c0
P1
L1
P1 = p1 R,
T0 = t0 R,
T1 = t1 R.
pT
1 Mp1 = 0,
pT
0 Mt0 = 0,
pT
1 Mt0 = 0.
d = p1 + t1 ,
, R.
167
B
T0
C
P0
P1
T1
D
(2.18)
2 T
2 T
= b Mb + d Md + 2b Md = 0.
b Mb bT Mb
)
(bT Md)2
bT Md
dT Md
.
= T
+ T
T
2
1,2
b Mb
(b Mb)
b Mb
t
Mt
t
1
1
1
1
0
1 .
0
0
0
0
0
1
168
T
If tT
0 Mt0 t1 Mt1 0 Eq. (2.19) is solvable and we can determine :
and : afterwards. The set of points that serves as connection
point is determined by the two conics through P0 and P1 that are
solutions of Eq. (2.18). Note, that the existence of a biarc cannot be
guaranteed in general. If both tangents are situated on dierent sides
of the quadric, there is no real point C = cR. This is the geometric
interpretation of the vanishing discriminant, see Eq. (2.19). If the
quadric Q is an oval one, there is no tangent in the interior of Q.
Thus, for oval quadrics we can always nd biarcs. Naturally, a biarc
is not determined uniquely. The set of all possible transition points is
contained in two dierent conics passing through P0 and P1 , see [56].
C(Cp,q,r ),
169
(2.20)
with
0
Q : xT Mx = 0, with X = xR P3 (R) and M =
0
0
0
1
0
0
0
0
1
0
0
0
.
0
1
170
(1 )e4
e2
( 1)e3
+
+
.
2
2
1+
(1 + ) 1 + 2
(1 + 2 )e2
((1 + 2) 2 2 1)e3
s() = e1
+
( 2 + 1) 2 2 2
( 2 + 1) 2 2 2
e4
.
( 2 + 1) 2 2 2
ford algebra. Therefore, let V = R3 and let Q = 0 1 0. The Clifford algebra has signature (p, q, r) = (3, 0, 0). The square of a general
grade-1 element a = a1 e1 + a2 e2 + a3 e3 is computed as a2 = a21 + a22 + a23 .
If the square is equal to 1 the element corresponds to a point on S 2 .
The points A = (0, 0, 1)T , B = (1, 0, 0)T , and C = (0, 1, 0)T that were
also used in the example above correspond to the elements a = e3 ,
b = e1 and c = e2 with a2 = b2 = c2 = 1. The chain cabc () dened by
these three vectors is given by
cabc () =
e1
( 1)e2
(1 )e3
+
+
.
2
2
1+
1+
1 + 2
171
(1 + 2 )e1
((1 + 2) 2 2 1)e2
s() =
+
( 2 + 1) 2 2 2
( 2 + 1) 2 2 2
e3
.
( 2 + 1) 2 2 2
172
1
0
T
3
Q : x Mx = 0, with X = xR P (R) and M =
0
0
0
1
0
0
0
0
1
0
0
0
.
0
2
e2
( 1)e3
1 2(1 )e4
cabc () = e1 +
+
+
.
1 + 2
1 + 2
2 1 + 2
( 2 1)e2
( 2 1 2 2 + )e3
cd () = e1 +
+ 2 2 2 2
+ 2 2 2 2
2e4
1
2 + 2 2 2 2
constructed
1 0 0
is determined as
cabc () =
( 1)e2
1 2(1 )e3
e1
+
+
.
1 + 2
1 + 2
2 1 + 2
173
( 2 1)e1
( 2 1 2 2 + )e2
cd () =
+
+ 2 2 2 2
+ 2 2 2 2
2e3
1
2 + 2 2 2 2
0
Q : xT Mx = 0, with X = xR P3 (R) and M =
0
0
0
1
0
0
0 0
0 0
.
1 0
0 1
(2 + 2) + 1 + 2
2 + 2 1 2
cabc () = e1 +
e2
e3 .
2 + 2 + 2 + 2
2 + 2 + 2 + 2
174
E
D
E
A
C
A
B
2(3 + 4 + 22 + 2 2 + 2 2)e1
4(2 2 + 3 + 2 + 2)e2
cd () =
(2 + 2 + 2)2
(2 + 2 + 2)2
2(2 2 3 + 4 + 22 + 2 2)e3
4(2 2 + 3)e4
+
.
(2 + 2 + 2)2
(2 + 2 + 2)2
2(12 + 8 2 2 + (4 2 8))e2
ce () = e1
2 + (2 2 4) 24 + 16 2
+
2 + (2 2 4) 24 + 16 2
2(36 + 24 2 2 + (6 2 12))e4
.
+
2 + (2 2 4) 24 + 16 2
175
Projective Setting The homogeneous Cliord algebra model is constructed with V = R4 as model for P3 (R) together with the quadratic
form Q with the coecient matrix
0 0 0
0 1 0
T
3
Q : x Mx = 0, with X = xR P (R) and M =
0 0 1
1
0 0
2
1
2
0
.
0
0
As example
we again choose three points on the paraboloid. Let A =
aR = (1, 2, 0, 2)T R, B = bR = (1, 1, 0, 1)T R, and C = cR = (1,0, 1, 1)T R
be three dierent points with algebra representation a = e1 + 2e2 +2e4 ,
b = e1 + e2 + e4 , c = e1 + e3 + e4 . We compute the chain cabc with Eq.
(2.8).
( + 8 + 6 2)e2
(1 + )e3
cabc () = e1 +
2
2
+ (2 2 + 2) + 6 + 4 2 + (2 2 + 2) + 6 + 4 2
(2 + (2 2 + 4) 12 8 2)e4
2 + (2 2 + 2) + 6 + 4 2
( + 21 + 15 2)e2
cd () = e1
2 + (3 2 + 4) 21 2 30
(2 + ( 2 + 1) 21 15 2)e3
+
2 + (3 2 + 4) 21 2 30
(2 + ( 2 + 2) + 30 + 21 2)e4
2 + (3 2 + 4) 21 2 30
176
ce () = e1 +
2
32 + (16 2 24) 96 2 + 136
+
2
32 + (16 2 24) 96 2 + 136
2
(3 + (18 2 28) + 340 240 2)e4
.
+
32 + (16 2 24) 96 2 + 136
O
I
I
,
O
where O denotes 3 3 zero matrix and I the 3 3 identity matrix. As example we take the same conic as in Ex. 2.5, but with
another parametrisation. We start with the Pl
ucker coordinates
A = aR = (0, 1, 0, 0, 0, 1)T R, B = bR = (1, 1, 1, 0, 0, 0)T R, and C = cR =
(0, 0, 1, 0, 1, 0)T R. Obviously, A, B, C M24 . The corresponding null
vectors are given by
a = e2 + e 6 ,
b = e1 + e2 + e3 ,
c = e 3 + e5 .
The chain dened by these three null vectors is computed with Eq.
(2.8)
cabc () =
2e1
(1 + )e2
(1 + )e3
( 1)e5
( 1)e6
+
+
+
.
1 + 2
1 + 2
1 + 2
1 + 2
1 + 2
177
s1
s2
B
D
C
The chain cd containing d that has contact with cabc is given by:
24(23+28)e1
28(12+5)e2
(704+3362 +371)e3
+
607+1040+3362 607+1040+3362
607+1040+3362
2
336e4
(1376+336 +827)e5
28(12+29)e6
+
+
.
607+1040+3362
607+1040+3362
607+1040+3362
cd () =
178
dual angle with respect to a xed line. In general, two such congruences with constant dual angles with respect to a xed line may have
zero, one, two, or innitely many oriented lines in common. The adequate ane model is the Cliord algebra C(0,3,0) (D) given by M = D3
0 0
ucker coordinates
1 0. As example we give the Pl
0 1
of four lines A = aR = (0, 1, 0, 0, 0, 1)T R, B = bR = ( 12 , 23 , 23 , 49 , 19 , 19 )T R,
C = cR = (0, 0, 1, 0, 1, 0)T R, D = dR = ( 35 , 0, 45 , 85 , 2, 65 )T R. The corres
1
and Q = 0
0
1 4
2 1
2 1
b = ( + )e1 +( + )e2 +( + )e3 ,
3 9
3 9
3 9
8
6
3
4
1
d = ( )e1 +( +2)e2 +( + )e3 .
5
5
5
5
5
We compute the chains cabc and cacd for a dual cross ratio with Eq.
(2.8).
(+1)(2 ++)
+
e1
2 +1+
(2 +1+)2
+1
(2+3 +2 1)
+
+
e2
2 +1+
(2 +1+)2
(+1) 2+3 1
+
e3 ,
2 +1+
(2 +1+)2
3(1)
6+3+42 122 483 +404 +4
e1
cacd (, ) =
1+52
(1+52 )2
(109+5152 +503 1)
1
+
+
e2
(1+52 )
(1+52 )2
(4+1)
(1)(303 312 +791)
+
+
e3 .
(1+52 )
(1+52 )2
cabc (, ) =
179
B
C
s2
s1
D
chain cd that is in contact with cabc can be computed with a dual cross
ratio + , since the scalar domain of the algebra is the commutative
ring of dual numbers. Both chains correspond to congruences of lines
with constant dual angle to an uniquely determined oriented line.
Hence, the two congruences of lines dene two cylinders that intersect
each other in a single point if the congruences of lines are in contact.
Furthermore both congruences of lines have a whole bundle of parallel
lines in common. This is obvious because we can translate the line
that belongs to both congruences and that passes through the point
of contact of the two cylinders along the axis of one of the cylinders
and it is still contained in both congruences. A parametrisation of
this congruence can be obtained with Eq. (2.20):
( + 3)
4 + 16 + 40 + 6 + 2 + 3
cd (, ) =
e1
5 + 3 + 2
(5 + 3 + 2 )2
282 5 + 70 + 50 + 2 + 73 + 4
+
+
e2
5 + 3 + 2
(5 + 3 + 2 )2
2
+ 3 + 4 2 + 52 + 3 + 12 + 30 + 3
+
+
e3 .
5 + 3 + 2
(5 + 3 + 2 )2
180
1 u2
u
1
u + v u + 1 ,
f1 (u, v) = 2
u +u+1
(u + 1)u
u
1
(1 u2 , u, u)T .
u2 + u + 1
4
u 7u3 28u2 68u40
1
f2 (u, v) = 2
u3 u2 20u50
(u +3u+5)2
u3 +6u2 +26u+30
15+14u+u3 +6u2
,
+v
5uu3 3u2
3
4
2
20+27u+6u +u +18u
u2
1
(1 u2 4u, u 10, 13u + 18)T .
+ 3u + 5
182
n 1
183
184
e23 i,
e31 j,
e12 k,
e1234 ,
e14 i,
e24 j,
e34 k.
(3.3)
185
(3.4)
SE(3) A = (a0 , a1 , a2 , a3 , c0 , c1 , c2 , c3 ) .
T
186
sin
cos
a
x
.
+
b
y
x1
,
x0
y=
x2
,
x0
z=
x3
,
x0
if x0 = 0.
T
A1 : x0 = x3 = 0,
187
188
gg = 1.
In
the projective representation the condition gg = 1 changes to
3
2
i=0 ai = 0 and further we have a0 c0 a1 c1 + a2 c2 a3 c3 = 0. The
action of a Spin group element can be written in terms of the sandwich
operator
gpg = (a20 + a21 + a22 + a23 )x0 e123
+ 2(a0 c1 a1 c0 a2 c3 a3 c2 )x0 + (a20 + a21 a22 a23 )x1
+ 2(a2 a1 a0 a3 )x2 + 2(a0 a2 + a3 a1 )x3 e234
+ 2(a0 c2 + a1 c3 a2 c0 a3 c1 )x0 + 2(a0 a3 + a2 a1 )x1
+ (a20 a21 + a22 a23 )x2 + 2(a3 a2 a0 a1 )x3 e134
+ (2a2 c1 + 2a1 c2 2a0 c3 2a3 c0 )x0 + 2(a3 a1 a0 a2 )x1
+ 2(a3 a2 + a0 a1 )x2 + (a20 a21 a22 + a23 )x3 e124 .
189
with
1
l
M=
m
n
0
a20 + a21 a22 a23
2(a0 a3 + a2 a1 )
2(a3 a1 a0 a2 )
0
2(a2 a1 a0 a3 )
a20 a21 + a22 a23
2(a3 a2 + a0 a1 )
0
2(a0 a2 + a3 a1 )
2(a3 a2 a0 a1 )
a20 a21 a22 + a23
and
= a20 + a21 + a22 + a23 ,
l = 2(a0 c1 a3 c2 a2 c3 a1 c0 ),
m = 2(a0 c2 + a1 c3 a2 c0 a3 c1 ),
n = 2(a2 c1 + a1 c2 a0 c3 a3 c0 ).
When
using homogeneous coordinates the norming condition becomes
3
2
i=0 ai = 0. Furthermore, the parameters (a0 , . . . , a3 , c0 , . . . , c3 ) have
to satisfy the condition (3.7). We get
ge0 = a0 e0 + a3 e12 + a2 e13 + c1 e14 + a1 e23 + c2 e24 + c3 e34 + c0 e1234 ,
ge23 = a1 e0 a2 e12 + a3 e13 c0 e14 + a0 e23 c3 e24 + c2 e34 + c1 e1234 ,
ge13 = a2 e0 + a1 e12 + a0 e13 c3 e14 a3 e23 + c0 e24 + c1 e34 c2 e1234 ,
ge12 = a3 e0 + a0 e12 a1 e13 c2 e14 + a2 e23 + c1 e24 c0 e34 + c3 e1234 ,
ge1234 = a1 e14 + a2 e24 a3 e34 + a0 e1234 ,
ge14 = a0 e14 a3 e24 a2 e34 + a1 e1234 ,
190
a0
a 1
a 2
+ a 3
G =
c0
c1
c
2
c3
a1
a0
a3
a2
c1
c0
c3
c2
a2
a3
a0
a1
c2
c3
c0
c1
a3
a2
a1
a0
c3
c2
c1
c0
0
0
0
0
a0
a1
a2
a3
0
0
0
0
a1
a0
a3
a2
0
0
0
0
a2
a3
a0
a1
0
0
0
.
a3
a2
a
1
a0
In order to write [G+ ] in the usual form, see Eq. (1.9) the isomorphism
(3.3) has to be used.
The condition that this element is a Spin group element reads now
gg = a20 + a21 = 1.
In this case we have just one equation in the scalar part. If we change
to homogeneous coordinates, the condition gets a20 + a21 = 0. Again we
have the equation of a pair of complex conjugate planes intersecting
in a real line. This results in the same image space as in the previous
section, see Eq. (3.5).
191
x0
a20 + a21
1
x1 =
2(a1 c0 + a0 c1 )
a20 + a21
x2
2(a1 c1 c0 a0 )
0
a20 a21
2a0 a1
0
x0
2a0 a1 x1 .
a20 a21
x2
(3.8)
In the rst line of Eq. (3.8) we see that the line at innity is xed under
these collineations. It is clear that all planar Euclidean displacements
are obtained in this way, because the Spin group is a double cover
of SE(2). The factor is added articially to guarantee that we just
look at matrices belonging to the Spin group. These two examples
show that the Spin group can be used to get the matrix group of the
direct isometries of a Cayley-Klein geometry, if it is representable as
homogeneous Cliord algebra model.
Collineations in the Image Space The Spin group elements induce
in the image space a group of transformations that x the pair of
complex conjugate planes, their intersection line and the two complex
conjugate points on this line. To get this group we have to look at the
matrix representation of the Spin group. Let g = a0 e0 + a1 e12 + c0 e23 +
c1 e13 Spin(2,0,1) . Calculating the matrix representation as described
in section 1.3.4 results in
a0 a1 0
0
a1 a 0
PGL(P2 (R), F) =
c1 c0 a0
c0 c1 a1
0
0
a0 , a1 , c 0 , c 1 R .
a1
a0
192
CK space
absolute Figure F
elliptic
Q13,0
hyperbolic
Q13,1
Euclidean
Q11,0 A1 Q02,0
pseudo-Euclidean
Q11,0 A1 Q02,1
quasi-elliptic
Q12,0 A0 Q1
1,0
1
0
quasi-hyperbolic
Q2,1 A Q1
1,0
totally isotropic Q11,0 A1 Q01,0 A0 Q1
1,0
Cliord alg.
C(3,0,0) , C(0,3,0)
C(2,1,0) , C(1,2,0)
C(2,0,1) , C(0,2,1)
C(1,1,1)
C(2,0,1) , C(0,2,1)
C(1,1,1)
C(1,0,2) , C(0,1,2)
193
gures F and the corresponding Cliord algebras for possible twodimensional Cayley-Klein spaces are given. Note that the points of the
Cayley-Klein geometry are always written in the grade-2 subspace of
the Cliord algebra except the quasi-elliptic and the quasi-hyperbolic
case. In these two cases we use the Poincare isomorphism to describe
these Cayley-Klein geometries as dual partners of the Euclidean
respectively, the pseudo-Euclidean Cayley-Klein geometry. Therefore,
points of the quasi-elliptic (quasi-hyperbolic) Cayley-Klein geometry
are represented as grade-1 elements in the algebra corresponding to the
Euclidean (pseudo-Euclidean) Cayley-Klein space. Table 3.2 shows
the kinematic image spaces of the seven planar Cayley-Klein spaces.
Note that every image space is again a Cayley-Klein space.
Table 3.2: Kinematic image spaces of two-dimensional Cayley-Klein geometries presented as three-dimensional Cayley-Klein spaces
with their absolute gures F
pre-image
CK image space
absolute Figure F
elliptic
elliptic
Q24,0
hyperbolic
hyperbolic idx. 1
Q24,2
Euclidean
quasi-elliptic
Q22,0 A1 Q02,0
pseudo-Euclidean quasi-hyperb. idx. 0
Q22,1 A1 Q02,1
quasi-elliptic
quasi-elliptic
Q22,0 A1 Q02,0
quasi-hyperbolic quasi-hyperb. idx. 0
Q22,1 A1 Q02,1
totally isotropic
totally isotropic
Q21,0 A2 Q12,0 A0 Q1
1,0
194
195
CK-space
elliptic
hyperbolic idx. 0
hyperbolic idx. 1
Euclidean
pseudo-Euclidean
quasi-elliptic
quasi-hyperb. idx. 0
double isotropic
C(p,q,r)
C(4,0,0) , C(0,4,0)
C(3,1,0) , C(1,3,0)
C(2,2,0)
C(3,0,1) , C(0,3,1)
C(2,1,1) , C(1,2,1)
C(2,0,2) , C(0,2,2)
C(1,1,2)
C(1,0,3) , C(0,1,3)
Q62
3
2
2
i=0 ai +
i=0 ci
3
3
2
2
i=0 ai
i=0 ci
a20 a21 a22 +a23 +c20 c21 c22 +c23
3
2
i=0 ai
a20 a21 a22 + a23
a20 + a23
a20 a23
a20
3
(3.9)
+ 2(c0 a0 a1 c1 c3 a3 + c2 a2 )e1234
+ 2(c7 c1 + c4 a0 c5 c2 c0 a5 + a6 a2 a1 a4 + c6 c3 a7 a3 )e1235
+ 2(c5 a1 + a6 c1 a5 a3 c2 a4 c7 a2 + c6 a0 + c0 a7 c4 c3 )e1245
+ 2(c4 c2 c6 a1 a5 a2 c3 a4 c0 a6 + c5 a0 + c7 a3 + a7 c1 )e1345
196
(3.10)
+ 2(c2 a2 c3 a3 + c0 a0 a1 c1 )e1234
+ 2(c4 a0 a7 a3 c7 c1 + c0 a5 a1 a4 + a6 a2 c6 c3 + c5 c2 )e1235
+ 2(c7 a2 + a6 c1 c5 a1 a5 a3 + c6 a0 c2 a4 + c4 c3 c0 a7 )e1245
+ 2(c0 a6 c3 a4 c4 c2 + c6 a1 c7 a3 + a7 c1 + c5 a0 a5 a2 )e1345
+ 2(a7 c2 + c4 c1 a5 a1 + c5 a3 c6 a2 c3 a6 c0 a4 + c7 a0 )e2345 .
9
2
i=1
with
N1 : c22 + c23 + a22 + a23 + a20 + c21 + a21 + c20 = 0,
Q1 : c0 a0 a1 c1 c3 a3 + c2 a2 = 0,
Q2 : c7 c1 + c4 a0 c5 c2 c0 a5 + a6 a2 a1 a4 + c6 c3 a7 a3 = 0,
Q3 : c5 a1 + a6 c1 a5 a3 c2 a4 c7 a2 + c6 a0 + c0 a7 c4 c3 = 0,
Q4 : c4 c2 c6 a1 a5 a2 c3 a4 c0 a6 + c5 a0 + c7 a3 + a7 c1 = 0,
Q5 : c6 a2 + c0 a4 c4 c1 c5 a3 + a7 c2 c3 a6 + c7 a0 a5 a1 = 0,
Q6 : c4 a0 a7 a3 c7 c1 + c0 a5 a1 a4 + a6 a2 c6 c3 + c5 c2 = 0,
Q7 : c7 a2 + a6 c1 c5 a1 a5 a3 + c6 a0 c2 a4 + c4 c3 c0 a7 = 0,
Q8 : c0 a6 c3 a4 c4 c2 + c6 a1 c7 a3 + a7 c1 + c5 a0 a5 a2 = 0,
Q9 : a7 c2 + c4 c1 a5 a1 + c5 a3 c6 a2 c3 a6 c0 a4 + c7 a0 = 0.
197
198
Now we use these four quadratic equations and formulate the Spin
group condition gg = 1,
gg = (a21 + c22 + c23 + c20 + c21 + a20 + a23 + a22 )e0
+ 2(a1 c1 c3 a3 + c0 a0 + c2 a2 )e1234
+ 2(c3 a6 + a7 c2 + c7 a0 a5 a1 )e2345
+ 2(c4 a0 a7 a3 + a6 a2 a1 a4 )e1235
+ 2(c6 a0 a5 a3 + a6 c1 c2 a4 )e1245
+ 2(a7 c1 a5 a2 c3 a4 + c5 a0 )e1345
and get ve more quadratic equations that have to vanish. All in all we
have the corresponding projective variety described by nine quadratic
equations and one exceptional quadric with equation
N1 : a20 + a21 + a22 + a23 + c20 + c21 + c22 + c23 = 0,
R1 : c6 a2 + c0 a4 c4 c1 c5 a3 = 0,
R2 : c7 c1 c6 c3 + c5 c2 + c0 a5 = 0,
199
R3 : c7 a2 + c0 a7 + c5 a1 c4 c3 = 0,
R4 : c0 a6 + c6 a1 c4 c2 c7 a3 = 0,
R5 : a1 c1 c3 a3 + c0 a0 + c2 a2 = 0,
R6 : c3 a6 + a7 c2 + c7 a0 a5 a1 = 0,
R7 : c4 a0 a7 a3 + a6 a2 a1 a4 = 0,
R8 : c6 a0 a5 a3 + a6 c1 c2 a4 = 0,
R9 : a7 c1 a5 a2 c3 a4 + c5 a0 = 0.
9
2
i=1
R2 =
1
(Q6 Q2 ),
2
R 5 = Q1 ,
R8 =
1
(Q3 + Q7 ),
2
1
(Q3 Q7 ),
2
1
R6 = (Q9 + Q5 ),
2
1
R9 = (Q8 + Q4 ).
2
R3 =
Conclusion
In chapter 1 we presented geometric models and representation that are
used today. Afterwards, Cliord algebras, the homogeneous Cliord
algebra model, and the conformal model were introduced. Moreover,
we presented a method to transfer general projective transformations
acting in three-dimensional projective space to elements of the homogeneous Cliord algebra model C(3,3,0) . All entities known from
line geometry occur naturally in this model and can be transformed
projectively by the application of the sandwich operator. It was
shown that the sandwich action of non-null vectors corresponds to
regular null polarities, i.e., correlations that are involutions as the
basic elements building up the group of regular projective transformations. The sandwich action of null vectors corresponds to singular null
polarities and we proved that every regular projective transformation
can be expressed as the product of six null polarities at the most.
Furthermore, we present a Cliord algebra model that serves as a
generalization of the conformal geometric algebra. A generalization
of inversions with respect to conics, quadrics, and even hyperquadrics
in any dimension is possible with the use of the sandwich operator.
Hyperquadrics in principal position are simply represented as grade-1
elements. Classical representations of groups are embedded in this
algebra naturally.
The basics of chain geometry were introduced in chapter 2. This theory was applied to Cliord algebras. We showed how the cross ratio
can be used to parametrize chains. Moreover, it was proven that the
connected components of the Pin- and Spin groups dene subspaces
of chain spaces. Therefore, every element contained in a chain that is
dened by three elements of the same connected component of the
Pin- or the Spin group is contained in the same connected component
of the Pin- or the Spin group. Of special interest was the question for
the cross ratio of dual unit quaternions that was answered in detail.
D. Klawitter, Clifford Algebras, DOI 10.1007/978-3-658-07618-4,
Springer Fachmedien Wiesbaden 2015
202
Conclusion
Index
absolute gure, 182
admissible pairs, 103
algebra, 102
Cliord, 26
degenerate, 26
geometric, 26
Benz planes, 116
biarc, 165
Blaschke cylinder, 117
Cayley-Klein
geometry, 182
space, 182, 193
center, 102
chain, 106
geometry, 106
space, 107
Cliord algebra
conjugation, 29
even part, 28
main involution, 30
odd part, 28
universal, 28
Cliord group, 33
contact
relation, 121
space, 121
cross ratio, 109
distance
relation, 101
space, 101
dual
number, 8
orthogonal matrix, 13
quaternion, 10
Rodriguez formula, 14
sphere, 136
vector, 9
general linear group, 102
geometric inner product null
space, 44
geometric outer product null
space, 44
Grassmann
algebra, 39
coordinates, 39, 40
group
Spin group, 35
horosphere, 43
ideal, 199
inner product null space, 32
Jordan-closed, 119
Jordan-system, 118, 119
204
Index
subalgebra, 118
versor, 28
trace, 197
transitive
2--transitive, 104
3--transitive, 104
List of Symbols
(P, )
(l, m)
(p, q, r)
V
C+
(p,q,r)
C(p,q,r)
(J )
(, )
(C(p,q,r) )
(R)
,
[]m
[A+ ] , [B ]
i, j, k
N (R)
C(R)
F
G
206
J
K
M
N1
O
R
V
A
C
D
Dn
En
H
Hd
NI(A)
NIG (A)
NO(A)
NOG (A)
P1 (J )
P1 (R)
Pn (R)
Ud
A
a
C(L, R)
g
I
Q
Q
GL(R, 2)
I
PGL(R, 2)
Q
SE(3)
List of Symbols
Jordan-system, p. 119
eld, p. 102
module, p. 77
null vectors of a Cliord algebra, p. 128
dual orthogonal matrix, p. 13
set of units of a ring R , p. 33
projective variety, p. 199
split-complex numbers, p. 117
complex numbers, p. 117
ring of dual numbers, p. 8
n-dimensional module over D, p. 9
n-dimensional Euclidean space, p. 7
skew eld of quaternions, p. 6
dual unit quaternions, p. 10
inner product null space of a k-blade A k V , p.
32
geometric
inner product null space of a k-blade
k
A
V , p. 44
outer product null space of a k-blade A k V , p.
31
geometric
outer product null space of a k-blade
k
A
V , p. 44
projective line over the Jordan-system J , p. 120
projective line over the ring R, p. 103
n-dimensional real projective space, p. 5
group of dual unit quaternions, p. 10
blade in C(p,q,r) , p. 31
grade-1 elements resp. vectors in C(p,q,r) , p. 30
set of chains in P1 (R), p. 106
element of the Cliord group, p. 33
pseudoscalar, p. 46
projective quadric, p. 115
non-double points of Q, p. 116
general linear group over a ring R, p. 102
identity matrix, p. 13
projective linear group over P1 (R), p. 104
quadratic form, p. 26
group of spatial displacements, p. 21
List of Symbols
SO(3)
T(2)
(, )
Pin(p,q,r)
(K, R, J )
(L, R)
Spin(p,q,r)
An1
cabc
cr(a, b, c, d)
e + , e
ei
J
Ln+1
1
N (a)
N26
O(p, q, r)
PI
Qnr,q
S26
2
SD
tr()
V
V3
v = v +
v
Vkd
k
207
group of special orthogonal 3 3 matrices, p. 6
additive group of planar translations, p. 92
bilinear form corresponding to a quadric, p. 16
dual angle, p. 14
parallel relation, p. 102
Pin group, p. 34
chain geometry over (K, R, J ), p. 120
chain geometry over the L-algebra R, p. 106
Spin group, p. 35
map of Blaschke and Gr
unwald, p. 186
Studys kinematic mapping, p. 21
absolute n1 -subspace, p. 182
chain dened by a, b, c, p. 129
cross ratio of a, b, c and d, p. 1
generators of CGA, p. 43
generators of a Cliord algebra, p. 26
Poincare duality, p. 41
Lies quadric (ndimensional), p. 74
norm of a, p. 33
singular quadric in P7 (R), p. 185
orthogonal group, p. 34
image of identity e SE(3) under Study mapping,
p. 23
quadric of rank r, index q and dimension n, p. 182
Studys quadric, p. 21
Studys sphere, p. 25
trace of an algebra element, p. 197
vector space, p. 26
the exceptional generator, p. 21
small bold Latin letters denote dual vectors, with
real part v and dual part v, p. 9
algebraic variety of dimension d and degree k, p.
17
grade-k subspace of the exterior algebra corresponding to V , p. 28
Acknowledgment
Foremost, I would like to express my deepest thanks to my two supervisors, Professor Dr. Gunter Weiss and Prof. Dr. Marco Hamann.
Their patience, encouragement, and immense knowledge were key motivations throughout my PhD. They persuasively conveyed an interest
in my work, and I am grateful to my advisers.
Prof. Dr. Gunter Weiss has been my supervisor through three and a
half years of geometric research. I am truly thankful for his steadfast
integrity, and seless dedication to both my personal and academic
development. I cannot think of a better supervisor to have. Prof. Dr.
Marco Hamann is a mentor, from whom I have learnt the vital skill
of disciplined critical thinking. His forensic scrutiny of my work has
been invaluable. He has always found the time to propose consistently
excellent improvements. I owe a great debt of gratitude to Prof. Dr.
Gunter Weiss and Prof. Dr. Marco Hamann.
I would like to thank PD Dr. Boris Odehnal for oering thorough
and excellent feedback on an earlier version of this thesis. In addition,
a thank you to Prof. Dr. Brehm. He gave feedback on research,
and suggested general improvements to my models. Furthermore, I
thank Prof. Dr. Andrea Blunck and Prof. Dr. Hans Havlicek. Both
answered me questions on chain geometry during email correspondences. A special mention for Prof. Dr. Daniel Lordick, for our fruitful
discussions.
I thank my parents Manuela Klawitter and Burghard Klawitter for
the encouragement during the last years.
And last but not least, I thank my girlfriend Sarah Kunte. Her patience and her emotional support did always motivate me to go further.
210
Acknowledgment
Bibliography
[1] Ablamowicz, R., and Fauser, B. Mathematics of CLIFFORD
- a Maple Package for Cliord and Grassmann Algebras. Advances
in Applied Cliord Algebras 15, 2 (2005), 157181.
[2] Ablamowicz, R., and Fauser, B. CLIFFORD/Bigebra, A
c
Maple Package for Cliord (Co)Algebra Computations: 19962011, RA&BF. Available at http://www.math.tntech.edu/rafal/
(2011).
[3] Ablamowicz, R., and Sobczyk, G. Lectures on Cliord
(geometric) algebras and applications. Birkh
auser, Boston, 2004.
212
Bibliography
ur
[10] Blunck, A. Chain spaces via Cliord algebras. Monatshefte f
Mathematik, 123 (1997), 93107.
[11] Blunck, A., and Herzer, A.
Einf
uhrung. Shaker, Aachen, 2005.
Kettengeometrien: Eine
Bibliography
213
214
Bibliography
[35] Hestenes, D. Old Wine in new Bottles - A new algebraic framework for computational geometry. E. Bayro-Corrochano & G.
Sobczyk (eds), Advances in Geometric Algebra with Applications
in Science and Engineering (2001), 114.
[36] Hestenes, D., and Sobczyk, G. Cliord Algebra to Geometric
Calculus: A unied Language for Mathematics and Physics. D.
Reidel, Dordrecht and Boston, 1984.
[37] Hotje, H. Zur Einbettung von Kettengeometrien in projektive
R
aume. Math. Zeitschr., 151 (1976), 517.
[38] Husty, M., Karger, A., Sachs, H., and Steinhilper, W.
Kinematik und Robotik. Springer, Berlin [u.a.], 1997.
[39] Karzel, H. Kinematische Algebren und ihre geometrischen
Ableitungen. Abhandlungen aus dem Mathematischen Seminar
der Universit
at Hamburg, 1 (1974), 158171.
[40] Klawitter, D. A Cliord algebraic Approach to Line Geometry.
Advances in Applied Cliord Algebras (2014).
[41] Klawitter, D., and Hagemann, M. Kinematic mappings
for Cayley-Klein geometries via Cliord algebras. Beitr
age zur
Algebra und Geometrie / Contributions to Algebra and Geometry
54, 2 (2013), 737761.
[42] Klein, F. Ueber Liniengeometrie und metrische Geometrie.
Mathematische Annalen, 5 (1872), 106126.
ber h
ohere Geo[43] Klein, F., and Blaschke, W. Vorlesungen u
metrie. Springer, Berlin, 1926.
[44] Kowol, G. Projektive Geometrie und Cayley-Klein Geometrien
der Ebene. Birkh
auser, Basel, 2009.
[45] Li, H., Hestenes, D., and Rockwood, A. A Universal Model
for Conformal Geometries of Euclidean, Spherical and DoubleHyperbolic Spaces. Geometric Computing with Cliord Algebras
(2001), 7799.
Bibliography
215
[46] Li, H., Hestenes, D., and Rockwood, A. Generalized Homogeneous Coordinates for Computational Geometry. Geometric
Computing with Cliord Algebras (2001), 2752.
[47] Li, H., and Zhang, L. Line geometry in terms of the null
geometric algebra over R(3, 3), and application to the inverse
singularity analysis of generalized stewart platforms. In Guide
to Geometric Algebra in Practice, L. Dorst and J. Lasenby, Eds.
Springer London, 2011, pp. 253272.
[48] Macdonald, A. A Survey of Geometric Algebra and Geometric
Calculus, 2008.
[49] McCarthy, J., and Ahlers, S. Dimensional synthesis of
robots using a double quaternion.
[50] McCarthy, J. M. An Introduction to Theoretical Kinematics.
MIT Press, Cambridge and Mass, 1990.
ller, E. Vorlesungen u
ber Darstellende Geometrie: Kon[51] Mu
struktive Behandlung der Regel
achen bearb. von Josef Leopold
Krames, vol. 3. Deuticke, Leipzig and Wien, 1931.
[52] Onishchik, A., and Sulanke, R. Projektive und CayleyKleinsche Geometrie. Springer, 2008.
[53] Perwass, C. B. U. Geometric Algebra with Applications in
Engineering. Springer, Berlin and Heidelberg, 2009.
[54] Pfurner, M. Analysis of spatial serial manipulators using
kinematic mapping. PhD thesis, University of Innsbruck, Austria,
Innsbruck.
[55] Porteous, I. R. Cliord Algebras and the Classical Groups.
Cambridge University Press, Cambridge, 2009.
[56] Pottmann, H., and Wallner, J. Computational Line Geometry. Springer, Berlin and New York, 2001.
[57] Ravani, B., and Sprott, K. Kinematic Generation of Ruled
Surfaces. Advances in Computational Mathematics, 17 (2002),
115133.
216
Bibliography