PGA(3,0,1) for Euclidean 3D Geometry

• PGA(3,0,1) is a 16-dimensional Clifford algebra that encodes Euclidean 3D geometry using homogeneous, projective coordinates.
• It employs both plane-based and point-based dual representations to unify projective and metric structures for points, lines, and planes.
• The algebra supports efficient operations—such as reflections, rotations, and translations—via sandwich operations for rigid-body kinematics.

@@@@1@@@@ (PGA) with signature $(3,0,1)$, denoted PGA$(3,0,1)$ or $P(\mathbb{R}^*_{3,0,1})$, is a 16-dimensional Clifford algebra constructed over a four-dimensional real vector space equipped with a symmetric bilinear form of signature $(3,0,1)$. This algebra is optimized for encoding the full range of Euclidean 3D geometry—points, lines, planes, distances, angles, rigid-body kinematics, and mechanics—using homogeneous (projective) coordinates. The algebra’s plane-based and point-based dual forms enable a unified, coordinate-free treatment of both projective geometry and Euclidean metric structure, making it a canonical model for applications requiring structure-preserving representations, direct support for dualities, and sandwich (versor) operations for the group of Euclidean isometries (Gunn, 2022, Gunn, 2014, Gunn, 2020, Brehmer et al., 2023).

1. Algebraic Structure and Metric Signature

Let $V \cong \mathbb{R}^4$ be a vector space with orthogonal basis $\{e_0, e_1, e_2, e_3\}$. The metric is specified as:

• $e_1^2 = e_2^2 = e_3^2 = +1$
• $e_0^2 = 0$ (degenerate)
• $e_i e_j + e_j e_i = 0$ for $i \neq j$.

Under this construction, $e_0$ encodes the ideal plane (the hyperplane at infinity in projective 3-space), providing a direct projectivization of Euclidean geometry. The total algebra $\mathcal{G} = \text{PGA}(3,0,1)$ has $2^4 = 16$ dimensions, naturally graded from scalars (grade 0) to the pseudoscalar (grade 4). The defining relations yield robust representations for all "flat" geometric primitives and their mutual incidences and metric relations (Gunn, 2014, Sobczyk, 2018).

The volume element (pseudoscalar) is $I = e_0 e_1 e_2 e_3$, satisfying $I^2 = 0$ due to the null square of $e_0$ (Gunn, 2022).

2. Representation of Geometric Entities

Projective geometric algebra natively encodes Euclidean geometric primitives of $\mathbb{R}^3$ as multivectors with specific grades:

Geometric entity	Plane-based: $\mathcal{G}$ (grade)	Point-based: $\mathcal{G}^*$ (grade)
Point	$p = p_0 e_0 + p_1 e_1 + p_2 e_2 + p_3 e_3\,\ (3)$	$P = P^0 e^0 + P^1 e^1 + P^2 e^2 + P^3 e^3\ (1)$
Line	$\ell = p \wedge q\,\ (2)$	$m = P \wedge Q\,\ (2)$
Plane	$a\,\ (1)$	$\pi = P \wedge Q \wedge R\ (3)$

• Points in $\mathcal{G}$ are intersections (meets) of three planes: $P = a \wedge b \wedge c$.
• Lines are intersections of two planes: $\ell = a \wedge b$.
• Planes are 1-vectors: $a = a_0 e_0 + a_1 e_1 + a_2 e_2 + a_3 e_3$.
• The dual algebra $\mathcal{G}^*$ ("point-based") reverses this identification by exchanging roles of points and planes (Gunn, 2022, Brehmer et al., 2023).

Canonical representatives utilize homogeneous or projective coordinates. For instance, the Euclidean point $(x, y, z) \in \mathbb{R}^3$ is represented (in the simplest convention) as $P = e_0 + x e_1 + y e_2 + z e_3$ (Gunn, 2020, Sobczyk, 2018).

3. Products, Duality, and the Role of Degeneracy

Exterior and Regressive Product

• The exterior (wedge) product $\wedge$ models the meet operation: intersection of geometric entities.
• The regressive (join) product $\vee$ computes the span. In degenerate algebras, join is implemented via duality:
  • $X \vee Y = J^{-1}(J(X) \wedge J(Y))$
  • $J : \mathcal{G} \to \mathcal{G}^*$ is the coordinate-free dual coordinate map, pairing $k$-blades with $(4-k)$-blades corresponding to the same geometric subspace (Gunn, 2022, Gunn, 2014).

Duality Maps

• The e-duality map $J$ interchanges the plane- and point-based models (double algebra duality), providing a coordinate-free, functorial identification essential for geometric duality in projective geometry.
• The Hodge map $H$ is a grade-reversing operator defined via combinatorial complement, but it is not invariant under change of basis and depends on the metric structure (Gunn, 2022).
• In degenerate metrics ($I^2=0$), the classical metric (pseudoscalar) polarity fails; all dualities must be defined combinatorially or via basis index complements.

