# Calculus/Points, paths, surfaces, and volumes

This chapter will provide an intuitive interpretation of vector calculus using simple concepts such as multi-points, multi-paths, multi-surfaces, and multi-volumes. Scalar fields will not be simply treated as a function $f:\mathbb {R} ^{3}\to \mathbb {R}$ that returns a number given an input point, and vector fields will not be simply treated as a function $\mathbf {F} :\mathbb {R} ^{3}\to \mathbb {R} ^{3}$ that returns a vector given an input point.

## Basic structures

The basic structures are multi-points, multi-paths, multi-surfaces, and multi-volumes.

### Multi-points

A point $\mathbf {q} _{0}$ is an arbitrary position. A "multi-point" is a set of point/weight pairs: $\mathbf {Q} =\{(\mathbf {q} _{1},w_{1}),(\mathbf {q} _{2},w_{2}),...,(\mathbf {q} _{k},w_{k})\}$ where $w_{i}$ is the "weight" that is assigned to point $\mathbf {q} _{i}$. Given two point/weight pairs $(\mathbf {q} ,w_{1})$ and $(\mathbf {q} ,w_{2})$ that cover the same point $\mathbf {q}$, the weights add up to get $(\mathbf {q} ,w_{1}+w_{2})$ which replaces $(\mathbf {q} ,w_{1})$ and $(\mathbf {q} ,w_{2})$. Any pair $(\mathbf {q} ,0)$ is removed. $\mathbf {Q}$ can consist of infinitely many points, and each point may have an infinitesimal weight.
An arbitrary point $\mathbf {q} _{0}$ can be described by the scalar field $\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})=\left\{{\begin{array}{cc}+\infty ^{3}&(\mathbf {q} =\mathbf {q} _{0})\\0&(\mathbf {q} \neq \mathbf {q} _{0})\end{array}}\right.$. This is the "Dirac delta function" centered on point $\mathbf {q} _{0}$. The $+\infty ^{3}$ is the inverse of an infinitely small volume that wraps point $\mathbf {q} _{0}$. To further explain this, let $\omega _{0}(\mathbf {q} _{0},v)$ be a tiny volume with volume $v$ that wraps point $\mathbf {q} _{0}$. $\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})$ can be approximated by $\Delta _{0}(\mathbf {q} ;\mathbf {q} _{0},v)=\left\{{\begin{array}{cc}1/v&(\mathbf {q} \in \omega _{0}(\mathbf {q} _{0},v))\\0&(\mathbf {q} \notin \omega _{0}(\mathbf {q} _{0},v))\end{array}}\right.$. A mass of 1 is being crammed into $\omega _{0}(\mathbf {q} _{0},v)$ yielding an infinitely high density. Since $\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})$ is essentially a density function, it brings with it the units $[{\text{length}}^{-3}]$.
Multi-point $\mathbf {Q} =\{(\mathbf {q} _{1},w_{1}),(\mathbf {q} _{2},w_{2}),...,(\mathbf {q} _{k},w_{k})\}$ can be described by the scalar field $\delta _{0}(\mathbf {q} ;\mathbf {Q} )=\sum _{i=1}^{k}w_{i}\delta _{0}(\mathbf {q} ;\mathbf {q} _{i})$. If $\mathbf {Q}$ consists of infinitely many points with each point having infinitesimal weight, then $\delta _{0}(\mathbf {q} ;\mathbf {Q} )$ is a density function.
In the image below, the multi-point in the left panel is converted to the scalar field in the center panel by averaging the point weight over each cell. The volume of each cell should be infinitesimal. The multi-point in the right panel corresponds to the same scalar field, and is in a more canonical form where oppositely weighted points have cancelled out.
The image below shows how a continuous scalar field $\rho :\mathbb {R} ^{3}\to \mathbb {R}$ can be generated as a collection of points. Consider position $\mathbf {q} _{0}$ and the infinitesimal volume $\omega _{0}(\mathbf {q} _{0},v)$ with volume $v$. The total point weight contained by $\omega _{0}(\mathbf {q} _{0},v)$ is $\iiint _{\mathbf {q} \in \omega _{0}(\mathbf {q} _{0},v)}\rho (\mathbf {q} )dV\approx v\cdot \rho (\mathbf {q} _{0})$. This weight of $v\cdot \rho (\mathbf {q} _{0})$ is then split up over an arbitrarily large number of points that are scattered over the volume $\omega _{0}(\mathbf {q} _{0},v)$.
In summary, a multi-point is denoted by a scalar field that quantifies the density at each point, and any scalar field that quantifies density at each point is best interpreted as a multi-point.

### Multi-paths

