The inverses of Jacobi elliptic functions possess an apparently-non-crucial property: they provide almost-everywhere-conformal maps on a hemisphere onto a torus and so, onto a parallelogram. Thus, they produce map projections on the sphere generalizing the famous quincuncial projection of Charles S. Peirce. Besides providing a general practical definition of -uncial map and proving that all the considered inverse elliptic functions are -uncial, we give operative handy formulas to calculate these maps. To the best of our knowledge, these useful formulas have not been all together published before, except for Pierce projection. We look forward to their numerical implementation. By the way, we also classify the resulting map projections according the number of singularities.
In the theory of functions of a complex variable, a function is called elliptic if it is meromorphic doubly periodic, cf. Lang [11]. These fundamental properties, arised from the addition theorems, were absolutely necessary to Abel [1] and Jacobi [9] in their foundational works on the subject. During the nineteenth century, other related properties of elliptic functions were discovered little by little. Here we concentrate our efforts in a property that proves to be a key ingredient to conceive geographical projections with great utility and gracefulness. We refer specifically to a problem addressed in a paper by Richelot [19] and in Durege’s well-known textbook [4]. A step forward was given by Peirce in his laconic paper [15], which was unriddled a little bit later by Pierpont [16].
In what follows we generalize Peirce projection through the notion of -uncial function and prove that all the inverse functions of Jacobi elliptic functions are -uncial with . The number denotes the number of singularities of the map. Within this framework, Peirce projection results to be the unique function of this type which has . Certainly, the adjective quincuncial, introduced by Peirce [15], comes from the Latin word quincuncialis, relating to the noun quincux, an arrangement of five objects in a square or rectangle, one at each corner and one in the middle. In simple words, quincux is the usual arrangement of the five-spot face on six-sided dice, playing cards, and dominoes. Quincux also denotes a type of church from the second Byzantine Golden Age: a five-domed temple based on the domed cross element. The idea is interesting to the geometer since it evokes an almost-everywhere-conformal transformation of a hemisphere onto a cuadrilateral with five singular points (one in the middle and one at each corner), a key feature of Peirce quincuncial projection.
The importance of the -uncial property is both theoretical and practical. On the one hand, being a characterizing property of the inverses of Jacobi elliptic functions, this property can be used to give an alternative definition of elliptic function, as far as the inverse of a function carries essentially the same information of the original function. This is particularly meaningful for elliptic functions because they are intrinsically meromorphic and the -uncial property is about the number of singular points. On the other hand, -uncial projections provide almost-everywhere conformal maps on the sphere onto a plane region and this is helpful to numerous applications. According to Peirce [15], the quincuncial map is useful for “meteorological, magnetological and other purposes”. Remarkably enough, he does not mention “cartographic purposes”. In general, cartographers indeed find the quincuncial projection interesting (cf. [12]), although completely impractical and this often happens to be the case with map projections based upon complex functions (and the associated computational intricateness). However, Peirce projection has shown to be advantageous because it tessellates the plane, that is, the hemispheres mosaic or tile indefinitely to cover the whole plane in an almost-everywhere-conformal manner. Being so, it faithfully displays nearly every locus on the sphere near to its neighboring loci (not all map projections share this feature). This explains the convenience of the map to earth magnetology and metereology, the primary interest of Peirce in 1879. In 1946, the U.S.A. Coast and Geodetic Survey [23] use the projection to draw a world map of air routes. In astronomy, Taylor and Bell [24] have graphed suitable quincuncial maps of the surface of the earth showing umbral limit lines for eclipses. And just as importantly, a group of biologists [17] have recently (2019) plotted on a Peirce chart the geographic distribution of Eurasian and American samples of genetic human male data in order to formulate hypotheses concerning the history “that shaped the present-day biological and cultural diversity” of the Americas. In three-dimensional computer graphic design, this mapping has also been applied to display spherical panoramas (e.g. [6]). Perhaps most remarkably, the projection serves as well for developing generalized longitudes and latitudes allowing to display the complete surface of a highly irregular body, such a comet. The literature on these appropriate maps for a comet surface is interesting and we refer the interested reader to [7, 8].
That an inverse of an elliptic function can supply a plane conformal representation of a sphere might be suggested by the theory of general Schwarz-Christoffel mappings, cf. Driscoll and Trefethen [3]. However, here we are interested in maps of Riemann surfaces with boundary [10, 5] that are conformal everywhere, except in a finite number of points. By the way, the striking 3D graphs of the complex modulus of elliptic functions exhibit certain cornettos about the poles and so, they insinuate us to consider the projections of these graphs onto their domains.
This paper is organized as follows. After reviewing some preparatory facts on Jacobi elliptic functions and setting up the basic notations in Section 2, we describe a general procedure to compute the values of Jacobi elliptic inverse functions in Section 3. Along the way, we introduce innovative convenient notations to write the complex modulus and argument of an arbitrary Jacobi elliptic function and, as a noteworthy outcome, we obtain a complete list of practical formulas for the inverses. In Section 4, we discuss the novel notion of -uncial map. Among other things, we prove that Peirce projection is just but a particular case of these maps; in fact, it is the only -uncial map with . Finally, in Section 5, we draw conclusions from the previous results and explain what we should expect from the graphical representations of these general mappings.
Jacobi elliptic functions
By their definition, elliptic functions map meromorphically a torus to the Riemann sphere . We furnish and with their usual conformal structures. That is to say, is the one-point compactification of the complex plane , provided with the atlas containing the stereographic projection; and , where is a lattice of rank 2, related to some fundamental parallelograms and to a fundamental pair of periods. Briefly, we endow both and with their customary structures of Riemann surfaces.
In his Fundamenta nova [9], Jacobi himself introduced a handful of elliptic functions as a sort of inverses to incomplete elliptic real integrals of the first kind
where is a parameter called elliptic modulus, cf. [13, 22]. is odd and increasing and its inverse function is the elliptic amplitude. This amplitude is likewise an odd increasing function, which increases by when increases by
is known as the complete elliptic integral of Legendre’s first kind with elliptic modulus . Given an elliptic modulus , its complementary modulus is defined by . The complete elliptic integral associated with the complementary modulus of is written .
Jacobi elliptic complex functions are the function composition of the amplitude followed by a trigonometric or a trigonometric-like function. For example, the sinus amplitudinis, sine of the amplitude or elliptic sine, is just but the function . After the improved notations of Gudermann and Glaisher, the twelve Jacobi elliptic functions are nowadays denoted by , where and are any of the letters (cosine), (sine), and (delta). The symbol stands for the unity. The basic Jacobi elliptic functions are
The remaining nine Jacobi elliptic functions are defined by
After several complicate attempts using only the addition formulas, Jacobi [9] discovered that the right setting for the study of these functions should be based on the notion of theta function. In terms of the modern Neville theta functions and , Jacobi elliptic functions can be certainly defined as
We refer the interested reader to the details in the classic book of Neville [14].
With the aid of theta functions, it is not hard to see that the periods, zeros and poles of Jacobi elliptic functions are those shown in Table 1. Each pair of integer values determines a unique zero and a unique pole. The periods, zeros and poles for the remaining Jacobi elliptic functions can be easily found from Table 1 and the definition of the functions.
Fundamental periods, zeros and poles of Jacobi elliptic functions.
Function
Fundamental periods
Zeros
Poles
cn
sn
dn
cs
cd
sd
Now, Jacobi elliptic functions defined on a fundamental parallelogram spanned by a pair of periods (or a torus, if we like) on the sphere are not injective. Such fundamental regions contains two zeros and two poles and, more generally, any of these functions has multiplicity two or it is “two-to-one”. Nonetheless, each fundamental parallelogram can be split into two smaller parallelograms in such a way that the restriction of the elliptic function to any of them is one-to-one. In the interior of such smaller parallelograms lie exactly one pole and one zero (singular points). For our purposes here, we will even divide each of such smaller parallelograms into two lesser parallelograms, each containing only one singular point (pole or zero). From now on, with no harm to the argumentation, we shall deal exclusively with these injective restrictions of Jacobi elliptic functions. Stated briefly, the domain of such restricted functions is a “small enough” open parallelogram and the range is an open hemisphere of the Riemann sphere. In addition, the only singular point lies always at the center of the parallelogram. Due to the symmetries introduced by the group generated by the fundamental periods, there are many possible choices of such open parallelograms.
In the next section we will make use of many famous identities of the Jacobi elliptic functions, or shortly, elliptic identities. Among them, the most important are perhaps the celebrated addition formulas such as
There are as well essential relations between squares of the functions or Pythagorean-like identities, for example,
For the inversion of odd elliptic functions we will resort to some, less-obvious, identities. They comprise the so-called real Jacobi transformations [9, 20], for instance,
Also, an important example of an imaginary Jacobi transformation will be given and used in the next section. By the way, the functions and their multiplicative inverses (reciprocal functions) are even. The remaining Jacobi elliptic functions are odd.
In the subsequent sections, we will use freely all the identities we could prove, or find in the literature. In particular, we apply those elliptic identities in Chapter 16 of Abramowitz and Stegun [2]. We have also found serviceable the 74 identities in Weisstein [25] and those in Reinhardt and Walker [18].
Inverse elliptic functions
Up to this point, we have denoted the elliptic functions by the symbols . Since from now it is necessary to specify the elliptic modulus , we will employ the notation . The idea behind the following procedure is to generalize the methods of Pierpont [16] and Solanilla, Oostra and Yañez [21], which elucidate the Peirce’s laconic paper [15]. Simply put, we extend a well-known formula for the function to an arbitrary Jacobi elliptic function .
In the first place, we need to introduce a new notation for the complex modulus and the principal argument of these functions. The new symbols consist of a principal (Greek) letter and four attributes ascribing qualities of the principal letter. Such principal letter are, as it is usual, (rho) and (theta), which is distinct from a theta function. The attributes are written counterclockwise beginning from the bottom right corner of the principal letter. They specify a Jacobi elliptic function, an (integer) exponent or power, an elliptic modulus and the point at which the function is evaluated. With reference to the common notations of Complex Analysis, they translate as
The exponent of the principal argument Arg is omitted because it will always be the unity.
Inverses of the even elliptic functions
By generalizing the method in Solanilla, Oostra and Yañez [21], we notice first that any inverse function can be expressed as a function of the complex moduli of some Jacobi elliptic functions.
.
Let be an even Jacobi elliptic function and write its inverse function as . Then, the real and imaginary parts of satisfy the conditions in Table 2, according to the respective inverse elliptic function.
Even functions: dependence of on the complex moduli of some elliptic functions
Function
Condition on
Condition on
Proof..
The function serves as example of the general procedure. By virtue of the addition formulas,
Now, for the involved functions, . Finally, again by an addition theorem,
∎
.
By these means, the problem of finding the complex inverses is reduced to the restrictions of the elliptic functions to the real axis. Accompanying relations hold for the other even inverse elliptic functions. The real inverses in the formulas for the imaginary part ‘’ of ‘’ have elliptic modulus , instead of .
In a second step, for each function in Table 2, the complex moduli necessary to compute and depend exclusively on the complex modulus and argument of . This makes the formulas in Table 2 useful.
.
The complex moduli required to compute the inverse functions of each even Jacobi elliptic function can be expressed in terms of the complex modulus and argument of its corresponding Jacobi elliptic function. The expressions accomplishing this goal are given in Table 3.
Even functions: the complex moduli in Table 2 are properly related to the elliptic function
Function
Required moduli
Proof..
Again for , we start with , i.e., . By multiplying this equation by its complex conjugate,
In like manner, or . Therefore,
∎
To sum up, we arrange together the previous results to obtain a set of appropriate formulas.
.
The inverses of the even Jacobi elliptic functions can be calculated by the real-variable formulas in Table 4.
Formulas for the inverse functions of the even Jacobi elliptic functions
Function
Real and imaginary parts
Inverses of the odd elliptic functions
In order to find the inverses of these functions we rely on some real and imaginary Jacobi transformations [9, 20] and so, some extra work must be made to do so. This method constitute one of the central contributions of this paper. We certainly simplify the early efforts of Richelot [19] and Durège [4]. As before, we first reduce the calculation of to some expressions involving real Jacobi elliptic functions and certain complex moduli.
.
Let be an odd Jacobi elliptic function with inverse . Then, and fulfill the conditions in Table 5.
Odd functions: dependence of on the complex moduli of some elliptic functions
Function
Condition on
Condition on
Proof..
The three inverses are respectivelly found by the following Jacobi transformations:
The inverse of sn is clear by the previous results on the even Jacobi functions. For the case of , we use in the formulas of , i.e.,
The treatment of is completely analogous. ∎
This time we have to be more careful to establish the right relations among the complex moduli of the elliptic functions.
.
The complex moduli needed to compute the inverse functions of the odd Jacobi elliptic functions depend on the complex modulus and argument of its respective Jacobi elliptic function. The formulas appear in Table 6.
Odd functions: the complex moduli in Table 5 are correctly related to the elliptic function
Function
Required moduli
Auxiliary moduli
, where
, where
Proof..
We consider each of the functions.
On the one hand, the Jacobi transformation for sn in polar form provides
On the other hand, the addition formula for sn together with the well-known values of the elliptic functions at yield
The auxiliary moduli , are easily obtained with the help of the Pythagorean-like identities.
First, the imaginary Jacobi transformation for sc is written as
Hence,
We start with
and employ the standard Pythagorean-like identities.
∎
We recapitulate the results for the odd functions in one result.
.
The real and imaginary parts of the inverses to the odd Jacobi elliptic functions are given by the formulas in Table 7.
Formulas for the inverse functions of the odd Jacobi elliptic functions
Function
Real and imaginary parts
Calculation of the inverse elliptic functions
For the sake of rendering usable the formulas en Tables 4 and 7, we must give procedures to compute the restrictions of the inverse Jacobi elliptic functions to the real line. There are at least two ways of accomplishing this purpose: either we can use the Jacobi original definitions or we can utilize certain integral formulas.
When we write Legendre elliptic integral of the first kind in the form , the amplitude is simply and we get at once
With a little more effort, by using elliptic identities, we derive
Many of the most popular programming languages for Mathematics come with a routine to approximate function and the usual trigonometric inverses. This fact makes this approach easy to implement.
Still and all, there are also direct integral formulas to calculate these inverse functions, v. [18]. They are presented, together with the former formulas, in Table 8. As long as the integral formulas in Table 8 are less known, their implementation may demand some extra work.
Two ways for computing the inverses of Jacobi elliptic functions
Jacobi definition
Direct integral formula
Domain
Stereographic projection
With regard to the conformal structure for an open hemisphere of the Riemann sphere , we need to introduce the stereographic projection in our formulas. Let denote the “North pole”. Once we parametrize the sphere by means of a geographic system, i.e., by local coordinates (latitude) and (longitude or azimuth), the translation permits to write the restricted stereographic projection , , under the suitable form (cf. [15, 21])
Conformality and singular points
For a particular Jacobi elliptic function pq, we have . Thus, in each of the preceding formulas in Tables 4 and 7 we must substitute
to get the correct values. In other words, we define the map projection associated with a specific Jacobi elliptic function as the composite function of the inverse Jacobi elliptic function following the stereographic projection . In our setting, these map projections are bijections applying a punctured open hemisphere onto a puncture open parallelogram where a Jacobi elliptic function is injective. Due to the azimuthal symmetry of the stereographic projection, the magnification of is given almost everywhere as
Without trouble,
and, happily, the derivatives of the basic Jacobi elliptic functions are widely known:
As a result,
.
Table 9 displays the magnification reciprocals of the map projections associated with Jacobi elliptic functions.
Magnification reciprocals of the map projections
Reciprocal of
Proof..
The expressions for the basic elliptic functions and dn follow forthrightly from Tables 3 and 6. So, we address the remaining functions.
We begin with the Pythagorean-like identities
Now, we rewrite these expressions in polar form, conjugate them and multiply each of them by its complex conjugate to obtain
This time the starting point is
Therefore,
The relations
produce
We notice the use of the formulas for the derivatives of Jacobi elliptic functions. ∎
One of the remarkable features of the map projections involved is the way they extend to the boundary of the punctured hemisphere. In the process, only a finite number of singularities is introduced.
.
The extension of each map projection – associated with a Jacobi elliptic function – to its domain boundary is conformal almost everywhere and its singularities are located according to the following rules
The “North pole” is mapped to the center of the parallelogram.
In the rest of the boundary of the hemisphere, i.e. on the “Equator”, an even number (0, 2 or 4) of singularities may appear. These singularities have order 2, where they do exist.
Proof..
The image of the “North pole” is evident. On the “Equator”, and so, . As the magnifications of map projections associated to the inverse elliptic functions and have a factor of the form , they blow up at . For the function , the corresponding factor is and the map projection blows up at . Now we seek other singularities for each :
The magnification has a factor . Hence, we have singularities at precisely when .
There are no other possible singular points. We summarize these results in Table 10. ∎
Singularities of the map projections on the “Equator” ()
Singularities independent of
Feasible singularities depending on
;
None
None
None
None
none
None
The functions obtained from Jacobi elliptic functions evaluated at the limiting elliptic moduli (whenever they exist and make sense) can be called improper. Jacobi elliptic functions defined in Section 2 above are all proper. By admitting improper Jacobi elliptic functions, we find new singularities:
The factor with leads to a singularity precisely when , .
The factor vanishes at if and only if , .
holds when together with , .
In what follows, we shall not refer to these improper cases.
-uncial maps
The foregoing findings can be put together into one all-embracing notion.
.
Let denote a positive integer. A map on the Riemann sphere into the complex plane is -uncial if it satisfies the following conditions:
The restriction of of to a punctured open hemisphere , being its pole, is everywhere conformal onto an open parallelogram , punctured at its center .
The extension of to the boundary of its domain is conformal except at points. is always a singular point. Each singularity has order 2.
When , an -uncial map is called uncial; when , teruncial; and when , quincuncial.
Corollary 1 implies at once the following consequence.
.
The valid map projections associated to the inverse functions of the Jacobi elliptic functions are -uncial with .
At long last, we can classify the maps projections associated to the proper inverse Jacobi elliptic functions according to the number of singularities:
or uncial. , . or teruncial. , ; , ; , ; , ; , , and , . or quincuncial. Peirce quincuncial projection:, .
.
Peirce quincuncial projection is the only quincuncial map projection obtained from a proper Jacobi elliptic function.
Fundamental rectangle of with .
Concluding remarks
The practical purpose of devising map projections based on Jacobi elliptic functions leads naturally to an important theoretical notion, namely, that of -uncial map. This is a property may be used to characterize (Jacobi) elliptic functions. Indeed, to what extend does this -uncial property determine a (Jacobi) elliptic function? Apropos this question, how can this approach to (Jacobi) elliptic functions be of use to describe the Riemann surfaces of the inverse (Jacobi) elliptic functions? We believe these two questions are not too challenging and can be answered within the framework of the theory of functions of a complex variable. Further, the formulas for in Tables 4 and 7 can easily be simplified to describe the images of the parallels of latitude and the meridians of the sphere.
More attractively for applications, does every map projection having a -uncial property yield a tessellation of the plane similar to Peirce’s tessellation in [15]? A preliminary analysis reveals that the answer is positive. In Fig. 11
This figure has been adapted (cutting, changing colors and drawing a small rectangle) from the file Sn-k-08.png in Wikimedia Commons. It is described as the “elliptic Jacobi function, sn, corresponding to , generated using a version of the domain coloring method”. It is published here under the Creative Commons Attribution-Share Alike 4.0 International license.
we have sketched the fundamental rectangle for sinus amplitudinis with elliptic modulus . The image of a hemisphere correspond to the area framed by the small rectangle (the half part of a region where the function is one-to-one). A domain coloring technique with contour lines of the modulus and the argument exposes the images of the lines of constant latitude and longitude on this rectangle. Then, a map of a hemisphere can be drawn in this small rectangle and consequently we obtain a tessellation similar to Peirce’s. It may be noted in passing that this map is teruncial, as it has been explained above.
Some other interesting questions arise from these results. They will be reported at a later date.
Footnotes
Acknowledgments
This research was partially funded by the Comité Central de Investigaciones, Universidad del Tolima, Ibagué, Colombia, grant number 60120. We thank the Facultad de Ciencias, Universidad del Tolima, for its logistic support to craft the manuscript.
Conflict of interest
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the analyses; in the writing of the manuscript, and in the decision to publish the results.
References
1.
AbelN.H., Recherches sur les fonctions elliptiques, Journal für die Reine und Angewandte Mathematik, Herausgeben vor Crelle2–3 (1827–1828), 263–388.
2.
AbramowitzM. and StegunI.A. (Eds.), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, 1972.
3.
DriscollT.A. and TrefethenL.N., Schwarz-Christoffel Mapping, Cambridge University Press, 2002.
4.
DurègeH., Theorie der elliptischen Funktionen, Druck und Verlag von B.G. Teubner, 1887.
5.
FarkasH.M. and KraI., Riemann Surfaces, Springer-Verlag, 1980.
6.
GermanD.M.d’AngeloP. et al., New methods to project panoramas for practical and aesthetic purposes, Computational Aesthetics in Graphics, Visualization, and Imaging (2007), 15–22.
7.
GriegerB., Quincuncial adaptive closed kohonen (QuACK) map for the irregularly shaped comet 67P/Churyumov-Gerasimenko, Astronomy & Physics, special issueA1 (2019).
8.
GriegerB. and VincentJ.-B., A closed self-organizing map of Chury, European Planetary Science Congress 2018, Vol. 12 EPSC2018-298 (2018).
9.
JacobiC.G.J., Fundamenta nova theoriae functionum ellipticarum, Sumptibus fratrum Bornträger, 1829.
10.
JostJ., Compact Riemann Surfaces. An Introduction to Contemporary Mathematics, Springer, 1997.
LeeL.P., Some conformal projections based on elliptic functions, Geographical Review55 (1965), 563–580.
13.
LegendreA.M., Traité des fonctions elliptiques et des intégrales eulériennes. Tome premier, Imprimerie de Huzard-Courcier, 1825.
14.
NevilleE.H., Jacobian Elliptic Functions, Oxford at the Clarendon Press, 1944.
15.
PeirceC.S., A quincuncial projection of the sphere, American Journal of Mathematics2 (1879), 394–396.
16.
PierpontJ., Note on the C.S. Peirce’s Paper on ‘A Qincuncial Projection of the Sphere’, American Journal of Mathematics18 (1896), 145–152.
17.
PinottiJ.BergströmA. et al., Y chromosome sequences reveal a short beringian standstill, rapid expansion, and early population structure of native american founders, Current Biology29 (2019), 149–157.
18.
ReinhardtW.P. and WalkerP.L., Digital Library of Mathematical Functions. Chapter 22. Jacobian Elliptic Functions (2019, April 21), Retrieved from https://dlmf.nist.gov/22.
19.
RichelotF.J., Darstellung einer beliebigen gegebenen Größe durch 𝑠𝑖𝑛𝑎𝑚(u+w,k), Journal für die Reine und Angewandte Mathematik45 (1853), 225–232.
20.
SolanillaL., Las transformaciones elípticas de Jacobi, Sello Editorial Universidad del Tolima, 2014.
21.
SolanillaL.OostraA. and YáñezJ.P., Peirce quincuncial projection, Revista Integración34 (2016), 23–38.
22.
SolanillaL.TamayoA.C. and ParejaG., Integrales elípticas con notas históricas, Sello Editorial Unversidad de Medellín, 2010.
23.
StanleyA.A., A quincuncial projection of the world, Surveying and Mapping6(1) (1946), 19.
24.
TaylorD.B. and BellS.A., Astronomical applications of the quincuncial map projection, Astronomy & Geophysics54 5 (2’13), 5.13–5.15.
25.
WeissteinE.W., Wolfram MathWorld. Jacobi Elliptic Functions (2019, April 21), Retrieved from http://mathworld.wolfram.com/JacobiEllipticFunctions.html.