## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Flat tori in three-dimensional space and convex integration

Edited by Yakov Eliashberg, Stanford University, Stanford, CA, and approved March 9, 2012 (received for review November 9, 2011)

## Abstract

It is well-known that the curvature tensor is an isometric invariant of *C*^{2} Riemannian manifolds. This invariant is at the origin of the rigidity observed in Riemannian geometry. In the mid 1950s, Nash amazed the world mathematical community by showing that this rigidity breaks down in regularity *C*^{1}. This unexpected flexibility has many paradoxical consequences, one of them is the existence of *C*^{1} isometric embeddings of flat tori into Euclidean three-dimensional space. In the 1970s and 1980s, M. Gromov, revisiting Nash’s results introduced convex integration theory offering a general framework to solve this type of geometric problems. In this research, we convert convex integration theory into an algorithm that produces isometric maps of flat tori. We provide an implementation of a convex integration process leading to images of an embedding of a flat torus. The resulting surface reveals a *C*^{1} fractal structure: Although the tangent plane is defined everywhere, the normal vector exhibits a fractal behavior. Isometric embeddings of flat tori may thus appear as a geometric occurrence of a structure that is simultaneously *C*^{1} and fractal. Beyond these results, our implementation demonstrates that convex integration, a theory still confined to specialists, can produce computationally tractable solutions of partial differential relations.

A geometric torus is a surface of revolution generated by revolving a circle in three-dimensional space about an axis coplanar with the circle. The standard parametrization of a geometric torus maps horizontal and vertical lines of a unit square to latitudes and meridians of the image surface. This unit square can also be seen as a torus; the top line is abstractly identified with the bottom line and so are the left and right sides. Because of its local Euclidean geometry, it is called a square flat torus. The standard parametrization now appears as a map from a square flat torus into the three-dimensional space having as its image a geometric torus. Although natural, this map distorts the distances: The lengths of latitudes vary whereas the lengths of the corresponding horizontal lines on the square remain constant.

It was a long-held belief that this defect could not be fixed. In other words, it was presumed that no isometric embedding of the square flat torus—a differentiable injective map that preserves distances—could exist into three-dimensional space. In the mid 1950s Nash (1) and Kuiper (2) amazed the world mathematical community by showing that such an embedding actually exists. However, their proof relies on an intricated construction that makes it difficult to analyze the properties of the isometric embedding. In particular, these atypical embeddings have never been visualized. One strong motivation for such a visualization is the unusual regularity of the embedding: A continuously differentiable map that cannot be enhanced to be twice continuously differentiable. As a consequence, the image surface is smooth enough to have a tangent plane everywhere, but not sufficient to admit extrinsic curvatures.

In the 1970s and 1980s, Gromov, revisiting the results of Nash and others such as Phillips, Smale, or Hirsch, extracted the underlying notion of their works: the *h*-principle (3, 4). This principle states that many partial differential relation problems reduce to purely topological questions. The raison d’être of this counterintuitive phenomenon was later brought to light by Eliashberg and Mishachev (5). To prove that the *h*-principle holds in many situations, Gromov introduced several powerful methods for solving partial differential relations. One of which, convex integration theory (5⇓–7), provides a quasi-constructive way to build sequences of embeddings converging toward isometric embeddings. Nevertheless, because of its broad purpose, this theory remains far too generic to allow for a precise description of the resulting map.

In this article, we convert convex integration theory into an explicit algorithm. We then provide an implementation leading to images of an embedded square flat torus in three-dimensional space. This visualization has led us in turn to discover a unique geometric structure. This structure, described in the corrugation theorem below, reveals a remarkable property: The normal vector exhibits a fractal behavior.

## General Strategy

The general strategy (1) starts with a strictly short embedding, i.e., an embedding of the square torus that strictly shrinks distances. To build an isometric embedding, this initial map is corrugated along the meridians with the purpose of increasing their length (Fig. 1). This corrugation is performed while keeping a strictly short map, achieving a smaller isometric default in the vertical direction. The isometric perturbation in the horizontal direction is also kept under control by choosing the number of oscillations sufficiently large. Corrugations are then applied repeatedly in various directions to produce a sequence of maps, reducing step-by-step the isometric default. Most importantly, the sequence of oscillation numbers can be chosen so that the limit map achieves a continuously differentiable isometry onto its image.