A simple path (also called a simple curve) $C$ is an oriented continuous curve that extends from a starting point $C(0)$ to an ending point $C(1)$. Intermediate points are indexed by $t\in [0,1]$ and are denoted by $C(t)$. A simple path should be continuous (no breaks), and may intersect or retrace itself. A "multi-path" is a set of simple-path/weight pairs: $\mathbf {C} =\{(C_{1},w_{1}),(C_{2},w_{2}),...,(C_{k},w_{k})\}$ where $w_{i}$ is the weight that is assigned to path $C_{i}$. Given two path/weight pairs $(C,w_{1})$ and $(C,w_{2})$ that cover the same path $C$, the weights add up to get $(C,w_{1}+w_{2})$ which replaces $(C,w_{1})$ and $(C,w_{2})$. Any pair $(C,0)$ is removed. In addition given two path/weight pairs $(C_{1},w)$ and $(C_{2},w)$ with the same weight $w$ and $C_{1}(1)=C_{2}(0)$, then $C_{1}$ and $C_{2}$ can be linked end-to-end to get the pair $(C_{1}+C_{2},w)$ which replaces $(C_{1},w)$ and $(C_{2},w)$. Assigning a path a negative weight effectively reverses its orientation: if $-C$ denotes path $C$ with the opposite orientation, then $(C,-w)$ is equivalent to $(-C,w)$. $\mathbf {C}$ can consist of infinitely many paths, and each path may have an infinitesimal weight.
An arbitrary curve $C$ can be described by the vector field $\delta _{1}(\mathbf {q} ;C)=\left\{{\begin{array}{cc}(+\infty ^{2}){\hat {\mathbf {n} }}(\mathbf {q} ;C)&(\mathbf {q} \in C)\\\mathbf {0} &(\mathbf {q} \notin C)\end{array}}\right.$. This is the "Dirac delta function" for the curve $C$. ${\hat {\mathbf {n} }}(\mathbf {q} ;C)$ is the unit length tangent vector to path $C$ at point $\mathbf {q} \in C$. ${\hat {\mathbf {n} }}(\mathbf {q} ;C)=\mathbf {0}$ if $\mathbf {q} \notin C$. If there are multiple tangent vectors due to $C$ intersecting itself, then ${\hat {\mathbf {n} }}(\mathbf {q} ;C)$ is the sum of these tangent vectors. The $+\infty ^{2}$ is the inverse of the cross-sectional area of an infinitely thin tube that encloses $C$. To further explain this, let $\omega _{1}(C,a)$ be a thin tube with cross-sectional area $a$ that encloses $C$. $\delta _{1}(\mathbf {q} ;C)$ can be approximated by $\Delta _{1}(\mathbf {q} ;C,a)=\left\{{\begin{array}{cc}(1/a){\hat {\mathbf {n} }}_{*}(\mathbf {q} ;C,a)&(\mathbf {q} \in \omega _{1}(C,a))\\\mathbf {0} &(\mathbf {q} \notin \omega _{1}(C,a))\end{array}}\right.$. ${\hat {\mathbf {n} }}_{*}(\mathbf {q} ;C,a)$ is the generalization of ${\hat {\mathbf {n} }}(\mathbf {q} ;C)$ to the tube $\omega _{1}(C,a)$. A path weight of 1 is being crammed into the cross-sectional area of $\omega _{1}(C,a)$ yielding an infinitely high path density. Since $\delta _{1}(\mathbf {q} ;C)$ is essentially a density over area, it brings with it the units $[{\text{length}}^{-2}]$.
The image to the right gives a depiction of the Dirac delta function for a simple curve. The vector field $\delta _{1}(\mathbf {q} ;C)$ is $\mathbf {0}$ everywhere outside of an infinitely thin tube that encloses the path. Inside the tube, the vectors are parallel to the path, and have a magnitude equal to the inverse of the cross-sectional area. The Dirac delta function is the limit as the tube becomes infinitely thin.
Multi-path $\mathbf {C} =\{(C_{1},w_{1}),(C_{2},w_{2}),...,(C_{k},w_{k})\}$ can be described by the vector field $\delta _{1}(\mathbf {q} ;\mathbf {C} )=\sum _{i=1}^{k}w_{i}\delta _{1}(\mathbf {q} ;C_{i})$. If $\mathbf {C}$ consists of infinitely many paths with each path having infinitesimal weight, then $\delta _{1}(\mathbf {q} ;\mathbf {C} )$ is a flow density function.
In the image below, the multi-path in the left panel is converted to the vector field in the center panel by computing the total displacement in each cell and averaging over the volume. The volume of each cell should be infinitesimal. The multi-path in the right panel corresponds to the same vector field, and is in a more canonical form where the individual paths do not cross each other.
In summary, a multi-path is denoted by a vector field that quantifies the path/flow density at each point, and any vector field that quantifies a flow density at each point (such as current density) is best interpreted as a multi-path. (Flow density is a vector that points in the net direction of a flow, and has a length equal to the flow rate per unit area through a surface that is perpendicular to the net flow.)

### Multi-surfaces

A simple surface $\sigma$ is an oriented continuous surface. A simple surface may intersect or fold back on itself. A "multi-surface" is a set of simple-surface/weight pairs: $\mathbf {S} =\{(\sigma _{1},w_{1}),(\sigma _{2},w_{2}),...,(\sigma _{k},w_{k})\}$ where $w_{i}$ is the weight that is assigned to surface $\sigma _{i}$. Given two surface/weight pairs $(\sigma ,w_{1})$ and $(\sigma ,w_{2})$ that cover the same surface $\sigma$, the weights add up to get $(\sigma ,w_{1}+w_{2})$ which replaces $(\sigma ,w_{1})$ and $(\sigma ,w_{2})$. Any pair $(\sigma ,0)$ is removed. In addition given two surface/weight pairs $(\sigma _{1},w)$ and $(\sigma _{2},w)$ with the same weight $w$, then $\sigma _{1}$ and $\sigma _{2}$ can be combined to get the pair $(\sigma _{1}+\sigma _{2},w)$ which replaces $(\sigma _{1},w)$ and $(\sigma _{2},w)$. Assigning a surface a negative weight effectively reverses its orientation: if $-\sigma$ denotes surface $\sigma$ with the opposite orientation, then $(\sigma ,-w)$ is equivalent to $(-\sigma ,w)$. $\mathbf {S}$ can consist of infinitely many surfaces, and each surface may have an infinitesimal weight.
An arbitrary surface $\sigma$ can be described by the vector field $\delta _{2}(\mathbf {q} ;\sigma )=\left\{{\begin{array}{cc}(+\infty ){\hat {\mathbf {n} }}(\mathbf {q} ;\sigma )&(\mathbf {q} \in \sigma )\\\mathbf {0} &(\mathbf {q} \notin \sigma )\end{array}}\right.$. This is the "Dirac delta function" for the surface $\sigma$. ${\hat {\mathbf {n} }}(\mathbf {q} ;\sigma )$ is the unit length normal vector to surface $\sigma$ at point $\mathbf {q} \in \sigma$. ${\hat {\mathbf {n} }}(\mathbf {q} ;\sigma )=\mathbf {0}$ if $\mathbf {q} \notin \sigma$. If there are multiple normal vectors due to $\sigma$ intersecting itself, then ${\hat {\mathbf {n} }}(\mathbf {q} ;\sigma )$ is the sum of these normal vectors. The $+\infty$ is the inverse of the thickness of an infinitely thin membrane that encloses $\sigma$. To further explain this, let $\omega _{2}(\sigma ,t)$ be a thin membrane with thickness $t$ that encloses $\sigma$. $\delta _{2}(\mathbf {q} ;\sigma )$ can be approximated by $\Delta _{2}(\mathbf {q} ;\sigma ,t)=\left\{{\begin{array}{cc}(1/t){\hat {\mathbf {n} }}_{*}(\mathbf {q} ;\sigma ,t)&(\mathbf {q} \in \omega _{2}(\sigma ,t))\\\mathbf {0} &(\mathbf {q} \notin \omega _{2}(\sigma ,t))\end{array}}\right.$. ${\hat {\mathbf {n} }}_{*}(\mathbf {q} ;\sigma ,t)$ is the generalization of ${\hat {\mathbf {n} }}(\mathbf {q} ;\sigma )$ to the membrane $\omega _{2}(\sigma ,t)$. A surface weight of 1 is being sandwiched into the thickness of $\omega _{2}(\sigma ,t)$ yielding an infinitely high surface density. Since $\delta _{2}(\mathbf {q} ;\sigma )$ is essentially a density over length, it brings with it the units $[{\text{length}}^{-1}]$.
Multi-surface $\mathbf {S} =\{(\sigma _{1},w_{1}),(\sigma _{2},w_{2}),...,(\sigma _{k},w_{k})\}$ can be described by the vector field $\delta _{2}(\mathbf {q} ;\mathbf {S} )=\sum _{i=1}^{k}w_{i}\delta _{2}(\mathbf {q} ;\sigma _{i})$. If $\mathbf {S}$ consists of infinitely many surfaces with each surface having infinitesimal weight, then $\delta _{2}(\mathbf {q} ;\mathbf {S} )$ is a rate-of-gain function.
In the image below, the multi-surface in the left panel is converted to the vector field in the center panel by computing the total surface in each cell and averaging over the volume. The volume of each cell should be infinitesimal. The multi-surface in the right panel corresponds to the same vector field, and is in a more canonical form where the individual surfaces do not intersect each other.
In summary, a multi-surface is denoted by a vector field that quantifies the rate of gain at each point. To describe the rate-of-gain, imagine that passing through a surface in the preferred direction gives "energy". The rate of gain is a vector that points in the direction that yields the greatest rate of energy increase per unit length, and has a length equal to the maximum rate of energy increase per unit length. Any vector field that quantifies a rate of gain at each point (such as a force field) is best interpreted as a multi-surface.

### Multi-volumes

A volume $\Omega$ is an arbitrary region of space. A "multi-volume" is a set of volume/weight pairs: $\mathbf {U} =\{(\Omega _{1},w_{1}),(\Omega _{2},w_{2}),...,(\Omega _{k},w_{k})\}$ where $w_{i}$ is the "weight" that is assigned to volume $\Omega _{i}$. Given two volume/weight pairs $(\Omega ,w_{1})$ and $(\Omega ,w_{2})$ that cover the same volume $\Omega$, the weights add up to get $(\Omega ,w_{1}+w_{2})$ which replaces $(\Omega ,w_{1})$ and $(\Omega ,w_{2})$. Any pair $(\Omega ,0)$ is removed. In addition given two volume/weight pairs $(\Omega _{1},w)$ and $(\Omega _{2},w)$ with the same weight $w$ and $\Omega _{1}\cap \Omega _{2}=\emptyset$, then $\Omega _{1}$ and $\Omega _{2}$ can be combined to get the pair $(\Omega _{1}\cup \Omega _{2},w)$ which replaces $(\Omega _{1},w)$ and $(\Omega _{2},w)$. $\mathbf {U}$ can consist of infinitely many volumes, and each volume may have an infinitesimal weight.
An arbitrary volume $\Omega$ can be described by the scalar field $\delta _{3}(\mathbf {q} ;\Omega )=\left\{{\begin{array}{cc}1&(\mathbf {q} \in \Omega )\\0&(\mathbf {q} \notin \Omega )\end{array}}\right.$. This is the "Dirac delta function" analog for volumes, and is essentially an indicator function that indicates whether or not a position is contained by $\Omega$ or not, 1 being yes and 0 being no. Since $\delta _{3}(\mathbf {q} ;\Omega )$ is simply an indicator function, it brings with it no units (it is dimensionless).
Multi-volume $\mathbf {U} =\{(\Omega _{1},w_{1}),(\Omega _{2},w_{2}),...,(\Omega _{k},w_{k})\}$ can be described by the scalar field $\delta _{3}(\mathbf {q} ;\mathbf {U} )=\sum _{i=1}^{k}w_{i}\delta _{3}(\mathbf {q} ;\Omega _{i})$. If $\mathbf {U}$ consists of infinitely many volumes with each volume having infinitesimal weight, then $\delta _{3}(\mathbf {q} ;\mathbf {U} )$ is a potential function.
In the image below, the multi-volume in the left panel is converted to the scalar field in the center panel by averaging the volume weight in each cell. The volume of each cell should be infinitesimal. The multi-volume in the right panel corresponds to the same scalar field, and is in a more canonical form where oppositely weighted volumes have cancelled out, and the remaining volume has diffused to fill each cell.
In summary, a multi-volume is denoted by a scalar field that quantifies a potential at each point, and any scalar field that quantifies a potential at each point is best interpreted as a multi-volume.

### At infinity

An important requirement is that all multi-points, multi-paths, multi-surfaces, and multi-volumes not extend to infinity. All structures can extend to an arbitrarily large range, as long as this range is not unbounded. Allowing the structures to extend to infinity will cause problems in the later discussions.
| Paths that extend to infinity are generally not allowed for most theorems related to vector calculus. | Surfaces that extend to infinity are generally not allowed for most theorems related to vector calculus. | Volumes that extend to infinity are generally not allowed for most theorems related to vector calculus. |

## Totals

These sections will describe the total weight of multi-points, the total displacement of multi-paths, the total surface of multi-surfaces, and the total volumes of multi-volumes.

### Total point weight

Given a multi-point $\mathbf {Q} =\{(\mathbf {q} _{1},w_{1}),(\mathbf {q} _{2},w_{2}),...,(\mathbf {q} _{k},w_{k})\}$, the total point weight is clearly $\sum _{i=1}^{k}w_{i}$. Given a scalar field $\rho$ that denotes a multi-point, the total weight of $\rho$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\rho (\mathbf {q} )dV$. Given a simple point $\mathbf {q} _{0}$, the total weight is 1 so $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})dV=1$.

### Total displacement

Given a simple path $C$ that starts at point $C(0)$ and ends at point $C(1)$, the total displacement generated by $C$ is $\int _{\mathbf {q} \in C}d\mathbf {q} =C(1)-C(0)$. This displacement is solely dependent on the endpoints as indicated by the top image to the right.
The displacement generated by a closed loop is $\mathbf {0}$.
Given a multi-path $\mathbf {C} =\{(C_{1},w_{1}),(C_{2},w_{2}),...,(C_{k},w_{k})\}$, the total displacement generated by $\mathbf {C}$ is $\sum _{i=1}^{k}w_{i}\int _{\mathbf {q} \in C_{i}}d\mathbf {q} =\sum _{i=1}^{k}w_{i}(C_{i}(1)-C_{i}(0))$.
Given a vector field $\mathbf {J}$ that denotes a multi-path, the total displacement generated by $\mathbf {J}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )dV$. Since the displacement generated by a simple path $C$ is $\int _{\mathbf {q} \in C}d\mathbf {q} =C(1)-C(0)$, it is the case that $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{1}(\mathbf {q} ;C)dV=\int _{\mathbf {q} \in C}d\mathbf {q} =C(1)-C(0)$.
One important observation from $\int _{\mathbf {q} \in C}d\mathbf {q} =\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{1}(\mathbf {q} ;C)dV$ is that given a path integral over path $C$, the differential $d\mathbf {q}$ is equal to $\delta _{1}(\mathbf {q} ;C)dV$ in a volume integral: $\int _{\mathbf {q} \in C}f(\mathbf {q} ,d\mathbf {q} )=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}f(\mathbf {q} ,\delta _{1}(\mathbf {q} ;C)dV)$ provided that function $f$ is linear in the second parameter. In the lower image to the right, the displacement differential $d\mathbf {q} ={\hat {\mathbf {n} }}\cdot \Delta l$ is equated to the volume differential $\left({\frac {\hat {\mathbf {n} }}{\Delta A}}\right)dV=\delta _{1}(\mathbf {q} ;C)dV$ by diffusing the path over an infintely thin cross-sectional area $\Delta A$. The integrand at points outside of the infinitely thin tube is 0: for all points $\mathbf {q} \notin C$, $f(\mathbf {q} ,\delta _{1}(\mathbf {q} ;C)dV)=f(\mathbf {q} ,\mathbf {0} )=0$.

### Total surface vector

Given an arbitrary oriented surface $\sigma$, its "counter-clockwise boundary", denoted by $\partial \sigma$, is the boundary of $\sigma$ whose orientation is determined in the following manner: Looking at $\sigma$ so that the preferred direction through $\sigma$ is oriented towards the viewer, the boundary $\partial \sigma$ wraps $\sigma$ in a counter-clockwise direction.
Given a flat surface as shown in the image to the right, this surface can be quantified by the "surface vector" which is a vector that is perpendicular (normal) to the surface in the preferred orientation, and has a length equal to the area of the surface. In the image to the right, a flat surface has an area of $A$ and is oriented to be perpendicular to unit-length normal vector ${\hat {\mathbf {n} }}$. The "surface vector" of this surface is $A\cdot {\hat {\mathbf {n} }}$.
Given a non-flat surface $\sigma$, the total surface vector of $\sigma$ is computed by summing the surface vectors of each infinitesimal portion of $\sigma$. The total surface vector is $\mathbf {S} =\iint _{\mathbf {q} \in \sigma }d\mathbf {S}$.
In a manner similar to how the total displacement of a path is solely a function of the endpoints, the total surface vector of a surface is solely a function of its counter-clockwise boundary. This is not intuitive, and will be explained in greater detail below using two approaches:

#### Generalizing from surfaces in 2D space

Below are shown two images related to surface vectors in 2D space. The image to the left shows surface vectors in 2D space. In 2 dimensions, surfaces are called 1D surfaces and are similar to paths. The boundary of a 1D surface consists of 2 points. The surface vector of a 1D surface segment is a 90 degree rotation of the segment and is oriented in the direction of the surface's orientation. The total surface vector of a 1D surface is the sum of all of the surface vectors of the individual components. For each component of the surface, the surface vector is a 90 degree rotation of the displacement that traverses the component, so the total surface vector is a 90 degree rotation of the displacement between the points that form the boundary of the surface. This proves that in two dimensions, the total surface vector depends only on the boundary of the 1D surface.
The image to the right extrudes the 1D surfaces in two dimensional space into 2D "ribbons" in 3 dimensional space. At the top a closed "ribbon" is shown. This "ribbon" is a surface that is always parallel to the vertical dimension, and whose boundary forms two identical loops that are vertically displaced from each other. The boundary loops are also perpendicular to the vertical dimension. The ribbon itself is partitioned into tiny rectangles whose height is equal to that of the ribbon. To the bottom left, a view of the same ribbon from the top down is shown. It can be seen that the length of the each surface vector is proportional to the length of the corresponding rectangular segment, since the heights are all uniform. To the bottom right, by rotating the surface vectors 90 degrees around the vertical dimension, the surface vectors now sum to $\mathbf {0}$, so the sum of the unrotated surface vectors is also $\mathbf {0}$.
| This image depicts how in 2 dimensions, the total surface vector of a 1D surface is a 90 degree rotation of the displacement between the two endpoints (the boundary of a 1D surface), and is therefore purely a function of the endpoints. In the left panel, a 1D surface is a sequence of black line segments, and the surface vectors of each segment are denoted by the dashed red arrows. Each surface vector is a 90 degree rotation of the displacement along the surface. The long grey line is the net displacement between the endpoints of the surface, and the dashed pink arrow is a 90 degree rotation of this net displacement. In the right panel, the pink arrow is shown as the sum of the dashed red arrow vectors, hence the "total surface" is purely a function of the 1D surface's endpoints. | This image demonstrates that the total surface vector of a surface that is a closed ribbon is 0. The top image shows a surface that is a closed ribbon where the width of the ribbon is constant, the width is always parallel to the vertical dimension, and the edge is always perpendicular to the vertical dimension. The surface is sub-divided into tiny rectangular portions, the surface vectors of which are shown. The lower-left image shows the same surface from a top down perspective. In the lower-right image, the surface vectors are all rotated 90 degrees counter-clockwise around the vertical dimension and clearly sum to 0. |
The fact that the total surface vector of a closed ribbon is $\mathbf {0}$ means that if relief is added to a surface without changing its boundaries, the total surface vector is conserved. The two left images below give examples of distorting the interior of a surface by hammering in relief. The vertical surfaces introduced by the relief are ribbons which contribute $\mathbf {0}$ to the total surface vector, while the horizontal surfaces are simply displaced vertically be the relief. The rightmost image below shows how the total surface vector is preserved if the "texture" of the surface at infinitesimal scales is converted from "steps" (a union of horizontal and vertical surfaces) to "smooth slopes" and vice versa. The surface formed from the red and green planes is a step, while the surface formed from the blue plane is a slope. These two surfaces can be seen to have equal total surface vectors from the right-angled triangle at the right side of the image.
| Adding elevation (depression in this image) or relief to a surface does not change the total surface vector. The red colored horizontal surfaces are clearly conserved, albeit at different elevations. The green colored vertical surfaces sum to 0 at each tier/elevation. | Adding elevation (depression in this image) or relief to a surface does not change the total surface vector. The horizontal surfaces are clearly conserved, albeit at different elevations. The green colored vertical surfaces sum to their initial value above the lower red surface, and sum to 0 beneath the lower red surface. | In this image there are two surfaces. The first surface is the union of the red and green planes, and the counter-clockwise boundary is shown by the thick black line. The second surface is the blue plane and the counter-clockwise boundary is shown by the dashed blue line. The surface vectors of the red, green, and blue planes are shown. The total surface vector of the first surface is the sum of the surface vectors of the red and green planes, and is equal to the surface vector of the blue plane. This all implies that the total surface vector of a sloped flat surface is unchanged by replacing the surface with its horizontal and vertical components (projections). |

#### Generalizing from displacement vectors

The total displacement along a simple oriented curve can be used to compute the net displacement in a specific direction. Given a simple oriented curve $C$ and an oriented straight line with the direction indicated by normal vector ${\hat {\mathbf {n} }}$, the total displacement $\Delta \mathbf {q}$ along $C$ can be used to compute the net displacement in the direction indicated by the line. This displacement is ${\hat {\mathbf {n} }}\cdot \Delta \mathbf {q}$, and depends only on the endpoints of the curve.
In a direct analogy, given a simple oriented surface $\sigma$ with counter-clockwise boundary $\partial \sigma$, and an oriented flat plane whose surface normal is ${\hat {\mathbf {n} }}$, a quantity of interest is the total signed area of $\sigma$ perpendicularly projected onto the plane. The signed area that is projected by a flat infinitesimal portion of $\sigma$ with surface vector $d\mathbf {S}$ is ${\hat {\mathbf {n} }}\cdot d\mathbf {S}$, and the total signed area is $\iint _{\mathbf {q} \in \sigma }{\hat {\mathbf {n} }}\cdot d\mathbf {S} ={\hat {\mathbf {n} }}\cdot \iint _{\mathbf {q} \in \sigma }d\mathbf {S} ={\hat {\mathbf {n} }}\cdot \mathbf {S}$ where $\mathbf {S}$ is the total surface vector of $\sigma$.
The total signed projected area ${\hat {\mathbf {n} }}\cdot \mathbf {S}$ onto the plane is purely a function of the boundary $\partial \sigma$, and does not depend on how $\sigma$ fills its boundary $\partial \sigma$. This is much more obvious and clearer than the claim that the total surface vector $\mathbf {S}$ is only a function of $\partial \sigma$: the area enclosed by a boundary in 2D space is purely a function of that boundary. Since the projected area is signed, "upside down" surfaces project negative area, and folds and overhangs cancel each other out.
Since ${\hat {\mathbf {n} }}\cdot \mathbf {S}$ is purely a function of $\partial \sigma$ for all choices of plane orientation ${\hat {\mathbf {n} }}$, then the total surface vector $\mathbf {S}$ is purely a function of $\partial \sigma$.
| Given an arbitrary oriented path, the total displacement covered by the perpendicularly projected path onto an oriented straight line does not depend on the placement of the interior points of the path. The displacement only depends on the endpoints. Since this is true no matter the choice of straight line, the total 3D displacement vector generated by an oriented curve is purely a function of its endpoints, and does not change if the interior points are moved. | The total signed area of the projection of an oriented surface onto an oriented flat plane depends only on the boundary and not on any of the interior points. The "shadow" does not change if the interior points are moved around. If the surface is deformed so that there is an "overhang" where some projected points fall outside of the projected boundary, such as in the example of the right, these points cancel out with the points on the opposite side (top or bottom) of the overhang. An upside down surface projects negative area, and in the example on the right, all negative projected area is cancelled out with the positive area projected by the upright surface on top of the overhang. |
* Computing the signed projected area of a flat surface onto a flat plane is equivalent to computing the signed projected length of the surface vector onto the line that is perpendicular to the plane.

#### Summary

The total surface vector generated by a closed surface is $\mathbf {0}$.
Given a multi-surface $\mathbf {S} =\{(\sigma _{1},w_{1}),(\sigma _{2},w_{2}),...,(\sigma _{k},w_{k})\}$ the total surface vector generated by $\mathbf {S}$ is $\sum _{i=1}^{k}w_{i}\iint _{\mathbf {q} \in \sigma _{i}}d\mathbf {S}$.
Given a vector field $\mathbf {F}$ that denotes a multi-surface, the total surface vector generated by $\mathbf {F}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )dV$. Since the surface vector generated by simple surface $\sigma$ is $\iint _{\mathbf {q} \in \sigma }d\mathbf {S}$, it is the case that $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{2}(\mathbf {q} ;\sigma )dV=\iint _{\mathbf {q} \in \sigma }d\mathbf {S}$. One important observation is that given a surface integral over $\sigma$, the differential $d\mathbf {S}$ is equal to $\delta _{2}(\mathbf {q} ;\sigma )dV$ in a volume integral: $\iint _{\mathbf {q} \in \sigma }f(\mathbf {q} ,d\mathbf {S} )=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}f(\mathbf {q} ,\delta _{2}(\mathbf {q} ;\sigma )dV)$ provided that function $f$ is a linear in the second parameter.

### Total volume

