# Geometric calculus

To be distinguished from vector calculus.

In mathematics, geometric calculus extends the geometric algebra to include differentiation and integration. The formalism is powerful and can be shown to encompass other mathematical theories including differential geometry and differential forms.[1]

## Differentiation

With a geometric algebra given, let a and b be vectors and let F(a) be a multivector-valued function. The directional derivative of F(a) along b is defined as

$b \cdot \nabla F(a) = \lim_{\epsilon \rightarrow 0}{\frac{F(a + \epsilon b) - F(a)}{\epsilon}}$

provided that the limit exists, where the limit is taken for scalar ε. This is similar to the usual definition of a directional derivative but extends it to functions that are not necessarily scalar-valued.

Next, choose a set of basis vectors $\{e_i\}$ and let $a=x^ie_i$, where we use the Einstein summation notation. This allows the geometric derivative to be treated as an operator

$\nabla=e^i\frac{\partial}{\partial x^i}=e^i\partial_i,$

which is independent of the choice of frame. This is similar to the usual definition of the gradient, but it, too, extends to functions that are not necessarily scalar-valued.

The standard order of operations for the geometric derivative is that it acts only on the function closest to its immediate right. Given two functions F and G, then for example we have

$\nabla FG = (\nabla F)G.$

### Product rule

Although the partial derivative exhibits a product rule, the geometric derivative only partially inherits this property. Consider two functions F and G:

\begin{align}\nabla(FG) &= e^i\partial_i(FG) \\ &= e^i((\partial_iF)G+F(\partial_iG)) \\ &= e^i(\partial_iF)G+e^iF(\partial_iG) \end{align}

Since the geometric product is not commutative with $e^iF \ne Fe^i$ in general, we cannot proceed further without new notation. A solution is to adopt the overdot notation, in which the scope of a geometric derivative with an overdot is the multivector-valued function sharing the same overdot. In this case, if we define

$\dot{\nabla}F\dot{G}=e^iF(\partial_iG),$

then the product rule for the geometric derivative is

$\nabla(FG) = \nabla FG+\dot{\nabla}F\dot{G}$

### Interior and exterior derivative

Let F be an r-grade multivector. Then we can define an additional pair of operators, the interior and exterior derivatives,

$\nabla \cdot F = \langle \nabla F \rangle_{r-1} = e^i \cdot \partial_i F$
$\nabla \wedge F = \langle \nabla F \rangle_{r+1} = e^i \wedge \partial_i F.$

In particular, if F is grade 1 (vector-valued function), then we can write

$\nabla F = \nabla \cdot F + \nabla \wedge F$

and identify the divergence and curl as

$\nabla \cdot F = \operatorname{div} F$
$\nabla \wedge F = I \, \operatorname{curl} F.$

Note, however, that these two operators are considerably weaker than the geometric derivative counterpart for several reasons. Neither the interior derivative operator nor the exterior derivative operator is invertible. Unlike the exterior product, the exterior derivative is not even associative.

## Integration

Let $\{e_1, \, ... \, e_n\}$ be a set of basis vectors that span an n-dimensional vector space. From geometric algebra, we interpret the pseudoscalar $e_1 \wedge e_2 \wedge\cdots\wedge e_n$ to be the signed volume of the n-parallelotope subtended by these basis vectors. If the basis vectors are orthonormal, then this is the unit pseudoscalar.

More generally, we may restrict ourselves to a subset of k of the basis vectors, where $1 \le k \le n$, to treat the length, area, or other general k-volume of a subspace in the overall n-dimensional vector space. We denote these selected basis vectors by $\{e_{i_1}, \, ... \, e_{i_k} \}$. A general k-volume of the k-parallelotope subtended by these basis vectors is the grade k multivector $e_{i_1} \wedge e_{i_2} \wedge\cdots\wedge e_{i_k}$.

Even more generally, we may consider a new set of vectors $\{x^{i_1}e_{i_1}, \, ... \, x^{i_k}e_{i_k} \}$ proportional to the k basis vectors, where each of the $\{x^{i_j}\}$ is a component that scales one of the basis vectors. We are free to choose components as infinitesimally small as we wish as long as they remain nonzero. Since the outer product of these terms can be interpreted as a k-volume, a natural way to define a measure is