### One-Dimensional Convex Integration.

The above general strategy leaves a considerable latitude in generating the corrugations. This great flexibility is one of the surprises of the Nash–Kuiper result because it produces a plethora of solutions to the isometric embedding problem. It is a remarkable fact that Gromov’s convex integration theory provides both the deep reason of the presence of corrugations and the analytic recipe to produce them. Nevertheless, this theory does not give preference to any particular corrugation and is not constructive in that respect. Here we refine the traditional analytic approach of corrugations, adding a geometric point of view.

A corrugation is primarily a one-dimensional process. It aims to produce, from an initial regular smooth curve (as usual denotes the three-dimensional Euclidean space), a new curve *f* whose speed is equal to a given function with . In the general framework of convex integration (5, 7), one starts with a one-parameter family of loops satisfying the isometric condition , for all , and the barycentric condition [1]This last condition expresses the derivative as the barycenter of the loop *h*(*t*,·). One then chooses the number *N* of oscillations of the corrugated map *f* and set [2]Here, *Nu* must be considered modulo . It appears that not only as desired, but *f* can also be made arbitrarily close to the initial curve *f*_{0}, see Fig. 2.

We have where denotes the *C*^{0} norm of a function .

Let *t*∈[0,1]. We put *n*≔[*Nt*] (the integer part of *Nt*) and set for 0 ≤ *j* ≤ *n* - 1 and . We write with *S*_{j}≔∫_{Ij}*h*(*v*,*Nv*)d*v* and . By the change of variables *u* = *Nv* - *j*, we get for each *j*∈[0,*n* - 1] It ensues that . The lemma then follows from the obvious inequalities and .

Here we set [3]where *e*^{iθ}≔ cos *θ* **t** + sin *θ* **n** with , is a smooth unit vector field normal to the initial curve and the function α is determined by the barycentric condition [**1**]. We claim that our convex integration formula captures the natural geometric notion of a corrugation. Indeed, if the initial curve *f*_{0} is planar then the signed curvature measure of the resulting curve is connected to the signed curvature measure *μ*_{0}≔*k*_{0}d*s* of the initial curve by the following simple formula Our corrugation thus modifies the curvature in the simplest way by sine and cosine terms with frequency *N*.

### Two-Dimensional Convex Integration.

The classical extension of convex integration to the two-dimensional case consists in applying the one-dimensional process to a one-parameter family of curves that foliates a two-dimensional domain. Given a strictly short smooth embedding of the square flat torus, a nowhere vanishing vector field , and a function , the aim is to produce a smooth map whose derivative in the direction *W* has the target norm *r*. The natural generalization of our one-dimensional process leads to the following formula: [4]where *φ*(*t*,*s*) denotes the point reached at time *s* by the flow of *W* issuing from *tV*. The vector *V* is chosen so that the line of initial conditions is a simple closed curve transverse to the flow. We also use the notation *e*^{iθ} = cos *θ* **t** + sin *θ* **n**, where **t** is the normalized derivative of *f*_{0} along *W* and **n** is a unit normal to the embedding *f*_{0}. Similarly as above, *θ*(*q*,*u*)≔*α*(*q*) cos 2*πNu*, α is determined by the barycentric condition [**1**] and . The resulting map *f* is formally defined over a cylinder. In general *f* does not descend to the flat square torus . This defect is rectified by adding a term that smoothly spreads out the gap preventing the map to be doubly periodic.

### Basis of the Embeddings Sequence.

