Boolean Combination of Circular Arcs using Orthogonal Spheres

An orthogonal sphere representation of arcs on spatial circles can be used to compactly perform Boolean combinations of such arcs. We formulate this using conformal geometric algebra, of which the oriented nature allows both minor and major arcs to be treated. Easily computable quantities discriminate the cases of relative positions. An application in the first stages of a problem in Discretizable Molecular Distance Geometry is included. We give a suggestion on how to extend this characterization by orthogonal spheres to the manifolds of arcs in the subsequent stages, using probabilistic eigenspheres of the distributions.


Introduction
The central problem in distance geometry [10] is the determination of a spatial configuration of points based on their given mutual distances. As a mathematical puzzle, it was originated by Menger [11]. The field has practical applications for detecting molecular structure of proteins, since nuclear magnetic resonance techniques exist to determine distances between neighboring atoms [1].
As a start to an algorithmic solution to this puzzle, simplifications have been made. It is for instance assumed that bond lengths and angles along the protein carbon backbone are known exactly (unit lengths, and 120 degree angles). Furthermore, in the Discretizable Molecular Distance Geometry Problem (DMDGP) [8], it is assumed that a few more distances between neighboring atoms along the chain are known as intervals [7]. For a set of This article is part of the Topical Collection on Proceedings of AGACSE 2018, IMECC-UNICAMP, Campinas, Brazil, edited by Sebastià Xambó-Descamps and Carlile Lavor. vertices of the molecular graph, one can then provide a distance matrix such as this one from [2] and [6]: where entry (i, j) is the distance d(i, j) (or a distance interval), between atom v i and atom v j and * denotes an unknown distance. The known exact distances restrain relative positions between atoms to spheres. Two such given distances then determine a circle, and three distances a possible point pair of admissible positions. The interval distances determine spherical shells of possible positions, which combine with the circles to restrain possible positions to arcs (for an example see Fig. 1). As one propagates from an initial triangle (v 1 , v 2 , v 3 ) along the molecular chain, the uncertainties combine to arcs of arcs of possible positions, etc.
Because of the prevalence of spheres and their intersections in this problem, [1,2] and [6] have addressed the problem by conformal geometric algebra (CGA). This is the geometric algebra of a vector space whose vectors can be interpreted as spheres, whose 2D subspaces can be interpreted as circles, and whose 3D subspaces are point pairs (the intersection of three spheres). Its elements can be given an orientation, which is convenient to denote the orientation of circles and arcs. CGA has a compact way to represent rotations and translations of arbitrary elements through spinors. In an example from [6], an interval distance of X 4 to X 1 constrains the possibilities for X 4 to two arcs Starting from the first few vertices with interval distances, the DMDGP problem is one of combining circular arcs in space. However, conformal geometric algebra has no natural element to represent such arcs, and so the authors of [2,6] used the motion representation provided by the spinors to generate discretely sampled points (in [2]) or continuously parametrized orbits (in [6]) as their arc representation.
In the present paper, we show that there is a way to represent arcs in CGA by means of orthogonal spheres to the circle they reside on. This leads to a compact specification of their relative positions, and to a compact computation of the Boolean combinations of those arcs (such as the required intersections). That solves the first few steps along the molecular backbone in DMDGP. We end the paper with a suggestion of how a similar characterization based on a decomposition into orthogonal spheres (from [4]) may be employed for the 'arcs of arcs', 'arcs of arcs of arcs' et cetera that one encounters further down the backbone, albeit in a probabilistic manner.
The paper starts with a very brief summary of conformal geometric algebra in Sect. 2. Then we introduce the orthospheres, both real and imaginary, that characterize a spatial arc in Sect. 3. The different relative positions of these arcs can be determined in terms of their real orthospheres and this informs the results of their Boolean operations. In Sect. 4, we specify these operations for both minor and major arcs, using the orientability properties of the conformal elements. In Sect. 5, we treat the first steps in the processing of the distance matrix above. We then indicate how this deterministic representation of arcs is a special case of the probabilistic representation of more general manifolds like 'arcs of arcs' that will be encountered further down the protein backbone, and how those might then be processed conformally.

