given a curve @xmath4 , we consider the geometric flow @xmath5    where @xmath6 is the curvature and @xmath7 the binormal component of the frenet - serret formulas @xmath8    the flow can be expressed as @xmath9    where @xmath10 is the usual cross - product , @xmath11 is the time , and @xmath12 is the arc - length parameter .
it appeared for the first time in 1906 in the ph .
d. thesis of da rios @xcite and was rederived in 1965 by arms and hama @xcite as an approximation of the dynamics of a vortex filament under the euler equations .
this model is usually known as the vortex filament equation ( vfe ) ; we refer the reader to @xcite and @xcite for an analysis and discussion of its limitations .
equations and are also known as the binormal equation and the localized induction approximation ( lia ) , respectively .
some of their explicit solutions are the line , the circle , and the helix .
since the tangent vector @xmath13 remains with constant length , we can assume that @xmath14 , for all @xmath11 .
differentiating , we get the so - called schrdinger map equation onto the sphere : @xmath15    which is a particular case of the landau - lifshitz equation for ferromagnetism @xcite . equation can be rewritten in a more geometric way as @xmath16 where @xmath17 is the covariant derivative , and @xmath18 is the complex structure of the sphere .
written in this way , can be generalized to more general definition domains and images , as in @xcite , where the hyperbolic plane @xmath19 was chosen as the target space .
a relevant step forward in the understanding of was given by hasimoto in @xcite , where a transformation ( see section [ s : hasimoto ] ) that relates and with the nonlinear schrdinger ( nls ) equation was introduced : @xmath20    equation is time reversible , i.e. , if @xmath2 is a solution , so is @xmath21 .
bearing in mind this fact , an important property of ( and , hence , of ) is that it has a one - parameter family of regular self - similar solutions that develop a corner - shaped singularity at finite time .
this was proved in @xcite for the euclidean case , and in @xcite for the hyperbolic case .
the self - similar solutions are written in the form @xmath22    as was proved by buttke in @xcite , @xmath23 is a smooth curve characterized by the geometric properties @xmath24    with @xmath25 being a constant .
therefore , the shape of @xmath23 is the one of a hairpin : it looks like a segment of a circle of radius @xmath26 close to the origin , that becomes two helicoidal curves with a bigger pitch for bigger @xmath27 , and that winds around two different lines @xcite .
buttke also studied numerically these self - similar solutions .
later on , in @xcite , a careful numerical study of those solutions for both the euclidean and the hyperbolic cases was done . on the one hand
, the authors reproduced numerically in @xcite the formation of the corner - shaped singularity , and , on the other hand , they started with a corner - shaped initial datum , recovering numerically the self - similar solutions ; in all cases , the correct choice of boundary conditions was shown to be vital .
furthermore , they gave numerical evidence that fractality phenomena appeared in ( see figure 2.6 in p. 1059 of @xcite ) , if , for instance , fixed boundary conditions were imposed on the tangent vector @xmath28 .
remark that the relationship between fractals and the schrdinger map equation is not new ; indeed , in @xcite , an aortic valve model was proposed , in order to study the apparent fractal character of the valve s fiber architecture . instead of @xmath2 ,
the authors wrote @xmath29 , where @xmath30 , which corresponds to our time , is such that the curves @xmath31 are the fibers . after imposing fixed boundary conditions at @xmath32
, they showed that the curves @xmath33 have a fractal character . besides
, the fractal dimension of those curves was calculated numerically in @xcite .
it is well known that this fractal behavior already occurs at the linear level , as can be easily seen when solving numerically the free schrdinger equation @xmath34    with an initial datum that is piecewise constant as the sign function in , say , @xmath35 $ ] , and then extended by periodicity . the trigonometric series associated to these problems
have been extensively studied and , as we will see later , the behavior at rational times modulo @xmath36 is completely different to that at irrational ones . in @xcite , berry and
goldberg made a far - reaching connection of this behavior with the so called talbot effect in optics . in @xcite and in @xcite ,
the authors used to model the talbot effect and showed that at rational times the solution can be obtained as a finite overlapping of translates of the initial datum , while at irrational times the images have a fractal profile .
also in @xcite , a conjecture was made about the dimensions of those fractals .
as far as we know , the only available results on the nonlinear setting are the recently ones obtained by erdoan and tzirakis @xcite ( see also @xcite ) in the cubic nonlinear schrdinger equation ( i.e. equation , with @xmath37 ) , for the initial data given by piecewise continuous functions .
this means that the data are `` almost '' in the sobolev class @xmath38 , that is subcritical with respect to the one determined by the scaling of the equation ( i.e. @xmath39 ) .
this typically implies that the nonlinear potential can be seen as an external force .
then , its contribution is computed using duhamel s integral , which has a regularizing effect . as a consequence ,
the fractal behavior is due to the linear part of the solution . at this point ,
the available results for the linear theory are used ( see @xcite and @xcite ) .
more recently , erdoan and tzirakis , in collaboration with chousionis , have also studied in @xcite the schrdinger map equation , but still with subcritical initial conditions .    in this paper , we are interested in filaments with corners .
this means to consider delta functions as the initial data for or piecewise continuous functions as the initial data for .
therefore , we work with data in critical spaces and the problem turns out to be much more involved .
the case of an infinite curve with just one corner and that is otherwise smooth has been studied in the sequence of papers @xcite by banica and vega .
one of the conclusions is that is in fact ill - posed for this type of initial conditions , so that is also ill - posed if it is understood in a classical sense .
however , after introducing appropriate function spaces and using the appropriate definition of solution , it is proved that and are well - posed for this type of initial condition .
more concretely , the velocity at the point where the corner is located is determined by the self - similar solution that at time zero has a corner with the same angle .
in fact , in @xcite it was proved that there exists just one self - similar solution with this property .
even if the solutions of for an initial datum with a corner are well understood , nothing had been done for more general initial data with several corners , in particular , polygons .
nevertheless , in two recent papers @xcite by jerrard and smets , they prove a global existence theorem that allows to consider such types of initial data .
moreover , they simulate numerically the evolution of the unit square at different times , suggesting that the solution could become again polygonal ; indeed , at a certain time , the square seems to reappear , but with the axes turned by a @xmath40-angle with respect to the initial ones . in the following pages , we will show not only that these observations are correct , but that , given a regular planar polygon with @xmath41 sides as the initial datum , another polygon is obtained at times of the form @xmath42 , @xmath43 .
in fact , we will give a complete description of this polygon ; in particular , and except for @xmath44 , we will see that it is a skew polygon with @xmath45 sides ( if @xmath46 odd ) , or @xmath47 sides ( if @xmath46 even ) .
one could wander what , if any , is the connection of these solutions with real fluids . of course , one can not expect anything at the quantitative level , because da rios and arms and hama s approximations are based first in a truncation of the biot - savart integral and then in a renormalization of the time variable .
however , using vfe as an approximation of the evolution of real filaments can not be ruled out at the qualitative level , even with data as singular as the ones we consider .
in fact , in @xcite , the similitude between the self - similar solutions and a flow traversing a triangular wing is pointed out ( see figure 1.1 in p. 1051 of that paper ) .    in the concrete case of the solutions considered in this article
, we could ask if the macroscopic effect of the turning axes mentioned above can be seen in the dynamics of real fluids .
it appears that this is something quite well documented for noncircular jets .
it was first observed when the nozzle has an elliptic shape and , later on , seen for shapes given by a rectangle , a square , and an equilateral triangle among others .
the interest of considering nozzles with corners started in the mid 80 s and has received considerable attention since then .
we refer the reader to the survey @xcite , where the papers related to the above observations can be found , together with a detailed analysis on how sensitive the axes - turning phenomenon is to the initial geometric conditions ; and several numerical experiments are also reviewed .
it has to be understood that , at the numerical level , the corners are always effectively rounded and , in real experiments , the corners are immediately smoothed out once the jet leaves the nozzle , due to viscosity effects .
this is consistent with the behavior of the self - similar solutions mentioned above , because the initial corner evolves into a smooth curve with the shape of a hairpin .
furthermore , the appearance of smooth polygonal structures with more angles has been appreciated for nozzles with the shape of an equilateral triangle ( see ( * ? ? ?
* figure 6 ) ) .
also in ( * ? ? ?
* figure 10 ) , while in @xcite it is said `` ... a consistent eightfold distribution pattern is also suggested ... '' for square nozzles .
similar results are obtained in the numerical simulations of @xcite ( see figure 8 in p. 9 , and p. 11 , where the @xmath48-degree `` flip - flop '' in the case of an equilateral triangle is emphasized ) .
we think it would be worth studying in more detail whether or not these more complicated polygonal structures appear in noncircular jet flows .
the structure of this paper is as follows . in section [ s : hasimoto ] , we review the hasimoto transformation for . in section [
s : xt ] , we gather the main theoretical arguments that support our numerical experiments in the form of a theorem .
we start recalling in section [ s : formulation ] some elementary geometric symmetries that regular polygons have and that are preserved by vfe . on the one hand
, those symmetries greatly simplify the numerical implementation and , on the other hand , they play a fundamental role in computing some important quantities that govern the dynamics of regular polygons .
then , in section [ s : galilean ] , we pay especial attention to a group of symmetries that , as far as we know , can only be visualized through the use of the hasimoto transformation and the nls equation .
they are the so - called galilean transformations , and leave invariant the set of solutions of vfe .
we proceed as in @xcite , where the solution of the nls equation that has the delta function as an initial condition is characterized by the fact that it leaves invariant all the galilean transformations ; this property determines the solution , except for one function of time that is easily computed by solving a simple ode and choosing @xmath49 in appropriately . in our case , we take a planar regular polygon with @xmath41 sides as the initial datum . then , as in @xcite , the solution of is determined except , again , by a function that depends on time ; but , in this case , it is very delicate how to specify this function and does not seem to be easy . our approach in this paper is to integrate the frenet - serret frame with @xmath50    where @xmath51 , which is in this case a distribution , is obtained from the hasimoto transformation .    in section [ s : rational ] , we evaluate @xmath52 at times @xmath53 that are rational multiples of @xmath54 , i.e. , @xmath42 , @xmath55 , because @xmath56 takes then the form , so the corresponding curve is a skew polygon with @xmath45 sides , for @xmath46 odd , and with @xmath47 sides , for @xmath46 even .
remark that the structure of the polygon is determined by the generalized quadratic gau sums @xmath57 appearing in ( see [ s : gausssums ] ) ; in particular , when @xmath46 is even , half of the sums are zero , which ultimately explains why the polygon has @xmath47 sides in that case , instead of @xmath45 .    the value @xmath58 is obtained in section [ s : recovering ] , by imposing the condition that the polygon has to be closed . as a consequence , the tangent vector @xmath28 is determined , except for a rotation that can depend on time . this rotation is computed in section [ s : rigid ] by using the extra symmetries that are available for regular polygons and that we have already mentioned .
the next step is to integrate @xmath28 , in order to obtain the corresponding skew polygon @xmath59 at rational times .
using again the symmetries , we are able to calculate all the necessary quantities , except for some vertical translation of @xmath59 that can be easily computed numerically ( although we do not see a simple way of doing it by using just theoretical arguments ) .    in section [ s : numerical ] , we propose a numerical method for and . in section [
s : symmetries ] , we explain carefully how to take advantage of the symmetries of @xmath59 and @xmath28 for the types of solutions considered , which leads to a dramatically improvement in the computational cost , as explained in section [ s : cost ] .    in section [ s : experiments ] , we perform numerical experiments , simulating the evolution of and for different numbers of initial sides @xmath41 .
concerning the vertical translation of @xmath59 , we give strong evidence that the center of mass of @xmath59 moves upward with a constant speed @xmath60 that depends on @xmath41 .
we compare the numerical results with the algebraic values from section [ s : xt ] , obtaining a complete agreement between both ( totally different ) approaches .
in section [ s : riemann ] , we study the evolution of @xmath2 , for @xmath61 . at this point , it is important to recall the invariances of the so - called jacobi theta function , closely related to : @xmath62    this function is precisely the solution of the free schrdinger equation with a one - periodic delta as the initial condition .
therefore , the galilean transformations still hold in this case ; moreover , they also imply the extra symmetry @xmath63    which , together with the periodicity property @xmath64 , generates the so - called unimodular subgroup .
in fact , in @xcite , all these symmetries are used to prove the existence of a continuum range of exponents for the hlder regularity of @xmath65,\ ] ]    whose real part is precisely riemann s non - differentiable function . more concretely , in @xcite , it was proved that riemann s non - differentiable function is a multifractal ( see also @xcite , where other relevant trigonometric sums are considered ) . in other words , it was proved that the set of times @xmath11 that have the same hlder exponent is a fractal with a dimension depending on the hlder exponent , in such a way that the conjecture stated by frisch and parisi in @xcite is fulfilled ( see also @xcite , for more details at this respect and the connection of this question with fully developed turbulence and intermittency ) .
we could ask if something similar to the properties of riemann s non - differentiable function also holds in our case .
it is very easy to find the analogous formulas to in out setting .
in fact , for a given number of initial sides @xmath41 , if we write @xmath66 , bearing in mind the symmetries of the problem , we conclude that @xmath3 is a planar curve , so we identify the plane where it lives with @xmath67 and define @xmath68    this curve , or , more precisely , @xmath69 , can be seen as a nonlinear version of . is @xmath70 a multifractal ? if the answer is positive , what is its spectrum of singularities in the sense of @xcite ?
we deem these two questions rather challenging from both a numerical and an analytical point of view .
the results of section [ s : riemann ] strongly suggest that the answer to the first question is positive .
in fact , after computing @xmath69 for different @xmath41 , and measuring the difference in the @xmath71-norm between it and an appropriately scaled and rotated version of @xmath72 , the convergence rate is rather strong .
we also give some numerical evidence that the hlder exponent of @xmath70 is @xmath73 for rational times @xmath53 ; nevertheless , how the constants depend of the denominator @xmath46 , which is a fundamental ingredient in the arguments in @xcite , is unclear .
this question deserves a much more detailed analysis that we plan to make in a forthcoming paper .    in section [ s : fractalst ]
, we study some interesting fractality phenomena associated to the tangent vector @xmath28 . in section [
s : conclusions ] , we draw the main conclusions , as well as some open questions that we postpone for the future . finally , to conclude this paper , we offer in [ s : gausssums ] a detailed study on the generalized quadratic gau sums .
a central point of this paper is the natural connection between - and the nls equation .
indeed , applying the hasimoto transformation @xcite to : @xmath74    the function @xmath75 satisfies the nls equation @xmath76    where @xmath49 is a certain real constant that depends on time . nevertheless , for our purposes
, it is not convenient to work with the torsion @xmath77 . instead
, we consider another version of the frenet - serret trihedron ; it is easy to check that all its possible generalizations have the form @xmath78    for some vectors @xmath79 and @xmath80 that form an orthonormal base with the tangent vector @xmath28 .
moreover , we can make , without loss of generality , one of the coefficients @xmath81 , @xmath82 , or @xmath83 equal zero .
if we choose @xmath84 , denoting @xmath85 and @xmath86 , we recover the habitual trihedron . on the other hand , in this paper ,
we choose @xmath87 , in order to avoid working with the torsion . in that case
, the hasimoto transformation takes the form @xmath88    where @xmath89    it is quite straightforward to check that this new definition of @xmath75 also satisfies .
the key point is to combine the basis vectors @xmath79 and @xmath80 into @xmath90 .
since the proof closely matches that in @xcite , we omit here the details .
the main idea is to work with , and , at a given @xmath11 , recover the curve @xmath2 and the tangent vector @xmath91 from @xmath52 by integrating , up to a rigid movement that we determine by the symmetries of the problem .
observe that , if we define @xmath92    and integrate for @xmath93 , the new solution of is @xmath94 , where @xmath95 , and @xmath96 , i.e. , @xmath28 does not change .
the aim of this paper is to understand the evolution of for polygonal initial data .
in fact , we limit ourselves to studying the simplest case , i.e. , that of a regular planar polygon with @xmath41 sides . more precisely , we prove the following theorem :    [ t : theorem ]    assume that there is a unique solution of the initial value problem @xmath0    where @xmath98 is a regular planar polygon with @xmath41 sides . then , at a time @xmath53 that is a rational multiple of @xmath54 , i.e. , @xmath99 , with @xmath100 , @xmath101 , @xmath43 , the solution is a skew polygon with @xmath45 sides ( if @xmath46 odd ) or @xmath47 sides ( if @xmath46 even ) . all the new sides have the same length ; and the angle @xmath102 between two adjacent sides is constant .
furthermore , the polygon at a time @xmath53 is the solution of the frenet - serret system @xmath103    where @xmath104 , and @xmath56 is the @xmath105-periodic function defined over the first period @xmath106 as @xmath107    with @xmath108    being a generalized quadratic gau sum ( see [ s : gausssums ] ) .
the rest of this section is devoted to the proof of this theorem .
moreover , we also conjecture that the angle @xmath102 is given by @xmath109    * remark . * from our point of view
, uniqueness is a very challenging problem for this type of equations , because of the critical regularity that we are assuming in the initial data .
as we say in the introduction , uniqueness has been established in a series of papers by banica and vega , when the initial datum has just one corner and is otherwise regular .
banica and vega prove that uniqueness holds for the binormal equation and for the schrdinger map equation , but the cubic nls equation is ill - posed as an initial value problem with initial datum being the dirac delta function . therefore , the work of erdoan and tzirakis @xcite can not be used in our case , and the only general result for the binormal flow for data as singular as ours is so far the one by jerrard and smets @xcite , where , on the other hand , they just prove existence .      since and
are invariant with respect to rotations , given a regular planar polygon with @xmath41 sides , we can assume without loss of generality that @xmath1 , and hence @xmath110 , live in the plane @xmath111 .
identifying @xmath111 with @xmath67 , and assuming without loss of generality a total length of @xmath36 , @xmath1 is the polygon parameterized by arc - length whose @xmath41 vertices , located at @xmath112 , @xmath113 , are @xmath114    and a point @xmath1 , for @xmath115 , lays in the segment that joints @xmath116 and @xmath117 .
observe that it is straightforward to extend @xmath1 periodically to the whole real with period @xmath36 .
then , since @xmath118 , the corresponding tangent vector is the periodic function with period @xmath36 such that @xmath119    the curve @xmath1 is continuous , while its tangent vector @xmath110 is only piecewise continuous .
in fact , @xmath1 can be regarded as a curve whose @xmath105-periodic curvature is a sum of dirac deltas : @xmath120    and the constant @xmath105 has been chosen in order for the integral of the curvature over an interval of length @xmath36 to be equal to @xmath36 .
at this point , it is necessary to understand the role played by symmetries in and for these initial data .
both equations are invariant with respect to rotations , i.e , given a rotation matrix @xmath121 , if @xmath122 and @xmath123 are , respectively , solutions of them , so are @xmath124 and @xmath125 ( in this paper , all the vectors are represented in column form ) .
therefore , if @xmath126 and @xmath127 , then , if the solution is unique , @xmath128 and @xmath129 , for all @xmath11 . in particular , since @xmath1 and @xmath110 , as defined in and , are invariant with respect to rotations of angle @xmath130 around the @xmath131-axis , for all @xmath132 , we conclude that also @xmath2 and @xmath91 are invariant with respect to that kind of rotations , for all @xmath11 .
more precisely , @xmath133    consequently , at a given @xmath11 , @xmath134 lay in the same orthogonal plane to the @xmath131-axis , for all @xmath135 .
this fact has important implications for implementing efficient numerical schemes , as we will see in section [ s : numerical ] .
besides invariance by rotations , is also mirror invariant : if @xmath59 is a solution of , so is @xmath136 , etc .
therefore , if @xmath137 , then that symmetry will be again preserved during the evolution , i.e , @xmath138 , for all @xmath11 , or , in other words , @xmath2 and @xmath139 will be symmetric with respect to the plane containing the @xmath131-axis and the @xmath140-axis . in our case , a regular polygon with @xmath41 sides like has @xmath141 of those mirror symmetries .
an important corollary that will be used later is that @xmath142 is a positive multiple of the vector @xmath143 .
one of the fundamental symmetries of the nls equation is the so - called galilean transformations , i.e. , if @xmath75 is a solution of , so is @xmath144    therefore , if we choose an initial datum such that @xmath145 , for all @xmath146 , i.e. , such that @xmath147 , for all @xmath146 , then , if the solution is unique , @xmath148 , for all @xmath149 . in this paper , this symmetry is enough for our purposes ; for the full symmetry group of the nls equation , see ( * ? ? ?
* section 17.2 ) .    in the case of a planar initial datum @xmath1 for
, the torsion is equal to zero and , from , @xmath150 is the curvature of @xmath1 .
hence , the @xmath150 corresponding to a regular polygon with @xmath41 sides is given by , i.e. , @xmath151    like , @xmath150 is @xmath105-periodic and , since is invariant with respect to space translations , @xmath52 is also @xmath105-periodic , for all @xmath152 .
moreover , it satisfies @xmath153 , for all @xmath132 .
therefore , the galilean transformation holds , i.e. , @xmath154 , for all @xmath132 , and for all @xmath152 . bearing in mind the previous arguments
, we calculate the fourier coefficients of @xmath52 : @xmath155ds      \cr & = \rme^{-\rmi ( mk)^2 t - \rmi m(j - k)(2mkt)}\hat\psi(j - k , t).\end{aligned}\ ] ]    this identity holds for all @xmath156 and @xmath157 .
in particular , evaluating both sides at @xmath158 , @xmath159 , so @xmath75 can be expressed as @xmath160    where @xmath161 is a constant that depends on time and that has to be chosen in such a way that the corresponding curve @xmath59 and tangent vector @xmath28 are periodic of period @xmath36 for a fixed @xmath11 .
remark that , from the arguments following , @xmath161 can be taken up to a constant with modulus one , so , in this paper , we assume without loss of generality that @xmath161 is real for any given @xmath11 .
when @xmath162 , we have trivially @xmath163 ; then , combining and , we recover the following well - known identity : @xmath164    observe that , if we solve the free schrdinger equation @xmath165 for , we get @xmath166    where @xmath167 is the well - known jacobi theta function .
equation is the same as , but with @xmath168 , for all @xmath152 .
let us evaluate at @xmath170 , where @xmath100 , @xmath101 , and we can suppose without loss of generality that @xmath43 .
then , bearing in mind the identity , @xmath171    where @xmath57 was defined in .
moreover , from , it follows that @xmath172    for a certain angle @xmath173 that depends on @xmath174 ( as well as , of course , @xmath175 and @xmath46 ) . therefore , @xmath52 in @xmath176 has evolved from @xmath41 dirac deltas at @xmath162 , to @xmath45 deltas at @xmath53 , for @xmath46 odd ; and to @xmath47 deltas at @xmath53 , for @xmath46 even .
furthermore , introducing into , and restricting ourselves to @xmath177 , i.e. , @xmath106 , we conclude that @xmath178    the coefficients multiplying the dirac deltas are in general not real , except for @xmath162 and @xmath179 .
therefore , @xmath56 does not correspond to a planar polygon , but to a skew polygon with @xmath45 sides , for @xmath46 odd ; and to a skew polygon with @xmath47 sides , for @xmath46 even .
moreover , @xmath180 has to be determined in such a way that the polygon is closed ; this is done in the next section , where we show that its rather involved choice is unique , concluding that @xmath52 is periodic in time , with period @xmath181 .
since the dirac deltas are equally spaced at @xmath182 , the length of the sides is the same .
observe that , if @xmath183 is odd , i.e. , if @xmath184 , then there is no vertex at @xmath61 .
for instance , if @xmath185 , @xmath44 , i.e. , @xmath186 , then @xmath187    in this case , the coefficients multiplying the dirac deltas are real , so @xmath188 is simply the curvature , and @xmath189 has been chosen in order for the integral of @xmath188 over @xmath190 to be @xmath36 .
bearing in mind the symmetries of the problem , we conclude that the corresponding polygon is a regular planar polygon with @xmath41 sides , which lives at a certain plane @xmath191 , but with the axes turned by @xmath192 degrees with respect to the initial datum , as r. jerrard and d. smets predicted in @xcite .      given a function @xmath193 ,
recovering the basis vectors @xmath28 , @xmath79 and @xmath80 from @xmath75 implies integrating . as seen in the previous section , at @xmath182
, @xmath52 is a sum of @xmath45 ( if @xmath46 odd ) or @xmath47 ( if @xmath46 even ) equally spaced dirac deltas in @xmath176 , that corresponds to a skew polygon @xmath194 of @xmath45 or @xmath47 sides . to integrate
, we have to understand the transition from one side of the polygon to the next one .
in order to do that , we reduce ourselves , without loss of generality , to a certain @xmath75 formed by one single dirac delta located at @xmath61 , i.e , @xmath195 , so we have to integrate @xmath196    this system corresponds to some basis vectors @xmath197 , @xmath198 and @xmath199 constant everywhere , except at @xmath61 , where there is a singularity .
let us suppose that , for @xmath200 , @xmath201 , @xmath202 , @xmath203 , so we want to calculate @xmath204 , @xmath205 , @xmath206 , such that , for @xmath207 , @xmath208 , @xmath209 , @xmath210 .
then , @xmath211      \cdot \left (      \begin{array}{c } { \mathbf t}(0 ^ -)^t      \cr \hline { \mathbf{e}_1}(0 ^ -)^t      \cr \hline { \mathbf{e}_2}(0 ^ -)^t      \end{array } \right )      \cr      & = \exp(\mathbf a)\cdot \left (      \begin{array}{c } { \mathbf t}(0 ^ -)^t      \cr \hline { \mathbf{e}_1}(0 ^ -)^t      \cr \hline { \mathbf{e}_2}(0 ^ -)^t      \end{array } \right ) .\end{aligned}\ ] ]    we can compute explicitly the matrix exponential . writing @xmath212 , then @xmath213c_\theta   s_\theta      \cr -s_\rho s_\theta & [ c_\rho - 1]c_\theta s_\theta & c_\rho s_\theta^2 + c_\theta^2 } , \end{aligned}\ ] ]    where @xmath214 , @xmath215 , @xmath216 , @xmath217 .
the matrix @xmath218 is a rotation matrix @xcite ; more precisely , it is the rotation matrix corresponding to a rotation about an axis @xmath219 by an angle @xmath102 .
coming back to the general form of @xmath75 , we have to integrate @xmath45 or @xmath47 times to obtain a closed skew , i.e. , non - planar polygon with @xmath45 or @xmath47 sides .
but , according to , in @xmath106 , @xmath220    hence , we conclude that , at @xmath221 , the angle @xmath102 between two adjacent sides is constant .
@xmath222    furthermore , the structure of the polygon is completely determined by the quantities @xmath173 appearing in the generalized quadratic gau sum , where @xmath223 .
let be @xmath224 the rotation matrix corresponding to @xmath225 .
if @xmath226 , @xmath224 is simply an identity matrix and can be ignored .
otherwise , @xmath224 is obtained after evaluating , at @xmath227 .
therefore , @xmath228    or , equivalently , @xmath229    where @xmath230 , and @xmath231 is periodic modulo @xmath46 , i.e. , @xmath232 .
this formula is valid for both @xmath46 odd and @xmath46 even , bearing in mind that half of the @xmath231 matrices are simply identity matrices if @xmath46 is even .
this explains why the polygon has @xmath233 sides , when @xmath46 is even ; equivalently , it could be regarded as having @xmath46 sides , but half of them being indistinguishable , because when the angle between two adjacent sides is zero , they merge into one single side .
in order for the polygon to be closed , we have to choose @xmath102 in in such a way that the basis vectors @xmath28 , @xmath79 and @xmath80 are periodic , i.e. , @xmath234    which is equivalent to imposing the condition @xmath235    let us define @xmath236    from , @xmath237 is an @xmath41-th root of the identity matrix .
moreover , it is also a rotation matrix that induces a rotation of @xmath105 degrees around a certain rotation axis .
therefore , we have to choose @xmath102 in order that any of the following properties is satisfied : @xmath238    were @xmath239 and @xmath240 denote the trace and the spectrum of @xmath237 , respectively .
we worked with the trace , because it is algebraically easier . in order to understand its structure
, we analyzed a few cases ( @xmath241 ) with the help of a symbolic manipulator .
the obtained results strongly suggest that , for any @xmath46 ( and for any @xmath175 coprime with it ) the only possible choice for @xmath242 is given by .
remark that if @xmath46 is odd , the angle @xmath102 between two adjacent sides is the same for @xmath46 and for @xmath243 .
although finding a universal proof that holds for any @xmath46 goes beyond the scope of this paper , we have checked it for a few more @xmath46 . moreover ,
this conjecture is in agreement with our numerical simulations , as we will see in section [ s : experiments ] .
therefore , we think that there is compelling evidence that is valid for all @xmath101 .    once we have found the correct choice of @xmath102 , we get from , @xmath244    it is straightforward to check that @xmath245 . notice also that @xmath246    finally , from and , we conclude that @xmath56 , @xmath247 , is given by .
this completes the proof of theorem [ t : theorem ] .      knowing the function @xmath56 , the basis vectors @xmath28 ( and , hence , the curve @xmath59 ) , @xmath79 and @xmath80 can be completely determined up to a rigid movement , which
, on the other hand , for an arbitrary polygon , can be pretty complicated to specify .
however , since our initial datum is a regular planar polygon , the symmetries of the solution are very advantageous .    the easiest way to understand and to implement numerically the correct rotation (
although not necessarily the simplest option for symbolic manipulation ) is to work with @xmath59 . to that aim , given a rational time @xmath53 , we compute first the associated rotation matrices @xmath224 .
then , by means of , we obtain up to a rotation the piecewise constant vectors @xmath28 , @xmath79 and @xmath80 , which we denote @xmath248 , @xmath249 and @xmath250 : @xmath251    and the case @xmath177 can be assigned any value , for instance the identity matrix .
again , it is not necessary to consider separately @xmath46 odd and @xmath46 even .
once @xmath248 has been obtained , @xmath252 , which is @xmath59 up to a rigid movement , is computed recursively : @xmath253    where , again , @xmath254 can be given any value , for instance @xmath255 .
to obtain the correct rotation for @xmath59 , we use the symmetries explained in section [ s : formulation ] .
in particular , we use the fact that , at any time @xmath11 , @xmath256 , @xmath113 , have to be coplanar and lay on a plane orthogonal to the @xmath131-axis ; and that @xmath257 is a positive multiple of @xmath143 . an efficient algorithm is as follows :    1 .
compute @xmath258 and @xmath259 .
2 .   compute @xmath260 .
3 .   compute the scalar product between @xmath261 and @xmath262 , i.e. , @xmath263 .
if @xmath264 , @xmath265 is the identity matrix . if not
, @xmath265 is the rotation matrix that induces a rotation of @xmath266 degrees around the axis given by the vector @xmath267 .
5 .   compute @xmath268 .
compute the rotation matrix @xmath269 that induces a rotation of @xmath270 degrees around the @xmath131-axis .
compute the sought rotation matrix , @xmath271 .
update @xmath272 and @xmath273 .
once we have computed the correct rotation of @xmath59 and @xmath28 , we still have to specify a translation for @xmath59 . again , bearing in mind the symmetries of the problem , we translate @xmath59 in such a way that its mass center , given by the mean of @xmath256 , @xmath274 , is in the @xmath131-axis .
thus , the position of @xmath59 is completely specified , up to a vertical translation .
in this section , we simulate numerically and taking respectively and as the initial data .
there have been a couple of papers devoted to reproducing numerically the self - similar solutions of the schrdinger map equation . in @xcite ,
a crank - nicholson scheme was consider , together with a finite - difference discretization of @xmath275 .
later on , in @xcite , both finite - difference discretizations and pseudo - spectral discretizations were considered to simulate and its equivalent on the hyperbolic plane . in particular , given the asymptotical structure of the self - similar solutions of , a truncated domain @xmath276 $ ] , @xmath277 , with a grid based on the chebyshev nodes was found to be very convenient to represent those solutions . nevertheless , due to the clustering of the chebyshev nodes , an explicit scheme to advance in time implied the undesirable restriction @xmath278 , where @xmath279 is the number of nodes . in order to solve
that , it was chosen to work with the stereographic projection of the tangent vector @xmath280 over  @xmath67 , @xmath281    which transforms into a nonlinear schrdinger equation : @xmath282    the advantage of this equation , as opposed to , is that the higher - order term @xmath283 can be treated implicitly , therefore eliminating or at least reducing significantly the restrictions on @xmath284 . indeed , working with was a very adequate choice for the purposes of @xcite , where a second - order semi - implicit backward differentiation formula was chosen .
this scheme is very stable because it imposes a very strong decay in the high frequency modes ; moreover , its low order is easily and effectively compensated by the use of an adaptive method , both in space and in time .
however , in this paper , working with does not seem to be such a good option .
indeed , unlike in @xcite , where an extremely high accuracy was required only for @xmath285 , we are now interested in the behaviour of @xmath28 at all times , for which a second - order scheme seems a very poor option .
after many unsuccessful attempts , we admit that we have been unable to find a good higher - order scheme for .
furthermore , the stereographic projection may introduce an artificial singularity at the points near the south pole of the sphere .
finally , we are also interested in the evolution of the curve @xmath59 , for which working with is of limited help . bearing in mind the previous arguments , we have opted to work directly with or , more precisely , with a combination of and : @xmath286    to which we have applied a fourth - order runge - kutta scheme in time :
@xmath287    the approximations @xmath288 and @xmath289 refer to the numerically obtained values of @xmath59 and @xmath28 at @xmath290 ; the last line guarantees that @xmath291 , for all @xmath292 . additionally , we can also project @xmath293 , @xmath294 , @xmath295 onto the unit sphere , but without significant improvement in the numerical results .
we have combined with a pseudo - spectral discretization directly in space .
more precisely , since we deal with periodic solutions in @xmath176 , we have simulated the evolution of the curve @xmath59 and of the tangent vector @xmath28 at @xmath279 equally spaced nodes @xmath296 , @xmath297 . in order to compute @xmath298 , @xmath275 , etc . , we remember that , given a periodical function @xmath299 evaluated at @xmath300 , its derivatives at @xmath300 can be spectrally approximated as @xmath301    where @xmath302    both and , which are respectively inverse and direct discrete fourier transforms ( dft ) of @xmath279 elements , can be computed efficiently by means of the fast fourier transform ( fft ) algorithm @xcite . moreover , since @xmath59 and @xmath28 are invariant with respect to rotations of @xmath105 degrees around the @xmath131-axis , as shown in , we can reduce the computational cost of and to a dft of @xmath303 elements . in what follows , we explain how to apply this idea to @xmath28 .
this is valid , with no change , for any vector satisfying , i.e. , @xmath59 , @xmath298 , @xmath275 , etc .
now , let us consider the first two components of @xmath28 .
denoting @xmath304 , and bearing in mind that @xmath305 , @xmath306    therefore , the non - zero @xmath307 are obtained by a dft of @xmath303 elements : @xmath308\rme^{-2\pi \rmi jk/(n / m)}.\ ] ]    for the third component of @xmath28 , bearing in mind that @xmath309 , we conclude @xmath310    so the non - zero @xmath311 are obtained again by a dft of @xmath303 elements : @xmath312    since @xmath313 is real , can be further simplified to a dft of @xmath314 elements .      combining with a pseudo - spectral discretization in space
, we obtain a scheme with a time - step restriction on @xmath284 that appears to be of the form @xmath315 , with @xmath316 . however , for a fixed @xmath41 , we only want to simulate the evolution of the curve @xmath59 and of the tangent vector @xmath28 for @xmath317 $ ] . denoting @xmath318 the biggest @xmath284 that makes stable , the smallest number of time - steps @xmath319 needed to reach @xmath320 satisfies @xmath321    from , together with and , it follows that , if the number of spatial nodes @xmath279 is of the form @xmath322 , with @xmath323 , then , the computational cost for simulating is exactly the same for any number of initial sides @xmath41 .
this important fact enables us to make consistent comparisons between different @xmath41 , because we are considering an equivalent spatial resolution , i.e , we are using the same number of points @xmath303 per initial side .
we take @xmath303 a power of two to take the greatest possible advantage of the fft . obviously , for @xmath303 constant , @xmath284 decreases as @xmath324 ; therefore , the most accurate results will be expected for the largest @xmath41 , as we will observe in our numerical experiments .    in the following section ,
we simulate @xmath59 and @xmath28 for initial data with different @xmath41 .
as we will see , despite its simplicity , gives surprisingly good results .
in section [ s : xt ] , in order to construct algebraically the evolution of a regular polygon at rational times , we have done some very strong assumptions , the most important one being uniqueness . however , as we will see in this section , the numerical experiments are in complete agreement with their theoretical predictions . to illustrate this , we have simulated the evolution of the curve @xmath59 and of the tangent vector @xmath28 by means of , taking regular polygons with different numbers of sides @xmath41 as the initial data , and making @xmath59 and @xmath28 evolve until @xmath320 .
the initial data @xmath1 and @xmath110 are given respectively by and ; in the case of @xmath59 , the non - vertex points are computed by linear interpolation .
we have taken @xmath279 equally spaced nodes @xmath325 , @xmath297 ; nevertheless , from the symmetries of @xmath59 and @xmath28 , as explained in the previous section , we only have to describe the evolution of @xmath59 and @xmath28 at the first @xmath303 nodes @xmath326 , @xmath327 .
we have divided the time - interval @xmath328 $ ] in @xmath319 equally spaced time steps of length @xmath329 .
therefore , during the simulation of , we obtain @xmath288 and @xmath289 , at @xmath330 , @xmath331 . in our experiments , bearing in mind , we have chosen @xmath332 , for @xmath333 .
remember that , in section [ s : xt ] , we had constructed algebraically the curve @xmath59 up to a vertical movement .
therefore , in order to completely specify @xmath59 at a given time , we would need to give the height of one point or , more conveniently , the height @xmath334 of the mass center , which is precisely the mean of all the values @xmath335 : @xmath336    we have observed that @xmath334 can be approximated with very high accuracy by means of a constant multiplied by @xmath11 ; more precisely , @xmath337    where @xmath338 is the mean speed . in table
[ t : speed ] , we give the maximum discrepancy between @xmath334 and its linear approximation @xmath339 , i.e. , @xmath340 .
we have considered different numbers of initial sides @xmath41 and different numbers of nodes @xmath279 ; the errors are very small and they seem to decrease as @xmath341 .
this gives very strong evidence that @xmath334 is linear of the form @xmath342 .
observe that @xmath343 , which is also offered in table [ t : speed ] for @xmath344 , grows with @xmath41 , tending to @xmath345 , as @xmath41 tends to infinity , i.e. as the initial polygon @xmath1 tends to a circle .
.@xmath340 , for different numbers of initial sides @xmath41 , and of nodes @xmath279 .
@xmath346 ; @xmath338 .
the values corresponding to @xmath343 have been calculated with @xmath344 . [ cols="^,^,^,^,^,^,^",options="header " , ]     on the other hand , comparing the algebraically constructed @xmath347 with the numerically obtained @xmath348 is more problematic .
for example , in figure [ f : tm3p1q3 ] , we have compared @xmath348 ( left ) with @xmath347 ( right ) , for an initial triangle , @xmath349 , at @xmath350 , @xmath344 .
the exact value of @xmath28 at that time is given by @xmath347 , which , by construction , is a piecewise constant function with exactly @xmath351 pieces , that we denote @xmath352 , @xmath353 , and whose explicit values , in this case , have an easy algebraic expression : @xmath3542 - 1 \cr -\sqrt{\sqrt[3]4 - 1 } \cr 1 - \sqrt[3]4 } , \quad { \mathbf t}_2 = \pmatrix { 1 \cr 0 \cr 0 } , \quad { \mathbf
t}_3 = \pmatrix { \sqrt[3]2 - 1 \cr \sqrt{\sqrt[3]4 - 1 } \cr \sqrt[3]4 - 1 } ; \ ] ]    @xmath355 , @xmath356 and @xmath357 are respectively @xmath358 , @xmath359 and @xmath360 rotated @xmath361 degrees around the @xmath131-axis ; @xmath362 , @xmath363 and @xmath364 are respectively @xmath358 , @xmath359 and @xmath360 rotated @xmath365 degrees around the @xmath131-axis . observe that @xmath348 is almost identical to @xmath347 , but clearly exhibits a gibbs - type phenomenon .
however , in @xmath348 , if we take the central values of each interval ( altogether 27 circles , indicated in figure [ f : tm3p1q3 ] with a black circle ) , and compare them componentwise with , we obtain a maximum error equal to @xmath366 .     versus @xmath347 , for @xmath349 , at @xmath350 .
@xmath367 appears in blue , @xmath368 in green , @xmath369 in red . in @xmath348 ,
the gibbs phenomenon is clearly visible .
the black circles denote the points chosen for the comparisons.,title="fig:",scaledwidth=50.0% ]   versus @xmath347 , for @xmath349 , at @xmath350 .
@xmath367 appears in blue , @xmath368 in green , @xmath369 in red . in @xmath348 ,
the gibbs phenomenon is clearly visible .
the black circles denote the points chosen for the comparisons.,title="fig:",scaledwidth=50.0% ]      from the mirror symmetries of the initial data , we conclude that the curve @xmath370 , @xmath371 , lives in a plane that contains the @xmath131-axis .
indeed , @xmath372 are the only values for which @xmath2 describe a planar curve . in
what follows , we will describe @xmath3 , although all said here is immediately applicable to any @xmath372 . since @xmath373 is planar , bearing in mind that @xmath374 , @xmath375 , we rotate @xmath3 clockwise @xmath376 degrees around the @xmath131-axis , until it lays on the plane @xmath377 , which we identify with @xmath67
. then @xmath3 becomes @xmath378    in figure [ f : xs0 ] , we have plotted the numerically obtained @xmath70 for @xmath379 . besides the conspicuous fractal character of the curves , which immediately reminds us of the works @xcite , their most striking feature ( see figure [ f : riemann ] ) is how much @xmath380 , i.e. , z(t ) without the vertical movement , resembles the graph of @xmath381.\ ] ]    , for @xmath379 .
we have plotted the @xmath70 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .,title="fig:",scaledwidth=33.33% ]   , for @xmath379 .
we have plotted the @xmath70 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .,title="fig:",scaledwidth=33.33% ]   , for @xmath379 .
we have plotted the @xmath70 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .,title="fig:",scaledwidth=33.33% ]    the function @xmath72 was conveniently used in @xcite , in order to study riemann s non - differentiable function , which is precisely the real part of @xmath72 , i.e. , @xmath384    this function is non - differentiable everywhere , except at those rational points @xmath385 , with @xmath175 and @xmath46 both odd . on the other hand , when constructing algebraically @xmath3 for a given t , we have observed that the only times at which @xmath3 has no corner are of the form @xmath386 , with @xmath175 odd and @xmath184 .
therefore , in order to compare @xmath72 and @xmath3 , we have to redefine slightly @xmath72 in : @xmath387,\ ] ]    i.e. , we multiply @xmath388 in by @xmath389 , in order to orientate it correctly ; then we change the period from @xmath390 $ ] to @xmath391 $ ] , and the orientation of the parameterization , by doing @xmath392 . on the other hand ,
given a number of initial sides @xmath41 , we define @xmath393.\ ] ]    in figure [ f : riemann ] , we have plotted @xmath72 ( red ) , as defined in , versus @xmath394 , for @xmath395 . since we have obtained numerically @xmath394 at @xmath396 points in
@xmath397 $ ] , we have evaluated @xmath72 precisely at those points .
although the figures are not identical , they are extremely similar . moreover ,
if we scale them accordingly , we observe that , the bigger @xmath41 grows , the more similar to @xmath72 is @xmath394 .
indeed , @xmath72 and @xmath398 , except for the scaling , are visually undistinguishable .
( red ) , as defined in , against @xmath394 ( blue ) , for @xmath395 . to compute @xmath394 , we took @xmath344 ; to compute @xmath72 , we took @xmath399 .
we have plotted the @xmath394 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .
we have computed algebraically and plotted @xmath72 at those times .
observe that , except for the scaling , @xmath72 and @xmath398 are visually undistinguishable.,title="fig:",scaledwidth=33.33% ]    ( red ) , as defined in , against @xmath394 ( blue ) , for @xmath395 . to compute @xmath394
, we took @xmath344 ; to compute @xmath72 , we took @xmath399 .
we have plotted the @xmath394 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .
we have computed algebraically and plotted @xmath72 at those times .
observe that , except for the scaling , @xmath72 and @xmath398 are visually undistinguishable.,title="fig:",scaledwidth=33.33% ]    ( red ) , as defined in , against @xmath394 ( blue ) , for @xmath395 . to compute @xmath394
, we took @xmath344 ; to compute @xmath72 , we took @xmath399 .
we have plotted the @xmath394 corresponding to all the @xmath382 points @xmath383 obtained during the execution of .
we have computed algebraically and plotted @xmath72 at those times .
observe that , except for the scaling , @xmath72 and @xmath398 are visually undistinguishable.,title="fig:",scaledwidth=33.33% ]    in order to give some numerical evidence that an adequately scaled @xmath394 converges to @xmath72 , we study @xmath400 , for certain @xmath401 and @xmath402 . after applying a least - square fitting ,
we choose @xmath403}{{\operatorname{mean}}(|z_m(t ) - { \operatorname{mean}}(z_m(t))|^2)}\right ) ,      \\
\mu_m = { \operatorname{mean}}(\phi(t ) ) - \lambda_m{\operatorname{mean}}(z_m(t ) ) , } \ ] ]    where @xmath404 denotes the mean over the @xmath319 + 1 points , etc .
figure [ f : convergence ] shows the absolute and relative discrepancies between @xmath72 and the scaled @xmath394 , i.e. , @xmath405 and @xmath406 , for @xmath407 .
bearing in mind that we are comparing @xmath408 different points one by one , the results give a pretty strong evidence of convergence as the number of initial sides @xmath41 tends to infinity , and strongly support the idea that @xmath409 is also a multifractal , for all @xmath410 .
( starred line ) and @xmath406 ( circled line ) , with @xmath411 and @xmath412 given by . bearing in mind that we are comparing @xmath408 points one by one , the convergence is pretty strong.,scaledwidth=50.0% ]    as we mentioned in the introduction , the proof that @xmath413 , as defined in , is a multifractal was done in @xcite . in this respect , it is fundamental the identity given for @xmath388 , as defined in , which states that @xmath414    with @xmath415 , @xmath416 , @xmath43 , @xmath417 .
this identity is proved in theorem 4.2 in @xcite .
of particular relevance is the fact that the hlder exponent is @xmath73 , i.e. , that @xmath418 lower - order terms .    in order to prove analytically that @xmath394 , is a multifractal , we would need an equivalent of .
we have made some numerical experiments ( see for example figure [ f : holder ] , for @xmath349 , @xmath185 , @xmath419 ) , and all of them give strong evidence that the hlder exponent of @xmath70 is @xmath73 for rational times , i.e. , that @xmath420 .     versus @xmath421 , for @xmath349 .
the asymptotically linear relation between both quantities as @xmath422 is evident.,scaledwidth=50.0% ]    nevertheless , how the constants depend on the denominator @xmath46 , something which is a fundamental ingredient in the arguments in @xcite and in @xcite , is unclear .
this question deserves a much more detailed analysis that we plan to make in a forthcoming paper .      in the previous subsection ,
we have given evidence of the multifractal character of the curve @xmath3 .
however , in @xcite , fractal - like phenomena were observed also on the tangent vector @xmath28 when imposing fixed boundary conditions . in the following lines
, we will show that similar observations are valid for periodic boundary conditions , too .
as we have seen , @xmath59 can be recovered algebraically up to a vertical movement , while @xmath28 can be completely recovered algebraically at rational times .
furthermore , @xmath347 is really a piecewise constant function , with no noise associated to the gibbs phenomenon .
this gives us a very powerful tool to progress in the understanding of @xmath28 ( and also of @xmath59 ) , avoiding numerical simulations at all .
a very interesting question is , given a rational time @xmath53 , with _ small _ @xmath46 , corresponding in @xmath59 to a polygon with @xmath45 or @xmath47 sides , and in @xmath28 to a piecewise continuous function with @xmath45 or @xmath47 jumps , what happens at a time @xmath425 , @xmath426 ?
let us take @xmath427 , with @xmath428 , in order that @xmath429 . assuming that @xmath430 , then , @xmath431 , so @xmath425 would correspond to a polygon with @xmath432 or @xmath433 sides .     and
@xmath347 , at @xmath434,title="fig:",scaledwidth=50.0% ]    and @xmath347 , at @xmath434,title="fig:",scaledwidth=50.0% ]    in figure [ f : xteps ] , we have have plotted @xmath435 and @xmath347 , for @xmath349 , at @xmath436 . while visually there is no difference whatsoever between @xmath435 at @xmath437 and @xmath435 at @xmath438 , we do not have a skew polygon with 6 sides , but a skew polygon with @xmath439 sides that closely resembles a polygon with 6 sides .
on the other hand , the corresponding @xmath347 are very different . indeed , in the plot of @xmath347 , we can clearly appreciate six spiral - like structures whose centers are precisely the six constant values of @xmath28 at @xmath440 .
nevertheless , it is important to underline that these structures are not really spirals , but the plot of the @xmath441 different values taking by @xmath28 , that closely resembles a curve with six spirals .
furthermore , these spirals remind us of the cornu spirals that appeared in @xcite . an open question that arises naturally is up to what extent the multiple - corner problem can be explained as a sum of several one - corner problems .     and
@xmath347 , at @xmath442 . , title="fig:",scaledwidth=50.0% ]    and @xmath347 , at @xmath442 . ,
title="fig:",scaledwidth=50.0% ]    another interesting question is what happens if we take a rational time @xmath53 with _
@xmath46 , such that there is no pair @xmath443 , with both @xmath444 and @xmath445 _ small_. in this case , the situation is very different . in figure
[ f : xteps2 ] , we have plotted @xmath435 and @xmath347 for @xmath349 , at @xmath442 . while the left - hand side is not so different from the left - hand side of figure [ f : xteps ] , the right - hand side is a set of @xmath446 points in @xmath447 that creates a spectacular fractality sensation with spiral - like structures at three or four different scales .
this can be better appreciated in figure [ f : tepszoom ] , where we have plotted the stereographic projection of @xmath28 onto @xmath67 .    in general , a rational time @xmath53 , with @xmath424 ,
can be regarded as an approximation of an irrational time . trying to fully understand the phenomena exhibited by @xmath91 at those times
is a challenging question that we postpone for the future .
in this paper , we have studied the evolution of and , for a regular planar polygonal with @xmath41 sides as the initial datum .
the algebraic calculations , backed by complete numerical simulations , prove ( under a certain uniqueness assumption ) that the curve @xmath2 is a polygon at times which are rational multiples of @xmath54 , i.e. , @xmath42 , @xmath43 , with the number of sides depending on @xmath46 , while the tangent vector @xmath91 is piecewise constant at those times . this could be seen as a nonlinear version of the talbot effect , in the spirit of the work @xcite by berry and goldberg .
we have obtained a striking connection between the curve @xmath3 and the so - called riemann s non - differentiable function . in @xcite ,
s. jaffard proved that this function is an example of multifractal whose spectrum of singularities satisfies the frisch - parisi conjecture .
although there is strong numerical evidence that @xmath3 is also a multifractal , an analytical proof seems to be challenging .
in fact , a first step in this direction is to give sense to our solutions from an analytical point of view .    in the future
, we also plan to extend these ideas to arbitrary polygons , and to do a more detailed study on fractality ; for such purposes , a more complete algebraically description of @xmath59 and @xmath28 is no doubt required
. finally , it would be interesting to prove .
we are indebted to the referees and editors for the thorough , constructive and helpful comments and suggestions .
we also wish to thank v. banica and f. chamizo for enlightening conversations .
this work was supported by the spanish ministry of economy and competiveness , with the project mtm2011 - 24054 , and by the basque government , with the project it641 - 13 .
the generalized quadratic gau sums are defined by @xmath448    for given integers
@xmath449 , with @xmath450 . from now on ,
we assume @xmath451 , and @xmath452 , which are the cases dealt with in this paper .
the value of these sums was calculated for the first time by gau for the case @xmath453 @xcite .
more precisely , given two integers @xmath454 such that @xmath452 , then @xmath455    where @xmath456 denotes the jacobi symbol . as we will see in the following lines , the generalized quadratic gau sums can be reduced to normal quadratic gau sums by completing the square . since they are multiplicative , i.e. , @xmath457    we can assume without loss of generality that @xmath458 is either odd or a power of two . if @xmath458 is odd , let us find a certain @xmath459 such that @xmath460 .
this @xmath459 is unique modulo @xmath458 and is precisely the inverse of @xmath461 in @xmath462 ; since @xmath463 , its existence is guaranteed by bzout s lemma , and it can be efficiently computed , for instance , by the extended euclidean algorithm .
then , we have @xmath464    if @xmath458 is a power of two , we first observe that @xmath465    therefore , we have two cases : @xmath466 and @xmath467 . if @xmath466 , then @xmath25 is odd , so @xmath468    on the other hand , if @xmath467 , @xmath469 , if @xmath470 is odd .
let us suppose @xmath470 is even . then , since @xmath25 is odd , we follow the previous reasoning to find a certain @xmath459 such that @xmath471 , so @xmath472    in general , given an arbitrary @xmath473 , we can factorize it as @xmath474 and use the multiplicative character of the generalized quadratic gau sums to calculate @xmath475 . for instance , if @xmath458 is even , but @xmath476 is odd , then @xmath477    i.e , @xmath478 , if @xmath470 is even , etc .    to conclude this appendix ,
let us mention explicitly the value of @xmath479 , deduced from the previous calculations , and which is of especial relevance in this paper : @xmath480