The iterated process leading to an isometric embedding requires to start with a strictly short embedding of the square flat torus The metric distortion induced by *f*_{init} is measured by a field of bilinear forms obtained as the pointwise difference: As usual, denotes the pullback of the Euclidean inner product by *f*_{init}. Notice that *f*_{init} is strictly short if and only if the isometric default Δ is a metric, i.e., a map from the square flat torus into the positive cone of inner products of the plane. The convexity of this cone implies the existence of linear forms of the plane *ℓ*_{1},…,*ℓ*_{S}, *S*≥3, such that for nonnegative functions *ρ*_{j}. By a convenient choice of the initial map *f*_{init} and of the *ℓ*_{j}s, the number *S* can be set to three. In practice we set the linear forms to the normalized duals of the following constant vector fields: where (*e*_{1},*e*_{2}) is the canonical basis of . For later use, we set Note that the parallelogram spanned by *U*(*j*) and *V*(*j*) is a fundamental domain for the action over . As an initial map we choose the standard parametrization of a geometric torus. It is easy to check that the range of the isometric default Δ lies inside the positive cone spanned by the *ℓ*_{j}⊗*ℓ*_{j}s, *j*∈{1,2,3}, whenever the sum of the minor and major radii of this geometric torus is strictly less than one (Fig. S1).

We define a sequence of metrics converging toward the Euclidean inner product, [5]with *δ*_{k} = 1 - *e*^{-γk} for some fixed *γ* > 0. We then construct a sequence of maps such that every *f*_{k} is quasi-isometric for *g*_{k}, i.e., . In other words, each *f*_{k}, seen as a map from the square flat torus to Euclidean three space, has an isometric default approximately equal to *e*^{-γk}Δ.

The map *f*_{k} is obtained from *f*_{k-1} by a succession of corrugations. Precisely, if *S*_{k} linear forms *ℓ*_{k,1},…*ℓ*_{k,Sk} are needed for the convex decomposition of the difference then *S*_{k} convex integrations will also be needed to (approximately) cancel every coefficient *ρ*_{k,j}, *j*∈{1,…,*S*_{k}}. As a key point of our implementation, we manage to set each number *S*_{k} to three and to keep unchanged our initial set of linear forms {*ℓ*_{k,1},*ℓ*_{k,2},*ℓ*_{k,3}} = {*ℓ*_{1},*ℓ*_{2},*ℓ*_{3}}. We therefore generate a sequence with an infinite succession of three terms blocks. Each map is the result of a two-dimensional convex integration process applied to the preceding term of the sequence. We eventually obtain the desired sequence of maps, setting *f*_{k}≔*f*_{k,3} for .

### Reduction of the Isometric Default.

The aim of each convex integration process is to reduce one of the coefficients *ρ*_{k,1}, *ρ*_{k,2}, or *ρ*_{k,3} without increasing the two others. This goal is achieved with a careful choice for the field of directions along which we apply the corrugations. Suppose we are given a map *f*_{k,0}≔*f*_{k-1,3} whose isometric default with respect to *g*_{k} lies inside the cone : [6](the *ρ*_{k,j}s being positive functions). We would like to build a map *f*_{k,1} with the requirement that its isometric default is roughly equal to the sum of the last two terms *ρ*_{k,2}*ℓ*_{2}⊗*ℓ*_{2} + *ρ*_{k,3}*ℓ*_{3}⊗*ℓ*_{3}. To this end we introduce the intermediary metric [7]and observe that the above requirement amounts to ask that *f*_{k,1} is quasi-isometric for *μ*_{k,1}. Although natural, it turns out that performing a two-dimensional convex integration along the constant vector field *U*(1) does not produce a quasi-isometric map for *μ*_{k,1}. Instead, we consider the following nonconstant vector field: where the scalar *ζ*_{k,1} is such that the field *W*_{k,1} is orthogonal to *V*(1) for the metric *μ*_{k,1}. With this choice the integral curves *φ*(*t*,.) of *W*_{k,1} issuing from the line of define a diffeomorphism . We now build a new map *F*_{k,1} by applying to *f*_{k,0} a two-dimensional convex integration (see Eq. **4**) along the integral curves *φ*(*t*,.), i.e., [8]The isometric condition in the direction *W*_{k,1} for the metric *μ*_{k,1} is . By differentiating [**8**] with respect to *s* we get , hence we must choose . Furthermore, the map *f*_{k,0} is strictly short for *μ*_{k,1} because We finally set *θ*(*q*,*u*)≔*α*(*q*) cos 2*πN*_{k,1}*u*, where *N*_{k,1} is the frequency of our corrugations.

