View source for Dupin cyclide

{{short description|Geometric inversion of a torus, cylinder or double cone}}
[[Image:Cyclide.png|thumb|A Dupin cyclide]]
In [[mathematics]], a '''Dupin cyclide''' or '''cyclide of Dupin''' is any [[Inversive geometry|geometric inversion]] of a [[standard torus]], [[Cylinder (geometry)|cylinder]] or [[cone|double cone]]. In particular, these latter are themselves examples of Dupin cyclides. They were discovered {{circa|1802}} by (and named after) [[Charles Dupin]], while he was still a student at the [[École polytechnique]] following [[Gaspard Monge]]'s lectures.<ref>{{Harvnb|O'Connor|Robertson|2000}}</ref> The key property of a Dupin cyclide is that it is a [[channel surface]] (envelope of a one-parameter family of spheres) in two different ways. This property means that Dupin cyclides are natural objects in [[Lie sphere geometry]].

Dupin cyclides are often simply known as ''cyclides'', but the latter term is also used to refer to a more general class of [[quartic surface]]s which are important in the theory of separation of variables for the [[Laplace equation]] in three dimensions.

Dupin cyclides were investigated not only by Dupin, but also by [[Arthur Cayley|A. Cayley]], [[James Clerk Maxwell|J.C. Maxwell]] and [[Mabel M. Young]].

Dupin cyclides are used in [[computer-aided design]] because cyclide patches have rational representations and are suitable for blending canal surfaces (cylinder, cones, tori, and others).

==Definitions and properties==
There are several equivalent definitions of Dupin cyclides. In <math>\R^3</math>, they can be defined as the images under any inversion of tori, cylinders and double cones. This shows that the class of Dupin cyclides is invariant under [[Möbius transformation|Möbius (or conformal) transformations]].
In complex space <math>\Complex^3</math> these three latter varieties can be mapped to one another by inversion, so Dupin cyclides can be defined as inversions of the torus (or the cylinder, or the double cone).

Since a standard torus is the orbit of a point under a two dimensional [[Abelian group|abelian]] [[subgroup]] of the Möbius group, it follows that the cyclides also are, and this provides a second way to define them.

A third property which characterizes Dupin cyclides is that their [[curvature line]]s are all circles (possibly through the [[point at infinity]]). Equivalently, the [[curvature sphere]]s, which are the spheres [[tangent]] to the surface with radii equal to the [[Multiplicative inverse|reciprocals]] of the [[principal curvature]]s at the point of tangency, are constant along the corresponding curvature lines: they are the tangent spheres containing the corresponding curvature lines as [[great circle]]s. Equivalently again, both sheets of the [[focal surface]] degenerate to conics.<ref>{{Harvnb|Hilbert|Cohn-Vossen|1999}}</ref> It follows that any Dupin cyclide is a [[channel surface]] (i.e., the envelope of a one-parameter family of spheres) in two different ways, and this gives another characterization.

The definition in terms of spheres shows that the class of Dupin cyclides is invariant under the larger group of all [[Lie sphere transformation]]s; any two Dupin cyclides are [[Lie sphere geometry|Lie-equivalent]]. They form (in some sense) the simplest class of Lie-invariant surfaces after the spheres, and are therefore particularly significant in [[Lie sphere geometry]].<ref>{{Harvnb|Cecil|1992}}</ref>