Consider a multi-volume $\mathbf {U} =\{(\Omega _{1},w_{1}),(\Omega _{2},w_{2}),...,(\Omega _{k},w_{k})\}$, where the volumes of $\Omega _{1},\Omega _{2},...,\Omega _{k}$ are respectively $V_{1},V_{2},...,V_{k}$, then the total volume of $\mathbf {U}$ is $\sum _{i=1}^{k}w_{i}V_{i}$. Each volume $V_{i}$ can be computed by $V_{i}=\iiint _{\mathbf {q} \in \Omega _{i}}dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{3}(\mathbf {q} ;\Omega _{i})dV$. The total volume of $\mathbf {U}$ is $V=\sum _{i=1}^{k}w_{i}V_{i}=\sum _{i=1}^{k}w_{i}\iiint _{\mathbf {q} \in \Omega _{i}}dV=\sum _{i=1}^{k}w_{i}\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{3}(\mathbf {q} ;\Omega _{i})dV$ $=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\left(\sum _{i=1}^{k}w_{i}\delta _{3}(\mathbf {q} ;\Omega _{i})\right)dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{3}(\mathbf {q} ;\mathbf {U} )dV$.
If a multi-volume $\mathbf {U}$ can be denoted by scalar field $U$, then the volume of $\mathbf {U}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U(\mathbf {q} )dV$.
Given an arbitrary volume $\Omega$, a volume integral over $\Omega$ can be converted to a volume integral over $\mathbb {R} ^{3}$ by replacing the differential $dV$ with $\delta _{3}(\mathbf {q} ;\Omega )dV$:
$\iiint _{\mathbf {q} \in \Omega }f(\mathbf {q} ,dV)=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}f(\mathbf {q} ,\delta _{3}(\mathbf {q} ;\Omega )dV)$ provided that $f$ is linear in the second parameter.

## Intersections

The union of two multi-points denoted by scalar fields $\rho _{1}$ and $\rho _{2}$ is simply $\rho _{1}+\rho _{2}$, and the same is true for the union of two multi-paths, the union of two multi-surfaces, and the union of two multi-volumes. The union of two structures with different types, such as a multi-point with a multi-path, is forbidden however.
Unions :
| structure | multi-point $\rho _{2}$ | multi-path $\mathbf {J} _{2}$ | multi-surface $\mathbf {F} _{2}$ | multi-volume $U_{2}$ |
| multi-point $\rho _{1}$ | multi-point $\rho _{1}+\rho _{2}$ | n/a | n/a | n/a |
| multi-path $\mathbf {J} _{1}$ | n/a | multi-path $\mathbf {J} _{1}+\mathbf {J} _{2}$ | n/a | n/a |
| multi-surface $\mathbf {F} _{1}$ | n/a | n/a | multi-surface $\mathbf {F} _{1}+\mathbf {F} _{2}$ | n/a |
| multi-volume $U_{1}$ | n/a | n/a | n/a | multi-volume $U_{1}+U_{2}$ |
The intersection on the other hand, is less trivial and can occur between structures of different types.

### Point-Volume intersections

When a point $\mathbf {q}$ with weight $w_{1}$ intersects a volume $\Omega$ with weight $w_{2}$, then the intersection is point $\mathbf {q}$ with weight $w_{1}w_{2}$. Given a multi-point and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple point with each simple volume. The image below gives an example of the intersection of a multi-point with a multi-volume.
Given a multi-point with scalar field $\rho$, and a multi-volume with scalar field $U$, then the intersection is a multi-point with scalar field $\rho U$.
The total intersection between a multi-point $\rho$ and a multi-volume $U$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\rho (\mathbf {q} )U(\mathbf {q} )dV$.
If $\rho$ denotes a simple point $\mathbf {q} _{0}$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})U(\mathbf {q} )dV=U(\mathbf {q} _{0})$.
If $U$ denotes a simple volume $\Omega$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\rho (\mathbf {q} )\delta _{3}(\mathbf {q} ;\Omega )dV=\iiint _{\mathbf {q} \in \Omega }\rho (\mathbf {q} )dV$.

### Path-Surface intersections

When a path $C$ with weight $w_{1}$ intersects a surface $\sigma$ with weight $w_{2}$ at point $\mathbf {q}$, then the intersection is point $\mathbf {q}$ with weight $\pm w_{1}w_{2}$. The weight is $+w_{1}w_{2}$ if $C$ passes through $\sigma$ in the direction in which $\sigma$ is oriented. The weight is $-w_{1}w_{2}$ if $C$ passes through $\sigma$ opposite the direction in which $\sigma$ is oriented. Given a multi-path and a multi-surface, the intersection is the sum of the pair-wise intersections of each simple path with each simple surface. The images below give examples of the intersections of a multi-path with a multi-surface.
| A 2D image showing the intersection of a multi-path (dark blue dashed curves) with a multi-surface (dark red solid curves). Positive intersection points (red circles) occur when a path intersects a surface in the preferred direction. Negative intersection points (teal circles) occur when a path intersects a surface in the opposite direction. The intersection is effectively a multi-point. | A 3D image showing the intersection of a simple path (red curve) with a simple surface (green surface with the counter-clockwise boundary highlighted). The positive intersection points are denoted by red "+" signs, and the negative intersection points are denoted by blue "-" signs. | The intersection between a multi-path shown as a blue tube with a multi-surface shown as layers of red sheets. Vector F is the flow density through the blue tube. Vector G is the surface density in the red sheets. The green parallelogram is a 2D projection of the volume of the intersection. The intersection points become more dilute as the angle theta increases, so the intersection point density is the dot product of F and G. |
In the image above to the far right, the multi-path is denoted by a vector field which has the value $\mathbf {F}$ inside the blue tube, and is $\mathbf {0}$ everywhere else. The multi-surface is denoted by a vector field which has the value $\mathbf {G}$ among the red sheets, and is $\mathbf {0}$ everywhere else. The total path weight in the blue tube is $|\mathbf {F} |\Delta A$. The total surface weight in the red sheets is $|\mathbf {G} |\Delta t$. The total weight of all the intersection points is $(|\mathbf {F} |\Delta A)(|\mathbf {G} |\Delta t)=|\mathbf {F} ||\mathbf {G} |\Delta A\Delta t$. The volume that the intersection points are evenly spread out in is $\Delta A\Delta t/\cos \theta$. The intersection point density is ${\frac {|\mathbf {F} ||\mathbf {G} |\Delta A\Delta t}{\Delta A\Delta t/\cos \theta }}=|\mathbf {F} ||\mathbf {G} |\cos \theta =\mathbf {F} \cdot \mathbf {G}$.
Given a multi-path with vector field $\mathbf {J}$, and a multi-surface with vector field $\mathbf {F}$, then the intersection is a multi-point with scalar field $\mathbf {J} \cdot \mathbf {F}$.
The total intersection between a multi-path $\mathbf {J}$ and a multi-surface $\mathbf {F}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} (\mathbf {q} )\cdot \mathbf {F} (\mathbf {q} ))dV$.
If $\mathbf {J}$ is a simple path $C$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\delta _{1}(\mathbf {q} ;C)\cdot \mathbf {F} (\mathbf {q} ))dV=\int _{\mathbf {q} \in C}\mathbf {F} (\mathbf {q} )\cdot d\mathbf {q}$.
If $\mathbf {F}$ is a simple surface $\sigma$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} (\mathbf {q} )\cdot \delta _{2}(\mathbf {q} ;\sigma ))dV=\iint _{\mathbf {q} \in \sigma }\mathbf {J} (\mathbf {q} )\cdot d\mathbf {S}$.

### Path-Volume intersections

When a path $C$ with weight $w_{1}$ intersects a volume $\Omega$ with weight $w_{2}$, then the intersection is path $C\cap \Omega$ with weight $w_{1}w_{2}$. Given a multi-path and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple path with each simple volume. The image below gives an example of the intersection of a multi-path with a multi-volume.
Given a multi-path with vector field $\mathbf {J}$, and a multi-volume with scalar field $U$, then the intersection is a multi-path with vector field $\mathbf {J} U$.
The total intersection between a multi-path $\mathbf {J}$ and a multi-volume $U$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )U(\mathbf {q} )dV$.
If $\mathbf {J}$ denotes a simple path $C$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{1}(\mathbf {q} ;C)U(\mathbf {q} )dV=\int _{\mathbf {q} \in C}U(\mathbf {q} )d\mathbf {q}$.
If $U$ denotes a simple volume $\Omega$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\delta _{3}(\mathbf {q} ;\Omega )dV=\iiint _{\mathbf {q} \in \Omega }\mathbf {J} (\mathbf {q} )dV$.

### Surface-Surface intersections

When a surface $\sigma _{1}$ with weight $w_{1}$ intersects a surface $\sigma _{2}$ with weight $w_{2}$, then the intersection is the path $\sigma _{1}\cap \sigma _{2}$ with weight $w_{1}w_{2}$. The orientation given to path $\sigma _{1}\cap \sigma _{2}$ is defined as follows: viewing the intersection where the surface normal vectors of $\sigma _{1}$ and $\sigma _{2}$ are oriented towards the viewer, the intersection path has $\sigma _{1}$ to its right, and $\sigma _{2}$ to its left. Put another way, the intersection path is oriented according to the "right-hand rule" where the surface normals of $\sigma _{1}$ are the "x" direction, and the surface normals of $\sigma _{2}$ are the "y" direction. The images below give examples of the intersections of a multi-surface with a multi-surface.
| A 3D image that shows the intersection of 2 surfaces. Surface 1 is blue and the normal vectors are oriented upwards. Surface 2 is red and the normal vectors are oriented to the right. The intersection is the black curve. The orientation of the intersection curve is determined via the right-hand rule with the surface normals of surface 1 as the "x" direction, and the surface normals of surface 2 as the "y" direction. | The intersection between two multi-surfaces. The first multi-surface is the layered blue sheets, and the second multi-surface is the layered red sheets. Vector F is the surface density in the blue sheets. Vector G is the surface density in the red sheets. The green parallelogram is a 2D cross-section of the prism that forms the intersection. The intersection paths become more dilute the further angle theta deviates from 90 degrees, so the intersection path density is the cross product of F and G. The intersection paths are also oriented out of the screen in this example. |
In the image above to the right, the first multi-surface is denoted by a vector field that has the value $\mathbf {F}$ among the blue sheets, and is $\mathbf {0}$ everywhere else. The second multi-surface is denoted by a vector field that has the value $\mathbf {G}$ among the red sheets, and is $\mathbf {0}$ everywhere else. The total surface weight in the blue sheets is $|\mathbf {F} |\Delta t_{1}$, and the total surface weight in the red sheets is $|\mathbf {G} |\Delta t_{2}$. The total weight of all the intersection paths is $(|\mathbf {F} |\Delta t_{1})(|\mathbf {G} |\Delta t_{2})=|\mathbf {F} ||\mathbf {G} |\Delta t_{1}\Delta t_{2}$. The cross-sectional area that the intersection paths are evenly spread out over is $\Delta t_{1}\Delta t_{2}/\sin \theta$. The intersection path density is ${\frac {|\mathbf {F} ||\mathbf {G} |\Delta t_{1}\Delta t_{2}}{\Delta t_{1}\Delta t_{2}/\sin \theta }}=|\mathbf {F} ||\mathbf {G} |\sin \theta =|\mathbf {F} \times \mathbf {G} |$. Lastly, it should be noted that the intersection paths are oriented out of the screen as per the right-hand rule.
Given a multi-surface with vector field $\mathbf {F} _{1}$, and a multi-surface with vector field $\mathbf {F} _{2}$, then the intersection is the multi-path with vector field $\mathbf {F} _{1}\times \mathbf {F} _{2}$.
The total intersection between multi-surface $\mathbf {F} _{1}$ and multi-surface $\mathbf {F} _{2}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {F} _{1}(\mathbf {q} )\times \mathbf {F} _{2}(\mathbf {q} ))dV$.
If $\mathbf {F} _{2}$ denotes a simple surface $\sigma$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {F} _{1}(\mathbf {q} )\times \delta _{2}(\mathbf {q} ;\sigma ))dV=\iint _{\mathbf {q} \in \sigma }\mathbf {F} _{1}(\mathbf {q} )\times d\mathbf {S}$.

### Surface-Volume intersections

When a surface $\sigma$ with weight $w_{1}$ intersects a volume $\Omega$ with weight $w_{2}$, then the intersection is surface $\sigma \cap \Omega$ with weight $w_{1}w_{2}$. Given a multi-surface and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple surface with each simple volume. The image below gives an example of the intersection of a multi-surface with a multi-volume.
Given a multi-surface with vector field $\mathbf {F}$, and a multi-volume with scalar field $U$, then the intersection is a multi-surface with vector field $\mathbf {F} U$.
The total intersection between a multi-surface $\mathbf {F}$ and a multi-volume $U$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )U(\mathbf {q} )dV$.
If $\mathbf {F}$ denotes a simple surface $\sigma$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{2}(\mathbf {q} ;\sigma )U(\mathbf {q} )dV=\iint _{\mathbf {q} \in \sigma }U(\mathbf {q} )d\mathbf {S}$.
If $U$ denotes a simple volume $\Omega$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\delta _{3}(\mathbf {q} ;\Omega )dV=\iiint _{\mathbf {q} \in \Omega }\mathbf {F} (\mathbf {q} )dV$.

### Volume-Volume intersections

When a volume $\Omega _{1}$ with weight $w_{1}$ intersects a volume $\Omega _{2}$ with weight $w_{2}$, then the intersection is the volume $\Omega _{1}\cap \Omega _{2}$ with weight $w_{1}w_{2}$. Given two multi-volumes, the intersection is the sum of the pair-wise intersections of each simple volume from the first multi-volume with each simple volume from the second multi-volume. The image below gives an example of the intersection between two multi-volumes.
Given a multi-volume with scalar field $U_{1}$, and a multi-volume with scalar field $U_{2}$, then the intersection is a multi-volume with scalar field $U_{1}U_{2}$.
The total intersection between multi-volume $U_{1}$ and multi-volume $U_{2}$ is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U_{1}(\mathbf {q} )U_{2}(\mathbf {q} )dV$.
If $U_{2}$ denotes a simple volume $\Omega$, then the total intersection is $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U_{1}(\mathbf {q} )\delta _{3}(\mathbf {q} ;\Omega )dV=\iiint _{\mathbf {q} \in \Omega }U_{1}(\mathbf {q} )dV$.

### Other intersections

Other types of intersections, such as Point-Point intersections, Point-Path intersections, Point-Surface intersections, and Path-Path intersections, are not considered since these kinds of intersections can occur only by design. For example, the probability that two randomly chosen points will intersect each other is 0, but if a point and a volume are randomly chosen, then the probability of the point landing in the volume is nonzero. Given two unrelated points, these two points will never land on each other, since a prior relationship has to exist for the points to coincide. Below is summarized the various types of intersections:
Intersections :
| structure | multi-point $\rho _{2}$ | multi-path $\mathbf {J} _{2}$ | multi-surface $\mathbf {F} _{2}$ | multi-volume $U_{2}$ |
| multi-point $\rho _{1}$ | n/a | n/a | n/a | multi-point $\rho _{1}U_{2}$ |
| multi-path $\mathbf {J} _{1}$ | n/a | n/a | multi-point $\mathbf {J} _{1}\cdot \mathbf {F} _{2}$ | multi-path $\mathbf {J} _{1}U_{2}$ |
| multi-surface $\mathbf {F} _{1}$ | n/a | multi-point $\mathbf {F} _{1}\cdot \mathbf {J} _{2}$ | multi-path $\mathbf {F} _{1}\times \mathbf {F} _{2}$ | multi-surface $\mathbf {F} _{1}U_{2}$ |
| multi-volume $U_{1}$ | multi-point $U_{1}\rho _{2}$ | multi-path $U_{1}\mathbf {J} _{2}$ | multi-surface $U_{1}\mathbf {F} _{2}$ | multi-volume $U_{1}U_{2}$ |

## Boundaries
### The endpoints of paths