Note that the map *F*_{k,1} is properly defined over a cylinder, but does not descend to the torus in general. We eventually glue the two cylinder boundaries with the following formula, leading to a map *f*_{k,1} defined over , [9]where *w*: [0,1] → [0,1] is a smooth *S*-shaped function satisfying *w*(0) = 0, *w*(1) = 1, and *w*^{(k)}(0) = *w*^{(k)}(1) = 0 for all .

To cancel the last two terms in [**6**], we apply two more corrugations in a similar way. For every *j*, the intermediary metric *μ*_{k,j} involves *f*_{k,j-1} and the *j*th coefficient of the isometric default . Notice that the three resulting maps *f*_{k,1}, *f*_{k,2}, and *f*_{k,3} are completely determined by their numbers of corrugations *N*_{k,1}, *N*_{k,2}, and *N*_{k,3}.

For *j*∈{1,2,3}, there exists a constant *C*_{k,j} independent of *N*_{k,j} (but depending on *f*_{k,j-1} and its derivatives) such that

The first point ensures that *f*_{k,j} is *C*^{0} close to *f*_{k,j-1}, whereas the second point keeps the increase of the differentials under control and the last point guarantees that *f*_{k,j} is quasi-isometric for *μ*_{k,j}.

We deduce from [**9**] that We then apply Lemma 1 to the right-hand side with *f*≔*F*_{k,j}∘*φ*(*t*,·) and *f*_{0}≔*f*_{k,j-1}∘*φ*(*t*,·) to obtain (*i*). For (*ii*), it is sufficient to bound ‖*df*_{k,j}(*X*) - *df*_{k,j-1}(*X*)‖_{C0} for *X* = *V*(*j*) and *X* = *W*_{k,j}. Because for some nonvanishing function *c*_{φ}, the norm ‖*df*_{k,j}(*V*(*j*)) - *df*_{k,j-1}(*V*(*j*))‖_{C0} is bounded by , up to a multiplicative constant. It is readily seen that results from a convex integration process applied to and Lemma 1 shows that . For *X* = *W*_{k,j}, we differentiate [**9**] with respect to *s* and obtain with Ψ≔*F*_{k,j} - *f*_{k,j-1}. We bound the second term of the right-hand side as for (*i*). Let be the Bessel function of 0 order. On the one hand, for every nonnegative α lower than the first positive root of *J*_{0} we have: 1 + *J*_{0}(*α*) - 2*J*_{0}(*α*) cos(*α*) ≤ 7[1 - *J*_{0}(*α*)]. On the other hand, , where . Because *df*_{k,j-1}(*W*_{k,j}) = *rJ*_{0}(*α*)**t** we obtain From [**7**], we deduce , hence (*ii*). Once the differential of Ψ is under control, (*iii*) reduces to a meticulous computation of the coefficients of in the basis (*W*_{k,j},*V*(*j*)).

## Corrugation Numbers

We make a repeated use of Theorem 1 to show that the map *f*_{k,3} is quasi-isometric for *g*_{k} and strictly short for *g*_{k+1}. Thereby, the whole process can be iterated. Moreover, the *C*^{0} control of the maps and the differentials, as provided by Theorem 1, allows in turn to control the *C*^{0} and *C*^{1} convergences of the sequence . With a suitable choice of the *N*_{k,j}s, this sequence can be made *C*^{1} converging, thus producing a *C*^{1} isometric map *f*_{∞} in the limit. By the *C*^{0} control of the sequence we also obtain a *C*^{0} density property: Given *ϵ* > 0, the *N*_{k,j}s can be chosen so that

### Our Isometric Constraint.

Compared to Nash’s and Kuiper’s proofs, we have an extra constraint. For each integer *k* the isometric default must lie inside the convex hull spanned by the *ℓ*_{j}⊗*ℓ*_{j}, *j*∈{1,2,3}. The reason for this constraint is to avoid the numerous local gluings required by Nash’s and Kuiper’s proofs. Using a single chart substantially simplifies the implementation. More importantly, keeping the same linear forms *ℓ*_{j} all through the process enlightens the recursive structure of the solution that was hidden in the previous methods. To deal with this constraint, we introduce some more notations. Let where the *ρ*_{j}s are, as above, the coefficients of the decomposition of Δ over the *ℓ*_{j}⊗*ℓ*_{j}s. We also denote by *err*_{k,j} the norm of the isometric default of *f*_{k,j}. By Theorem 1, this number can be made as small as we want provided that the number of corrugations *N*_{k,j} is chosen large enough.