The definition also means that a Dupin cyclide is the envelope of the one-parameter family of spheres tangent to three given mutually tangent spheres. It follows that it is tangent to infinitely many [[Soddy's hexlet]] configurations of spheres.

== Parametric and implicit representation ==
: '''(CS):''' A Dupin cyclide can be represented in two ways as the [[Envelope (mathematics)|envelope]] of a one parametric pencil of spheres, i.e. it is a [[canal surface]] with two [[Directrix (generatrix)|directrices]]. The pair of directrices are [[focal conics]] and consists either of an ellipse and a hyperbola or of two parabolas. In the first case one defines the cyclide as ''elliptic'', in the second case as ''parabolic''. In both cases the conics are contained in two mutually orthogonal planes. In extreme cases (if the ellipse is a circle) the hyperbola degenerates to a line and the cyclide is a torus of revolution.

A further special property of a cyclide is:
: '''(CL):''' Any [[curvature line]] of a Dupin cyclide is a ''circle''.

=== Elliptic cyclides ===
An elliptic cyclide can be represented parametrically by the following formulas (see section [[#Cyclide as channel surface|Cyclide as channel surface]]):
[[File:Zyklide-el-hy-def.svg|thumb|upright=1.3|Meanings of the design parameters  <math>a,b,c,d</math>:<br>
<math>d</math> is the radius of the generating sphere at the co-vertices of the ellipse<br>
The two circles in the x-z-plane with centers <math>(\pm a,0,0)</math> have radii <math>d\mp c</math>.<br>
Here: <math>a=1,\; b=0.98 \to c= 0.199 </math> and <math> d=0.3</math>]]
:<math> x=\frac{d(c-a\cos u\cos v)+b^2\cos u}{a-c\cos u \cos v} \ ,</math>
:<math> y=\frac{b\sin u (a-d\cos v)}{a-c\cos u \cos v} \ ,</math>
:<math> z=\frac{b\sin v (c \cos u-d)}{a-c\cos u \cos v} \ ,</math>
:<math>0\le u,v <2\pi \ .</math>
The numbers <math>a,b,c,d</math> are the semi major and semi minor axes and <math>c</math> the linear eccentricity of the ellipse: 
:<math>\frac{x^2}{a^2}+\frac{y^2}{b^2}=1, z=0\ .</math> 
The hyperbola <math>\frac{x^2}{c^2}-\frac{z^2}{b^2}=1, y=0</math>
is the focal conic to the ellipse. That means: The foci/vertices of the ellipse are the vertices/foci of the hyperbola. The two conics form the two degenerated [[focal surface]]s of the cyclide.

<math>d</math> can be considered as the average radius of the generating spheres.

For <math>u=const</math> , <math>v=const </math> respectively one gets the curvature lines (circles) of the surface.

The corresponding [[Implicit surface|implicit representation]] is:
:<math>(x^2+y^2+z^2+b^2-d^2)^2-4(ax-cd)^2-4b^2y^2=0 \ .</math>

In case of  <math>a=b</math> one gets <math>c=0</math>, i. e. the ellipse is a circle and the hyperbola degenerates to a line. The corresponding cyclides are tori of revolution.

{| class="wikitable" style="background-color: white; margin: 1em auto 1em auto; text-align: center;"
! colspan="6" style="background-color:white;"     | (ellipt.) Dupin cyclides for designparameters a,b,c,d
|-

| style="padding-left:3em; padding-right:3em" | <math> d=0</math>

| <math>0<d<c</math>
| <math>d=c</math>
| <math>c<d</math>
| <math> d=a </math>
|<math> a<d</math>

|-
|[[File:Zyklide-0.svg|150px]]
| [[File:Zyklide-11.svg|150px]]
|[[File:Zyklide-1.svg|150px]]
|[[File:Zyklide-2.svg|150px]]
|[[File:Zyklide-4.svg|150px]]
|[[File:Zyklide-3.svg|150px]]

|-
|symm. horn cyclide
|horn cyclide
|horn cyclide
|ring cyclide
|ring cyclide
|spindle cyclide
|}

More intuitive design parameters are the intersections of the cyclide with the x-axis. See section [[#Cyclide through 4 points on the x-axis|Cyclide through 4 points on the x-axis]].

=== Parabolic cyclides ===
A parabolic cyclide can be represented by the following parametric representation (see section [[#Cyclide as channel surface|Cyclide as channel surface]]):
[[File:Zyklide-parab-1.svg|thumb|upright=1.2|parabolic cyclide with its directrices (focal parabolas)]]
:<math> x=\frac{p}{2}\, \frac{2v^2+k(1-u^2-v^2)}{1+u^2+v^2} \ , </math>
:<math> y=pu\, \frac{v^2+k}{1+u^2+v^2} \ , </math>
:<math> z=pv\, \frac{1+u^2-k}{1+u^2+v^2} \ , </math>
:<math> -\infty<u,v<\infty \ . </math>
The number <math>p</math> determines the shape of both the parabolas, which are focal conics: 
:<math>y^2=p^2-2px, \ z=0\ </math> and <math>\ z^2=2px, \ y=0\ .</math>

<math>k</math> determines the relation of the diameters of the two holes (see diagram). <math>k=0.5</math> means: both diameters are equal. For the diagram is <math>k=0.7</math>.

A corresponding implicit representation is
:<math>\left(x + \left(\frac{k}{2} - 1\right)p \right)\left(x^2+y^2+z^2 - \frac{k^2p^2}{4}\right) + pz^2 = 0 \ .</math>

{| class="wikitable" style="background-color: white; margin: 1em auto 1em auto; text-align: center;"
! colspan="3" style="background-color:white;"     | parabolic Dupin cyclides for designparameters p=1, k
|-

| style="padding-left:3em; padding-right:3em" |<math>k=0.5</math>
|<math>k=1</math>
|<math>k=1.5</math>

|-
|[[File:Zyklide-p-0.svg|200px]]
|[[File:Zyklide-p-1.svg|200px]]
|[[File:Zyklide-p-2.svg|200px]]

|-
|ring cyclide
|horn cyclide
|horn cyclide
|}

''Remark'': By displaying the circles there appear gaps which are caused by the necessary restriction of the parameters <math>u,v</math>.

== Cyclide as channel surface ==

[[File:Zyklide-kanalfl-1-2.svg|thumb|upright=1|Dupin cyclide as channel surface ([[Envelope (mathematics)|envelope]] of a family of spheres)]]

There are two ways to generate an elliptic Dupin cyclide as a [[channel surface]]. The first one uses an ellipse as directrix, the second one a hyperbola:<ref>W. Blaschke: ''Analytische Geometrie'', Springer-Verlag, 2013, {{ISBN|303486812X}}, S. 115</ref>

=== Ellipse as directrix ===

In the x-y-plane the directrix is the ellipse with equation
:<math>\frac{x^2}{a^2}+\frac{y^2}{b^2}=1\ ,\ z=0\quad </math> and <math>\ a>b</math>. 
It has the parametric representation
:<math>\ x=a\cos\varphi\; ,\ y=b\sin \varphi \;, \ z=0\ .</math>
<math>a</math> is the semi major and <math>b</math> the semi minor axis.
<math>c</math> is the linear eccentricity of the ellipse. Hence: <math>c^2=a^2-b^2</math>.
The radii of the generating spheres are
:<math>r(\varphi)=d-c\cos\varphi \ .</math>
<math>d</math> is a design parameter. It can be seen as the average of the radii of the spheres. In case of <math>c=0</math> the ellipse is a circle and the cyclide a torus of revolution with <math>d</math> the radius of the generating circle (generatrix).

In the diagram: <math>\; a=1,\; b=0.99,\ d=0.25\; </math>.

[[File:Zyklide-kanalfl-3.svg|thumb|upright=1|Maxwell: Property of the foci of the directrix ellipse. The ellipse is the [[equidistant set]] of the blue and purple circles.]]
=== Maxwell property ===
The following simple relation between the actual sphere center (ellipse point) and the corresponding sphere radius is due to Maxwell:<ref>mentioned in W. Boehm: ''On Cyclides in Geometric Modeling''.Computer Aided Geometric Design 7, 1990, p. 243–255.</ref>

* The difference/sum of the sphere's radius and the distance of the sphere's center (ellipse point) from one (but fixed) of the foci is constant. 

;Proof
The foci of the ellipse <math>\ E=(a\cos\varphi,b\sin\varphi,0)\ </math> are  <math>\ F_i=(\pm c,0,0)\ </math>. If one chooses <math>\ F_1=(c,0,0)\ </math> and calculates the distance <math>\ |EF_1|\ </math>, one gets <math>\ |EF_1|=a-c\cos\varphi\ </math>. Together with the radius of the actual sphere (see above) one gets <math>\ |EF_1|-r=a-d\ </math>.
<br>
Choosing the other focus yields: <math>\ |EF_2|+r=a+d\ . </math> 

Hence:

In the x-y-plane the envelopes of the circles of the spheres are two circles with the foci of the ellipse as centers and the radii <math>a\pm d</math> (see diagram).

=== Cyclide through 4 points on the x-axis ===

[[File:Zyklide-maxwell-abcd.svg|thumb|upright=1.4|Determining the design parameters a,b,c,d belonging to the given reals <math>x_1,x_2,x_3,x_4</math>]]
[[File:Zyklide-x1234.svg|thumb|upright=1.2|
top: ring cyclide with <math>x_1=3,x_2=1,x_3=-3,x_4=-7</math><br>
middle: horn cyclide with <math>x_1=1,x_2=3,x_3=-3,x_4=-7</math><br>
bottom: spindel cyclide with <math>x_1=3,x_2=-3,x_3=1,x_4=-7</math>]]

The Maxwell-property gives reason for determining a ring cyclide by prescribing its intersections with the x-axis:

''Given:'' Four points <math>x_1> x_2>x_3>x_4</math> on the x-axis (see diagram).

''Wanted:'' Center <math>m_0</math>, semiaxes <math>a,b</math>, linear eccentricity <math>c</math> and foci of the directrix ellipse  and the parameter <math>d</math> of the corresponding ring cyclide.

From the Maxwell-property one derives
:<math>2(a+d)=x_1-x_4 \; ,\quad \ 2(a-d)=x_2-x_3 \ .</math> 
Solving for <math>a,c</math> yields
:<math> a=\frac{1}{4}\left(x_1+x_2-x_3-x_4\right),\quad</math> <math> d=\frac{1}{4}\left(x_1-x_2+x_3-x_4\right)\ .</math>
The foci (on the x-axis) are  
:<math>f_1=\frac{1}{2}(x_2+x_3),\quad f_2=\frac{1}{2}(x_1+x_4)</math> and hence
: <math>c=\frac{1}{4}\left(-x_1+x_2+x_3-x_4\right)\; ,\quad \ b=\sqrt{a^2-c^2} \ </math> 
The center of the focal conics (ellipse and hyperbola) has the x-coordinate
:<math>m_0=\frac{1}{4}\left(x_1+x_2+x_3+x_4\right)\ .</math>

If one wants to display the cyclide with help of the parametric representation  above one has to consider the shift <math>m_0</math> of the center !

;Meaning of the order of the numbers <math>x_1,x_2,x_3,x_4</math>:
(The calculation above presumes <math>x_1> x_2>x_3>x_4</math>, see diagram.)

'''(H)''' Swapping <math>x_1,x_2</math> generates a horn cyclide.<br>
'''(S)''' Swapping <math>x_2,x_3</math>, generates a spindel cyclide.<br>
'''(H1)''' For <math>x_1=x_2=2</math> one gets a 1-horn cyclide.<br>
'''(R)''' For <math>x_2=x_3=0</math> one gets a ring cyclide touching itself at the origin.

=== Parallel surfaces ===

[[File:Zyklide-parallel-flaeche.svg|thumb|upright=1.2|Parallel surfaces of a ring cyclide with parameters <math>a=1,b=0.98, \ d=0.30, 0.45, 0.60</math>]]
By increasing or decreasing parameter <math>d</math>, such that the type does not change, one gets parallel surfaces (similar to [[parallel curve]]s) of the same type (see diagram).

=== Hyperbola as directrix ===

The second way to generate the ring cyclide as channel surface uses the focal hyperbola as directrix. It has the equation
[[File:Zyklide-kanalfl-hyp-1-2.svg|thumb|upright=2|Cyclide with two touching spheres with centers on the directrix hyperbola]]
:<math>\frac{x^2}{c^2}-\frac{z^2}{b^2}=1\ ,\ y=0\ .</math>
In this case the spheres touch the cyclide from outside at the second family of circles (curvature lines). To each arm of the hyperbola belongs a subfamily of circles. The spheres of one family enclose the cyclide (in diagram: purple). Spheres of the other family are touched from outside by the cyclide (blue).

Parametric representation of the hyperbola:
:<math>(\pm c\cosh\psi,0,b\sinh\psi) \ . </math>
The radii of the corresponding spheres are 
:<math>R(\psi)=a\cosh\psi\mp d\ .</math>

In case of a torus (<math>c=0</math>) the hyperbola degenerates into the axis of the torus.
[[File:Zyklide-maxwell-hyp.svg|thumb|upright=1|Maxwell-property for the hyperbola case. The hyperbola arm <math>H_+</math> is an equidistant curve of both the grey circles.]]
===Maxwell-property (hyperbola case)===
The foci of the hyperbola <math>\;H_\pm(\psi)=(\pm c\cosh\psi, 0, b\sinh\psi)\;</math> are <math>\;F_i=(\pm a,0,0)\;</math>. The distance of hyperbola point <math>\;H_+(\psi)=(c\cosh\psi,0,b\sinh\psi)\;</math> to the focus <math>F_1</math> is <math>\; |H_+F_1|=a\cosh\psi-c\; </math> and together with the sphere radius <math>R_+(\psi)=a\cosh\psi-d</math> one gets <math>\;|H_+F_1|-R_+=d-c\;</math>. Analogously one gets <math>\;|H_+F_2|-R_+=d+c\;</math> . For a point on the second arm of the hyperbola one derives the equations: <math>\;R_--|H_-F_1|=d-c\; ,\ R_--|H_-F_2|=d+c\; .</math>

Hence: 

In the x-z-plane the circles of the spheres with centers <math>H_\pm(\psi)</math> and radii <math>R_\pm(\psi)</math> have the two circles (in diagram grey) with centers <math>F_{1/2}=(\pm a,0,0)</math> and radii <math>d\mp c</math> as envelopes.
[[File:Zyklide-punkt-konstr.svg|thumb|upright=1.2|Point construction]]
=== Derivation of the parametric representation ===
==== Elliptic cyclide ====

The ellipse and hyperbola (focal conics) are the degenerated focal surfaces of the elliptic cyclide. For any pair <math>E, H</math> of points of the ellipse and hyperbola the following is true (because of the definition of a focal surface):
:1) The line <math>\overline{EH}</math> is a normal of the cyclide and
:2) the corresponding point <math>P</math> of the cyclide divides the chord <math>EH</math> with relation <math>r:R</math> (see diagram).

From the parametric representation of the focal conics and the radii of the spheres
:Ellipse: <math>\quad\ \ \ E(u)=(a\cos u,b\sin u,0),\quad r(u)=d-c\cos u \ \ ,</math>
:Hyperbola: <math>\quad H(v)=(\frac{c}{\cos v}, 0, b\tan v), \quad R(u)=\frac{a}{\cos v}-d\ .</math>

one gets the corresponding point <math>P</math> of the cyclide (see diagram):
:<math> \ P(u,v)=\frac{R(v)}{r(u)+R(v)}\;E(u) +\frac{r(u)}{r(u)+R(v)}\;H(v)\ .</math>
(For the unusual but convenient parametric representation of the hyperbola: see  [[hyperbola#Parametric equations|hyperbola]].)

Calculation in detail leads to the parametric representation of the elliptic cyclide given [[#Parametric and implicit representation|above]].

If one uses the parametric representation given in the article on channel surfaces, then, in general, only one family of parametric curves consists of circles.

==== Parabolic cyclide ====
[[File:Zyklide-parab-canal.svg|thumb|upright=1.2|Generation of a parabolic cyclide as channel surface]]
The derivation of the parametric representation for the parabolic case runs analogously:

With the parametric representations of the focal parabolas (degenerated focal surfaces) and the radii of the spheres:
:<math>P_1(u)= \left(\frac p 2 (1-u^2),pu,0\right), \quad r_1=\frac p 2 (1-k +u^2)\; </math>
:<math>P_2(v)= \left(\frac p 2 v^2,0,pv\right), \qquad \qquad r_2=\frac p 2 (k+v^2) \ </math>

one gets 
:<math> \ P(u,v)=\frac{r_2(v)}{r_1(u)+r_2(v)}\;P_1(u) +\frac{r_1(u)}{r_1(u)+r_2(v)}\;P_2(v) </math>
which provides the parametric representation above of a parabolic cyclide.

== Dupin cyclides and geometric inversions ==

An advantage for investigations of cyclides is the property:
: '''(I):''' Any Dupin cyclide is the image either of a [[right circular cylinder]] or a [[Cone|right circular double cone]] or a [[Torus|torus of revolution]] by an [[Inversive geometry|inversion]] (reflection at a sphere).

The inversion at the sphere with equation <math>x^2+y^2+z^2=R^2</math> can be described analytically by:
: <math> (x,y,z) \rightarrow \frac{R^2\cdot(x,y,z)}{x^2+y^2+z^2} \ .</math>
The most important properties of an inversion at a sphere are:
# Spheres and circles are mapped on the same objects.
# Planes and lines containing the origin (center of inversion) are mapped on themselves.
# Planes and lines ''not'' containing the origin are mapped on spheres or circles passing the origin.
# An inversion is [[Involution (mathematics)|involutory]] (identical with the inverse mapping).
# An inversion preserves angles.

One can map arbitrary surfaces by an inversion. The formulas above give in any case parametric or implicit representations of the image surface, if the surfaces are given parametrically or implicitly. In case of a parametric surface one gets:
:<math>(x(u,v),y(u,v),z(u,v)) \rightarrow \frac{R^2\cdot(x(u,v),y(u,v),z(u,v))}{x(u,v)^2+y(u,v)^2+z(u,v)^2} \ .</math>
[[File:Zyklide-iv-zyl1.svg|upright=1|thumb|ring cyclide generated by an inversion of a cylinder at a sphere (magenta)]]
[[File:Zyklide-iv-zyl2.svg|upright=1|thumb|parabolic ring cyclide generated by an inversion of a cylinder containing the origin]]
[[File:Zyklide-iv-keg1.svg|upright=1|thumb|horn cyclide generated by an inversion of a cone]]
[[File:Zyklide-iv-tor1.svg|upright=1|thumb|ring cyclide generated by an inversion of a torus]]

But: Only in case of right circular cylinders and cones and tori of revolution one gets Dupin cyclides and vice versa.

=== Example cylinder ===
a) Because lines, which do not contain the origin, are mapped by an inversion at a sphere (in picture: magenta) on circles containing the origin the image of the cylinder is a ring cyclide with mutually touching circles at the origin. As the images of the line segments, shown in the picture, there appear on line circle segments as images. The spheres which touch the cylinder on the inner side are mapped on a first pencil of spheres which generate the cyclide as a canal surface. The images of the tangent planes of the cylinder become the second pencil of spheres touching the cyclide. The latter ones pass through the origin.<br />
b) The second example inverses a cylinder that contains the origin. Lines passing the origin are mapped onto themselves. Hence the surface is unbounded and a parabolic cyclide.

=== Example cone ===
The lines generating the cone are mapped on circles, which intersect at the origin and the image of the cone's vertex. The image of the cone is a double horn cyclide. The picture shows the images of the line segments (of the cone), which are circles segments, actually.

=== Example torus ===
Both the pencils of circles on the torus (shown in the picture) are mapped on the corresponding pencils of circles on the cyclide. In case of a self-intersecting torus one would get a spindle cyclide.

;Villarceau circles
[[File:Zyklide-torus-vill.svg|thumb|upright=0.8|Ring cyclide with Villarceau-circles]]
Because Dupin ring-cyclides can be seen as images of tori via suitable inversions and an inversion maps a circle onto a circle or line, the images of the [[Villarceau circle]]s form further two families of circles on a cyclide (see diagram). 

;Determining the designparameters <math>a,b,c,d</math>
The formula of the inversion of a parametric surface (see above) provides a parametric representation of a cyclide (as inversion of a torus) with circles as parametric curves. But the points of a parametric net are not well distributed. So it is better to calculate the design parameters <math>a,b,c,d</math> and to use the parametric representation above:

[[File:Zyklide-tor-invers-abcd.svg|thumb|upright=1.5|Cyclide (blue) as image by an inversion of a torus (black) at the unit sphere (red)]]
''Given:'' A torus, which is shifted out of the standard position along the x-axis. 
Let be <math>x_1, x_2,x_3,x_4</math> the intersections of the torus with the x-axis (see diagram). All not zero. Otherwise the inversion of the torus would not be a ring-cyclide.<br>
''Wanted:'' semi-axes <math>a,b</math> and linear eccentricity <math>c</math> of the ellipse (directrix) and parameter <math>d</math> of the ring-cyclide, which is the image of the torus under the inversion at the unitsphere.

The inversion maps <math>x_i</math> onto <math>\tfrac{1}{x_i}</math>, which are the x-coordinates of 4 points of the ring-cyclide (see diagram). From section [[#Cyclide through 4 points on the x-axis|Cyclide through 4 points on the x-axis]] one gets 
:<math> a=\frac{1}{4}\left(\frac{1}{x_1}+\frac{1}{x_2}-\frac{1}{x_3}-\frac{1}{x_4}\right),\quad</math> <math> d=\frac{1}{4}\left(-\frac{1}{x_1}+\frac{1}{x_2}-\frac{1}{x_3}+\frac{1}{x_4}\right)\ </math> and
:<math>c=\frac{1}{4}\left(\frac{1}{x_1}-\frac{1}{x_2}-\frac{1}{x_3}+\frac{1}{x_4}\right)\; ,\quad \ b=\sqrt{a^2-c^2} \ </math> 
The center of the focal conics has the x-ccordinate
:<math>m_0=\frac{1}{4}\left(\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\frac{1}{x_4}\right)\ .</math>

==Separation of variables==
Dupin cyclides are a special case of a more general notion of a cyclide, which is a natural extension of the notion of a [[quadric surface]]. Whereas a quadric can be described as the zero-set of second order polynomial in Cartesian coordinates (''x''<sub>1</sub>,''x''<sub>2</sub>,''x''<sub>3</sub>), a cyclide is given by the zero-set of a second order polynomial in (''x''<sub>1</sub>,''x''<sub>2</sub>,''x''<sub>3</sub>,''r''<sup>2</sup>), where
''r''<sup>2</sup>=''x''<sub>1</sub><sup>2</sup>+''x''<sub>2</sub><sup>2</sup>+''x''<sub>3</sub><sup>2</sup>. Thus it is a quartic surface in Cartesian coordinates, with an equation of the form:
:<math>
A r^4 + \sum_{i=1}^3 P_i  x_i r^2 + \sum_{i,j=1}^3 Q_{ij}  x_i  x_j + \sum_{i=1}^3 R_i  x_i + B = 0
</math>
where ''Q'' is a 3x3 matrix, ''P'' and ''R'' are a 3-dimensional [[vector (geometric)|vectors]], and ''A'' and ''B'' are constants.<ref>{{Harvnb|Miller|1977}}</ref>

Families of cyclides give rise to various cyclidic coordinate geometries.

In [[Maxime Bôcher]]'s 1891 dissertation, ''Ueber die Reihenentwickelungen der Potentialtheorie'', it was shown that the [[Laplace equation]] in three variables can be solved using separation of variables in 17 conformally distinct quadric and cyclidic coordinate geometries.  Many other cyclidic geometries can be obtained by studying R-separation of variables for the Laplace equation.<ref>{{Harvnb|Moon|Spencer|1961}}</ref>

==See also==
* [[:de:Inversion (Geometrie) |Inversion of curves and surfaces (German)]]
== External links ==
* [http://www.mathcurve.com/surfaces.gb/cycliddedupin/cyclidededupin.shtml  Dupin cyclides on MathCurve]

==Notes==
{{reflist}}

==References==
* {{citation | last = Cecil | first = Thomas E. | title = Lie sphere geometry | publisher = Universitext, Springer-Verlag|place= New York | year = 1992|isbn =978-0-387-97747-8}}.
* {{citation|last=Eisenhart|first= Luther P.|chapter=§133 Cyclides of Dupin|title= A Treatise on the Differential Geometry of Curves and Surfaces|place=  New York|publisher= Dover|pages=312–314|year= 1960}}.
* {{citation | title = Geometry and the Imagination | author1-link = David Hilbert|first1=David|last1=Hilbert|first2=Stephan |last2=Cohn-Vossen |authorlink2=Stephan Cohn-Vossen| year = 1999 | publisher = American Mathematical Society | isbn= 0-8218-1998-4}}.
* {{citation | title = Field Theory Handbook: including coordinate systems, differential equations, and their solutions | first1=Parry|last1= Moon |first2=Domina Eberle|last2= Spencer| year = 1961 | publisher = Springer | isbn=0-387-02732-7}}.
* {{citation| last1=O'Connor|first1= John J.|first2=Edmund F.|last2= Robertson|chapter-url=http://www-groups.dcs.st-and.ac.uk/~history/Biographies/Dupin.html|chapter=Pierre Charles François Dupin|title=[[MacTutor History of Mathematics archive]]|year=2000}}.
* {{citation|last=Pinkall|first=Ulrich|author-link=Ulrich Pinkall|chapter=§3.3 Cyclides of Dupin|title=Mathematical Models from the Collections of Universities and Museums|title-link= Mathematical Models: From the collections of universities and museums |editor=G. Fischer|place= Braunschweig, Germany|pages= 28–30|year= 1986|publisher=Vieweg}}.
* {{citation | first = Willard | last = Miller | year = 1977 | title = Symmetry and Separation of Variables}}.
* A. Cayley (1873) "On the cyclide", ''Quarterly Journal of Pure and Applied Mathematics'' 12: p.&nbsp;148–163.
* V. Chandru, D. Dutta, C.M. Hoffmann (1989) "On the geometry of Dupin cyclides", ''The Visual Computer.'' (5), p.&nbsp;277–290.
* C. Dupin (1822) ''Applications de Geometrie et de Mechanique.'' Bachelier, Paris.
* F. Klein, W. Blaschke (1926) ''Vorlesungen Über Höhere Geometrie.'' Springer-Verlag, {{ISBN|978-3-642-98494-5}}, p.&nbsp;56.
* J. C. Maxwell (1868) "On the cyclide", [[Quarterly Journal of Pure and Applied Mathematics]] 9: p.&nbsp;111–126.
* M. J. Pratt (1989) ''Cyclide Blending in Solid Modelling.'' In: Wolfgang Strasser, Hans-Peter Seidel (Hrsg.): ''Theory and Practice in Geometric Modelling.'' Springer-Verlag, {{ISBN|0-387-51472-4}}, p.&nbsp;235.
* Y. L. Srinivas, V. Kumar, D. Dutta (1996) "Surface design using cyclide patches", ''Computer-Aided Design'' 28(4): 263–276.
* [[Mabel M. Young]] (1916) "Dupin's cyclide as a self-dual surface", [[American Journal of Mathematics]] 38(3): 269–286

==External links==
{{commons category|Dupin cyclide}}
* {{Mathworld|Cyclide|Cyclide}}
* E. Berberich, M. Kerber: [https://web.archive.org/web/20150924142245/http://people.mpi-inf.mpg.de/~mkerber/bk-arrangements-eurocg08.pdf ''Arrangements on Surfaces of Genus One: Tori and Dupin Cyclides.'']

[[Category:Surfaces]]
[[Category:Eponyms in geometry]]