An Algebra of Spheres
Our solution to the arc intersection is cast in terms of conformal geometric algebra. That algebra is very well suited to deal with n-D spheres, their orthogonality and intersection. 1 A sphere in 3D Euclidean space, with center at Euclidean position vector c and radius ρ, is represented homogeneously by a vector in a space with two additional dimensions (orthogonal to the three Euclidean dimensions), which are here denoted by the additional basis vectors o and ∞. This augmented space is given a metric for which the Euclidean part (denoted by bold font) retains its usual dot product. With that metric, the dot product between two sphere representatives is: L. Dorst Adv. Appl. Clifford Algebras Spheres that intersect orthogonally thus have a zero dot product. As a special case, a point is a zero radius sphere and it has zero dot product with any sphere it lies on; we denote the normalized vector representative of a point at Euclidean location c by CGA(c), as in Eq. (1). The dot product c 1 · c 2 between two normalized points c 1 and c 2 gives (minus one half times) their squared Euclidean distance d 2 (c 1 , c 2 ). 'Normalized' means that the homogeneous representation weight of the sphere (or point) is chosen so that the coefficient of o equals 1; this implies that −∞ · σ = 1. The additional basis vectors o and ∞ can be interpreted as the point at the (arbitrary) origin, and the point at infinity, respectively. Planes are also included in this model, as elements of the form π = n + δ ∞; it is easy to verify that a point represented by vector x = CGA(x) = o + x + 1 2 x 2 ∞ is on a plane iff x · π = 0, since this is equivalent to x · n = δ. Since ∞ · π = 0, a plane is actually just a sphere containing the point at infinity.
The representation of spheres can be dualized, through division by the volume element of the (3 + 2)-dimensional space, which converts a vector into an 4-dimensional hyperplane. This dualization is denoted by a * , as in X → X * ; its inverse dualization is denoted by X → X − * . In 3D, X − * = −X * . This duality swaps a sphere from its vector representation by center and radius to the four-vector description as the element that orthogonally hits 4 points (this was dubbed the 'plunge' operation in [5], Chapter 15). Dualization of the 2D element σ 1 ∧ σ 2 representing the meet of two spheres is to the 3D element (σ 1 ∧ σ 2 ) * , representing the common circle directly, in a form that can alternatively be written as the plunge of three points p, q and r on the circle: p ∧ q ∧ r. Dualization raises the required techniques in the representational space from those of Grassmann algebra (spanning vectors) to those of Clifford algebra (allowing dot products and division by vectors and subspaces). Using the Clifford geometric product, the square of a normalized sphere is the square of its radius, and the square of a plane is the squared norm of its normal vector. The dualization moreover permits conversion between dot and outer products: which allows flexible switching between 'meet' and 'plunge' of elements and their duals. The intersection operation is closed: it always gives a result in the algebra, which moreover can be interpreted geometrically. That is how we can use the intersection of spheres (and planes) to produce circles, lines, point pairs, and more. When spheres do not intersect in a real circle, they always still intersect in an imaginary circle (or, exceptionally, a 2D tangent element at a point). In the conformal model, this is simply an element of which the square (employing the geometric product) is negative (or zero in the exceptional case). Complex numbers are not required for this, since only squared radii occur in the CGA representation. In fact, analyzing the dot product of the (3 + 2)-D representational space, we find that it is R 4,1 , a space with a Minkowski metric in which an orthogonal basis can be given with 4 vectors with a positive square, and 1 with a negative square. Geometric points are represented by null vectors of this space.
This brief summary of the essence of the conformal model is the bare minimum required to interpret the equations in this paper. If you have not seen this before, the compact but expressive results of the paper may inspire further reading. An accessible introduction at the geometrical level is the tutorial book [5], which comes with interactive software to familiarize yourself with this powerful language for interacting spheres. All equations in this paper can be easily typed in, and are automatically visualized, producing graphics like Fig. 2.

Characterization of Circular Arcs in Space
An arc is a segment of a circle; that circle can be given a natural orientation relative to the orientation of the plane it resides in. The arc can be characterized by the carrier circle and the oriented point pair of its end points. The orientation of the point pair then determines whether it represents the minor or the major arc.
In the context of the distance geometry puzzle, we need to be able to combine arcs on the same circle: determining whether or not they overlap, and if they do, find a smaller arc representing their intersection. That is rather involved to execute in the direct point pair representation; this led the authors of [9] to represent the arc by sampled points, and later in [2,6] to introduce rotors whose motion is the arc. In the present paper, we show that arc combination becomes easily tractable when we factorize the point pair in terms of spheres; and subsequent computations are easier when those spheres are orthogonal to the circle, and to each other. So 'orthogonal factorization' will be our approach to associate unique spheres with arcs.