If then lies inside .

We want to show that *ρ*_{j}(*D*) > 0 for *j*∈{1,2,3}. Let . Because *D* = (*δ*_{k+1} - *δ*_{k})Δ + *B*, we have by linearity of the decomposition coefficient *ρ*_{j}: [10]In particular, the condition ‖*ρ*_{j}(*B*)‖ < (*δ*_{k+1} - *δ*_{k})*ρ*_{min}(Δ) implies *ρ*_{j}(*D*) > 0. Now, it follows by some easy linear algebra that and a computation shows that ‖*B*‖_{C0} ≤ 3(*err*_{1} + *err*_{2} + *err*_{3}).

We are now in position to choose the corrugation numbers, and doing so, to settle a complete description of the sequence .

### Choice of the Corrugation Numbers.

Let *ϵ* > 0. At each step, we choose the corrugation number *N*_{k,j} large enough so that the following three conditions hold (*j*∈{1,2,3}):

Here we have put, similarly as above, . The first condition ensures the *C*^{0} closeness of *f*_{∞} to *f*_{init}. Thanks to Lemma 2, the third condition implies that the isometric default lies inside the cone . It can be shown to also imply that the intermediary bilinear forms *μ*_{k,2} and *μ*_{k,3} are metrics, an essential property to apply the convex integration process to *f*_{k,1} and *f*_{k,2}. Finally, we can prove the *C*^{1} convergence of the resulting sequence with the help of the second condition.

Note that at each step, the map *f*_{k,j} is ensured to be a *C*^{1} embedding if *N*_{k,j} is chosen large enough. This property follows from the two conditions (*i*) and (*ii*), because a *C*^{1} immersion, which is *C*^{1} close to a *C*^{1} embedding, must be an embedding.

*C*^{1} Fractal Structure

The recursive definition of the sequence paves the way for a geometric understanding of its limit. Because the resulting embedding is *C*^{1} and not *C*^{2}, this geometry consists merely of the behavior of its tangent planes or, equivalently, of the properties of its Gauss map. We denote by **v**_{k,j} the normalized derivative of *f*_{k,j} in the direction *V*(*j*) and by **n**_{k,j} the unit normal to *f*_{k,j}. We also set . Obviously, there exists a matrix such that Here, (**a** **b** **c**)^{t} stands for the transpose of the matrix with column vectors **a**, **b**, and **c**. We call a corrugation matrix because it encodes the effect of one corrugation on the map *f*_{k,j-1}. Despite its natural and simple definition, the corrugation matrix has intricate coefficients with integro-differential expressions. The situation is further complicated by some technicalities such as the elaborated direction field of the corrugation or the final stitching of the map used to descend to the torus. Remarkably, all these difficulties vanish when considering the dominant terms of the two parts of a specific splitting of .

[Corrugation Theorem] The matrix can be expressed as the product of two orthogonal matrices where and and where is the norm of the isometric default, *β*_{j} is the angle between *V*(*j*) and *V*(*j* + 1), and, as above, *θ*_{k,j}(*p*,*u*) = *α*_{k,j}(*p*) cos 2*πN*_{k,j}*u*.

The matrix maps to where **t**_{k,j-1} is the normalized derivative of *f*_{k,j-1} in the direction *W*_{k,j} (Fig. S2). This last vector field converges toward *U*(*j*) when the isometric default tends to zero. Hence, reduces to a rotation matrix of the tangent plane that maps *V*(*j* - 1) to *V*(*j*). The matrix accounts for the corrugation along the flow lines. From the proof of Theorem 1 (*ii*) we have . Therefore, modulo , the transversal effect of a corrugation is not visible. In other words, a corrugation reduces at this scale to a purely one-dimensional phenomenon. Hence the simple expression of the dominant part of this matrix. Notice also that Theorem 1 (*i*) implies that the perturbations induced by the stitching are not visible as well.