Given a simple path $C$ that starts at point $C(0)$ and ends at point $C(1)$, the "endpoints" of $C$ is the multi-point $\{(C(0),+1),(C(1),-1)\}$ that consists of the starting point with a weight of +1, and the final point with a weight of -1. While $C$ is denoted by the vector field $\delta _{1}(\mathbf {q} ;C)$, the endpoints are denoted by the scalar field $\delta _{0}(\mathbf {q} ;C(0))-\delta _{0}(\mathbf {q} ;C(1))$. The image below gives several examples of simple paths and their associated endpoints.
Given a multi-path $\mathbf {C} =\{(C_{1},w_{1}),(C_{2},w_{2}),...,(C_{k},w_{k})\}$, the endpoints of $\mathbf {C}$ is the multi-point $\{(C_{1}(0),+1),(C_{1}(1),-1),(C_{2}(0),+1),(C_{2}(1),-1),...,(C_{k}(0),+1),(C_{k}(1),-1)\}$.
Given a vector field $\mathbf {J}$ that denotes a multi-path, the multi-point that denotes the endpoints of $\mathbf {J}$ is denoted by scalar field $\nabla \cdot \mathbf {J}$. The scalar field $\nabla \cdot \mathbf {J}$ evaluated at point $\mathbf {q}$ is denoted by $\nabla \cdot \mathbf {J} (\mathbf {q} )$, $(\nabla \cdot \mathbf {J} )(\mathbf {q} )$ or $\nabla \cdot \mathbf {J} |_{\mathbf {q} }$.
The requirement that no path extends to infinity means that every starting point is paired with a final point, and therefore the total weight of all of the endpoints together is 0: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \mathbf {J} (\mathbf {q} ))dV=0$.
The similarity of the notation $\nabla \cdot \mathbf {J}$ to the intersection of multi-path $\mathbf {J}$ with multi-surface $\mathbf {F}$, denoted by $\mathbf {F} \cdot \mathbf {J}$, makes sense if $\nabla$ is interpreted as the "surface of reality". A starting point forms when a path pokes into reality, and a final point forms when a path pokes out of reality.
In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. On the right is a simple path $\mathbf {F}$, along with its endpoints $\nabla \cdot \mathbf {F}$. On the left $\mathbf {F} _{\text{ext}}$ is an extension of $\mathbf {F}$ that is behind the "veil" of surface $\mathbf {G} _{\nabla }$. $\mathbf {F} _{\text{ext}}$ pokes out of and into $\mathbf {G} _{\nabla }$ at points consistent with the endpoints of $\mathbf {F}$: i.e. $\mathbf {G} _{\nabla }\cdot \mathbf {F} _{\text{ext}}=\nabla \cdot \mathbf {F}$.

### The counter-clockwise oriented boundaries of surfaces

Given an oriented surface $\sigma$, the "counter-clockwise oriented boundary" of $\sigma$ is a path $\partial \sigma$ that traces the boundary of $\sigma$ in a counter-clockwise direction. The counter-clockwise direction is better described as follows: While located on the boundary, the counter-clockwise direction is the "forwards" direction when the surface normal vectors point "up" and the surface itself is on the "left". The image below gives several examples of oriented surfaces and their counter-clockwise boundaries. Note in particular the 4th panel that shows that the orientation around a hole in the surface appears to be clockwise.
Given a multi-surface $\mathbf {S} =\{(\sigma _{1},w_{1}),(\sigma _{2},w_{2}),...,(\sigma _{k},w_{k})\}$, the counter-clockwise boundary of $\mathbf {S}$ is the multi-path $\{(\partial \sigma _{1},w_{1}),(\partial \sigma _{2},w_{2}),...,(\partial \sigma _{k},w_{k})\}$.
Given a vector field $\mathbf {F}$ that denotes a multi-surface, the multi-path that denotes the counter-clockwise boundary of $\mathbf {F}$ is denoted by vector field $\nabla \times \mathbf {F}$. The vector field $\nabla \times \mathbf {F}$ evaluated at point $\mathbf {q}$ is denoted by $\nabla \times \mathbf {F} (\mathbf {q} )$, $(\nabla \times \mathbf {F} )(\mathbf {q} )$, or $\nabla \times \mathbf {F} |_{\mathbf {q} }$.
The requirement that no surface weight extends to infinity means that all counter-clockwise boundaries form closed loops, and therefore the total displacement of the total counter-clockwise boundary is $\mathbf {0}$: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} (\mathbf {q} ))dV=\mathbf {0}$.
It is also important to note that the counter-clockwise boundary has no endpoints: $\nabla \cdot (\nabla \times \mathbf {F} )=0$.
The similarity of the notation $\nabla \times \mathbf {F} _{2}$ to the intersection of multi-surface $\mathbf {F} _{1}$ with multi-surface $\mathbf {F} _{2}$, denoted by $\mathbf {F} _{1}\times \mathbf {F} _{2}$, again makes sense if $\nabla$ is interpreted as the "surface of reality". An edge is formed when a surface "slices" into reality.
In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. On the right is a simple surface $\mathbf {F}$, along with its counter-clockwise boundary $\nabla \times \mathbf {F}$. On the left $\mathbf {F} _{\text{ext}}$ is an extension of $\mathbf {F}$ that is behind the "veil" of surface $\mathbf {G} _{\nabla }$. $\mathbf {F} _{\text{ext}}$ slices into $\mathbf {G} _{\nabla }$ at curves consistent with the boundary of $\mathbf {F}$: i.e. $\mathbf {G} _{\nabla }\times \mathbf {F} _{\text{ext}}=\nabla \times \mathbf {F}$.

### The inwards-oriented surfaces of volumes

Given a volume $\Omega$, the "inwards oriented surface" of $\Omega$ is a surface $\partial \Omega$ that wraps the volume $\Omega$ with the surface normals pointing inwards. The image below gives several examples of volumes and their inwards oriented surfaces.
Given a multi-volume $\mathbf {U} =\{(\Omega _{1},w_{1}),(\Omega _{2},w_{2}),...,(\Omega _{k},w_{k})\}$, the inwards oriented surface of $\mathbf {U}$ is the multi-surface $\{(\partial \Omega _{1},w_{1}),(\partial \Omega _{2},w_{2}),...,(\partial \Omega _{k},w_{k})\}$.
Given a scalar field $U$ that denotes a multi-volume, the multi-surface that denotes the inwards oriented surface of $U$ is denoted by vector field $\nabla U$. The vector field $\nabla U$ evaluated at point $\mathbf {q}$ is denoted by $\nabla U(\mathbf {q} )$, $(\nabla U)(\mathbf {q} )$, or $\nabla U|_{\mathbf {q} }$.
The requirement that no volume weight extends to infinity means that all inwards oriented surfaces form closed surfaces, and therefore the total surface vector of the total inwards oriented surface is $\mathbf {0}$: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla U(\mathbf {q} ))dV=\mathbf {0}$.
It is also important to note that the inwards oriented surface has no boundary: $\nabla \times (\nabla U)=\mathbf {0}$.
The similarity of the notation $\nabla U$ to the intersection of multi-surface $\mathbf {F}$ with multi-volume $U$, denoted by $\mathbf {F} U$, again makes sense if $\nabla$ is interpreted as the "surface of reality". A surface is formed from the surface of reality when the volume "pushes" into reality.
In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. The image is a 2D cross-section for simplicity. On the right is a simple volume $U$, along with its inwards oriented surface $\nabla U$. On the left $U_{\text{ext}}$ is an extension of $U$ that is behind the "veil" of surface $\mathbf {G} _{\nabla }$. $U_{\text{ext}}$ pushes through $\mathbf {G} _{\nabla }$ at surfaces consistent with the surface of $U$: i.e. $\mathbf {G} _{\nabla }U_{\text{ext}}=\nabla U$.

### Closed loops and closed surfaces

A simple path is "closed" or a "loop" if its starting and final points are the same, so the total endpoints is 0 since the weights of the starting and final points cancel out. More generally, a multi-path $\mathbf {J}$ is "closed" or a "multi-loop" if $\nabla \cdot \mathbf {J} =0$. As previously noted, the counter-clockwise boundary of a surface is closed.
A simple surface is "closed" or a "bubble" if it has no boundary. More generally, a multi-surface $\mathbf {F}$ is "closed" or a "multi-bubble" if $\nabla \times \mathbf {F} =\mathbf {0}$. As previously noted, the inwards oriented surface of a volume is closed.
It is clear that the total displacement present in a closed multi-path is $\mathbf {0}$: $\nabla \cdot \mathbf {J} =0\implies \iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} dV=\mathbf {0}$, and it is also clear that the total surface vector of a closed multi-surface is also $\mathbf {0}$: $\nabla \times \mathbf {F} =\mathbf {0} \implies \iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} dV=\mathbf {0}$.
Given both a simple loop and a simple bubble, the total point weight of all intersection points is 0: every time the loop enters the bubble, it must also leave the bubble, and the weights of these two intersection points cancel out. More generally, given a closed multi-path $\mathbf {J}$ and a closed multi-surface $\mathbf {F}$, then the total intersection point weight is 0: $(\nabla \cdot \mathbf {J} =0\;{\text{and}}\;\nabla \times \mathbf {F} =\mathbf {0} )\implies \iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} \cdot \mathbf {F} )dV=0$.
The above identity gives rise to the following observations:
* The total intersection point weight of a multi-loop and a multi-surface is purely a function of the multi-loop and the multi-surface's counter-clockwise boundary: the interior of the multi-surface does not matter. If $\nabla \cdot \mathbf {J} =0$ and $\nabla \times \mathbf {F} _{1}=\nabla \times \mathbf {F} _{2}$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} \cdot \mathbf {F} _{1})dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} \cdot \mathbf {F} _{2})dV$.
* The total intersection point weight of a multi-path and a multi-bubble is purely a function of the multi-bubble and the multi-path's endpoints: the interior of the multi-path does not matter. If $\nabla \times \mathbf {F} =\mathbf {0}$ and $\nabla \cdot \mathbf {J} _{1}=\nabla \cdot \mathbf {J} _{2}$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} _{1}\cdot \mathbf {F} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\mathbf {J} _{2}\cdot \mathbf {F} )dV$.
The inwards oriented surface of a volume is closed. Conversely, given a closed surface, there exists a volume that "fills" the surface. More generally, given a multi-bubble $\mathbf {F}$, there exists a multi-volume $U$ for which $\mathbf {F}$ is the inwards oriented multi-surface of $U$: $\nabla \times \mathbf {F} =\mathbf {0} \implies \exists U:\nabla U=\mathbf {F}$. This multi-volume is referred to as the "scalar potential" of $\mathbf {F}$. The requirement that volumes cannot extend to infinity means that $U$ is unique.
The counter-clockwise oriented boundary of a surface is closed. Conversely, given a loop, there exists a surface that "fills" the loop. More generally, given a multi-loop $\mathbf {J}$, there exists a multi-surface $\mathbf {F}$ for which $\mathbf {J}$ is the counter-clockwise boundary of $\mathbf {F}$: $\nabla \cdot \mathbf {J} =0\implies \exists \mathbf {F} :\nabla \times \mathbf {F} =\mathbf {J}$. This multi-surface is referred to as the "vector potential" of $\mathbf {J}$. Even with the requirement that surfaces cannot extend to infinity, $\mathbf {F}$ is not unique without additional restrictions.

## Coordinate Systems