Primitives and the Plücker Relation

• Grade-2 elements (bivectors) $L \in \mathcal{G}^2$ encode lines, with their six Plücker coordinates $L_{ij}$ arranged as $$L = L_{01} e_0 \wedge e_1 + \cdots + L_{12} e_1 \wedge e_2$$. The geometric Plücker condition $L \wedge L = 0$ characterizes lines as simple bivectors and defines the Grassmannian $G(2,4)$ as a quadric in $\mathbb{P}^5$ (Sobczyk, 2018).

4. Metric Structure and Euclidean Motions

Despite the degenerate metric, PGA$(3,0,1)$ faithfully encodes:

• Distances: The norm $\|P \vee Q\|$ of the join of two normalized points yields the Euclidean distance $\|P - Q\|$.
• Angles: For normalized planes $a, b$, $\cos\phi = a \cdot b$, and for lines, $\cos\theta = (L_1 \cdot L_2) / (\|L_1\|\|L_2\|)$ (Gunn, 2014, Gunn, 2020).
• Reflections: Reflections in a (normalized) plane $a$ are performed as $X' = a X a$.
• Motors: Rotations and translations (the generators of $SE(3)$) are constructed via sandwiched exponentials of bivectors:
  • Rotation: $R(\theta) = e^{-\frac{\theta}{2} B}$, $X' = R X R^{-1}$;
  • Translation: $T(t) = 1 - \frac{t}{2} m e_0$, $X' = T X T^{-1}$;
  • Screw motions: $S = e^{\frac{t}{2}(B + p B I)}$ (Gunn, 2014, Gunn, 2020, Brehmer et al., 2023).

The even subalgebra isomorphic to dual quaternions provides a direct computational parallel to biquaternion approaches in kinematics.

5. Implementation Aspects and Duality-Neutral Software

A robust software strategy for PGA$(3,0,1)$ must accommodate both duality frameworks without bias:

• Each multivector is equipped with a Boolean flag isDual. Dualization (operator $\mathfrak{D}$) flips this bit and applies the combinatorial complement to coefficients.
• The same dual operator realizes $J$ and $J^{-1}$ in e-duality mode, or $H$ and $H^{-1}$ in Hodge mode.
• The regressive product is written uniformly as $$X \vee Y = \mathfrak{D}(\mathfrak{D}(X) \wedge \mathfrak{D}(Y))$$.
• Algebraic operations interrogate isDual to select the multiplication context, promoting or signaling type mismatches as needed (Gunn, 2022).

This approach ensures transparent handling of degenerate pseudoscalars and supports a duality-neutral computational environment aligned with the natural geometric structure of projective space.

6. Broader Applications and Comparative Models

PGA$(3,0,1)$ supports the entirety of flat Euclidean geometry, including:

• Robust, coordinate-free manipulations of points, lines, and planes
• Seamless integration of meets/joins, handling both incident and parallel cases with the same formulas
• Kinematics and rigid body mechanics: The algebra encodes the Lie algebra $\mathfrak{so}(3,0,1)$, naturally supports screw theory, and embeds the Newton–Euler equations for rigid bodies
• Dual quaternion algebra and automatic differentiation: The even subalgebra covers dual quaternions, and the scalar-pseudoscalar sector provides dual numbers for automatic differentiation (Gunn, 2014, Gunn, 2020).

In comparison to conformal geometric algebra (CGA, signature $(4,1,0)$), PGA$(3,0,1)$ realizes all required structure-preserving operations for flat geometries while remaining of minimal dimension. Spheres and circles are not represented as blades in PGA, but for purely flat entities, the feature sets are isomorphic (Gunn, 2014). In machine learning, the Geometric Algebra Transformer (GATr) leverages PGA's 16-dimensional representation to encode geometric objects and operators, achieving $E(3)$-equivariance in deep learning tasks (Brehmer et al., 2023).

7. Geometric and Algebraic Significance

PGA$(3,0,1)$ serves as a canonical model for Euclidean 3-space in both theoretical and applied domains. Its algebraic foundation directly implements projective geometry’s principle that all subspaces can be uniformly represented and manipulated. The algebra’s handling of the degenerate pseudoscalar circumvents pitfalls associated with metric-polarity methods, instead providing coordinate-free, basis-invariant duality through the $J$ map. This feature underpins rigorous computational geometry, automatic symbolic manipulation of joins and meets, and the direct embedding of classical screw theory and rigid-body dynamics (Gunn, 2022, Gunn, 2014, Gunn, 2020).

By integrating geometric, metric, and transformation structure within a single algebraic model, PGA$(3,0,1)$ offers a minimal, structure-preserving, and computationally efficient framework for both classical geometric tasks and emerging applications in robotics, computer graphics, and machine learning (Brehmer et al., 2023).