Geometric gradient

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.

## Differentiation

With a geometric algebra given, let

a

{\displaystyle a}

 and

b

{\displaystyle b}

 be vectors and let

F
(
a
)

{\displaystyle F(a)}

 be a multivector-valued function. The directional derivative of

F
(
a
)

{\displaystyle F(a)}

 along

b

{\displaystyle b}

 is defined as

∇

b

F
(
a
# )

lim

ϵ
→
0

F
(
a
ϵ
b
)
−
F
(
a
)

ϵ

,

{\displaystyle \nabla \_{b}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

ϵ

{\displaystyle \epsilon}

. 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

}

{\displaystyle \{e\_i\}}

 and consider the operators, denoted

∂

i

{\displaystyle \partial\_{i}}

, that perform directional derivatives in the directions of

e
i

{\displaystyle e\_i}

:

∂
i

:
F
↦
(
x
↦

∇

e
i

F
(
x
)
)
.

{\displaystyle \partial\_i : F \mapsto (x\mapsto \nabla\_{e\_i} F(x)).}

Then, using the Einstein summation notation, consider the operator:

e
i

∂
i

,

{\displaystyle e^i\partial\_i,}

which means

F
↦

e

i

∂

i

F
,

{\displaystyle F\mapsto e^{i}\partial \_{i}F,}

or, more verbosely:

F
↦
(
x
↦

e

i

∇

e

i

F
(
x
)
)
.

{\displaystyle F\mapsto (x\mapsto e^{i}\nabla \_{e\_{i}}F(x)).}

It can be shown that this operator is independent of the choice of frame, and can thus be used to define the *geometric derivative*:

# ∇

e
i

∂
i

.

{\displaystyle \nabla = e^i\partial\_i.}

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

It can be shown that the directional derivative is linear regarding its direction, that is:

∇

α
a
β
b

=
α

∇

a

β

∇

b

.

{\displaystyle \nabla \_{\alpha a+\beta b}=\alpha \nabla \_{a}+\beta \nabla \_{b}.}

From this follows that the directional derivative is the inner product of its direction by the geometric derivative. All needs to be observed is that the direction

a

{\displaystyle a}

 can be written

# a
(
a
⋅

e
i

)

e
i

{\displaystyle a = (a\cdot e^i) e\_i}

, so that:

∇
a

=

∇

(
a
⋅

e
i

)

e
i

=
(
a
⋅

e
i

)

∇

e
i

=
a
⋅
(

e
i

∇

e
i

# )
a
⋅
∇
.

{\displaystyle \nabla\_a = \nabla\_{(a\cdot e^i)e\_i} = (a\cdot e^i)\nabla\_{e\_i} = a\cdot(e^i\nabla\_{e^i}) = a\cdot \nabla.}

For this reason,

∇
a

F
(
x
)

{\displaystyle \nabla\_a F(x)}

 is often noted

a
⋅
∇
F
(
x
)

{\displaystyle a\cdot \nabla F(x)}

.

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

{\displaystyle F}

 and

G

{\displaystyle G}

, then for example we have

∇
F
# G
(
∇
F
)
G
.

{\displaystyle \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

{\displaystyle F}

 and

G

{\displaystyle G}

:

∇
(
F
G
)

=

e

i

∂

i

(
F
G
)

=

e

i

(
(

∂

i

F
)
G
F
(

∂

i

G
)
)

=

e

i

(

∂

i

F
)
G

e

i

F
(

∂

i

G
)
.

{\displaystyle {\begin{aligned}\nabla (FG)&=e^{i}\partial \_{i}(FG)\\&=e^{i}((\partial \_{i}F)G+F(\partial \_{i}G))\\&=e^{i}(\partial \_{i}F)G+e^{i}F(\partial \_{i}G).\end{aligned}}}

Since the geometric product is not commutative with

e

i

F
≠
F

e

i

{\displaystyle e^{i}F\neq 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

∇
˙

F

G
˙

=

e

i

F
(

∂

i

G
)
,

{\displaystyle {\dot {\nabla }}F{\dot {G}}=e^{i}F(\partial \_{i}G),}

then the product rule for the geometric derivative is

∇
(
F
G
# )
∇
F
G

∇
˙

F

G
˙

.

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

### Interior and exterior derivative

Let

F

{\displaystyle F}

 be an

r

{\displaystyle r}

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

∇
⋅
# F
⟨
∇
F

⟩

r
−
1

=

e

i

⋅

∂

i

F
,

{\displaystyle \nabla \cdot F=\langle \nabla F\rangle \_{r-1}=e^{i}\cdot \partial \_{i}F,}

∇
∧
# F
⟨
∇
F

⟩

r
1

=

e

i

∧

∂

i

F
.

{\displaystyle \nabla \wedge F=\langle \nabla F\rangle \_{r+1}=e^{i}\wedge \partial \_{i}F.}

In particular, if

F

{\displaystyle F}

 is grade 1 (vector-valued function), then we can write

∇
# F
∇
⋅
F
∇
∧
F

{\displaystyle \nabla F=\nabla \cdot F+\nabla \wedge F}

and identify the divergence and curl as

∇
⋅
# F
div
⁡
F
,

{\displaystyle \nabla \cdot F=\operatorname {div} F,}

∇
∧
# F
I

curl
⁡
F
.

{\displaystyle \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.

## Integration

Let

{

e
1

,
…
,

e
n

}

{\displaystyle \{e\_1,\ldots,e\_n\}}

 be a set of basis vectors that span an

n

{\displaystyle n}

-dimensional vector space. From geometric algebra, we interpret the pseudoscalar

e

1

∧

e

2

∧
⋯
∧

e

n

{\displaystyle e\_{1}\wedge e\_{2}\wedge \cdots \wedge e\_{n}}

 to be the signed volume of the

n

{\displaystyle 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

{\displaystyle k}

 of the basis vectors, where

1
≤
k
≤
n

{\displaystyle 1 \leq k \leq n }

, to treat the length, area, or other general

k

{\displaystyle k}

-volume of a subspace in the overall

n

{\displaystyle n}

-dimensional vector space. We denote these selected basis vectors by

{

e

i

1

,
…
,

e

i

k

}

{\displaystyle \{e\_{i\_{1}},\ldots ,e\_{i\_{k}}\}}

. A general

k

{\displaystyle k}

-volume of the

k

{\displaystyle k}

-parallelotope subtended by these basis vectors is the grade

k

{\displaystyle k}

 multivector

e

i

1

∧

e

i

2

∧
⋯
∧

e

i

k

{\displaystyle 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

}

{\displaystyle \{x^{i\_{1}}e\_{i\_{1}},\ldots ,x^{i\_{k}}e\_{i\_{k}}\}}

 proportional to the

k

{\displaystyle k}

 basis vectors, where each of the

{

x

i

j

}

{\displaystyle \{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

{\displaystyle k}

-volume, a natural way to define a measure is

d

k

X

=

(

d

x

i

1

e

i

1

)

∧

(

d

x

i

2

e

i

2

)

∧
⋯
∧

(

d

x

i

k

e

i

k

)

=

(

e

i

1

∧

e

i

2

∧
⋯
∧

e

i

k

)

d

x

i

1

d

x

i

2

⋯
d

x

i

k

.

{\displaystyle {\begin{aligned}d^{k}X&=\left(dx^{i\_{1}}e\_{i\_{1}}\right)\wedge \left(dx^{i\_{2}}e\_{i\_{2}}\right)\wedge \cdots \wedge \left(dx^{i\_{k}}e\_{i\_{k}}\right)\\&=\left(e\_{i\_{1}}\wedge e\_{i\_{2}}\wedge \cdots \wedge e\_{i\_{k}}\right)dx^{i\_{1}}dx^{i\_{2}}\cdots dx^{i\_{k}}.\end{aligned}}}

The measure is therefore always proportional to the unit pseudoscalar of a

k

{\displaystyle 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:

∫

V

F
(
x
)

d

k

# X

∫

V

F
(
x
)

(

e

i

1

∧

e

i

2

∧
⋯
∧

e

i

k

)

d

x

i

1

d

x

i

2

⋯
d

x

i

k

.

{\displaystyle \int \_{V}F(x)\,d^{k}X=\int \_{V}F(x)\left(e\_{i\_{1}}\wedge e\_{i\_{2}}\wedge \cdots \wedge e\_{i\_{k}}\right)dx^{i\_{1}}dx^{i\_{2}}\cdots dx^{i\_{k}}.}

More formally, consider some directed volume

V

{\displaystyle V}

 of the subspace. We may divide this volume into a sum of simplices. Let

{

x
i

}

{\displaystyle \{ x\_i \}}

 be the coordinates of the vertices. At each vertex we assign a measure

Δ

U

i

(
x
)

{\displaystyle \Delta U\_{i}(x)}

 as the average measure of the simplices sharing the vertex. Then the integral of

F
(
x
)

{\displaystyle F(x)}

 with respect to

U
(
x
)

{\displaystyle U(x)}

 over this volume is obtained in the limit of finer partitioning of the volume into smaller simplices:

∫

V

F

d
# U

lim

n
→
∞

∑

# i
1

n

F
(

x

i

)

Δ

U

i

(
x
)
.

{\displaystyle \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

L

(
A
;
x
)

{\displaystyle {\mathsf {L}}(A;x)}

 be a multivector-valued function of

r

{\displaystyle r}

-grade input

A

{\displaystyle A}

 and general position

x

{\displaystyle x}

, linear in its first argument. Then the fundamental theorem of geometric calculus relates the integral of a derivative over the volume

V

{\displaystyle V}

 to the integral over its boundary:

∫

V

L

˙

(

∇
˙

d
X
;
x

)

=

∮

∂
V

L

(
d
S
;
x
)
.

{\displaystyle \int \_{V}{\dot {\mathsf {L}}}\left({\dot {\nabla }}dX;x\right)=\oint \_{\partial V}{\mathsf {L}}(dS;x).}

As an example, let

L

(
A
;
x
# )
⟨
F
(
x
)
A

I

−
1

⟩

{\displaystyle {\mathsf {L}}(A;x)=\langle F(x)AI^{-1}\rangle }

 for a vector-valued function

F
(
x
)

{\displaystyle F(x)}

 and a (

n
−
1

{\displaystyle n-1}

)-grade multivector

A

{\displaystyle A}

. We find that

∫

V

L

˙

(

∇
˙

d
X
;
x

)

=

∫

V

⟨

F
˙

(
x
)

∇
˙

d
X

I

−
1

⟩

=

∫

V

⟨

F
˙

(
x
)

∇
˙

d
X

⟩

=

∫

V

∇
⋅
F
(
x
)

d
X

.

Likewise,

∮

∂
V

L

(
d
S
;
x
)

=

∮

∂
V

⟨
F
(
x
)

d
S

I

−
1

⟩

=

∮

∂
V

⟨
F
(
x
)

n
^

d
S

⟩

=

∮

∂
V

F
(
x
)
⋅

n
^

d
S

.

Thus we recover the divergence theorem,

∫

V

∇
⋅
F
(
x
)

d
X

=

∮

∂
V

F
(
x
)
⋅

n
^

d
S

.

## Covariant derivative

A sufficiently smooth

k

{\displaystyle k}

-surface in an

n

{\displaystyle n}

-dimensional space is deemed a manifold. To each point on the manifold, we may attach a

k

{\displaystyle k}

-blade

B

{\displaystyle B}

 that is tangent to the manifold. Locally,

B

{\displaystyle B}

 acts as a pseudoscalar of the

k

{\displaystyle k}

-dimensional space. This blade defines a projection of vectors onto the manifold:

P

B

(
A
# )
(
A
⋅

B

−
1

)
B
.

{\displaystyle {\mathcal {P}}\_{B}(A)=(A\cdot B^{-1})B.}

Just as the geometric derivative

∇

{\displaystyle \nabla }

 is defined over the entire

n

{\displaystyle n}

-dimensional space, we may wish to define an *intrinsic derivative*

∂

{\displaystyle \partial}

, locally defined on the manifold:

∂
# F

P

B

(
∇
)
F
.

{\displaystyle \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

P

B

(
∇
F
)

{\displaystyle {\mathcal {P}}\_{B}(\nabla F)}

, which necessarily does lie in the tangent space.)

If

a

{\displaystyle a}

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

a
⋅
∂
# F
a
⋅
∇
F
.

{\displaystyle a\cdot \partial F=a\cdot \nabla F.}

Although this operation is perfectly valid, it is not always useful because

∂
F

{\displaystyle \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
⋅
D
# F

P

B

(
a
⋅
∂
F
# )

P

B

(
a
⋅

P

B

(
∇
F
)
)
.

{\displaystyle 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
⋅
∂
# F

P

B

(
a
⋅
∂
F
)

P

B

⊥

(
a
⋅
∂
F
)
,

{\displaystyle 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

S

(
a
)

{\displaystyle {\mathsf {S}}(a)}

, which satisfies

F
×

S

(
a
# )

P

B

⊥

(
a
⋅
∂
F
)
,

{\displaystyle F\times {\mathsf {S}}(a)={\mathcal {P}}\_{B}^{\perp }(a\cdot \partial F),}

where

×

{\displaystyle \times}

 is the commutator product. In a local coordinate basis

{

e
i

}

{\displaystyle \{e\_i\}}

 spanning the tangent surface, the shape tensor is given by

S

(
a
# )

e

i

∧

P

B

⊥

(
a
⋅
∂

e

i

)
.

{\displaystyle {\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
⋅
D
,

b
⋅
D
]
# F
−
(

S

(
a
)
×

S

(
b
)
)
×
F
.

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

Clearly the term

S

(
a
)
×

S

(
b
)

{\displaystyle {\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:

R

(
a
∧
b
# )
−

P

B

(

S

(
a
)
×

S

(
b
)
)
.

{\displaystyle {\mathsf {R}}(a\wedge b)=-{\mathcal {P}}\_{B}({\mathsf {S}}(a)\times {\mathsf {S}}(b)).}

Lastly, if

F

{\displaystyle F}

 is of grade

r

{\displaystyle r}

, then we can define interior and exterior covariant derivatives as

D
⋅
# F
⟨
D
F

⟩

r
−
1

,

{\displaystyle D\cdot F=\langle DF\rangle \_{r-1},}

D
∧
# F
⟨
D
F

⟩

r
1

,

{\displaystyle D\wedge F=\langle DF\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

}

{\displaystyle \{e\_i\}}

. We can associate the components of a metric tensor, the Christoffel symbols, and the Riemann curvature tensor as follows:

g

i
j

=

e

i

⋅

e

j

,

{\displaystyle g\_{ij}=e\_{i}\cdot e\_{j},}

Γ

i
j

k

=
(

e

i

⋅
D

e

j

)
⋅

e

k

,

{\displaystyle \Gamma \_{ij}^{k}=(e\_{i}\cdot De\_{j})\cdot e^{k},}

R

i
j
k
l

=
(

R

(

e

i

∧

e

j

)
⋅

e

k

)
⋅

e

l

.

{\displaystyle 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 (

x
1

,
…
,

x
n

{\displaystyle x^1,\ldots,x^n}

), the coordinate differentials

d

x

1

{\displaystyle dx^{1}}

, ...,

d

x

n

{\displaystyle dx^{n}}

 form a basic set of one-forms within the coordinate chart. Given a multi-index

i

1

,
…
,

i

k

{\displaystyle i\_{1},\ldots ,i\_{k}}

 with

1
≤

i

p

≤
n

{\displaystyle 1\leq i\_{p}\leq n}

 for

1
≤
p
≤
k

{\displaystyle 1\leq p\leq k}

, we can define a

k

{\displaystyle k}

-form

# ω

f

I

d

x

I

=

f

i

1

,

i

2

⋯

i

k

d

x

i

1

∧
d

x

i

2

∧
⋯
∧
d

x

i

k

.

{\displaystyle \omega =f\_{I}\,dx^{I}=f\_{i\_{1},i\_{2}\cdots i\_{k}}\,dx^{i\_{1}}\wedge dx^{i\_{2}}\wedge \cdots \wedge dx^{i\_{k}}.}

We can alternatively introduce a

k

{\displaystyle k}

-grade multivector

A

{\displaystyle A}

 as

# A

f

i

1

,

i

2

⋯

i

k

e

i

1

∧

e

i

2

∧
⋯
∧

e

i

k

{\displaystyle 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

d

k

X

=

(

d

x

i

1

e

i

1

)

∧

(

d

x

i

2

e

i

2

)

∧
⋯
∧

(

d

x

i

k

e

i

k

)

=

(

e

i

1

∧

e

i

2

∧
⋯
∧

e

i

k

)

d

x

i

1

d

x

i

2

⋯
d

x

i

k

.

{\displaystyle {\begin{aligned}d^{k}X&=\left(dx^{i\_{1}}e\_{i\_{1}}\right)\wedge \left(dx^{i\_{2}}e\_{i\_{2}}\right)\wedge \cdots \wedge \left(dx^{i\_{k}}e\_{i\_{k}}\right)\\&=\left(e\_{i\_{1}}\wedge e\_{i\_{2}}\wedge \cdots \wedge e\_{i\_{k}}\right)dx^{i\_{1}}dx^{i\_{2}}\cdots dx^{i\_{k}}.\end{aligned}}}

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

ω
≅

A

†

⋅

d

k

# X
A
⋅

(

d

k

X

)

†

,

{\displaystyle \omega \cong A^{\dagger }\cdot d^{k}X=A\cdot \left(d^{k}X\right)^{\dagger },}

its derivative

d
ω
≅
(
D
∧
A

)

†

⋅

d

k
1

# X
(
D
∧
A
)
⋅

(

d

k
1

X

)

†

,

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

and its Hodge dual

⋆
ω
≅
(

I

−
1

A

)

†

⋅

d

k

X
,

{\displaystyle \star \omega \cong (I^{-1}A)^{\dagger }\cdot d^{k}X,}

embed the theory of differential forms within geometric calculus.