This section will describe how to compute various quantities such as intersections, endpoints, boundaries, and surfaces given a curvilinear coordinate system.
Let the curvilinear coordinate system be arbitrary. Let the 3 coordinates that index all points be $c_{1},c_{2},c_{3}$. Coordinates will be denoted by the triple $(c_{1},c_{2},c_{3})$.
The following notation will be used in the following discussions:
* Given an arbitrary expression $f:\{1,2,3\}\to \mathbb {R}$ that assigns a real number to each index $i=1,2,3$, then $(i;f(i))$ will denote the triple $(f(1),f(2),f(3))$.
* Given index variables $i,j\in \{1,2,3\}$, the expression $\mathbf {1} (i=j)$ equals 1 if $i=j$ and 0 if otherwise.
* Given an arbitrary expression $f:\{1,2,3\}\to \mathbb {R}$ that assigns a real number to each index $i=1,2,3$, then $\sum _{i}f(i)$ will denote the sum $f(1)+f(2)+f(3)$.
* Given an index variable $i\in \{1,2,3\}$, $i+1$ will rotate $i$ forwards by 1, and $i+2$ will rotate $i$ forwards by 2. In essence, $i+1=\left\{{\begin{array}{cc}i+1&(i=1,2)\\1&(i=3)\end{array}}\right.$ and $i+2=\left\{{\begin{array}{cc}3&(i=1)\\i-1&(i=2,3)\end{array}}\right.$.
Start with an arbitrary coordinate $(c'_{1},c'_{2},c'_{3})=(j;c'_{j})$ and infinitesimal differences $\Delta c_{1}$, $\Delta c_{2}$, and $\Delta c_{3}$. The following 3 paths, 3 surfaces, and volume will be associated with point $(j;c'_{j})$:
* For each $i\in \{1,2,3\}$ there exists an infinitely short path $C_{i}((j;c'_{j}))$ starting from point $(j;c'_{j})$ and ending on point $(j;c'_{j}+\Delta c_{i}\mathbf {1} (j=i))$ along the curve defined by $c'_{i}\leq c_{i}<c'_{i}+\Delta c_{i}$, $c_{i+1}=c'_{i+1}$ and $c_{i+2}=c'_{i+2}$. The displacement covered by $C_{i}((j;c'_{j}))$ is approximately $\Delta c_{i}\cdot l_{i}((j;c'_{j}))\cdot {\hat {\mathbf {a} }}_{i}((j;c'_{j}))$ where ${\hat {\mathbf {a} }}_{i}((j;c'_{j}))$ is a unit length vector that is parallel to the displacement between points $(j;c'_{j})$ and $(j;c'_{j}+\Delta c_{i}\mathbf {1} (j=i))$, and $\Delta c_{i}\cdot l_{i}((j;c'_{j}))$ is the length of the displacement. Note that the length of the displacement is proportional to $\Delta c_{i}$, with $l_{i}((j;c'_{j}))$ being the constant of proportionality. The set of vectors $\{{\hat {\mathbf {a} }}_{1}((j;c'_{j})),{\hat {\mathbf {a} }}_{2}((j;c'_{j})),{\hat {\mathbf {a} }}_{3}((j;c'_{j}))\}$ is the set of displacement basis vectors.
* For each $i\in \{1,2,3\}$ there exists an infinitely small surface $\sigma _{i}((j;c'_{j}))$ that is defined by the following: $c_{i}=c'_{i}$, $c'_{i+1}\leq c_{i+1}<c'_{i+1}+\Delta c_{i+1}$, and $c'_{i+2}\leq c_{i+2}<c'_{i+2}+\Delta c_{i+2}$. The orientation of $\sigma _{i}((j;c'_{j}))$ is in the direction of increasing $c_{i}$. The surface vector of $\sigma _{i}((j;c'_{j}))$ is approximately $\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}((j;c'_{j}))\cdot {\hat {\mathbf {a} }}^{i}((j;c'_{j}))$ where ${\hat {\mathbf {a} }}^{i}((j;c'_{j}))$ is a unit length vector that is perpendicular to $\sigma _{i}((j;c'_{j}))$, and $\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}((j;c'_{j}))$ is the area of $\sigma _{i}((j;c'_{j}))$. Note that the area of $\sigma _{i}((j;c'_{j}))$ is proportional to $\Delta c_{i+1}\Delta c_{i+2}$, with $A_{i}((j;c'_{j}))$ being the constant of proportionality. The set of vectors $\{{\hat {\mathbf {a} }}^{1}((j;c'_{j})),{\hat {\mathbf {a} }}^{2}((j;c'_{j})),{\hat {\mathbf {a} }}^{3}((j;c'_{j}))\}$ is the set of surface basis vectors.
* There is an infinitely small volume $\Omega ((j;c'_{j}))$ defined by $c'_{1}\leq c_{1}<c'_{1}+\Delta c_{1}$, $c'_{2}\leq c_{2}<c'_{2}+\Delta c_{2}$, and $c'_{3}\leq c_{3}<c'_{3}+\Delta c_{3}$. $\Omega ((j;c'_{j}))$ has a shape that is approximately that of a parallelepiped. The volume of $\Omega ((j;c'_{j}))$ is approximately $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V((j;c'_{j}))$. Note that the volume of $\Omega ((j;c'_{j}))$ is proportional to $\Delta c_{1}\Delta c_{2}\Delta c_{3}$, with $V((j;c'_{j}))$ being the constant of proportionality.
It is important to note that:
* $(i;c_{i})\in \Omega ((j;c'_{j}))$ if and only if $c'_{1}\leq c_{1}<c'_{1}+\Delta c_{1}$, $c'_{2}\leq c_{2}<c'_{2}+\Delta c_{2}$, and $c'_{3}\leq c_{3}<c'_{3}+\Delta c_{3}$ (note the strictness of the upper bounds).
* For all $i\in \{1,2,3\}$, $C_{i}((j;c_{j}))\subseteq \Omega ((j;c'_{j}))$ if and only if $c_{i}=c'_{i}$, $c'_{i+1}\leq c_{i+1}<c'_{i+1}+\Delta c_{i+1}$, and $c'_{i+2}\leq c_{i+2}<c'_{i+2}+\Delta c_{i+2}$ (note the strictness of the upper bounds).
* For all $i\in \{1,2,3\}$, $\sigma _{i}((j;c_{j}))\subseteq \Omega ((j;c'_{j}))$ if and only if $c'_{i}\leq c_{i}<c'_{i}+\Delta c_{i}$ (note the strictness of the upper bound), $c_{i+1}=c'_{i+1}$, and $c_{i+2}=c'_{i+2}$.
Converting between multi-points, multi-paths, multi-surfaces, and multi-volumes and their respective scalar fields and vector fields proceeds as follows:
This conversion is performed by subdividing space into discrete volumes or cells. Infinitesimal differences $\Delta c_{1}$, $\Delta c_{2}$, and $\Delta c_{3}$ are chosen, and a lattice consisting of the points $(j;k_{j}\Delta c_{j})$ where $(j;k_{j})$ is an arbitrary triple of integers is generated. The cell indexed by $(j;k_{j})$ consists of the point $(j;k_{j}\Delta c_{j})$, the paths $C_{i}((j;k_{j}\Delta c_{j}))$ for each $i\in \{1,2,3\}$, the surfaces $\sigma _{i}((j;k_{j}\Delta c_{j}))$ for each $i\in \{1,2,3\}$, and the volume $\Omega ((j;k_{j}\Delta c_{j}))$. All points $(j;c_{j})$ where $k_{i}\Delta c_{i}\leq c_{i}<(k_{i}+1)\Delta c_{i}$ for all $i\in \{1,2,3\}$ "belong" to the cell indexed by $(j;k_{j})$ (note that the upper bounds are excluded). Given an arbitrary point $(j;c_{j})$, the cell that contains $(j;c_{j})$ is indexed by $(j;k_{j})=\left(j;\left\lfloor {\frac {c_{j}}{\Delta c_{j}}}\right\rfloor \right)$. The point $(j;c'_{j})=(j;k_{j}\Delta c_{j})$ is the vertex that the cell is associated with.
A multi-point, multi-path, multi-surface, or multi-volume is converted to a scalar field or vector field by computing the total point weight, displacement, surface vector, or volume contained by each cell and then averaging over the cell's volume.
A scalar-field $\rho$ is converted to a multi-point by doing the following for each cell $(j;k_{j})$. First compute the total point weight contained inside the cell: $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}\rho (\mathbf {q} )dV\approx \rho ((j;k_{j}\Delta c_{j}))V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next assign this weight to the point $(j;k_{j}\Delta c_{j})$.
A vector-field $\mathbf {J} =\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}$ is converted to a multi-path by doing the following for each cell $(j;k_{j})$. First compute the total displacement contained inside the cell: $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}\mathbf {J} (\mathbf {q} )dV\approx \left(\sum _{i}J_{i}((j;k_{j}\Delta c_{j})){\hat {\mathbf {a} }}_{i}((j;k_{j}\Delta c_{j}))\right)V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next separate this total displacement into components according to the basis ${\hat {\mathbf {a} }}_{1}$, ${\hat {\mathbf {a} }}_{2}$, and ${\hat {\mathbf {a} }}_{3}$: for each $i\in \{1,2,3\}$ the coefficient of ${\hat {\mathbf {a} }}_{i}$ is $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}J_{i}(\mathbf {q} )dV\approx J_{i}((j;k_{j}\Delta c_{j}))V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next for each $i\in \{1,2,3\}$, divide the coefficient of ${\hat {\mathbf {a} }}_{i}$ by the length of $C_{i}((j;k_{j}\Delta c_{j}))$, which results in approximately $J_{i}((j;k_{j}\Delta c_{j})){\frac {V((j;k_{j}\Delta c_{j}))}{l_{i}((j;k_{j}\Delta c_{j}))}}\Delta c_{i+1}\Delta c_{i+2}$, and assign this weight to $C_{i}((j;k_{j}\Delta c_{j}))$.
A vector-field $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$ is converted to a multi-surface by doing the following for each cell $(j;k_{j})$. First compute the total surface vector contained inside the cell: $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}\mathbf {F} (\mathbf {q} )dV\approx \left(\sum _{i}F_{i}((j;k_{j}\Delta c_{j})){\hat {\mathbf {a} }}^{i}((j;k_{j}\Delta c_{j}))\right)V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next separate this total surface vector into components according to the basis ${\hat {\mathbf {a} }}^{1}$, ${\hat {\mathbf {a} }}^{2}$, and ${\hat {\mathbf {a} }}^{3}$: for each $i\in \{1,2,3\}$ the coefficient of ${\hat {\mathbf {a} }}^{i}$ is $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}F_{i}(\mathbf {q} )dV\approx F_{i}((j;k_{j}\Delta c_{j}))V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next for each $i\in \{1,2,3\}$, divide the coefficient of ${\hat {\mathbf {a} }}^{i}$ by the area of $\sigma _{i}((j;k_{j}\Delta c_{j}))$, which results in approximately $F_{i}((j;k_{j}\Delta c_{j})){\frac {V((j;k_{j}\Delta c_{j}))}{A_{i}((j;k_{j}\Delta c_{j}))}}\Delta c_{i}$, and assign this weight to $\sigma _{i}((j;k_{j}\Delta c_{j}))$.
A scalar-field $U$ is converted to a multi-volume by doing the following for each cell $(j;k_{j})$. First compute the total volume contained inside the cell: $\iiint _{\mathbf {q} \in \Omega ((j;k_{j}\Delta c_{j}))}U(\mathbf {q} )dV\approx U((j;k_{j}\Delta c_{j}))V((j;k_{j}\Delta c_{j}))\Delta c_{1}\Delta c_{2}\Delta c_{3}$. Next divide this weight by the volume of $\Omega ((j;k_{j}\Delta c_{j}))$, which results in approximately $U((j;k_{j}\Delta c_{j}))$, and assign this weight to $\Omega ((j;k_{j}\Delta c_{j}))$.

### Computing various intersections

Computing the intersection of any structure with a multi-volume is trivial matter: Simply multiply the scalar of vector field by the scalar field that denotes the multi-volume. When both structures are denoted by vector fields however, computing the intersection is far less trivial.

#### Computing path-surface intersections

To save space, the notation $(j;c_{j})$ and $(j;k_{j}\Delta c_{j})$ will be omitted from the various terms.
Given a multi-path $\mathbf {C}$ denoted by vector field $\mathbf {J} =\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}$, and a multi-surface $\mathbf {S}$ denoted by vector field $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$, the scalar field that denotes the intersection can be computed as follows:
The following computations applies to each cell:
For each $i\in \{1,2,3\}$, the weight assigned to $C_{i}$ by $\mathbf {C}$ is computed as follows: $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot J_{i}$ is the ${\hat {\mathbf {a} }}_{i}$ component of the total displacement contained by the current cell. Computing the weight assigned to $C_{i}$ requires that this displacement be spread over the length of $C_{i}$: ${\frac {\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot J_{i}}{\Delta c_{i}\cdot l_{i}}}={\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}$.
For each $i\in \{1,2,3\}$, the weight assigned to $\sigma _{i}$ by $\mathbf {S}$ is computed as follows: $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}$ is the ${\hat {\mathbf {a} }}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $\sigma _{i}$ requires that this surface vector be spread over the area of $\sigma _{i}$: ${\frac {\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}}{\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}}}={\frac {V}{A_{i}}}\cdot \Delta c_{i}\cdot F_{i}$.
The intersection between $C_{i}$ and $\sigma _{i}$ is the current lattice point with weight $\left({\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}\right)\left({\frac {V}{A_{i}}}\cdot \Delta c_{i}\cdot F_{i}\right)$$={\frac {V^{2}}{l_{i}A_{i}}}\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot J_{i}F_{i}$.
Aside from the intersections between $C_{i}$ and $\sigma _{i}$ for each cell and $i\in \{1,2,3\}$, no other intersections occur. The total weight of the intersection at the vertex of the current cell is $\sum _{i}{\frac {V^{2}}{l_{i}A_{i}}}\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot J_{i}F_{i}$$=V^{2}\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}\sum _{i}{\frac {1}{l_{i}A_{i}}}\cdot J_{i}F_{i}$.
The value of $\mathbf {J} \cdot \mathbf {F}$ at the current cell is approximately ${\frac {1}{\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V}}\cdot V^{2}\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}\sum _{i}{\frac {1}{l_{i}A_{i}}}\cdot J_{i}F_{i}$$=V\sum _{i}{\frac {1}{l_{i}A_{i}}}\cdot J_{i}F_{i}$. The coefficient of ${\frac {1}{\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V}}$ exists to spread the point weight over the current cell.
Therefore $\mathbf {J} \cdot \mathbf {F} =V\sum _{i}{\frac {1}{l_{i}A_{i}}}\cdot J_{i}F_{i}$. Note that $\mathbf {J}$ is expressed using the displacement basis vectors, while $\mathbf {F}$ is expressed using the surface basis vectors.

#### Computing surface-surface intersections

To save space, the notation $(j;c_{j})$ and $(j;k_{j}\Delta c_{j})$ will be omitted from the various terms.
Given a multi-surface $\mathbf {S} _{1}$ denoted by vector field $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$, and a multi-surface $\mathbf {S} _{2}$ denoted by vector field $\mathbf {G} =\sum _{i}G_{i}{\hat {\mathbf {a} }}^{i}$, the vector field that denotes the intersection can be computed as follows:
The following computations applies to each cell:
For each $i\in \{1,2,3\}$, the weight assigned to $\sigma _{i}$ by $\mathbf {S} _{1}$ is computed as follows: $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}$ is the ${\hat {\mathbf {a} }}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $\sigma _{i}$ requires that this surface vector be spread over the area of $\sigma _{i}$: ${\frac {\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}}{\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}}}={\frac {V}{A_{i}}}\cdot \Delta c_{i}\cdot F_{i}$. Similarly, the weight assigned to $\sigma _{i}$ by $\mathbf {S} _{2}$ is ${\frac {V}{A_{i}}}\cdot \Delta c_{i}\cdot G_{i}$.
The intersection between $\sigma _{i+1}$ and $\sigma _{i+2}$ is path $C_{i}$ with weight $\left({\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}\right)\left({\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot G_{i+2}\right)$$={\frac {V^{2}}{A_{i+1}A_{i+2}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot F_{i+1}G_{i+2}$. Conversely, the intersection between $\sigma _{i+2}$ and $\sigma _{i+1}$ is path $C_{i}$ with weight $-{\frac {V^{2}}{A_{i+2}A_{i+1}}}\cdot \Delta c_{i+2}\Delta c_{i+1}\cdot F_{i+2}G_{i+1}$.
Aside from the intersections between $\sigma _{i+1}$ and $\sigma _{i+2}$, and the intersections between $\sigma _{i+2}$ and $\sigma _{i+1}$, for each cell and $i\in \{1,2,3\}$, no other intersections occur. For each $i\in \{1,2,3\}$, the total weight assigned to $C_{i}$ is ${\frac {V^{2}}{A_{i+1}A_{i+2}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot F_{i+1}G_{i+2}-{\frac {V^{2}}{A_{i+2}A_{i+1}}}\cdot \Delta c_{i+2}\Delta c_{i+1}\cdot F_{i+2}G_{i+1}$ $={\frac {V^{2}}{A_{i+1}A_{i+2}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot (F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$.
The value of $\mathbf {F} \times \mathbf {G}$ at the current cell is approximately $\sum _{i}{\frac {l_{i}\cdot \Delta c_{i}\cdot {\hat {\mathbf {a} }}_{i}}{V\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}}}\cdot {\frac {V^{2}}{A_{i+1}A_{i+2}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot (F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$ $=\sum _{i}{\frac {Vl_{i}}{A_{i+1}A_{i+2}}}\cdot (F_{i+1}G_{i+2}-F_{i+2}G_{i+1}){\hat {\mathbf {a} }}_{i}$. The coefficient of ${\frac {l_{i}\cdot \Delta c_{i}\cdot {\hat {\mathbf {a} }}_{i}}{V\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}}}$ exists to spread the displacement of each path over the current cell.
Therefore $\mathbf {F} \times \mathbf {G} =\sum _{i}{\frac {Vl_{i}}{A_{i+1}A_{i+2}}}\cdot (F_{i+1}G_{i+2}-F_{i+2}G_{i+1}){\hat {\mathbf {a} }}_{i}$. Note that both $\mathbf {F}$ and $\mathbf {G}$ are both expressed using the surface basis vectors, but $\mathbf {F} \times \mathbf {G}$ is using the displacement basis vectors.

### Computing the endpoints of paths

To save space, the notation $(j;k_{j})$, $(j;c_{j})$, and $(j;k_{j}\Delta c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i\in \{1,2,3\}$, the notation $[Q]_{-i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$. This cell will be referred to as the $-i$ neighbor of the current cell.
Given a multi-path $\mathbf {C}$ denoted by vector field $\mathbf {J} =\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}$, the scalar field that denotes the endpoints can be computed as follows:
The following computations apply to each cell:
For each $i\in \{1,2,3\}$, the weight assigned to $C_{i}$ by $\mathbf {C}$ is computed as follows: $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot J_{i}$ is the ${\hat {\mathbf {a} }}_{i}$ component of the total displacement contained by the current cell. Computing the weight assigned to $C_{i}$ requires that this displacement be spread over the length of $C_{i}$: ${\frac {\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot J_{i}}{\Delta c_{i}\cdot l_{i}}}={\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}$.
For each $i\in \{1,2,3\}$, path $C_{i}$ contributes a weight of $+{\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}$ to the lattice point of the current cell, and path $[C_{i}]_{-i}$ contributes a weight of $-\left[{\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}\right]_{-i}$ to the lattice point of the current cell.
The total weight assigned to the lattice point of the current cell is $\sum _{i}\left(+{\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}-\left[{\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}\right]_{-i}\right)$ $\approx \sum _{i}\Delta c_{i}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}\cdot \Delta c_{i+1}\Delta c_{i+2}\cdot J_{i}\right)$ $=\Delta c_{1}\Delta c_{2}\Delta c_{3}\sum _{i}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}\cdot J_{i}\right)$.
Spreading the weight assigned to the current lattice point over the volume of the current cell gives: $\nabla \cdot \mathbf {J} ={\frac {1}{\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V}}\cdot \Delta c_{1}\Delta c_{2}\Delta c_{3}\sum _{i}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}\cdot J_{i}\right)$ $={\frac {1}{V}}\sum _{i}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}\cdot J_{i}\right)$.
Therefore: $\nabla \cdot \mathbf {J} ={\frac {1}{V}}\sum _{i}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}\cdot J_{i}\right)$. Note that $\mathbf {J}$ is expressed using the displacement basis vectors.