\begin{align}\mathrm{d}^kX &= \left(\mathrm{d}x^{i_1} e_{i_1}\right) \wedge \left(\mathrm{d}x^{i_2}e_{i_2}\right) \wedge\cdots\wedge \left(\mathrm{d}x^{i_k}e_{i_k}\right) \\ &= \left( e_{i_1}\wedge e_{i_2}\wedge\cdots\wedge e_{i_k} \right) \mathrm{d}x^{i_1} \mathrm{d}x^{i_2} \cdots \mathrm{d}x^{i_k}\end{align}

The measure is therefore always proportional to the unit pseudoscalar of a k-dimensional subspace of the vector space. Compare the Riemannian volume form in the theory of differential forms. The integral is taken with respect to this measure:

$\int_V F(x) \mathrm{d}^kX = \int_V F(x) \left( e_{i_1}\wedge e_{i_2}\wedge\cdots\wedge e_{i_k} \right) \mathrm{d}x^{i_1} \mathrm{d}x^{i_2} \cdots \mathrm{d}x^{i_k}$

More formally, consider some directed volume V of the subspace. We may divide this volume into a sum of simplices. Let $\{x_i\}$ be the coordinates of the vertices. At each vertex we assign a measure $\Delta U_i(x)$ as the average measure of the simplices sharing the vertex. Then the integral of F(x) with respect to U(x) over this volume is obtained in the limit of finer partitioning of the volume into smaller simplices:

$\int_V F dU = \lim_{n \rightarrow \infty} \sum_{i=1}^n F(x_i) \Delta U_i(x).$

### Fundamental theorem of geometric calculus

The reason for defining the geometric derivative and integral as above is that they allow a strong generalization of Stokes' theorem. Let $\mathsf{L}(A;x)$ be a multivector-valued function of r-grade input A and general position x, linear in its first argument. Then the fundamental theorem of geometric calculus relates the integral of a derivative over the volume V to the integral over its boundary:

 $\int_V \dot{\mathsf{L}} \left(\dot{\nabla} dX;x \right) = \oint_{\partial V} \mathsf{L} (dS;x)$

As an example, let $\mathsf{L}(A;x)=\langle F(x) A I^{-1} \rangle$ for a vector-valued function F(x) and a (n-1)-grade multivector A. We find that

\begin{align}\int_V \dot{\mathsf{L}} \left(\dot{\nabla} dX;x \right) &= \int_V \langle\dot{F}(x)\dot{\nabla} dX I^{-1} \rangle \\ &= \int_V \langle\dot{F}(x)\dot{\nabla} |dX| \rangle \\ &= \int_V \nabla \cdot F(x) |dX| . \end{align}

and likewise

\begin{align}\oint_{\partial V} \mathsf{L} (dS;x) &= \oint_{\partial V} \langle F(x) dS I^{-1} \rangle \\ &= \oint_{\partial V} \langle F(x) \hat{n} |dS| \rangle \\ &= \oint_{\partial V} F(x) \cdot \hat{n} |dS| \end{align}

Thus we recover the divergence theorem,

$\int_V \nabla \cdot F(x) |dX| = \oint_{\partial V} F(x) \cdot \hat{n} |dS|.$

## Covariant derivative

A sufficiently smooth k-surface in an n-dimensional space is deemed a manifold. To each point on the manifold, we may attach a k-blade B that is tangent to the manifold. Locally, B acts as a pseudoscalar of the k-dimensional space. This blade defines a projection of vectors onto the manifold:

$\mathcal{P}_B (A) = (A \cdot B^{-1}) B$

Just as the geometric derivative $\nabla$ is defined over the entire n-dimensional space, we may wish to define an intrinsic derivative $\partial$, locally defined on the manifold:

$\partial F = \mathcal{P}_B (\nabla )F$

(Note: The right hand side of the above may not lie in the tangent space to the manifold. Therefore it is not the same as $\mathcal{P}_B (\nabla F)$, which necessarily does lie in the tangent space.)

If a is a vector tangent to the manifold, then indeed both the geometric derivative and intrinsic derivative give the same directional derivative:

$a \cdot \partial F = a \cdot \nabla F$

Although this operation is perfectly valid, it is not always useful because $\partial F$ itself is not necessarily on the manifold. Therefore we define the covariant derivative to be the forced projection of the intrinsic derivative back onto the manifold:

$a \cdot DF = \mathcal{P}_B (a \cdot \partial F) = \mathcal{P}_B (a \cdot \mathcal{P}_B (\nabla F))$

Since any general multivector can be expressed as a sum of a projection and a rejection, in this case