The Gauss map **n**_{∞} of the limit embedding can be expressed very simply by means of the corrugation matrices: The Corrugation Theorem gives the key to understand this infinite product. It shows that asymptotically the terms of this product resemble each other, only the amplitudes *α*_{k,j}, the frequencies *N*_{k,j}, and the directions are changing. In particular, the Gauss map *n*_{∞} shows an asymptotic self-similarity: The accumulation of corrugations creates a fractal structure.

It should be noted that there is a clear formal similarity between the infinite product defining **n**_{∞} and, in a one-dimensional setting, the well-known Riesz products, where and are two given sequences. It is a fact (8) that an exponential growth of *N*_{k}, known as Hadamard’s lacunary condition, results in a fractional Hausdorff dimension of the Riesz measure *n*(*x*)*dx*. A similar result for the normal **n**_{∞} of the embedding of the flat square torus seems hard to obtain. It is likely that the graph of the Gauss map **n**_{∞} has Hausdorff dimension strictly larger than two. Yet, because the limit map is a continuously differentiable isometry onto its image, the embedded flat torus has Hausdorff dimension two.

## Implementation of the Convex Integration Theory

The above convex integration process provides us with an algorithmic solution to the isometric embedding problem for square flat tori. This algorithm has for initial data three numbers , *ϵ* > 0, *γ* > 0, and a map for which the isometric default Δ is lying inside the cone . From *f*_{init} a finite sequence of maps (*f*_{k,j})_{k∈{1,...,K},j∈{1,2,3}} is iteratively constructed. Each map *f*_{k,j} is built from the map *f*_{k,j-1} by first applying the convex integration formula [**8**] to obtain an intermediary map *F*_{k,j}. The gluing formula [**9**] is further applied to *F*_{k,j} resulting in the composition *f*_{k,j}∘*φ*, where *φ* is the flow of the vector field *W*_{k,j}. We finally get *f*_{k,j} by composing with the inverse map of the flow. The number *γ* rules the amplitude of the isometric default of each *f*_{k,3} via the formula [**5**]. Formulas [**8**] and [**9**] are completely explicit except for the corrugation number *N*_{k,j} which is to be determined so that *f*_{k,j} fulfills the postconditions expressed in the above choice of the corrugation numbers. Because there is no available formula, we obtain by binary search the smallest integer *N*_{k,j} that satisfies these postconditions. The algorithm stops when the map *f*_{K,3} is constructed. Note that Theorem 1 insures that the algorithm terminates. The map *f*_{K,3} satisfies ‖*f*_{K,3} - *f*_{init}‖_{C0} ≤ *ϵ* and its isometric default is less than . Moreover, the limit *f*_{∞} of the *f*_{K,3}s is a *C*^{1} isometric map and . Therefore, the algorithm produces an approximation *f*_{K,3} of a solution of the underdetermined partial differential system for isometric maps, and this approximation is *C*^{0} close to the initial map *f*_{init}.