### Computing the counter-clockwise boundaries of surfaces

To save space, the notation $(j;k_{j})$, $(j;c_{j})$, and $(j;k_{j}\Delta c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i\in \{1,2,3\}$, the notation $[Q]_{-i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$. This cell will be referred to as the $-i$ neighbor of the current cell.
Given a multi-surface $\mathbf {S}$ denoted by vector field $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$, the vector field that denotes the counter-clockwise boundary can be computed as follows:
The following computations apply to each cell:
For each $i\in \{1,2,3\}$, the weight assigned to $\sigma _{i}$ by $\mathbf {S}$ is computed as follows: $\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}$ is the ${\hat {\mathbf {a} }}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $\sigma _{i}$ requires that this surface vector be spread over the area of $\sigma _{i}$: ${\frac {\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V\cdot F_{i}}{\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}}}={\frac {V}{A_{i}}}\cdot \Delta c_{i}\cdot F_{i}$.
For each $i\in \{1,2,3\}$, surfaces that contain path $C_{i}$ as part of their boundary include $\sigma _{i+1}$, $[\sigma _{i+1}]_{-(i+2)}$, $\sigma _{i+2}$, and $[\sigma _{i+2}]_{-(i+1)}$. $C_{i}$ receives a mass of $-{\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}$ from $\sigma _{i+1}$; a mass of $+\left[{\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}\right]_{-(i+2)}$ from $[\sigma _{i+1}]_{-(i+2)}$; a mass of $+{\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot F_{i+2}$ from $\sigma _{i+2}$; and a mass of $-\left[{\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot F_{i+2}\right]_{-(i+1)}$ from $[\sigma _{i+2}]_{-(i+1)}$. The total mass assigned to $C_{i}$ is $-{\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}+\left[{\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}\right]_{-(i+2)}+{\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot F_{i+2}-\left[{\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot F_{i+2}\right]_{-(i+1)}$ $\approx -\Delta c_{i+2}{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}\cdot \Delta c_{i+1}\cdot F_{i+1}\right)+\Delta c_{i+1}\cdot {\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}\cdot \Delta c_{i+2}\cdot F_{i+2}\right)$ $=\Delta c_{i+1}\Delta c_{i+2}\left({\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}\cdot F_{i+2}\right)-{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}\cdot F_{i+1}\right)\right)$.
Spreading the displacement generated by each $C_{i}$ over the volume of the current cell gives: $\nabla \times \mathbf {F} =\sum _{i}{\frac {\Delta c_{i}\cdot l_{i}\cdot {\hat {\mathbf {a} }}_{i}}{\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V}}\cdot \Delta c_{i+1}\Delta c_{i+2}\left({\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}\cdot F_{i+2}\right)-{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}\cdot F_{i+1}\right)\right)$ $=\sum _{i}{\frac {l_{i}}{V}}\left({\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}\cdot F_{i+2}\right)-{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}\cdot F_{i+1}\right)\right){\hat {\mathbf {a} }}_{i}$.
Therefore: $\nabla \times \mathbf {F} =\sum _{i}{\frac {l_{i}}{V}}\left({\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}\cdot F_{i+2}\right)-{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}\cdot F_{i+1}\right)\right){\hat {\mathbf {a} }}_{i}$. Note that $\mathbf {F}$ is expressed using the surface basis vectors, but $\nabla \times \mathbf {F}$ is using the displacement basis vectors.

### Computing the inwards-oriented surfaces of volumes

To save space, the notation $(j;k_{j})$, $(j;c_{j})$, and $(j;k_{j}\Delta c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i\in \{1,2,3\}$, the notation $[Q]_{-i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$. This cell will be referred to as the $-i$ neighbor of the current cell.
Given a multi-volume $\mathbf {U}$ denoted by scalar field $U$, the vector field that denotes the inwards-oriented surface can be computed as follows:
The following computations apply to each cell:
The cell's volume $\Omega$ has the weight $U$.
For each $i\in \{1,2,3\}$, surface $\sigma _{i}$ receives a weight of $U$ from the current cell, and a weight of $-[U]_{-i}$ from the $-i$ neighbor of the current cell. The total weight is simply $U-[U]_{-i}\approx \Delta c_{i}{\frac {\partial U}{\partial c_{i}}}$. Spreading the surface vector generated by each $\sigma _{i}$ over the volume of the current cell gives: $\nabla U=\sum _{i}{\frac {\Delta c_{i+1}\Delta c_{i+2}\cdot A_{i}\cdot {\hat {\mathbf {a} }}^{i}}{\Delta c_{1}\Delta c_{2}\Delta c_{3}\cdot V}}\cdot \Delta c_{i}{\frac {\partial U}{\partial c_{i}}}$ $=\sum _{i}{\frac {A_{i}}{V}}{\frac {\partial U}{\partial c_{i}}}{\hat {\mathbf {a} }}^{i}$.
Therefore: $\nabla U=\sum _{i}{\frac {A_{i}}{V}}{\frac {\partial U}{\partial c_{i}}}{\hat {\mathbf {a} }}^{i}$. Note that $\nabla U$ uses the surface basis vectors.

### Summary

* Given multi-path $\mathbf {J} =\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}$ and multi-surface $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$, the intersection of $\mathbf {J}$ with $\mathbf {F}$ is multi-point $\mathbf {J} \cdot \mathbf {F} =\sum _{i}{\frac {V}{l_{i}A_{i}}}J_{i}F_{i}$.
* Given multi-surfaces $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$ and $\mathbf {G} =\sum _{i}G_{i}{\hat {\mathbf {a} }}^{i}$, the intersection of $\mathbf {F}$ with $\mathbf {G}$ is multi-path $\mathbf {F} \times \mathbf {G} =\sum _{i}{\frac {l_{i}V}{A_{i+1}A_{i+2}}}(F_{i+1}G_{i+2}-F_{i+2}G_{i+1}){\hat {\mathbf {a} }}_{i}$.
* Given multi-path $\mathbf {J} =\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}$, the endpoints of $\mathbf {J}$ is multi-point $\nabla \cdot \mathbf {J} =\sum _{i}{\frac {1}{V}}{\frac {\partial }{\partial c_{i}}}\left({\frac {V}{l_{i}}}J_{i}\right)$.
* Given multi-surface $\mathbf {F} =\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}$, the counter-clockwise boundary of $\mathbf {F}$ is multi-path $\nabla \times \mathbf {F} =\sum _{i}{\frac {l_{i}}{V}}\left({\frac {\partial }{\partial c_{i+1}}}\left({\frac {V}{A_{i+2}}}F_{i+2}\right)-{\frac {\partial }{\partial c_{i+2}}}\left({\frac {V}{A_{i+1}}}F_{i+1}\right)\right){\hat {\mathbf {a} }}_{i}$.
* Given multi-volume $U$, the inwards-oriented surface of $U$ is multi-surface $\nabla U=\sum _{i}{\frac {A_{i}}{V}}{\frac {\partial U}{\partial c_{i}}}{\hat {\mathbf {a} }}^{i}$.

#### Orthogonal coordinate systems

In the special case where the displacement basis vectors $\{{\hat {\mathbf {a} }}_{1},{\hat {\mathbf {a} }}_{2},{\hat {\mathbf {a} }}_{3}\}$ are all mutually orthogonal (perpendicular), then:
* The surface basis vectors are identical to the displacement basis vectors: $\forall i\in \{1,2,3\}:{\hat {\mathbf {a} }}^{i}={\hat {\mathbf {a} }}_{i}$.
* For each $i\in \{1,2,3\}$, $A_{i}=l_{i+1}l_{i+2}$.
* $V=l_{1}l_{2}l_{3}$.
The above formulas simplify to:
* $\left(\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}\right)\cdot \left(\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}\right)=\sum _{i}J_{i}F_{i}$.
* $\left(\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}\right)\times \left(\sum _{i}G_{i}{\hat {\mathbf {a} }}^{i}\right)=\sum _{i}(F_{i+1}G_{i+2}-F_{i+2}G_{i+1}){\hat {\mathbf {a} }}_{i}$.
* $\nabla \cdot \left(\sum _{i}J_{i}{\hat {\mathbf {a} }}_{i}\right)=\sum _{i}{\frac {1}{l_{1}l_{2}l_{3}}}{\frac {\partial }{\partial c_{i}}}\left(l_{i+1}l_{i+2}J_{i}\right)$.
* $\nabla \times \left(\sum _{i}F_{i}{\hat {\mathbf {a} }}^{i}\right)=\sum _{i}{\frac {1}{l_{i+1}l_{i+2}}}\left({\frac {\partial }{\partial c_{i+1}}}(l_{i+2}F_{i+2})-{\frac {\partial }{\partial c_{i+2}}}\left(l_{i+1}F_{i+1}\right)\right){\hat {\mathbf {a} }}_{i}$.
* $\nabla U=\sum _{i}{\frac {1}{l_{i}}}{\frac {\partial U}{\partial c_{i}}}{\hat {\mathbf {a} }}^{i}$.
For Cartesian coordinates, $c_{1}=x$, $c_{2}=y$, $c_{3}=z$, and ${\hat {\mathbf {a} }}^{1}={\hat {\mathbf {a} }}_{1}={\hat {\mathbf {x} }}$, ${\hat {\mathbf {a} }}^{2}={\hat {\mathbf {a} }}_{2}={\hat {\mathbf {y} }}$, ${\hat {\mathbf {a} }}^{3}={\hat {\mathbf {a} }}_{3}={\hat {\mathbf {z} }}$, and $l_{1}=1$, $l_{2}=1$, $l_{3}=1$. Therefore:
* $(J_{x}{\hat {\mathbf {x} }}+J_{y}{\hat {\mathbf {y} }}+J_{z}{\hat {\mathbf {z} }})\cdot (F_{x}{\hat {\mathbf {x} }}+F_{y}{\hat {\mathbf {y} }}+F_{z}{\hat {\mathbf {z} }})=J_{x}F_{x}+J_{y}F_{y}+J_{z}F_{z}$.
* $(F_{x}{\hat {\mathbf {x} }}+F_{y}{\hat {\mathbf {y} }}+F_{z}{\hat {\mathbf {z} }})\times (G_{x}{\hat {\mathbf {x} }}+G_{y}{\hat {\mathbf {y} }}+G_{z}{\hat {\mathbf {z} }})=(F_{y}G_{z}-F_{z}G_{y}){\hat {\mathbf {x} }}+(F_{z}G_{x}-F_{x}G_{z}){\hat {\mathbf {y} }}+(F_{x}G_{y}-F_{y}G_{x}){\hat {\mathbf {z} }}$.
* $\nabla \cdot (J_{x}{\hat {\mathbf {x} }}+J_{y}{\hat {\mathbf {y} }}+J_{z}{\hat {\mathbf {z} }})={\frac {\partial J_{x}}{\partial x}}+{\frac {\partial J_{y}}{\partial y}}+{\frac {\partial J_{z}}{\partial z}}$.
* $\nabla \times (F_{x}{\hat {\mathbf {x} }}+F_{y}{\hat {\mathbf {y} }}+F_{z}{\hat {\mathbf {z} }})=\left({\frac {\partial F_{z}}{\partial y}}-{\frac {\partial F_{y}}{\partial z}}\right){\hat {\mathbf {x} }}+\left({\frac {\partial F_{x}}{\partial z}}-{\frac {\partial F_{z}}{\partial x}}\right){\hat {\mathbf {y} }}+\left({\frac {\partial F_{y}}{\partial x}}-{\frac {\partial F_{x}}{\partial y}}\right){\hat {\mathbf {z} }}$.
* $\nabla U={\frac {\partial U}{\partial x}}{\hat {\mathbf {x} }}+{\frac {\partial U}{\partial y}}{\hat {\mathbf {y} }}+{\frac {\partial U}{\partial z}}{\hat {\mathbf {z} }}$.
For cylindrical coordinates, $c_{1}=\rho$, $c_{2}=\phi$, $c_{3}=z$, and ${\hat {\mathbf {a} }}^{1}={\hat {\mathbf {a} }}_{1}={\hat {\mathbf {\rho } }}$, ${\hat {\mathbf {a} }}^{2}={\hat {\mathbf {a} }}_{2}={\hat {\mathbf {\phi } }}$, ${\hat {\mathbf {a} }}^{3}={\hat {\mathbf {a} }}_{3}={\hat {\mathbf {z} }}$, and $l_{1}=1$, $l_{2}=\rho$, $l_{3}=1$. Therefore:
* $(J_{\rho }{\hat {\mathbf {\rho } }}+J_{\phi }{\hat {\mathbf {\phi } }}+J_{z}{\hat {\mathbf {z} }})\cdot (F_{\rho }{\hat {\mathbf {\rho } }}+F_{\phi }{\hat {\mathbf {\phi } }}+F_{z}{\hat {\mathbf {z} }})=J_{\rho }F_{\rho }+J_{\phi }F_{\phi }+J_{z}F_{z}$.
* $(F_{\rho }{\hat {\mathbf {\rho } }}+F_{\phi }{\hat {\mathbf {\phi } }}+F_{z}{\hat {\mathbf {z} }})\times (F_{\rho }{\hat {\mathbf {\rho } }}+F_{\phi }{\hat {\mathbf {\phi } }}+F_{z}{\hat {\mathbf {z} }})=(F_{\phi }G_{z}-F_{z}G_{\phi }){\hat {\mathbf {\rho } }}+(F_{z}G_{\rho }-F_{\rho }G_{z}){\hat {\mathbf {\phi } }}+(F_{\rho }G_{\phi }-F_{\phi }G_{\rho }){\hat {\mathbf {z} }}$.
* $\nabla \cdot (J_{\rho }{\hat {\mathbf {\rho } }}+J_{\phi }{\hat {\mathbf {\phi } }}+J_{z}{\hat {\mathbf {z} }})={\frac {1}{\rho }}\left({\frac {\partial }{\partial \rho }}(\rho J_{\rho })+{\frac {\partial F_{\phi }}{\partial \phi }}+{\frac {\partial }{\partial z}}(\rho F_{z})\right)$.
* $\nabla \times (F_{\rho }{\hat {\mathbf {\rho } }}+F_{\phi }{\hat {\mathbf {\phi } }}+F_{z}{\hat {\mathbf {z} }})={\frac {1}{\rho }}\left({\frac {\partial F_{z}}{\partial \phi }}-{\frac {\partial }{\partial z}}(\rho F_{\phi })\right){\hat {\mathbf {\rho } }}+\left({\frac {\partial F_{\rho }}{\partial z}}-{\frac {\partial F_{z}}{\partial \rho }}\right){\hat {\mathbf {\phi } }}+{\frac {1}{\rho }}\left({\frac {\partial }{\partial \rho }}(\rho F_{\phi })-{\frac {\partial F_{\rho }}{\partial \phi }}\right){\hat {\mathbf {z} }}$.
* $\nabla U={\frac {\partial U}{\partial \rho }}{\hat {\mathbf {\rho } }}+{\frac {1}{\rho }}{\frac {\partial U}{\partial \phi }}{\hat {\mathbf {\phi } }}+{\frac {\partial U}{\partial z}}{\hat {\mathbf {z} }}$.
For spherical coordinates, $c_{1}=r$, $c_{2}=\theta$, $c_{3}=\phi$, and ${\hat {\mathbf {a} }}^{1}={\hat {\mathbf {a} }}_{1}={\hat {\mathbf {r} }}$, ${\hat {\mathbf {a} }}^{2}={\hat {\mathbf {a} }}_{2}={\hat {\mathbf {\theta } }}$, ${\hat {\mathbf {a} }}^{3}={\hat {\mathbf {a} }}_{3}={\hat {\mathbf {\phi } }}$, and $l_{1}=1$, $l_{2}=r$, $l_{3}=r\sin \theta$. Therefore:
* $(J_{r}{\hat {\mathbf {r} }}+J_{\theta }{\hat {\mathbf {\theta } }}+J_{\phi }{\hat {\mathbf {\phi } }})\cdot (F_{r}{\hat {\mathbf {r} }}+F_{\theta }{\hat {\mathbf {\theta } }}+F_{\phi }{\hat {\mathbf {\phi } }})=J_{r}F_{r}+J_{\theta }F_{\theta }+J_{\phi }F_{\phi }$.
* $(F_{r}{\hat {\mathbf {r} }}+F_{\theta }{\hat {\mathbf {\theta } }}+F_{\phi }{\hat {\mathbf {\phi } }})\times (G_{r}{\hat {\mathbf {r} }}+G_{\theta }{\hat {\mathbf {\theta } }}+G_{\phi }{\hat {\mathbf {\phi } }})=(F_{\theta }G_{\phi }-F_{\phi }G_{\theta }){\hat {\mathbf {r} }}+(F_{\phi }G_{r}-F_{r}G_{\phi }){\hat {\mathbf {\theta } }}+(F_{r}G_{\theta }-F_{\theta }G_{r}){\hat {\mathbf {\phi } }}$.
* $\nabla \cdot (J_{r}{\hat {\mathbf {r} }}+J_{\theta }{\hat {\mathbf {\theta } }}+J_{\phi }{\hat {\mathbf {\phi } }})={\frac {1}{r^{2}\sin \theta }}\left({\frac {\partial }{\partial r}}(r^{2}\sin \theta J_{r})+{\frac {\partial }{\partial \theta }}(r\sin \theta J_{\theta })+{\frac {\partial }{\partial \phi }}(rF_{\phi })\right)$.
* $\nabla \times (F_{r}{\hat {\mathbf {r} }}+F_{\theta }{\hat {\mathbf {\theta } }}+F_{\phi }{\hat {\mathbf {\phi } }})={\frac {1}{r^{2}\sin \theta }}\left({\frac {\partial }{\partial \theta }}(r\sin \theta F_{\phi })-{\frac {\partial }{\partial \phi }}(rF_{\theta })\right){\hat {\mathbf {r} }}+{\frac {1}{r\sin \theta }}\left({\frac {\partial F_{r}}{\partial \phi }}-{\frac {\partial }{\partial r}}(r\sin \theta F_{\phi })\right){\hat {\mathbf {\theta } }}+{\frac {1}{r}}\left({\frac {\partial }{\partial r}}(rF_{\theta })-{\frac {\partial F_{r}}{\partial \theta }}\right){\hat {\mathbf {\phi } }}$.
* $\nabla U={\frac {\partial U}{\partial r}}{\hat {\mathbf {r} }}+{\frac {1}{r}}{\frac {\partial U}{\partial \theta }}{\hat {\mathbf {\theta } }}+{\frac {1}{r\sin \theta }}{\frac {\partial U}{\partial \phi }}{\hat {\mathbf {\phi } }}$.