$a \cdot \partial F = \mathcal{P}_B (a \cdot \partial F) + \mathcal{P}_B^{\perp} (a \cdot \partial F),$

we introduce a new function, the shape tensor $\mathsf{S}(a)$, which satisfies

$F \times \mathsf{S}(a) = \mathcal{P}_B^{\perp} (a \cdot \partial F),$

where $\times$ is the commutator product. In a local coordinate basis $\{e_i\}$ spanning the tangent surface, the shape tensor is given by

$\mathsf{S}(a) = e^i \wedge \mathcal{P}_B^{\perp} (a \cdot \partial e_i).$

Importantly, on a general manifold, the covariant derivative does not commute. In particular, the commutator is related to the shape tensor by

$[a \cdot D, \, b \cdot D]F=-(\mathsf{S}(a) \times \mathsf{S}(b)) \times F.$

Clearly the term $\mathsf{S}(a) \times \mathsf{S}(b)$ is of interest. However it, like the intrinsic derivative, is not necessarily on the manifold. Therefore we can define the Riemann tensor to be the projection back onto the manifold:

$\mathsf{R}(a \wedge b)=-\mathcal{P}_B (\mathsf{S}(a) \times \mathsf{S}(b)).$

Lastly, if F is of grade r, then we can define interior and exterior covariant derivatives as

$D \cdot F = \langle DF \rangle_{r-1}$
$D \wedge F = \langle D F \rangle_{r+1},$

and likewise for the intrinsic derivative.

## Relation to differential geometry

On a manifold, locally we may assign a tangent surface spanned by a set of basis vectors $\{e_i\}$. We can associate the components of a metric tensor, the Christoffel symbols, and the Riemann tensor as follows:

$g_{ij}=e_i \cdot e_j$
$\Gamma^k_{ij}=(e_i \cdot De_j) \cdot e^k$
$R_{ijkl}=(\mathsf{R}(e_i \wedge e_j) \cdot e_k) \cdot e_l$

These relations embed the theory of differential geometry within geometric calculus.

## Relation to differential forms

In a local coordinate system (x1, ..., xn), the coordinate differentials dx1, ..., dxn form a basic set of one-forms within the coordinate chart. Given a multi-index i1, ..., ik with 1 ≤ ipn for 1 ≤ pk, we can define a k-form

$\omega = f_I\mathrm{d}x^I=f_{i_1,i_2\cdots i_k}\mathrm{d}x^{i_1}\wedge \mathrm{d}x^{i_2}\wedge\cdots\wedge \mathrm{d}x^{i_k}.$

We can alternatively introduce a k-grade multivector A as

$A = f_{i_1,i_2\cdots i_k}e^{i_1}\wedge e^{i_2}\wedge\cdots\wedge e^{i_k}$

and a measure

\begin{align}\mathrm{d}^kX &= \left(\mathrm{d}x^{i_1} e_{i_1}\right) \wedge \left(\mathrm{d}x^{i_2}e_{i_2}\right) \wedge\cdots\wedge \left(\mathrm{d}x^{i_k}e_{i_k}\right) \\ &= \left( e_{i_1}\wedge e_{i_2}\wedge\cdots\wedge e_{i_k} \right) \mathrm{d}x^{i_1} \mathrm{d}x^{i_2} \cdots \mathrm{d}x^{i_k}\end{align}.

Apart from a subtle difference in meaning for the exterior product with respect to differential forms versus the exterior product with respect to vectors (indeed one should note that in the former the increments are covectors, whereas in the latter they represent scalars), we see the correspondences of the differential form

$\omega \cong A^{\dagger} \cdot \mathrm{d}^kX = A \cdot \left(\mathrm{d}^kX \right)^{\dagger},$

its derivative

$\mathrm{d}\omega \cong (D \wedge A)^{\dagger} \cdot \mathrm{d}^{k+1}X = (D \wedge A) \cdot \left(\mathrm{d}^{k+1}X \right)^{\dagger},$

and its Hodge dual

$\star\omega \cong (I^{-1} A)^{\dagger} \cdot \mathrm{d}^kX,$

embed the theory of differential forms within geometric calculus.

## History

Following is a diagram summarizing the history of geometric calculus.

History of geometric calculus.

## References

1. ^ David Hestenes, Garrett Sobczyk: Clifford Algebra to Geometric Calculus, a Unified Language for mathematics and Physics (Dordrecht/Boston:G.Reidel Publ.Co., 1984, ISBN 90-277-2561-6