We have implemented the above algorithm to obtain a visualization of a flat torus shown in Fig. 3. The initial map *f*_{init} is the standard parametrization of the torus of revolution with minor and major radii, respectively, and . Each map of the sequence (*f*_{k,j})_{k∈{1,...,K},j∈{1,2,3}} is encoded by a *n* × *n* grid whose node *i*_{1}, *i*_{2} contains the coordinates of *f*_{k,j}(*i*_{1}/*n*,*i*_{2}/*n*). Flows and integrals are common numerical operations for which we have used Hairer’s solver (9) based on DOPRI5 for nonstiff differential equations. To invert the flow *φ* of *W*_{k,j} we take advantage of the fact that the *U*(*j*) component of *W*_{k,j} is constant. The line of initial conditions is thus carried parallel to itself along the flow. It follows that the points of a uniform sampling of the flow are covered by *n* lines running parallel to the initial conditions. We now observe that the same set of lines also covers the nodes of the *n* × *n* uniform grid over . This last observation leads to a simple linear time algorithm for inverting the flow. It is worth noting that the integrand in Eq. **8** essentially depends upon the first order derivatives of *f*_{k,j-1} and upon the corrugation frequency *N*_{k,j}. The derivatives are accurately estimated with a finite difference formula of order four. Regarding the corrugation frequency, we have observed a rapid growth as from the four first values of *N*_{k,j}. For instance, for *γ* = 0.1, we have obtained the following: In fact, computing these values is already a challenge because a reasonable resolution of 10 grid samples per period of corrugation would require a grid with nodes. We have restricted the computations to a small neighborhood of the origin to obtain the above values. These values are therefore lower bounds with respect to the postconditions for the choice of *N*_{k,j}. However, these postconditions are only sufficient and we were able to reduce these numbers for the four first steps to, respectively, 12, 80, 500, and 9,000 after a number of trials over the entire grid mesh (see Figs. 1 and 3). The pointwise displacement between the last map *f*_{2,1} and the limit isometric map could hardly be detected as the amplitudes of the next corrugations decrease dramatically. Further corrugations would thus not be visible to the naked eye. To illustrate the metric improvement we have compared the lengths of a collection of meridians, parallels, and diagonals on the flat torus with the lengths of their images by *f*_{2,1}. The length of any curve in the collection differs by at most 10.2% with the length of its image. By contrast, the deviation reaches 80% when the standard torus *f*_{init} is taken in place of *f*_{2,1}.

In practice, calculations were performed on a 8-core cpu (3.16 GHz) with 32 GB of RAM and parallelized C++ code. We used a 10,000^{2} grid mesh (*n* = 10,000) for the three first corrugations and refined the grid to 2 milliards nodes for the last corrugation. Because of memory limitations, the last mesh was divided into 33 pieces. We then had to render each piece with a ray tracer software and to combine the resulting images into a single one as in Fig. 3. The computation of the final mesh took approximately 2 h. Two extra days were needed for the final image rendering with the ray tracer software Yafaray (10).

## Conclusion and Perspective

Convex integration is a major theoretical tool for solving underdetermined systems (4). After a substantial simplification, we obtained the implementation of a convex integration process, providing images of a flat torus. This visualization led us in turn to discover a geometric structure that combines the usual differentiability found in Riemannian geometry with the self-similarity of fractal objects. This *C*^{1} fractal structure is captured by an infinite product of corrugation matrices. In some way, these corrugations constitute an efficient and natural answer to the curvature obstruction (11) observed in the introduction, leading to an atypical solution. A similar process occurs in weak solutions of the Euler equation (12) and could be present in other natural phenomena (13, 14).

Despite its high power, convex integration theory remains relatively unknown to nonspecialists (15). We hope that our implementation will help to popularize this technique and will open a door to applications ranging from other isometric immersions, such as hyperbolic compact spaces, to solutions of underdetermined systems of nonlinear partial differential equations. Convex integration theory could emerge in a near future as a major operating tool in a very large spectrum of applied sciences.

## Acknowledgments

We thank Jean-Pierre Kahane for pointing out Riesz products, Thierry Dumont for helpful discussions on mathematical softwares, Damien Rohmer for valuable advice on graphics rendering and the Calcul Intensif/Modélisation/Expérimentation Numérique et Technologique (CIMENT) project for providing access to their computing platform. Research partially supported by Région Rhône-Alpes [Projet Blanc-Créativité-Innovation (CIBLE)] and Centre National de la Recherche Scientifique [Projet Exploratoire Premier Soutien (PEPS)].

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. E-mail: vincent.borrelli{at}math.univ-lyon1.fr.

Author contributions: V.B., S.J., F.L., and B.T. performed research and wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1118478109/-/DCSupplemental.

## References

- ↵
- ↵
- Kuiper N

*C*^{1}-isometric imbeddings. Indag Math 17:545–556. - ↵
- Gromov M

- ↵
- Gromov M

- ↵
- Eliashberg Y,
- Mishachev N

- ↵
- Geiges H

- ↵
- Spring D

- ↵
- ↵
- Hairer E,
- Norsett SP,
- Wanner G

- ↵
- ↵
- Han Q,
- Hong JX

- ↵
- Villani C

- ↵
- ↵
- Liang H,
- Mahadevan L

- ↵
- Berger M

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Physical Sciences
- Mathematics