## Intersection Boundaries
### The endpoints of intersections

Many identities related to vector calculus can be derived from examining the endpoints of path-volume intersections and surface-surface intersections.

#### The endpoints of path-volume intersections

Start with a multi-path $\mathbf {C}$, denoted by vector field $\mathbf {J}$, and a multi-volume $\mathbf {U}$, denoted by scalar field $U$. The intersection $\mathbf {C} \cap \mathbf {U}$ is denoted by vector field $\mathbf {J} U$.
Any time a path $C$ with weight $w_{1}$ starts in a volume $\Omega$ with weight $w_{2}$, the intersection $\mathbf {C} \cap \mathbf {U}$ gains an endpoint at the starting point of $C$ with weight $w_{1}w_{2}$. Any time a path $C$ with weight $w_{1}$ finishes in a volume $\Omega$ with weight $w_{2}$, the intersection $\mathbf {C} \cap \mathbf {U}$ gains an endpoint at the finishing point of $C$ with weight $-w_{1}w_{2}$. The endpoints for $\mathbf {C} \cap \mathbf {U}$ that are generated when paths from $\mathbf {C}$ start or finish in volumes from $\mathbf {U}$ is the intersection of the endpoints of $\mathbf {C}$ with multi-volume $\mathbf {U}$. This contributes the term $(\nabla \cdot \mathbf {J} )U$ to $\nabla \cdot (\mathbf {J} U)$.
Any time a path $C$ with weight $w_{1}$ enters a volume $\Omega$ with weight $w_{2}$, the intersection $\mathbf {C} \cap \mathbf {U}$ gains an endpoint at the point of entry with weight $w_{1}w_{2}$. Any time a path $C$ with weight $w_{1}$ leaves a volume $\Omega$ with weight $w_{2}$, the intersection $\mathbf {C} \cap \mathbf {U}$ gains an endpoint at the point of exit with weight $-w_{1}w_{2}$. The endpoints for $\mathbf {C} \cap \mathbf {U}$ that are generated when paths from $\mathbf {C}$ enter or exit volumes from $\mathbf {U}$ is the intersection of multi-path $\mathbf {C}$ with the inwards oriented multi-surface of $\mathbf {U}$. This contributes the term $\mathbf {J} \cdot (\nabla U)$ to $\nabla \cdot (\mathbf {J} U)$.
The total endpoints of $\mathbf {C} \cap \mathbf {U}$ are: $\nabla \cdot (\mathbf {J} U)=(\nabla \cdot \mathbf {J} )U+\mathbf {J} \cdot (\nabla U)$. In essence, the endpoints of $\mathbf {C} \cap \mathbf {U}$ are the endpoints of $\mathbf {C}$ that are contained in $\mathbf {U}$, plus the points at which paths from $\mathbf {C}$ enter or exit volumes from $\mathbf {U}$. This is depicted in the images on the right.
From the identity $\nabla \cdot (\mathbf {J} U)=(\nabla \cdot \mathbf {J} )U+\mathbf {J} \cdot (\nabla U)$, counting the total point weight gives: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \cdot (\mathbf {J} (\mathbf {q} )U(\mathbf {q} ))dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \mathbf {J} (\mathbf {q} ))U(\mathbf {q} )dV+\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\cdot (\nabla U(\mathbf {q} ))dV$. For the endpoints of a multi-path, every starting point must be paired with a finishing point so the total point weight of the endpoints of a multi-path is 0. $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \cdot (\mathbf {J} (\mathbf {q} )U(\mathbf {q} ))dV=0$ so hence $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \mathbf {J} (\mathbf {q} ))U(\mathbf {q} )dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\cdot (\nabla U(\mathbf {q} ))dV$. The total intersection between the endpoints of multi-path $\mathbf {C}$ and multi-volume $\mathbf {U}$ is the negative of the total intersection between $\mathbf {C}$ and the inwards oriented surface of $\mathbf {U}$.
If $\mathbf {J}$ denotes a simple path $C$ that starts at point $\mathbf {q} _{0}$ and ends at point $\mathbf {q} _{1}$, then the above integral identity becomes:
$\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \delta _{1}(\mathbf {q} ;C))U(\mathbf {q} )dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{1}(\mathbf {q} ;C)\cdot (\nabla U(\mathbf {q} ))dV$ $\iff \iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\delta _{0}(\mathbf {q} ;\mathbf {q} _{0})-\delta _{0}(\mathbf {q} ;\mathbf {q} _{1}))U(\mathbf {q} )dV=-\int _{\mathbf {q} \in C}(\nabla U(\mathbf {q} ))\cdot d\mathbf {q}$ $\iff U(\mathbf {q} _{0})-U(\mathbf {q} _{1})=-\int _{\mathbf {q} \in C}(\nabla U(\mathbf {q} ))\cdot d\mathbf {q}$ $\iff \int _{\mathbf {q} \in C}(\nabla U(\mathbf {q} ))\cdot d\mathbf {q} =U(\mathbf {q} _{1})-U(\mathbf {q} _{0})$ This is known as the gradient theorem.
If $U$ denotes a simple volume $\Omega$ with a outwards oriented surface $\sigma$, then the integral identity becomes:
$\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \mathbf {J} (\mathbf {q} ))\delta _{3}(\mathbf {q} ;\Omega )dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\cdot (\nabla \delta _{3}(\mathbf {q} ;\Omega ))dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla \cdot \mathbf {J} (\mathbf {q} ))dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\cdot (-\delta _{2}(\mathbf {q} ;\sigma ))dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla \cdot \mathbf {J} (\mathbf {q} ))dV=\iint _{\mathbf {q} \in \sigma }\mathbf {J} (\mathbf {q} )\cdot d\mathbf {S}$ This is known as Gauss's divergence theorem.
In summary:
* Given a multi-path denoted by vector field $\mathbf {J}$, and a multi-volume denoted by scalar field $U$, then the endpoints of the intersection are: $\nabla \cdot (\mathbf {J} U)=(\nabla \cdot \mathbf {J} )U+\mathbf {J} \cdot (\nabla U)$.
* Given a multi-path denoted by vector field $\mathbf {J}$, and a multi-volume denoted by scalar field $U$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \cdot \mathbf {J} (\mathbf {q} ))U(\mathbf {q} )dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {J} (\mathbf {q} )\cdot (\nabla U(\mathbf {q} ))dV$.
* Given a simple path $C$ that starts at point $\mathbf {q} _{0}$ and ends at point $\mathbf {q} _{1}$, and a multi-volume denoted by scalar field $U$, then $\int _{\mathbf {q} \in C}(\nabla U(\mathbf {q} ))\cdot d\mathbf {q} =U(\mathbf {q} _{1})-U(\mathbf {q} _{0})$. This is the gradient theorem.
* Given a multi-path denoted by vector field $\mathbf {J}$, and a simple volume $\Omega$ with outwards oriented surface $\sigma$, then $\iiint _{\mathbf {q} \in \Omega }(\nabla \cdot \mathbf {J} (\mathbf {q} ))dV=\iint _{\mathbf {q} \in \sigma }\mathbf {J} (\mathbf {q} )\cdot d\mathbf {S}$. This is Gauss's divergence theorem.

#### The endpoints of surface-surface intersections