Support and Surround of a Circle
First, the circle K on which an arc is defined is typically determined through the intersection (meet) of two dual spheres σ i and σ j , as K = (σ j ∧ σ i ) − * (in the protein example, the spheres are determined by the two previous exact distances from the current atom v i to the previous atoms v i−1 and v i−2 ; this represents the known angles and distances along the protein backbone). The value of K is given in this form by the data. When we need to do any computations with K, it is better to characterize this given K not as the meet of the two original specific non-orthogonal spheres, but as an orthogonal factorization of a plane and a sphere. This was treated in [5], Chapter 15, and plane and sphere were called 'support' and 'surround' of the circle, respectively. It is slightly easier to compute with duals, so we define κ = K * = σ j ∧ σ i . The (dual) circle κ can be then factorized through: Consequently κ = σ ⊥ π and π · σ ⊥ = 0 (and of course ∞ · π = 0), so this indeed achieves the orthogonal factorization of the dual circle κ by the (dual) plane π and the (dual) orthogonal sphere σ ⊥ .

The Real Orthosphere of an Arc
When an additional interval distance to a third atom is known, the circle determined by the two first distances is cut by the spheres corresponding to the upper and lower bound. Admissible solutions should be inside the upper bound and outside the lower bound, which produces spatial arcs of possible positions.
Algebraically, an arc of the circle K is determined by a point pair P = p − ∧ p + constructed from its end points p − and p + (with p − · κ = p + · κ = 0, since they are on the circle K dually represented as κ). We can factorize this point pair on K (or more conveniently its dual 3-blade) by means of the geometric product in an orthogonal manner. We expose two ways of doing so, the first more natural for the processing of arcs on the same circle, the second (in Sect. 3.4) an extension which will generalize in Sect. 5.3 to the representation of 'arcs of arcs' etc.
The circle on which the arc resides is characterized as the intersection of its surround σ ⊥ and its support plane π. For the arc, we compute the unique sphere σ orthogonal to the support π that cuts the surround σ ⊥ perpendicularly to produce the point pair P (see Fig. 2). So we desire, with an unspecified proportionality factor: (We will denote proportionality by ∼ rather than ∝ since the latter's resemblance to an α-symbol makes equations hard to parse visually.) The solution is that σ is positively proportional to −P · K: The derivation of Eq. (2) is: which is (positively) proportional to the above (since below we show that in K −1 = K/K 2 , the factor K 2 is non-negative). Let us call σ the 'orthosphere' of the circular arc. We will use the sign of the point pair P to distinguish minor and major arcs later on; we therefore need to be careful about the sign of its representing orthosphere σ as well, so as not to lose that interpretation. To be specific, we take the circle K as characterized positively proportional to p − ∧ r ∧ p + , for a point pair p − ∧ p + with a distance smaller than the diameter of K, and a normalized point r between the normalized points p − and p + when traveling along the circle in its positive orientation direction. Then: so as defined in Eq. (2), the sphere σ has non-negative weight. Equality to zero only occurs when the point pair forms a diameter of κ (which is clearly the ambiguous case), since only then do p − , r and p + form a right-angled triangle (for which Pythagoras' theorem cancels the second factor).
It makes sense, from dimensional analysis, to define σ = −P/K or even use its inverse 1/σ as the characterization of the orthosphere. Since those differ from the above by positive factors (because , this is permitted without changing our sign conventions. We have not found a fully satisfactory set of dimensional scaling conventions for our objects; but as long as we use scaling invariant interpretations of their geometrical parameters (provided in [5], Table 14.1), this is not a major concern. Instead, we typically normalize any sphere s through dividing by −∞ · s, to force its coefficient of o to be equal to 1.
We can rewrite the orthosphere in terms of the dual circle κ = K * as: This is a simple rewriting: But σ can also be expressed in other ways, revealing different aspects of its geometry: In the final rewriting, σ ⊥ ∧P is the circle through P , perpendicular to σ ⊥ , and σ is thus proportional to its (dual) surround. In the penultimate step, ∞ ∧ K is the support plane, and ∞ ∧ σ ⊥ ∧ P is another form of that same plane. Equivalence is easily established by looking at its dual: Spatially, the situation is as in Fig. 2.

Minor and Major Arcs by Oriented Point Pairs on Oriented Circles
The previous treatment and depictions used arcs that consist of points between the endpoints along the shortest (minor) circular arc; we could also wish to characterize major arcs that run along the longest (major) connection between the points.
When the points are united in a point pair by the outer product, they naturally acquire an order; and reversing that order changes the sign. Letting an arc run always from the first point in a point pair to the second, along the orientation of the circle, we can therefore encode both types of arcs. We establish a unique point order associated with a point pair element P using the point pair decomposition from [5] pg. 427. A point pair P can be decomposed into its constituent normalized points p − and p + by: Then P ∼ p − ∧ p + , modulo a positive scalar factor. For a point pair P whose points reside on an oriented circle K, we have fixed the orientation sign for the associated orthosphere by relating them through σ = −P · K = −P K. As a consequence, we can retrieve a properly oriented point pair from its orthosphere: or, if one wants to avoid the inverse, P is positively proportional to −σ · K. Let us verify that the signs of the minor and major arcs indeed correspond to orthospheres with positive and negative weight, respectively. We express the circle K by means of the two elements of the point pair, plus a third: K = p − ∧ r ∧ p + . If r is inside the orthosphere of p − ∧ p + , this orients the circle from p − to p + as we desire. If the point r has negative weight, the orientation of K is that of p − ∧ p + ∧ r, so traversed in the order p − , then p + then r. Let us now conform the sign of the weight of σ = c − 1 2 ρ 2 ∞: . So indeed when r is inside the orthosphere of p − ∧ p + , the proportionality factor between σ · K and (−P ) is positive. And when r is outside that orthosphere, and hence used to produce a major arc for the same point pair P , then σ changes its sign.

The Imaginary Orthosphere
In the 3D space we are considering, one can make many choices for an orthogonal basis of its 5D conformal geometric algebra. It makes sense to pick a basis adapted to the objects one wants to describe. For our circular arc, we have already introduced three mutually orthogonal conformal vectors: the plane π of the circle, its surrounding sphere σ ⊥ , and the orthosphere σ. Another obvious vector is the plane of symmetry of the arc, which is π ≡ ∞·P = p + − p − (the midplane of the two points). These four conformal vectors have positive squares; the fifth orthogonal vector should therefore be an imaginary sphere σ, see [15]. We show below that it is proportional to: Orthogonality with the other vectors is easily established. When drawn in the support plane of the circle K (see Fig. 3), the sphere σ is seen to contain the chord of the arc as its diameter. And indeed, the squared radius is σ 2 = 1 2 (p + · p − ) = − d(p − , p + )/2 2 , so that the associated diameter of σ is equal to the chord length of the arc. The imaginary orthosphere σ is algebraically orthogonal to the real orthosphere σ in the manner of imaginary spheres: σ is intersected by σ in a great circle on σ (see [5], Chapter 14) and it also cuts the other real basis vectors π, σ ⊥ and π equatorially. A simple derivation of the expression Eq. (4) for σ is to write the volume element of the representational space as π ∧ σ ⊥ ∧ π ∧ σ ∧ σ = π σ ⊥ π σ σ (due to orthogonality, all outer products are geometric products). This gives σ ∼ (π σ ⊥ π σ) * = (π κ σ) * ∼ (∞ · P ) P − * * = (∞ · P ) P ∼ (∞ · P )/P . In Sect. 3.3, we related the sign of σ to that of the point pair, so that we could represent minor and major arcs; let us do the same for σ. The definition Eq. (4) is insensitive to the sign (and even to the weight) of P . We can fix a sign explicitly by requiring that (∞ · σ) (∞ · σ) > 0, or implicitly by picking the appropriate order of terms in defining the pseudoscalar of the whole space in the derivation above.
An alternative rewriting of σ is: Deriving this form is done through:  Figure 3. A sketch of the plane π of the circle κ containing a minor arc. The thick red arc on the red circle κ can be characterized by the dashed green imaginary sphere σ related to the inner chord, and by the solid blue sphere σ related to the outer chord (the real orthosphere). The two spheres σ and σ are mutually algebraically orthogonal, and both algebraically orthogonal to the red circle κ. Since that makes them orthogonal to the carrier plane as well, we can show both of them cut through an equator in this plane π, giving this 2D drawing (color figure online) so that π can be factored out of numerator and denominator. Geometrically, σ∧σ ⊥ is an imaginary point pair (dual to the real circular intersection of σ ⊥ and σ); then σ is its (dual) surround. In a visualization, these two spheres σ and σ appear to intersect in the real circle (σ ∧ σ ⊥ ) * , even though the imaginary sphere σ does not contain any real points. (A more exact way to state this is to first convert σ into a real sphere with the same center and radius, which is done by σ → σ (σ ∧ ∞).) Note that the imaginary orthosphere σ is actually independent of the circle K, a fact which is better seen in Eq. (4) than in Eq. (5). It is therefore less descriptive of the arc than the real orthosphere σ, which does involve the carrier circle K.  Figure 3 illustrates that the orthospheres σ and σ to the circular arc have a radius that is simply related to basic properties of the arc. Let us first define the (inner) chord as usual, and the outer chord as the broken segment composed of the two tangents to the circular arc at its end points. Then the squared radius of the normalized σ, computed as σ 2 , is minus the squared size of half the inner chord of the arc (which is ρ/ sin(φ/2) for an arc of φ).

The Chords of an Arc
The squared radius of the normalized σ, computed as σ 2 , is the squared size of half of the outer chord (which is ρ/ tan(φ/2)). The centers of those spheres are the location of the center of the inner chord, and the vertex of the outer chord. Two orthospheres σ and σ thus encode the inner and outer chord to the arc, though they are not specific on the plane of the arc. The surround sphere of the arc's circle K (but not its support plane) can be determined from σ and σ. We observe (see below) that its imaginary equator is K = (σ ∧ σ) * , and then use that to determine the sphere as the surround of the equator: σ ⊥ = K/(∞ ∧ K). We show the equatorial nature of K = (σ ∧ σ) * by computing its squared radius from its dual κ as −(κ κ)/(∞ · κ) 2 according to [5], Table 14.1. First, defining ρ 2 = σ 2 and ρ 2 = σ 2 , Then Hence the reciprocal of the squared radius −r 2 of the imaginary circle K satisfies: 1 r 2 + 1 ρ 2 = 1 ρ 2 . This is the 'reciprocal Pythagorean theorem', the relationship in a right triangle between the magnitudes of the sides r and ρ and its altitude ρ on the hypotenuse [12]. Inspecting Fig. 3, this implies that r 2 = σ 2 ⊥ , and therefore r 2 is indeed equal to the squared radius of the circle of the arc.

Establishing Relative Positions of Arcs
The reason for using orthospheres to the circle K is that the overlapping and containment conditions of the spheres completely coincide with those of the arcs they represent, simply because they intersect the circle K containing the arcs perpendicularly; non-orthospheres (such as the surround of the end points) would not provide this equivalence. In this section, we first treat the case of 'inner arcs' (minor arcs) of the point pairs to focus our explanation; we can then easily expand this to major arcs by incorporating the point pair signs afterwards.
The meet of two real orthospheres σ 1 and σ 2 is a circle (like that of any two spheres).

L. Dorst
Adv. Appl. Clifford Algebras  • When this circle is real (so its square is positive and the square of its dual is negative), there is a maximal smaller sphere orthogonal to K that is contained in both. This is the orthosphere representing the intersection of the arcs of σ 1 and σ 2 (see Fig. 4, top row). • When this circle is imaginary (so its square is negative and the square of its dual is positive), the real spheres do not really meet. This could be because their arcs are not overlapping, or because one arc is fully contained in the other. In the former case, we can eliminate this incidence of arcs from consideration; in the latter, the smallest orthosphere should remain. The discrimination between these two cases is the sign of the inner product of the orthospheres: it is positive for containment, and negative for exclusion (see Fig. 4, bottom row). • When the common circle squares to zero, it has degenerated to a point (actually a tangent 2-vector). The arcs are now touching; one could return the common point, or report lack of overlap (see Fig. 4, middle row).
The dual representation of the intersection circle is σ 1 ∧ σ 2 , and we thus can specify the cases simply in terms of the signs of (σ 1 ∧ σ 2 ) 2 and of σ 1 ·σ 2 . The correspondence may be established by expressing those quantities in terms of the radii ρ 1 and ρ 2 , plus the distance between the centers δ ≡ d(c 1 , c 2 ), for normalized orthospheres σ 1 = c 1 − 1 2 ρ 2 1 ∞ and σ 2 = c 2 − 1 2 ρ 2 2 ∞: The former is zero when δ = |ρ 1 ± ρ 2 |, so at tangency (with equal-radius spheres at zero distance as special case). The latter is zero when δ 2 = ρ 2 1 + ρ 2 2 , which is when the real spheres intersect perpendicularly (but for the classification of arc intersections, the discrimination which the dot product offers on the acuteness of the angle of intersection is irrelevant). The signs for the situations on either side of zero are easily established. When containment occurs, one can compare the squared radii σ 2 1 and σ 2 2 to decide which of the two orthospheres is contained in the other.

The Four Tangent Orthospheres
The cases as distinguished in Fig. 4 will produce new arcs with new orthospheres as a result of arc intersection. These are usually composed from endpoints of the existing arcs. We will show how this is done in Sect. 4.3; let us decide on a consistent labeling beforehand, so that we can refer to those new potential orthospheres.
We will have to form new arcs representing the intersection from the endpoints of the input arcs. We visualize the cases that occur, independently of the circle K they occur on, by a simple conformal trick: we perform an inversion in a well chosen sphere that turns the circle K into a line. Any sphere σ I composed as σ I = κ − · (κ + ∧ ∞) from an anti-polar point pair κ − ∧ κ + of K will do; the anti-polar point pair can be made through intersection of K with any (non-carrier) plane through the center of K.
When K becomes a straight line K = −σ I K/σ I under this inversion, the two orthospheres of course remain orthospheres to this line (since spherical inversion is conformal). Now along this line K , introduce 1-D coordinates in its direction e. Let the centers and radii of the transformed orthospheres σ p and σ q be p, q, ρ p and ρ q , respectively. (These can be determined as the Euclidean parts of −σ p ∞/σ p and −σ q ∞/σ q , and from ρ 2 p = (σ p ) 2 and ρ 2 q = (σ q ) 2 .) Then the four tangent spheres have centers at: c ±± = 1 2 (p ± ρ p e) + 1 2 (q ± ρ q e), (where the index signs are in the same order as in the equation, so select the corresponding top or bottom sign). The radii of the new orthospheres are, respectively: Now construct the resulting normalized orthospheres as σ ±± = CGA(c ±± ) − 1 2 ρ 2 ±± ∞, (where CGA(·) produces the conformal point vector at a Euclidean location). Then just re-invert the situation relative to σ I to retrieve the four desired orthospheres tangent to the original σ p and σ q . Besides the original arcs, four additional arcs are thus part of our vocabulary of solutions to the arc intersection problem.
In actual practice, one does not compute the spheres σ ±± from this inversion and re-inversion, but simply from the conformal points p − and p + derived directly from P (or from P = −σ p /K) by decomposing it according to Eq. (3), and similarly for the points related to the point pair Q = −σ q /K. Then for a chosen sphere, say σ +− , those points are recomposed to a point pair P +− = p + ∧ q − , and that is made into the desired orthosphere relative to the original circle K as σ +− = −P +− · K. The inversion above was just done to easily depict the relationships between the labelled orthospheres, in terms which are invariant under that conformal transformation.
Incidentally, the inversion in σ I makes the imaginary orthosphere σ of a point pair look coincident with the corresponding real orthosphere σ, since both the real and imaginary spheres are orthogonal to the same 'flat' element (the line K ) that contains their common point pair.

Boolean Operations on Arcs
Even though the sign of an orthosphere σ was related to the arc orientation in Sect. 3.3, we note that the classification of the relative situation of the orthospheres as in Fig. 4 mostly involves squares, and is therefore not too dependent on the signs. Only when one orthosphere represents a minor arc, and the other a major arc, should we modify the label of the classification of the columns (since it is based on the dot product).
The orthosphere we should return as a result of a Boolean combination of arcs should have the correct sign, to represent the resulting arc appropriately. We characterize the results by the signs of the point pairs (and hence of the orthospheres relative to the circle K). And while we are at it, we might as well not only state the result of intersection of the point pairs, but also of their union (even thought we do not really need that in the DMDGP problem).
With the labeling as in Fig. 5, we can specify what the outcome should be of the various Boolean operations on the arcs represented by the orthospheres ±σ p and ±σ q . In Fig. 6, we list the outcomes following the order in Fig. 4 (from bottom left upwards, then down at the right). Usually, the outcome is a single arc; but when the original arcs are completely separated or contained, the result of their Boolean combination may be a pair of arcs.
In the situations where σ p and σ q are tangent, only two extra orthospheres are generated, one of which degenerates to a point; this is indicated in those columns as '(pt)'.
Clearly if one would label p and q conversely, all index signs also swap, so the difference between σ −+ and σ +− is not independent of this order,  Figure 5. After inversion, the circle K becomes the blue line K , the orthospheres σ p and σ q become orthospheres σ p and σ q to that line (in red). The parameters of the four blue tangent orthospheres are read off easily (the direction vector e of the line K is taken to the right, so a label like +− can be read as: connects axis point to the right (+) of p with that to the left (−) of q). There are three main situations: the red orthospheres intersect, one is contained in the other, or they are fully separate (color figure online)