Start with multi-surface $\mathbf {S} _{1}$, denoted by vector field $\mathbf {F} _{1}$, and a second multi-surface $\mathbf {S} _{2}$, denoted by vector field $\mathbf {F} _{2}$. The intersection $\mathbf {S} _{1}\cap \mathbf {S} _{2}$ is denoted by vector field $\mathbf {F} _{1}\times \mathbf {F} _{2}$.
Consider a surface $\sigma _{1}$ with weight $w_{1}$ from $\mathbf {S} _{1}$, and a surface $\sigma _{2}$ with weight $w_{2}$ from $\mathbf {S} _{2}$. Let $\partial \sigma _{1}$ denote the counter-clockwise boundary of $\sigma _{1}$, and let $\partial \sigma _{2}$ denote the counter-clockwise boundary of $\sigma _{2}$. There are 4 scenarios regarding the endpoints of $\sigma _{1}\cap \sigma _{2}$:
* When the $\partial \sigma _{1}$ intersects $\sigma _{2}$ in the preferred direction, the intersection point $\partial \sigma _{1}\cap \sigma _{2}$ has a weight of $+w_{1}w_{2}$, and an endpoint with weight $+w_{1}w_{2}$ (starting point) for $\sigma _{1}\cap \sigma _{2}$ forms at $\partial \sigma _{1}\cap \sigma _{2}$.
* When the $\partial \sigma _{1}$ intersects $\sigma _{2}$ in the opposite direction, the intersection point $\partial \sigma _{1}\cap \sigma _{2}$ has a weight of $-w_{1}w_{2}$, and an endpoint with weight $-w_{1}w_{2}$ (finishing point) for $\sigma _{1}\cap \sigma _{2}$ forms at $\partial \sigma _{1}\cap \sigma _{2}$.
* When the $\partial \sigma _{2}$ intersects $\sigma _{1}$ in the preferred direction, the intersection point $\sigma _{1}\cap \partial \sigma _{2}$ has a weight of $+w_{1}w_{2}$, and an endpoint with weight $-w_{1}w_{2}$ (finishing point) for $\sigma _{1}\cap \sigma _{2}$ forms at $\sigma _{1}\cap \partial \sigma _{2}$.
* When the $\partial \sigma _{2}$ intersects $\sigma _{1}$ in the opposite direction, the intersection point $\sigma _{1}\cap \partial \sigma _{2}$ has a weight of $-w_{1}w_{2}$, and an endpoint with weight $+w_{1}w_{2}$ (starting point) for $\sigma _{1}\cap \sigma _{2}$ forms at $\sigma _{1}\cap \partial \sigma _{2}$.
It can be seen that the intersection $\partial \sigma _{1}\cap \sigma _{2}$ forms endpoints for $\sigma _{1}\cap \sigma _{2}$ with the correct polarity, and that the intersection $\sigma _{1}\cap \partial \sigma _{2}$ forms endpoints for $\sigma _{1}\cap \sigma _{2}$ with the opposite polarity. This can be observed in the image of the right. This implies that the endpoints of $\mathbf {S} _{1}\cap \mathbf {S} _{2}$ are: $\nabla \cdot (\mathbf {F} _{1}\times \mathbf {F} _{2})=(\nabla \times \mathbf {F} _{1})\cdot \mathbf {F} _{2}-\mathbf {F} _{1}\cdot (\nabla \times \mathbf {F} _{2})$.
From the identity $\nabla \cdot (\mathbf {F} _{1}\times \mathbf {F} _{2})=(\nabla \times \mathbf {F} _{1})\cdot \mathbf {F} _{2}-\mathbf {F} _{1}\cdot (\nabla \times \mathbf {F} _{2})$, counting the total point weight gives: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \cdot (\mathbf {F} _{1}(\mathbf {q} )\times \mathbf {F} _{2}(\mathbf {q} ))dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot \mathbf {F} _{2}(\mathbf {q} )dV-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} _{1}(\mathbf {q} )\cdot (\nabla \times \mathbf {F} _{2}(\mathbf {q} ))dV$. For the endpoints of a multi-path, every starting point must be paired with a finishing point so the total point weight of the endpoints of a multi-path is 0. $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \cdot (\mathbf {F} _{1}(\mathbf {q} )\times \mathbf {F} _{2}(\mathbf {q} ))dV=0$ so hence $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot \mathbf {F} _{2}(\mathbf {q} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} _{1}(\mathbf {q} )\cdot (\nabla \times \mathbf {F} _{2}(\mathbf {q} ))dV$. The total intersection of the counter-clockwise boundary of multi-surface $\mathbf {S} _{1}$ with multi-surface $\mathbf {S} _{2}$ is the total intersection of the counter-clockwise boundary of $\mathbf {S} _{2}$ with $\mathbf {S} _{1}$. This is illustrated by the image on the right.
If $\mathbf {F} _{2}$ denotes a simple surface $\sigma$ with a counter-clockwise boundary $\partial \sigma$, then the above integral identity becomes:
$\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot \delta _{2}(\mathbf {q} ;\sigma )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} _{1}(\mathbf {q} )\cdot (\nabla \times \delta _{2}(\mathbf {q} ;\sigma ))dV$ $\iff \iint _{\mathbf {q} \in \sigma }(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot d\mathbf {S} =\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} _{1}(\mathbf {q} )\cdot \delta _{1}(\mathbf {q} ;\partial \sigma )dV$ $\iff \iint _{\mathbf {q} \in \sigma }(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot d\mathbf {S} =\int _{\mathbf {q} \in \partial \sigma }\mathbf {F} _{1}(\mathbf {q} )\cdot d\mathbf {q}$
This is known as Stokes' theorem.
In summary:
* Given two multi-surfaces denoted by vector fields $\mathbf {F} _{1}$ and $\mathbf {F} _{2}$, then the endpoints of the intersection are: $\nabla \cdot (\mathbf {F} _{1}\times \mathbf {F} _{2})=(\nabla \times \mathbf {F} _{1})\cdot \mathbf {F} _{2}-\mathbf {F} _{1}\cdot (\nabla \times \mathbf {F} _{2})$.
* Given two multi-surfaces denoted by vector fields $\mathbf {F} _{1}$ and $\mathbf {F} _{2}$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot \mathbf {F} _{2}(\mathbf {q} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} _{1}(\mathbf {q} )\cdot (\nabla \times \mathbf {F} _{2}(\mathbf {q} ))dV$.
* Given a multi-surface denoted by vector field $\mathbf {F} _{1}$ and a simple surface $\sigma$ with counter-clockwise oriented boundary $\partial \sigma$, then $\iint _{\mathbf {q} \in \sigma }(\nabla \times \mathbf {F} _{1}(\mathbf {q} ))\cdot d\mathbf {S} =\int _{\mathbf {q} \in \partial \sigma }\mathbf {F} _{1}(\mathbf {q} )\cdot d\mathbf {q}$. This is Stokes' theorem.

### The boundaries of intersections

In addition to the identities derived from examining the endpoints of intersections, some more identities can be derived by examining the counter-clockwise boundaries of surfaces that result from intersections.

#### The counter-clockwise boundary of surface-volume intersections

Start with a multi-surface $\mathbf {S}$, denoted by vector field $\mathbf {F}$, and a multi-volume $\mathbf {U}$, denoted by scalar field $U$. The intersection $\mathbf {S} \cap U$ is denoted by vector field $\mathbf {F} U$.
Consider a surface $\sigma$ with weight $w_{1}$ from $\mathbf {S}$, and a volume $\Omega$ with weight $w_{2}$ from $\mathbf {U}$. Let $\partial \sigma$ denote the counter-clockwise boundary of $\sigma$, and let $\partial \Omega$ denote the inwards oriented surface of $\Omega$. There are two sources for the counter-clockwise boundary of $\sigma \cap \Omega$. Any time $\partial \sigma$ intersects $\Omega$, the intersection $\partial \sigma \cap \Omega$ contributes to the boundary of $\sigma \cap \Omega$. When $\partial \sigma$ leaves $\Omega$, the boundary of $\sigma \cap \Omega$ cannot follow, and instead must trace along the surface of $\Omega$ while remaining in the surface $\sigma$ as indicated in the image to the right. The boundary of the total intersection $\mathbf {S} \cap \mathbf {U}$, denoted by $\nabla \times (\mathbf {F} U)$, consists of two parts: the intersection of the boundary of $\mathbf {S}$ with $\mathbf {U}$, denoted by $(\nabla \times \mathbf {F} )U$, and the intersection of the inwards-oriented surface of $\mathbf {U}$ with $\mathbf {S}$, denoted by $(\nabla U)\times \mathbf {F} =-\mathbf {F} \times (\nabla U)$. Therefore: $\nabla \times (\mathbf {F} U)=(\nabla \times \mathbf {F} )U+(\nabla U)\times \mathbf {F} =(\nabla \times \mathbf {F} )U-\mathbf {F} \times (\nabla U)$.
From the identity $\nabla \times (\mathbf {F} U)=(\nabla \times \mathbf {F} )U-\mathbf {F} \times (\nabla U)$, computing the total displacement gives: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \times (\mathbf {F} (\mathbf {q} )U(\mathbf {q} ))dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} (\mathbf {q} ))U(\mathbf {q} )dV-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\times (\nabla U(\mathbf {q} ))dV$. The counter-clockwise boundary of a multi-surface is a closed multi-loop, and the total displacement generated by a loop is $\mathbf {0}$. $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla \times (\mathbf {F} (\mathbf {q} )U(\mathbf {q} ))dV=\mathbf {0}$ so hence $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} (\mathbf {q} ))U(\mathbf {q} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\times (\nabla U(\mathbf {q} ))dV$. The total intersection of a the boundary of multi-surface $\mathbf {S}$ with multi-volume $\mathbf {U}$ is the total intersection of $\mathbf {S}$ with the surface of $\mathbf {U}$.
If $\mathbf {F}$ denotes a simple surface $\sigma$ with counter-clockwise boundary $\partial \sigma$, then the above integral identity becomes:
$\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \delta _{2}(\mathbf {q} ;\sigma ))U(\mathbf {q} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{2}(\mathbf {q} ;\sigma )\times (\nabla U(\mathbf {q} ))dV$ $\iff \iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{1}(\mathbf {q} ;\partial \sigma )U(\mathbf {q} )dV=\iint _{\mathbf {q} \in \sigma }d\mathbf {S} \times (\nabla U(\mathbf {q} ))$ $\iff \int _{\mathbf {q} \in \partial \sigma }U(\mathbf {q} )d\mathbf {q} =-\iint _{\mathbf {q} \in \sigma }(\nabla U(\mathbf {q} ))\times d\mathbf {S}$
If $\Omega$ denotes a simple volume with outwards oriented surface $\sigma$, then the integral identity becomes:
$\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} (\mathbf {q} ))\delta _{3}(\mathbf {q} ;\Omega )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\times (\nabla \delta _{3}(\mathbf {q} ;\Omega ))dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla \times \mathbf {F} (\mathbf {q} ))dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\times (-\delta _{2}(\mathbf {q} ;\sigma ))dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla \times \mathbf {F} (\mathbf {q} ))dV=-\iint _{\mathbf {q} \in \sigma }\mathbf {F} (\mathbf {q} )\times d\mathbf {S}$
In summary:
* Given a multi-surface denoted by vector field $\mathbf {F}$, and a multi-volume denoted by scalar field $U$, then the counter-clockwise boundary of the intersection is: $\nabla \times (\mathbf {F} U)=(\nabla \times \mathbf {F} )U-\mathbf {F} \times (\nabla U)$.
* Given a multi-surface denoted by vector field $\mathbf {F}$, and a multi-volume denoted by scalar field $U$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \times \mathbf {F} (\mathbf {q} ))U(\mathbf {q} )dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\mathbf {F} (\mathbf {q} )\times (\nabla U(\mathbf {q} ))dV$
* Given a simple surface $\sigma$ with counter-clockwise boundary $\partial \sigma$, and a multi-volume denoted by scalar field $U$, then $\int _{\mathbf {q} \in \partial \sigma }U(\mathbf {q} )d\mathbf {q} =-\iint _{\mathbf {q} \in \sigma }(\nabla U(\mathbf {q} ))\times d\mathbf {S}$.
* Given a multi-surface denoted by vector field $\mathbf {F}$, and a simple volume $\Omega$ with outwards oriented surface $\sigma$, then $\iiint _{\mathbf {q} \in \Omega }(\nabla \times \mathbf {F} (\mathbf {q} ))dV=-\iint _{\mathbf {q} \in \sigma }\mathbf {F} (\mathbf {q} )\times d\mathbf {S}$.

### The surfaces of intersections

Some more identities can be derived by examining the surfaces of volumes that result from intersections.

#### The inwards-oriented surface of volume-volume intersections

Start with a multi-volume $\mathbf {U} _{1}$, denoted by scalar field $U_{1}$, and a second multi-volume $\mathbf {U} _{2}$, denoted by scalar field $U_{2}$. The intersection $\mathbf {U} _{1}\cap \mathbf {U} _{2}$ is denoted by the scalar field $U_{1}U_{2}$.
Consider a volume $\Omega _{1}$ with weight $w_{1}$ from $\mathbf {U} _{1}$, and a volume $\Omega _{2}$ with weight $w_{2}$ from $\mathbf {U} _{2}$. Let $\sigma _{1}$ denote the inwards-oriented surface of $\Omega _{1}$, and let $\sigma _{2}$ denote the inwards-oriented surface of $\Omega _{2}$. There are two parts to the inwards-oriented surface of the intersection $\Omega _{1}\cap \Omega _{2}$, as shown in the image to the right. Part of the surface of $\Omega _{1}\cap \Omega _{2}$ consists of the portion of $\sigma _{2}$ that is contained by $\Omega _{1}$, which contributes the term $U_{1}(\nabla U_{2})$ to $\nabla (U_{1}U_{2})$. The other part of the surface of $\Omega _{1}\cap \Omega _{2}$ consists of the portion of $\sigma _{1}$ that is contained by $\Omega _{2}$, which contributes the term $(\nabla U_{1})U_{2}$ to $\nabla (U_{1}U_{2})$. Therefore the total surface of $\mathbf {U} _{1}\cap \mathbf {U} _{2}$ is $\nabla (U_{1}U_{2})=U_{1}(\nabla U_{2})+(\nabla U_{1})U_{2}$.
From the identity $\nabla (U_{1}U_{2})=U_{1}(\nabla U_{2})+(\nabla U_{1})U_{2}$, computing the total surface vector gives: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla (U_{1}(\mathbf {q} )U_{2}(\mathbf {q} ))dV=\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U_{1}(\mathbf {q} )(\nabla U_{2}(\mathbf {q} ))dV+\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla U_{1}(\mathbf {q} ))U_{2}(\mathbf {q} )dV$. The inwards-oriented surface of a multi-volume is a closed multi-surface, and the total surface vector of a closed surface is $\mathbf {0}$. $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\nabla (U_{1}(\mathbf {q} )U_{2}(\mathbf {q} ))dV=\mathbf {0}$ so hence $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U_{1}(\mathbf {q} )(\nabla U_{2}(\mathbf {q} ))dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla U_{1}(\mathbf {q} ))U_{2}(\mathbf {q} )dV$. The total surface vector of the intersection of multi-volume $\mathbf {U} _{1}$ with the inwards oriented surface of multi-volume $\mathbf {U} _{2}$ is the opposite of the total surface vector of the intersection of the inwards oriented surface of $\mathbf {U} _{1}$ with $\mathbf {U} _{2}$.
If $U_{1}$ denotes a simple volume $\Omega$ with outwards oriented surface $\sigma$, then the above integral identity becomes: $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}\delta _{3}(\mathbf {q} ;\Omega )(\nabla U_{2}(\mathbf {q} ))dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla \delta _{3}(\mathbf {q} ;\Omega ))U_{2}(\mathbf {q} )dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla U_{2}(\mathbf {q} ))dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(-\delta _{2}(\mathbf {q} ;\sigma ))U_{2}(\mathbf {q} )dV$ $\iff \iiint _{\mathbf {q} \in \Omega }(\nabla U_{2}(\mathbf {q} ))dV=\iint _{\mathbf {q} \in \sigma }U_{2}(\mathbf {q} )d\mathbf {S}$
In summary:
* Given two multi-volumes denoted by scalar fields $U_{1}$ and $U_{2}$, then the inwards-oriented surface of the intersection is: $\nabla (U_{1}U_{2})=U_{1}(\nabla U_{2})+(\nabla U_{1})U_{2}$.
* Given two multi-volumes denoted by scalar fields $U_{1}$ and $U_{2}$, then $\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}U_{1}(\mathbf {q} )(\nabla U_{2}(\mathbf {q} ))dV=-\iiint _{\mathbf {q} \in \mathbb {R} ^{3}}(\nabla U_{1}(\mathbf {q} ))U_{2}(\mathbf {q} )dV$.
* Given a simple volume $\Omega$ with outwards oriented surface $\sigma$ and a multi-volume denoted by scalar field $U_{2}$, then $\iiint _{\mathbf {q} \in \Omega }(\nabla U_{2}(\mathbf {q} ))dV=\iint _{\mathbf {q} \in \sigma }U_{2}(\mathbf {q} )d\mathbf {S}$.

### Summary

The tables below summarizes the results of the previous sections:
The endpoints, boundaries, and surfaces of intersections :
| structure 1 | structure 2 | intersection | endpoints, boundary, or surface |
| multi-path $\mathbf {J} _{1}$ | multi-volume $U_{2}$ | multi-path $\mathbf {J} _{1}U_{2}$ | multi-point $\nabla \cdot (\mathbf {J} _{1}U_{2})=(\nabla \cdot \mathbf {J} _{1})U_{2}+\mathbf {J} _{1}\cdot (\nabla U_{2})$ |
| multi-surface $\mathbf {F} _{1}$ | multi-surface $\mathbf {F} _{2}$ | multi-path $\mathbf {F} _{1}\times \mathbf {F} _{2}$ | multi-point $\nabla \cdot (\mathbf {F} _{1}\times \mathbf {F} _{2})=(\nabla \times \mathbf {F} _{1})\cdot \mathbf {F} _{2}-\mathbf {F} _{1}\cdot (\nabla \times \mathbf {F} _{2})$ |
| multi-surface $\mathbf {F} _{1}$ | multi-volume $U_{2}$ | multi-surface $\mathbf {F} _{1}U_{2}$ | multi-path $\nabla \times (\mathbf {F} _{1}U_{2})=(\nabla \times \mathbf {F} _{1})U_{2}-\mathbf {F} _{1}\times (\nabla U_{2})$ |
| multi-volume $U_{1}$ | multi-volume $U_{2}$ | multi-volume $U_{1}U_{2}$ | multi-surface $\nabla (U_{1}U_{2})=(\nabla U_{1})U_{2}+U_{1}(\nabla U_{2})$ |
Integral identities :
| Simple structure | Multi-structure | Integral identity | Identity name |
| simple path $C$, with starting point $C(0)$ and final point $C(1)$ | multi-volume $U$ | $\int _{\mathbf {q} \in C}(\nabla U)\cdot d\mathbf {q} =U(C(1))-U(C(0))$ | the gradient theorem |
| simple volume $\Omega$, with outwards-oriented surface $\sigma$ | multi-path $\mathbf {J}$ | $\iiint _{\mathbf {q} \in \Omega }(\nabla \cdot \mathbf {J} )dV=\iint _{\mathbf {q} \in \sigma }\mathbf {J} \cdot d\mathbf {S}$ | Gauss's divergence theorem |
| simple surface $\sigma$ with counter-clockwise oriented boundary $\partial \sigma$ | multi-surface $\mathbf {F}$ | $\iint _{\mathbf {q} \in \sigma }(\nabla \times \mathbf {F} )\cdot d\mathbf {S} =\int _{\mathbf {q} \in \partial \sigma }\mathbf {F} \cdot d\mathbf {q}$ | Stokes' theorem |
| simple surface $\sigma$ with counter-clockwise oriented boundary $\partial \sigma$ | multi-volume $U$ | $\iint _{\mathbf {q} \in \sigma }(\nabla U)\times d\mathbf {S} =-\int _{\mathbf {q} \in \partial \sigma }Ud\mathbf {q}$ | unnamed |
| simple volume $\Omega$, with outwards-oriented surface $\sigma$ | multi-surface $\mathbf {F}$ | $\iiint _{\mathbf {q} \in \Omega }(\nabla \times \mathbf {F} )dV=-\iint _{\mathbf {q} \in \sigma }\mathbf {F} \times d\mathbf {S}$ | unnamed |
| simple volume $\Omega$, with outwards-oriented surface $\sigma$ | multi-volume $U$ | $\iiint _{\mathbf {q} \in \Omega }(\nabla U)dV=\iint _{\mathbf {q} \in \sigma }Ud\mathbf {S}$ | unnamed |