Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: flow rate measurements of liquid metals are required in various technological processes ranging from the cooling of nuclear reactors to the dosing and casting of molten metals.@xcite electromagnetic flowmeters are essential in the diagnostics and automatic control of such processes . a variety of electromagnetic flowmeters have been developed starting from the late 1940s and described by shercliff.@xcite the standard approach is to determine the flow rate by measuring the potential difference induced between a pair of electrodes by a flow of conducting liquid in the magnetic field.@xcite this approach is now well developed and works reliably for common liquids like water,@xcite but not so for liquid metals . major problem in molten metals , especially at elevated temperatures , is the electrode corrosion and other interfacial effects , which can cause a spurious potential difference between the electrodes . the electrode problem is avoided by contactless eddy - current flowmeters , which determine the flow rate by sensing the flow - induced perturbation in an applied magnetic field.@xcite the main problem with this type of flowmeters is the weak field perturbation which may be caused not only by the flow . we showed recently that the flow - induced phase shift of ac magnetic field is more reliable for flow rate measurements than the amplitude perturbation.@xcite    another contactless techniques for flow rate measurements in liquid metals is the so - called magnetic flywheel invented by shercliff,@xcite who prescribes _ a `` plurality '' _ of permanent magnets distributed equidistantly along the circumference of a disk , which is mounted on an axle and placed close to a tube carrying the liquid metal flow . the eddy currents induced by the flow across the magnetic field interact with the magnets by entraining them , which makes the disk rotate with a rate proportional to that of the flow . this type of flowmeter , described also in the textbook by shercliff@xcite and extensively used by bucenieks,@xcite was recently successfully reembodied under the name of the lorentz force velocimetry ( lfv).@xcite    recently , we suggested an alternative and much more compact design of such a flowmeter , which conversely to shercliff s flywheel uses just a single magnet mounted on the axle it can freely rotate around and magnetized perpendicularly to it.@xcite we also introduced a basic mathematical model and presented first experimental implementation of this type of flowmeter.@xcite when such a magnet is placed properly at a tube with the liquid metal flow , it starts to revolve similarly to shercliff s flywheel . but in contrast to the latter , which is driven by the electromagnetic force acting on separate magnets , the single magnet is set into rotation only by the torque . this driving torque is due to the eddy currents induced by the flow across the magnetic field . as the magnet starts to rotate , additional eddy currents are induced , which brake the rotation . an equilibrium rotation rate is attained when the braking torque balances the driving one , and this rate depends only on the flow velocity and the flowmeter arrangement , whereas it is independent of the electromagnetic torque itself . thus , the equilibrium rotation rate is affected neither by the magnet strength nor by the electrical conductivity of the liquid metal provided that the friction on the magnet is negligible . this a major advantage of the single - magnet rotary flowmeter over the lfv approach , which relies on direct force measurements.@xcite    in this paper , we present an extended theory of the single - magnet rotary flowmeter and compare it with experimental results . two limiting cases of long and short magnets are analyzed using linear - dipole and single - dipole approximations . we obtain simple analytic solutions for the force and torque on slowly moving and rotating magnets due to eddy currents in the layers of infinite horizontal extent and arbitrary depth . this allows us to find the equilibrium rotation rate of the magnet at which the torques due to the translation and rotation balance each other . we also consider an active approach , where the force on the magnet is used to control its rotation rate so that the resulting force vanishes . this rotation rate , similarly to the equilibrium one , is proportional to the layer velocity and independent of its conductivity and the magnet strength . the torque on a magnetic dipole rotating about an axis normal to a thin sheet has been calculated by smythe using an original receding image method.@xcite reitz uses this method to calculate the lift and drag forces on the coils of various geometries moving with constant velocity above a conducting thin plate.@xcite the lift force , which at high speeds approaches the force between the coil and its image located directly below it , varies as the velocity squared in the low - speed limit considered in this paper . the drag force is found to vary inversely and directly with the velocity at low and high speeds , respectively . palmer finds analytical expressions for the eddy current forces on a circular current loop moving with a constant velocity parallel to a thin conducting sheet.@xcite the force on a rectangular coil moving above a conducting slab has been calculated numerically by reitz and davis using the fourier transform method.@xcite the same problem for the magnetic dipole of arbitrary orientation placed next to a thin slowly moving slab is addressed by kirpo _ @xcite the force and torque on a transversely oriented dipole above a slowly moving plane layer of arbitrary thickness has been found analytically in the context of the lfv.@xcite fast computation of forces on moving magnets are of interested also for the eddy current force testing techniques.@xcite    this paper is organized as follows . the following section presents two simple mathematical models of the single - magnet rotary flowmeter , which are used to calculate analytically the force and torque on the magnets moving and rotating slowly above a layer of infinite lateral extent . the limits of long and short magnets , which are approximated by linear and point dipoles , are considered in secs . [ sub:2d ] and [ sub:3d ] , respectively . section [ sec : exp ] presents the flowmeter implementation details and test results . the paper is concluded by a summary in sec . [ sec : con ] . schematic view of the single - magnet rotary flowmeter with the @xmath0-axis directed out of the plane of figure . ]    consider a horizontally unbounded planar layer of electrical conductivity @xmath1 occupying the lower half - space and moving as a solid body with a constant velocity @xmath2 parallel to a permanent magnet placed at a distance @xmath3 above its surface and rotating with a constant angular velocity @xmath4 around an axis parallel to the surface and perpendicular to @xmath5 velocities are assumed sufficiently low for the magnetic field of induced currents to be negligible compared to the field of the magnet . the origin of cartesian coordinates is set at the center of the magnet with the @xmath6 @xmath0 and @xmath7 axis directed along @xmath8 @xmath9 and downward normally to the surface , respectively , as shown in fig . [ cap : sketch ] . in the following , two limiting cases will be considered in which the magnet will be assumed either much longer or much shorter than @xmath10      we start with a long cylinder magnetized perpendicularly to its axis about which it can freely rotate.@xcite in this case , the magnetic field is approximated by that of a two - dimensional ( linear ) dipole with the vector potential @xmath11 which has only the @xmath0-component @xmath12 where @xmath13 is the vacuum permeability and @xmath14 is the linear dipole moment , which is perpendicular to @xmath15 and directed at the angle @xmath16 from the positive @xmath17 axis ; @xmath18 is the radius vector from the magnet axis , and @xmath19 and @xmath20 are the cylindrical radius and the polar angle in the cylindrical coordinates around the @xmath0-axis . the magnetic field of linear dipole is @xmath21 eddy currents are induced by two effects : the translation of the layer and the temporal variation of the magnetic field due to its rotation . the latter vanishes in the co - rotating frame of reference , where the magnetic field is stationary , while the layer appears to move with the resulting velocity @xmath22 which contains also an apparent rotational motion of the layer @xmath23 opposite to that of the magnet . the density of eddy currents is given by ohm s law for a moving medium @xmath24 where @xmath25 is the electric potential . in this case , no electric potential is induced because the e.m.f . , @xmath26 is both solenoidal and tangential to the surface . if the induced magnetic field is negligible as originally assumed , eddy currents can be represented as the superposition @xmath27 where @xmath28 and @xmath29 are the currents induced by the translation and rotation , respectively . it is important to note that in the approximation under consideration with a fixed magnetic field distribution , eddy currents are determined only by the instantaneous velocities . besides that eddy currents are coplanar to the surface and , thus , mutually independent over the depth of the layer consequently , a layer of finite thickness may be represented as a semi - infinite one with zero velocity at @xmath30 where @xmath31 is the distance of the lower boundary of finite - thickness layer from the magnet . this , in turn , is equivalent to the superposition of two semi - infinite layers with the second layer at @xmath32 moving oppositely to the first one at @xmath33 so that the resulting velocity vanishes at @xmath34 in the following , this approach allows us to construct the solution for a finite - thickness layer by taking the difference of two half - space solutions , which is subsequently denoted by @xmath35_{h_{1}}^{h_{2}},$ ] where @xmath36 stands either for the force or the torque due to the eddy currents in a half - space . moreover , by the same arguments , this approach can easily be extended to @xmath7-dependent velocity distributions @xmath37 , for which general solution can be constructed as a superposition of solutions for thin layers moving with various velocities given by @xmath38 where @xmath39_{h_{1}}^{h_{2}}/\left[z\right]_{h_{1}}^{h_{2}}$ ] for @xmath40 and @xmath41    the linear force density experienced by an infinitely long magnet due to the layer translation , which according to the momentum conservation law is opposite to that acting upon the layer itself , can be written as @xmath42 where the integral is taken over the @xmath43-cross - section of the layer . the @xmath0-component of force is absent due to the @xmath44 reflection symmetry . in the low - speed limit under consideration , when force varies linearly with the velocity , there is also no @xmath7-component of force . this is the case because according to the linearity @xmath45 when @xmath46 while the latter transformation is equivalent to the rotation of the coordinate system by @xmath47 around the @xmath7-axis , which leaves the @xmath7-components invariant . in polar coordinates with the surface defined by @xmath48 we obtain @xmath49 the linear torque density , which because of the aforementioned symmetries has only the @xmath0-component , can be found as @xmath50 where @xmath51 the linear force density due to the rotation , defined by @xmath52 is found as @xmath53 the linear torque density due to the rotation , as that due to the translation above , has only the @xmath0-component @xmath54 which is not defined for a semi - infinite layer because the last integral diverges . nevertheless , expression ( [ eq : m1-gen ] ) can still be evaluated for the layer of finite depth by substituting the infinite limit in the last integral by @xmath55 which results in @xmath56_{h_{2}}^{h_{1}}.\label{eq : m1-h2}\ ] ]    the solution above becomes unbounded as @xmath57 and , thus , inapplicable to thick layers . this implies that for a half - space layer the induced magnetic field can not be neglected however slow the rotation of magnet . the induced magnetic field is related to eddy currents by ampere s law @xmath58 which combined with expressions ( [ eq : b2d ] ) and ( [ eq : ohm ] ) leads to @xmath59 as suggested by the external magnetic field ( [ eq : a2d ] ) , we search for the vector potential in the complex form @xmath60.\label{eq : a - cmp}\ ] ] then eq . ( [ eq : a ] ) for the complex amplitude @xmath61 takes the form @xmath62 where the prime stands for the derivative with respect @xmath63 the general solution of eq . ( [ eq : a - rot ] ) can be written as @xmath64 where @xmath65 is an unknown constant to be determined by matching the induced and externally imposed magnetic fields ; @xmath66 is the modified bessel function of the second kind,@xcite and @xmath67 is the skin depth due to the rotation . at distances @xmath68 smaller than the skin depth @xmath69 @xmath70 is expected to approach the imposed field @xmath71 which yields @xmath72 ) and @xmath73 substituting expressions ( [ eq : a - cmp ] ) and ( [ eq : a - hat ] ) into integral ( [ eq : m1-gen ] ) , the torque on the magnet can be represented as @xmath74 where @xmath75r\,\d r$ ] and the asterisk denotes the complex conjugate . using eq . ( [ eq : a - rot ] ) , after some algebra we obtain @xmath76=\left(\frac{\mu_{0}}{2\pi}\right)^{2}\re\left[k_{0}((1+i)r / d)k_{1}((1-i)r / d)(1+i)r / d\right].\ ] ] for the low rotation rates that satisfy @xmath77 , we obtain @xmath78=\left(\frac{\mu_{0}}{2\pi}\right)^{2}(\ln(2d / r)-\gamma),\ ] ] where @xmath79 is euler s constant . finally , taking into account that @xmath80 and @xmath81 we obtain@xmath82 where @xmath83 is the effective skin depth , which owing to the low velocities under consideration is supposed to be large relative to @xmath84 . now we can use the results above to find the magnet rotation rate depending on the velocity of layer . first , if the magnet rotates steadily without a significant friction , the torques due to the translation and rotation are at equilibrium : @xmath85_{h_{2}}^{h_{1}}=0,$ ] which yields @xmath86 note that this equilibrium rotation rate @xmath87 depends neither on the magnet strength nor on the layer conductivity unless @xmath88 in the latter case , the skin effect becomes important and , thus , @xmath31 has to be substituted by @xmath89 in the expression above . the velocity of layer can be determined also in another way by measuring the force on the magnet , as in the lfv , to control the magnet rotation rate so that the resulting force vanishes , @xmath90_{h_{2}}^{h_{1}}=0.$ ] this results in@xmath91 which again depends linearly on the layer velocity , but does not depend on its conductivity or the magnet strength .      in the other limiting case of a short magnet , when the distance to the surface is large or at least comparable to the size of magnet , the latter can be considered as a dipole with the scalar magnetic potential @xmath92 where @xmath93 is the dipole moment and @xmath94 is the fundamental solution of laplace s equation , which satisfies @xmath95 with the dirac delta function on the r.h.s . in the following , we use simplified notation @xmath96 where @xmath97 is the unit vector . then the dipole magnetic field is given by @xmath98 where the last term is added to ensure the solenoidality of @xmath99 also at @xmath100 the associated dipole current distribution is @xmath101 which easily leads to the classical expressions for the force and torque used later on . in the following , we will be using also the spherical and cylindrical coordinates associated with the cartesian ones in the usual way . as in the previous section , we change to the co - rotating frame of reference and consider the eddy currents due to translation and rotation separately .      for the translation with @xmath8 the charge conservation @xmath102 applied to eq . ( [ eq : ohm ] ) results in laplace s equation for @xmath103@xmath104 at the surface @xmath105 the normal component of electric current vanishes : @xmath106 in order to find the induced electric potential , firstly , it is important to notice that @xmath99 being a free - space magnetic field satisfies the laplace equation itself . consequently , the cartesian components of @xmath99 satisfy this equation , too . secondly , if @xmath107 in bc ( [ bc : phi0 ] ) satisfies laplace s equation , then @xmath108 satisfies not only bc ( [ bc : phi0 ] ) but also eq . ( [ eq : phi0 ] ) because the integration along a straight line , similarly to the differentiation , are interchangeable with the laplacian . by the same argument , we can interchange integration and differentiation in expression ( [ sol : phi0-gen ] ) , which yields @xmath109 where @xmath110 and @xmath111 is the zeroth degree associated legendre function of the second kind;@xcite @xmath112 is the spherical polar angle from the positive @xmath7-axis ( see fig . [ cap : sketch ] ) and @xmath68 is the corresponding cylindrical radius . the `` constant '' of integration @xmath113 which similarly to the first term @xmath114 is supposed to be axisymmetric and also to satisfy the laplace equation , is chosen to regularize @xmath115 at @xmath116 by removing the logarithmic singularity @xmath117 for @xmath118 this results in @xmath119 and@xmath120 where @xmath121 is the spherical radius . for a transversal dipole @xmath122 considered also by thess _ _ at al.__@xcite , solution ( [ sol : phi0-h ] ) reduces to @xmath123 where the last term represents the magnetic potential of a dipole aligned with the @xmath0-axis . for a general dipole orientation , expression ( [ eq : h ] ) substituted into solution ( [ sol : phi0-h ] ) after some algebra yields@xmath124,\label{sol : phi0}\ ] ] where @xmath20 is the azimuthal angle from the positive @xmath17-axis in the @xmath125-plane .    the electric potential distribution ( [ sol : phi0 ] ) allows us to calculate the force acting upon the magnet , which , as noted above , is opposite to that acting upon the layer , i.e. , @xmath126 where the integral of the lorentz force density @xmath127 is taken over the layer volume @xmath128 for the longitudinal force component , we obtain@xmath129+\int_{h}^{\infty}(\partial_{xm}^{2}g)^{2}\,\d z\right\ } \,\d s\nonumber \\   & = & \frac{\mu_{0}^{2}\sigma v}{512\pi h^{3}}(3m_{x}^{2}+m_{y}^{2}+4m_{z}^{2}),\label{eq : fx}\end{aligned}\ ] ] where the first integral is taken over the surface @xmath130 at @xmath131 and can be swapped with the second one over the layer depth . for a transversal dipole @xmath132 , the expression above coincides with that of thess _ _ et al.__@xcite as well as with the result of reitz@xcite for a thin sheet in the limit of a slowly moving dipole . as seen , the force is the strongest on a transversal dipole and reduces on longitudinal and spanwise dipoles by factors of @xmath133 and @xmath134 , respectively . thus , the force ( [ eq : fx ] ) , in contrast to force ( [ eq : f0-gen ] ) for a long magnet , varies with the dipole orientation in the @xmath43-plane . on a horizontally inclined dipole @xmath135 , there is also a spanwise force component@xmath136 but there is no vertical force whatever the dipole orientation:@xmath137 this is consistent with the result of reitz@xcite stating that in the low - speed limit the lift force is proportional to the velocity squared , while the approximation under consideration takes into account only the part of force proportional to the velocity . alternatively , the force and torque acting on dipole can be found using the associated current distribution ( [ eq : j ] ) and the induced magnetic field @xmath138 which lead straightforwardly to the classical expressions@xcite @xmath139 where the integrals are taken over the space @xmath140 above the layer . this is a bit longer but algebraically more straightforward approach , which will be pursued in the following .    in order to find the induced magnetic field @xmath138 it is important to notice that solution ( [ sol : phi0-gen ] ) satisfies condition ( [ bc : phi0 ] ) not only at @xmath131 but at any @xmath141 thus , the @xmath7-component of the current is absent not only at the surface but throughout the whole layer . then @xmath102 and @xmath142 imply @xmath143 where @xmath144 is the electric stream function , whose isolines coincide with the eddy current lines . substituting this expression into eq . ( [ eq : ohm ] ) and taking the @xmath7-component of the _ curl _ of the resulting equation , we obtain@xmath145 note that in contrast to the electric potential in eq . ( [ eq : phi0 ] ) , no boundary conditions are required for @xmath144 at the surface . this is because , firstly , eq . ( [ eq : psi - gen ] ) contains no derivatives in @xmath7 and , secondly , the absence of the @xmath7-component of current at the surface is explicitly ensured by expression ( [ eq : j - psi ] ) .    for translational motion with @xmath8 the r.h.s of eq . ( [ eq : psi - gen ] ) takes the form @xmath146 which by the same arguments as above satisfies laplace s equation and equals to @xmath147 with @xmath148 which satisfies the same equation and , thus , represents the solution of eq . ( [ eq : psi - gen ] ) . for a transversal dipole @xmath149 ) , the general solution ( [ eq : psi0-h ] ) simplifies to @xmath150 where the last term represents the magnetic potential of a longitudinal dipole . thus , in this case , the eddy current lines coincide with the isolines of electric potential ( [ sol : phi0 ] ) rotated by @xmath151 about the @xmath7-axis . for a longitudinal dipole @xmath152 the comparison of solutions ( [ eq : psi0-h ] ) and ( [ sol : phi0-h ] ) shows that the eddy current lines and the electric potential isolines are swapped with the corresponding distributions induced by a spanwise dipole @xmath153 and rotated by @xmath151 about the @xmath7-axis . for an arbitrarily oriented dipole , expression ( [ eq : h ] ) substituted into the general solution ( [ eq : psi0-h ] ) after some algebra yields@xmath154.\label{sol : psi0}\ ] ] figure [ fig : mov - xyz ] shows the isolines of the electric potential ( [ sol : phi0 ] ) and those of @xmath155 which represent the eddy current lines , in the @xmath125-plane @xmath156 for three basic dipole orientations along the @xmath17- , @xmath0- and @xmath7-axis . note that the patterns are self - similar in the @xmath125-plane with the characteristic length scaling directly with @xmath141 namely , these distributions are functions of spatial angles only , while according to ( [ sol : phi0 ] ) and ( [ sol : psi0 ] ) the magnitude of both @xmath103 and @xmath157 falls off as @xmath158    isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) . , title="fig:"](-30,0)(a)isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) . , title="fig:"](-30,0)(b )    isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) . , title="fig:"](-30,0)(c )    in order to satisfy the solenoidality condition @xmath159 the induced magnetic field is sought as @xmath160 where @xmath161 is the vector potential , which is supposed to satisfy the coulomb gauge @xmath162 then ampere s law leads to @xmath163 which , in turn , results in @xmath164 substituting expression ( [ eq : j - psi ] ) into the last integral , after some algebra we obtain @xmath165 where @xmath166 is the same as used by thess _ _ et al.__@xcite the previous expression implies that @xmath167 similarly to its source @xmath168 has no @xmath7-component . then the induced magnetic field can be written as @xmath169 where @xmath170 represents the scalar magnetic potential , which completely defines @xmath171 outside the layer , where @xmath172 further , using the tensor notation with the einstein summation convention , expressions ( [ eq : f ] ) and ( [ eq : m ] ) can be represented in terms of @xmath173 as @xmath174 where @xmath175 is the anti - symmetric tensor and the subscript after the semicolon denotes the differential with respect to the corresponding coordinate with indices @xmath176 @xmath177 and @xmath178 standing for the @xmath6@xmath0 and @xmath179 directions , respectively . taking into account the symmetry of @xmath180 with respect to the interchange of observation and integration points , the derivatives of @xmath173 at the dipole location @xmath181 in expressions ( [ eq : fi ] ) and ( [ eq : mi ] ) are found as @xmath182 where @xmath183 integrals ( [ eq : x - ij ] ) and ( [ eq : x - ijk ] ) can be evaluated analytically using , for example , the computer algebra system mathematica,@xcite which also allows us to carry out all other analytical transformations . in such a way , we firstly verify that eq . ( [ eq : fi ] ) with expression ( [ eq : psi0-h ] ) substituted into integral ( [ eq : x - ijk ] ) indeed reproduces previous results ( [ eq : fx ] ) , ( [ eq : fy ] ) and ( [ eq : fz ] ) . secondly , eq . ( [ eq : mi ] ) with the same expression for @xmath157 substituted in integral ( [ eq : x - ij ] ) results in@xmath184 which for a transversal dipole again coincides with the results of thess _ _ et al.__@xcite      for a solid - body rotation with the velocity @xmath185 which appears in the co - rotating frame of reference , the charge conservation @xmath102 applied to eq . ( [ eq : ohm ] ) results in@xmath186 while the vanishing of the normal current component at the surface @xmath131 requires@xmath187 in order to solve this problem , firstly , it is important to notice that since @xmath188 satisfies the laplace equation , @xmath189 is a particular solution to eq . ( [ eq : phi1 ] ) . comparing expressions ( [ eq : phi1-b ] ) and ( [ sol : phi0 ] ) shows that @xmath190 where @xmath103 is given by solution ( [ sol : phi0-h ] ) . searching for the solution as @xmath191 reduces eq . ( [ eq : phi1 ] ) to the laplace equation for @xmath192 while bc ( [ bc : phi1 ] ) takes the form @xmath193 by the usual arguments , we obtain @xmath194 where @xmath195 and @xmath196 is the first order associated legendre function of the second kind.@xcite using the expressions above , after some transformations we obtain@xmath197.\label{sol : phi1}\end{aligned}\ ] ] which is the electric potential in the co - rotating frame of reference . instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(a)instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(b )    instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(c )    as before , condition ( [ bc : phi1 ] ) is satisfied not only at @xmath131 but any @xmath7 and , thus , the @xmath7-component of the current vanishes throughout the layer . consequently , the electric current can again be expressed as ( [ eq : j - psi ] ) with @xmath200 for which eq . ( [ eq : psi - gen ] ) takes the form@xmath201 since the r.h.s . of eq . ( [ eq : psi1 ] ) satisfies the laplace equation , after some algebra we obtain@xmath202 the eddy current lines in the @xmath125-plane , which are the isolines of @xmath200 along with the isolines of the electric potential ( [ sol : phi1 ] ) are shown in fig . [ fig : rot - xyz ] . again , the patterns are self - similar and scale directly with the distance @xmath7 from the dipole , while the magnitude of both @xmath199 and @xmath198 falls off as @xmath203 in contrast to the translation considered above , the eddy currents induced by rotation are time - dependent and vary periodically as the dipole orientation changes from the @xmath17- to the @xmath7-axis , which are shown in figs . [ fig : rot - xyz](a ) and ( c ) , respectively . only the electric potential , but no current , is induced by a solid - body rotation when the dipole is aligned with the axis of rotation ( @xmath0 ) . in his case , when the magnetic field is symmetric about the axis of rotation , the induced e.m.f . is irrotational and , thus , compensated by the electric potential gradient , which equals to @xmath204 where @xmath205 is the cylindrical radius from the axis of rotation , @xmath206 is the azimuthal unity vector , and @xmath207 is the azimuthal component of the vector potential , which can be used to describe a general axially symmetric poloidal magnetic field . for the dipole field , we have@xmath208 where @xmath209 is the associated dipole current ( [ eq : j ] ) and @xmath210 substituting this into expression ( [ eq : v1xb ] ) and equating it to the gradient of the electric potential @xmath211 we obtain @xmath212 which is the free - space electric potential in the co - rotating frame of reference shown in fig . [ fig : rot - xyz](b ) and coinciding with ( [ sol : phi1 ] ) when @xmath213 this potential vanishes in the laboratory frame of reference , where the magnetic field is invariant with respect to the rotation around the symmetry axis . further , using expression ( [ sol : psi1 ] ) for @xmath199 in integrals ( [ eq : x - ij ] ) and ( [ eq : x - ijk ] ) , which can be evaluated together with eqs . ( [ eq : fi ] ) and ( [ eq : mi ] ) in the same way as for the translation in the previous section , we obtain @xmath214    the first point to note is that both the force and torque vanishes when dipole is aligned with the axis of rotation @xmath215 which , as discussed above , is due to the absence of eddy currents in this case . secondly , the longitudinal @xmath216 component of force ( [ eq : f1 ] ) , in contrast to that due to the translation given by expression ( [ eq : fx ] ) , is independent of the dipole orientation in the @xmath43-plane and depends only on the magnitude of the dipole moment in this plane : @xmath217 thirdly , there is also a non - zero transversal @xmath218 force , which in contrast to the longitudinal one is purely oscillatory and varies periodically with the dipole orientation in the @xmath43-plane as @xmath219 where @xmath16 is the poloidal angle of the dipole orientation in the @xmath43-plane from the positive @xmath17-axis , which is the same as that for the linear dipole in expression ( [ eq : a2d ] ) ( see fig . [ cap : sketch ] ) . moreover , when dipole is inclined to the axis of rotation ( @xmath220 also a spanwise @xmath221 force appears , which similarly to the transversal one alternates periodically around zero as @xmath222 with the dipole rotation . as seen from expression ( [ eq : m1 ] ) , inclination also gives rise to alternating torque components around both the @xmath17- and @xmath7-axis . the @xmath0-component of torque , which is negative and , thus , opposing the rotation , has not only a constant but also an alternating part , which varies with the double frequency of rotation as @xmath223    instantaneous rotation rate ( [ eq : omgt - trq ] ) versus time ( solid ) and orientation ( dashed ) for a negligible inertia of the magnet.,scaledwidth=50.0% ]    thus , in contrast to the long magnet in @xmath224[sub:2d ] , there is no steady rotation rate , which could balance the @xmath0-component of the constant torque ( @xmath225 ) due to the translation . the variation of the magnet rotation rate in response to oscillatory torque is constrained by the inertia of magnet . for the oscillatory part of the rotation rate @xmath226 to be negligible compared to the mean one @xmath227 we need @xmath228 where @xmath229 is the moment of inertia of magnet . if this condition is satisfied , the oscillatory component may be neglected in the balance of torques ( [ eq : m0 ] ) and ( [ eq : m1 ] ) around the axis of rotation , which is the @xmath0-axis . for a layer of finite thickness , this results in @xmath230 the expression above differs only by a factor of @xmath231 from result ( [ eq : omg - frc ] ) for a long magnet in the case of a vanishing longitudinal force . in the opposite limit of a negligible inertia , the torque balance yields@xmath232 this instantaneous rotation rate @xmath233 which is shown in fig . [ fig : omg - t ] versus both the time @xmath234 and orientation @xmath235 alternates between @xmath236 and @xmath237 of the mean value ( [ eq : omgb - trq ] ) @xmath238 the last relation follows from the integration of eq . ( [ eq : omgt - trq ] ) as @xmath239 which also defines @xmath87 parametrically versus @xmath234 in fig . [ fig : omg - t ] . note that the oscillations of the rotation rate do not affect its mean value @xmath240 which is defined by ( [ eq : omgb - trq ] ) independently of condition ( [ eq : omeg - osc ] ) . alternatively , when the magnet rotation rate is actively controlled so that to balance the longitudinal force ( [ eq : fx ] ) with the @xmath17-component of force ( [ eq : f1 ] ) , the instantaneous rotation rate is @xmath241 where the temporal mean value @xmath242 follows from the solution of differential equation ( [ eq : omgt - frc ] ) , which can be written as @xmath243 single - magnet rotary flowmeter with a _ smco_-type magnet ( 1 ) of @xmath244 and @xmath245 in diameter and height , respectively , installed at the side of a rectangular stainless steel duct ( 2 ) of @xmath246 in cross - section carrying the liquid sodium flow . the rotation rate is measured with an inductive magnetic proximity sensor ( 3 ) and also monitored with a photoelectric barrier sensor ( 4 ) . , scaledwidth=50.0% ]    the laboratory model of the single - magnet flowmeter shown in fig . [ exp - pic ] consists of a cylindrical _ smco_-type permanent magnet with diameter @xmath247 and height @xmath248 , magnetized perpendicularly to its axis with @xmath249 surface induction , which holds up to @xmath250 the magnet is mounted on a stainless steel axle held by a housing made of `` macor '' ceramic , which ensures a low mechanical friction and can withstand up to @xmath251    the flowmeter was tested on a sodium loop using about 90 liters of molten _ na _ at @xmath252 with the electrical conductivity of @xmath253 the flow was driven by a linear electromagnetic pump , which provided a maximum velocity of @xmath254 at the flow rate of @xmath255 in a stainless steel duct of @xmath246 cross - section .    the frequency of magnet rotation was determined with an inductive magnetic proximity sensor ( sick mm12 - 60aps - zu0 ) , which was fixed at the distance of @xmath256 perpendicularly to the side surface of magnet at its midheight as shown in fig . [ exp - pic ] . the sensor briefly switched off and then on again as the magnetic field along the axis of the sensor changed its direction . the frequency of sensor transition from the off to the on state , which happens twice per revolution , was measured with keithley 2000 multimeter using the reciprocal frequency counting techniques with the measurement gate time set to @xmath257 the frequency was also monitored with tektronix tds 210 oscilloscope . after changing the pump current , the flow was allowed to develop for @xmath258 and then six measurements were taken with the intervals of approximately @xmath259 the standard deviation in the measured rotation frequencies was typically a few tenths of @xmath260    the rotation rates measured at two gap widths between the magnet and the duct are shown in fig . [ fig : exp](a ) depending on the average sodium velocity , which was determined using ultrasonic doppler velocimetry@xcite with the resolution of @xmath261@xcite although the rotation rate is seen to increase as the magnet is approached to the duct and to vary nearly linearly with the flow velocity in agreement with the theory , the linear fit shows a certain zero offset . obviously , there is some opposing force , which has to be overcome for the rotation to start . this requires the flow velocity of at least @xmath262 first , such an opposing force is caused by the dry friction in bearings . second , the magnet is oriented by the earth s magnetic field , which has a stronger effect than the friction and also has to be overcome for the rotation to start . third , even stronger effect is caused by the stainless steel walls of the duct , which are weakly ferromagnetic.@xcite the magnet , when approached to the duct , was observed to turn with the dipole moment perpendicular to the wall . it is also important to note that the magnet , when turned by hand away from the equilibrium orientation , returned to it monotonically without oscillations . this implies the inertia of magnet to be small relative to the electromagnetic drag force . thus , the inertia can not significantly contribute to the overcoming of the orienting torque due to the wall magnetization . however , when magnet rotates slowly , the orienting torque causes the rotation rate to fluctuate , which shows up as an increased scatter in the measured values seen in fig . [ fig : exp ] at low flow velocities .    when the zero offset is removed and the remaining rotation rate is multiplied by the distance @xmath263 between the magnet axis and the liquid metal , where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall , we obtain the relative rotation rate @xmath266 which represents a rescaled slope coefficient for fig . [ fig : exp](a ) and is plotted in fig . [ fig : exp](b ) versus the velocity . the short - magnet solution ( [ eq : omgb - trq ] ) in two limiting cases of a semi - infinite @xmath267 and a thin ( @xmath268 layer yields @xmath269 and @xmath270 respectively , which are smaller than the measured values @xmath271 seen in fig . [ fig : exp](b ) . on the other hand , the long - magnet solution ( [ eq : omg - trq ] ) , which appears more adequate for the experimental setup , yields @xmath272 for @xmath273 such quantitative differences between the experiment and theory are not surprising given the simplifications underlying the latter . first , theoretical model does not take into account the finite width of the duct . but this alone could hardly explain the high rotation rate of the magnet at which its field travels faster than the layer . this apparently being the case in the experiment implies the presence of significant velocity gradients , which are also ignored by the theory , but could be taken into account as outlined in sec . [ sub:2d ] . note that strong velocity gradients can be caused by the magnet itself , which due to its large size and strength may act as a magnetic obstacle partially blocking the flow and so increasing its local velocity .    measured ( a ) and rescaled ( b ) magnet rotation rates along with the best linear fits versus the average velocity for @xmath274 and @xmath275 gaps between the magnet and the duct with the distance between the magnet axis and liquid metal @xmath276 where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall.,title="fig:",scaledwidth=50.0%](-30,30)(a)measured ( a ) and rescaled ( b ) magnet rotation rates along with the best linear fits versus the average velocity for @xmath274 and @xmath275 gaps between the magnet and the duct with the distance between the magnet axis and liquid metal @xmath276 where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall.,title="fig:",scaledwidth=50.0%](-30,30)(b ) we have presented a theory of single - magnet rotary flowmeter for two limiting cases of long and short magnets , which were modeled as linear and single dipoles , respectively . simple analytical solutions were obtained for the force and torque on slowly translating and rotating magnets due to eddy currents in layers of arbitrary depth and infinite horizontal extent . the velocity was assumed to be constant and the motion so slow that the induced magnetic field could be neglected . the latter assumption was not applicable to the long magnet rotating above a conducting half - space . in this case , to obtain a finite braking torque , the skin effect due to the induced magnetic field had to be taken into account however slow the rotation . for a single dipole of arbitrary orientation , compact analytical solutions were obtained in terms of both the electric potential and stream function induced by the layer translation and rotation in the co - rotating frame of reference . the electric stream function was used further to find the scalar potential of induced magnetic field at the dipole location , which resulted in simple expressions for the force and torque on the dipole . eventually , we found the equilibrium rotation rate at which the driving torque due to the layer translation is balanced by the braking torque due to the magnet rotation . an alternative approach was also considered , where the force on the magnet could be used to control its rotation rate so that the resulting force vanishes . in either case , the resulting rotation rate is directly proportional to the layer velocity and inversely proportional to the distance between the magnet and the liquid metal . these results were found in a qualitative agreement with the measurements on the liquid sodium flow . a more accurate quantitative agreement with experiment is limited due to the substantial approximations underlying the theoretical model , which neglects the finite lateral extension of the layer as well as the spatial and temporal variations of the velocity distribution .    in conclusion , note that the resulting rotation rate is independent of the magnet strength and the electrical conductivity of the liquid metal provided that the mechanical friction or other external effects are negligible compared to the driving torque . this is the main advantage of rotary flowmeter over the lfv.@xcite              x. fu , l. hu , k. m. lee , j. zou , x. d. ruan , and h. y. yang , dry calibration of electromagnetic flowmeters based on numerical models combining multiple physical phenomena ( multiphysics ) , _ j. appl . phys . _ * 108 * , 083908 - 11 ( 2010 ) . i. bucenieks , modelling of rotary inductive electromagnetic flowmeter for liquid metal flow control , in : _ proceedings of the 5th int . symposium on magnetic suspension technology , _ 204208 ( dresden , germany , 2005 ) . g. gerbeth , j. priede , i. bucenieks , a. bojarevics , j. gelfgat , and d. buchenau , verfahren und anordnung zur kontaklosen messung des durchflusses elektrisch leitfhiger medien , _ german patent application _ de 102007046881a1 ( 2009 ) . j. priede , d. buchenau , and g. gerbeth , force - free and contactless sensor for electromagnetic flowrate measurements , in : _ pamir conference _ , 815820 ( presqule de giens , france , 2008 ) ; _ magnetohydrodynamics _ * 45 * , 451458 ( 2009 ) . s. eckert , d. buchenau , g. gerbeth , f. stefani and f .- weiss , some recent developments in the field of measuring techniques and instrumentation for liquid metal flows , _ j. nucl . tech . _ * 48 * , 490498 ( 2011 ) .
Abstract:

Answer: we present a theory of single - magnet flowmeter for liquid metals and compare it with experimental results . the flowmeter consists of a freely rotating permanent magnet , which is magnetized perpendicularly to the axle it is mounted on . when such a magnet is placed close to a tube carrying liquid metal flow , it rotates so that the driving torque due to the eddy currents induced by the flow is balanced by the braking torque induced by the rotation itself . the equilibrium rotation rate , which varies directly with the flow velocity and inversely with the distance between the magnet and the layer , is affected neither by the electrical conductivity of the metal nor by the magnet strength . we obtain simple analytical solutions for the force and torque on slowly moving and rotating magnets due to eddy currents in a layer of infinite horizontal extent . the predicted equilibrium rotation rates qualitatively agree with the magnet rotation rate measured on a liquid sodium flow in stainless steel duct .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the present paper is a sequel of a previous work@xcite concerning the evaluation of saddle - point integrals of the form @xmath0\,\mathrm{d}s , \label{sd.1}\ ] ] where @xmath1 is a suitable integration path in the complex @xmath2-plane , @xmath3 and @xmath4 are functions which , for simplicity , will be assumed to be nonsingular , and @xmath5 will be intended as a `` large '' ( in modulus ) complex parameter . as is well known , the numerical evaluation of integrals of the kind in eq . ( [ sd.1 ] ) is customarily required for solving several classes of physical problems , occurring in optics , quantum mechanics , statistical physics , fluid mechanics , and so on . in optics , the evaluation of several diffraction integrals is customarily carried out asymptotically by identifying the parameter @xmath5 as the wavenumber of the radiation . in quantum mechanics , the same role is played by the inverse of the planck s constant , while in fluid mechanics by the reynold s number@xcite .    in the stationary phase treatment of diffraction integrals the values of the associated complex wavefield are asymptotically evaluated by taking the contributions coming from the stationary points of @xmath4 , each of them associated to a  ray " in the corresponding geometrical picture . of particular importance is the birth and the death , as the spatial parameters ruling the  phase " function @xmath6 vary , of  evanescent " rays across sets of codimension 1 , named  stokes sets " @xcite . the @xmath7- , or weniger , transformation@xcite ( wt for short ) , is particularly efficient for resumming the factorial divergent asymptotic series well away from stokes sets , as well as sets where two or more saddles are symmetrically placed in the complex singulant space@xcite . unfortunately , as with other resummation techniques@xcite , the wt fails to perform across stokes sets . the reason of such a failure stems to the extreme `` specialization '' of the transformation itself , which requires , for a successful resummation , an alternating sign pattern of the sequence of the single terms of the series@xcite . several methods have been conceived for resumming nonalternating , slowly convergent or divergent , sequences @xcite , some of them being based on the serial combination of various resummation techniques@xcite . for the class of saddle - point integrals in eq . ( [ sd.1 ] ) , the marriage between hyperasymptotics@xcite ( h for short ) and the wt @xcite , generating the so - called h - wt ( which stands for hyperasymptotic - weniger transformation ) , allows the wt to successfully operate also across stokes sets .    basically , the h - wt consists in the sequential application , to the integral in eq . ( [ sd.1 ] ) , of a classical hyperasymptotic treatment , as described in ref . @xcite , followed by the action of the wt on all asymptotic divergent series generated by h. in particular , the results obtained have shown how the 1st - order h - wt , for which only the first - stage of h is anticipated to the wt , is able to provide relative errors several orders of magnitude smaller than those achievable via the use of full hyperasymptotic treatments and with considerably lighter computational complexity and effort . a key aspect is that , differently from h , the first truncation operated on the starting asymptotic series has not to be an optimal , in the sense of superasymptotics ( i.e. , at the least term ) one , but rather the corresponding truncation order , say @xmath8 , must be used as a free parameter for the subsequent application of the wt .    a question which was not addressed in ref . @xcite , but mentioned only in the last sentence , is whether wt and h can be combined to higher orders in h , and if so , how the accuracy improves with order . the present paper is aimed at giving a first answer to such a question . we shall limit our analysis only to the second stage of h. further increasing of the h - wt order would be achievable along the same guidelines outlined here . on the other hand , it should be noted how , on increasing the order of h , the number of asymptotic series associated to the corresponding remainder that have to be resummed exponentially grows for topologies involving more than two saddles and , at the same time , the number of free parameters ( i.e. , the truncation orders at each h step ) linearly increases . accordingly , from a mere computational viewpoint , it is mandatory to find a compromise between the h - wt order and the computational effort . some limits of the 1st - order h - wt have already been emphasized in ref . @xcite , where asymptotic evaluations of saddle - point integrals for  small " values of the asymptotic parameter were considered . in such cases , in fact , to increase the parameter @xmath8 of the 1st - order h - wt does not necessarily lead to an improvement of the reached accuracy , but often results in the opposite , i.e. , a worsening of it . it is just the above scenario that we are interested in when the 2nd - level h - wt will be developed . numerical experiments will be carried out on the class of saddle - point integrals already considered in the numerical sections of the first paper @xcite . moreover , asymptotic evaluations of the so - called swallowtail diffraction catastrophe @xcite will be proposed as a new numerical experiment . the swallowtail function is defined via eq . ( [ sd.1 ] ) with @xmath4 being a 5th - order polynomial with respect to @xmath2 , thus involving a four saddle network . we will present a study of the accuracy achievable via h - wt asymptotic evaluations of the swallowtail diffraction catastrophe for points placed at the stokes set , following the prescriptions by berry and howls@xcite . in doing so , we will find that the corresponding asymptotic expanding coefficients can be expressed in closed - form terms . for any practical implementation of the h - wt , a key role is played by the numerical evaluation of the corresponding hyperterminants@xcite which are defined through suitable multiple integrals . for the lowest - order hyperterminant the exact analytical expression is available from literature@xcite , but unfortunately this is not true for higher - order hyperterminants , included those involved in the 2nd - level h - wt . in the present paper we solve the problem of the 2nd - level hyperterminant exact evaluation for a particular , but very important , choice of the hyperterminant parameters , which often occurs in the implementation of h for evaluating a wide class of saddle - point integrals . up to our knowledge , this is a new result which also provides an interesting connection of such hypeterminants to the meijer - g functions @xcite . moreover , although the closed - form evaluation of 2nd - level hypeterminants for arbitrary choices of their parameters seems to remain an open problem , in the present paper we find a semi - analytical representation which turns out to be suitable for numerical calculations via standard integration packages . similarly as done in ref . @xcite , one of our aims is to keep the paper reasonably self - consistent . accordingly , in the next section a brief review of h , up to the 2nd - level , is given . as far as the wt is concerned , we believe that what is contained in ref . @xcite , together with the extensive bibliography , should be enough also for a nonexpert reader . for this reason , we do not repeat it in the present paper . for simplicity , we shall refer to the asymptotic evaluation of saddle - point integrals of the type in eq . ( [ sd.1 ] ) where the set of saddle points of @xmath4 will be denoted @xmath9 and the integration path @xmath1 will be thought of as the union of a finite number of steepest descent arcs each of them , say @xmath10 , passing through the contributive saddle point @xmath11 , which will be supposed to be a simple one . accordingly , the quantity @xmath12 can generally be written as @xmath13\,\mathrm{d}s= \displaystyle\sum_{n \in \mathcal{s}'}\,\mathcal{i}^{(n)}(k ) , \label{sdreview.5}\ ] ] where @xmath14 denotes the subset of @xmath9 containing all the contributive saddles , and @xmath15\,\mathrm{d}s . \label{sdreview.5.1}\ ] ] the last integral can be written as@xcite @xmath16 where @xmath17 , and where @xmath18 can _ formally _ be written through the following asymptotic series expansion : @xmath19 the expanding coefficients @xmath20 being expressed via the integral representation @xcite @xmath21^{r+1/2}}\,\mathrm{d}s , \label{sdreview.5.0}\ ] ] where the subscript @xmath22 denotes a small positive loop around the saddle @xmath11 . h starts by writing eq . ( [ sdreview.4 ] ) in the form @xmath23 where @xmath8 represents a positive integer and @xmath24 , denotes the corresponding remainder which , due to the diverging character of the asymptotic series , turns out to be a _ diverging _ quantity too . h is based on a couple of fundamental results , found via a nontrivial analysis in ref . the first is that the value of the expanding coefficients @xmath20 at the saddle @xmath11 is strictly related to the values of the expanding coefficients @xmath25 at all those saddles , say @xmath26 , which are _ adjacent _ to @xmath11 , via the following formal resurgence relation@xcite : @xmath27 where @xmath28 denotes the set containing the indexes pertinent to all saddles adjacent to @xmath11 , the quantities @xmath29 , called _ singulants _ , are defined by @xmath30 and the binary quantities @xmath31 are obtained through a topological rule@xcite . the other fundamental tool of h is the following integral representation of the remainder @xmath32@xcite : @xmath33 equations ( [ sdreview.4])-([remainder.2 ] ) allow hyperasymptotic expansions for the saddle integral in eq . ( [ sdreview.5.1 ] ) to be built up in principle to any order@xcite . for instance , the direct substitution of eq . ( [ sdreview.4 ] ) into eq . ( [ remainder.2 ] ) leads to the 1st - stage hyperasymptotic expansion @xcite , @xmath34 where the function @xmath35 , called _ hyperterminant _ of order 1@xcite , is defined through the integral @xmath36 where , in order for it to converge , @xmath37 . moreover , it can be shown that@xcite @xmath38 where @xmath39 denotes the exponential integral function @xcite , while @xmath40 equals 1 if @xmath41 and zero otherwise . the presence of the term containing @xmath40 has to be ascribed to the evaluation of the integral in eq . ( [ remainder.3.1 ] ) , when @xmath41 , in the cauchy principal - value sense . equation ( [ istage.complete ] ) represents the first hyperasymptotic stage , at which the divergence of the asymptotic series in eq . ( [ sdreview.4 ] ) is led back to the presence of adjacent saddles @xcite . furthermore , the asymptotic series in eq . ( [ istage.complete ] ) are only formal , since for @xmath42 the terminant @xmath43 diverges . in ref . @xcite eq . ( [ istage.complete ] ) was taken as the starting point for introducing the h - wt . in particular , instead of using the wt directly on the single terms of the series in eq . ( [ sdreview.4 ] ) , it is employed for resumming the asymptotic series associated to all saddles @xmath44 , with @xmath45 , which appear in eq . ( [ istage.complete ] ) . of course , due to the fact that @xmath46 in eq . ( [ istage.complete ] ) , it is mandatory that @xmath8 be left as a free parameter , in order for the wt to be able in decoding the above asymptotic series . the 2nd - level h can be derived by truncating each of the asymptotic series in eq . ( [ istage.complete ] ) at an order , say @xmath47 , and by generating , for _ each _ adjacent saddle @xmath44 , with @xmath48 , a list of asymptotic series associated to all saddles , say @xmath49 , such that @xmath50 . in appendix [ iistage ] , only for the reader convenience , the derivation of the 2nd - level hyperasymptotic expansion of the integral in eq . ( [ sd.1 ] ) is briefly recalled according to the formalism of ref . in particular , it is found that @xmath51 where @xmath52 , the hyperterminant of order 2 , is now defined through the double integral @xmath53 similarly as we done in ref . @xcite , eq . ( [ iistage.complete ] ) can be used to give estimates of @xmath18 , as functions of the two ( free ) parameters @xmath8 and @xmath47 , by resumming , via the wt , all asymptotic series generated at the second stage of h which are inside the double sum with respect to @xmath54 and @xmath55 . the numerical evaluation of the hyperterminants represents a step of fundamental importance for any practical implementation of the h - wt algorithm . unfortunately , differently from the lowest - order h - wt , for which the corresponding hyperterminants are achievable via the closed - form expression in eq . ( [ remainder.3.1.1 ] ) , there are no analytical expressions available for higher - order hyperterminants . in a series of important papers , olde daalhuis  @xcite addressed the general problem of the hyperterminants evaluation , up to arbitrary precisions , through the use of convergent series representations based on hypergeometric functions . however , for the particular case of the 2nd - level hyperterminant , it seems that some new , at least up to our knowledge , results could be established .    from eq . ( [ terminant2 ] ) where , in order for it to converge , @xmath37 and @xmath56 , the hyperterminant can be written as @xmath57 and , by formally expanding the factor @xmath58 as a geometric series , after some algebra takes the form @xmath59 where @xmath60 denotes the hypergeometric function@xcite . although the series in eq . ( [ k2.0.2 ] ) is divergent , it can be decoded via borel summation , i.e. , by replacing the term @xmath61 by its integral representation , i.e. , @xmath62 which , once substituted into eq . ( [ k2.0.2 ] ) , leads to @xmath63 although , as we shall in a moment , it is possible to express the series inside the last equation through a closed form , it is better to carry out the evaluations for the case @xmath64 and @xmath65 separately .    on letting into eq . ( [ k2.0.3.1 ] ) @xmath64 we have @xmath66 the integral in eq . ( [ k2.0.3.2 ] ) can be evaluated by using the representation of the hypergeometric function given by formula 9.34.7 of ref . in particular , it turns out that @xmath67 where @xmath68 denotes the meijer function@xcite . finally , by using formulas 9.31.5 , 7.813.1 , and 9.31.2 of @xcite , after some algebra it is found that @xmath69 equation  ( [ terminant2.1 ] ) represents one of the main results of the present paper . as we shall see in the numerical section , in applying the 2nd - level h - wt the evaluation of the hyperterminants @xmath52 is often required for @xmath64 . this happens whenever the contributive saddle @xmath70 turns out to be adjacent to itself after two hyperasymptotic stages , i.e. , when @xmath71 into eq . ( [ iistage.complete])@xcite .    for @xmath65 , the series inside the integral in eq . ( [ k2.0.2 ] ) can still be expressed through a closed form , although the subsequent integral unfortunately not . however , a semi - analytical expression , which turns out to be suitable for being evaluated via standard numerical integration packages can be derived . in doing this , eq . ( [ k2.0.2 ] ) is first rewritten as @xmath72 , \end{array } \label{gn1.1}\ ] ] where @xmath73 so that the task is to evaluate the series in eq . ( [ gn1.2 ] ) for @xmath65 . on substituting from eq . ( [ k2.0.3 ] ) into eq . ( [ gn1.2 ] ) we have @xmath74 where use has been made of the relation[see ref .  @xcite p.  347 ] @xmath75 finally , on writing eq .  ( [ gn1.3 ] ) as @xmath76 where @xmath77 denotes the pochhammer symbol , formula 6.7.1.8 of ref . @xcite gives at once @xmath78 notice that , when @xmath79>0 $ ] , the function @xmath80 can also be evaluated through the alternative form @xmath81}{p}\ , f(1,1;n+1;1-\gamma p ) . \end{array } \label{gn1.7}\ ] ] although it seems that the above expressions can not be further simplified , the numerical evaluation of the function @xmath80 can be done with high accuracies by using standard integration packages . finally , it should be stressed that , for @xmath41 , the evaluation of the double integral defining @xmath52 has to be done , with respect to the @xmath82-variable , in the cauchy principal value sense , in order to overcome the singularity placed at @xmath83 . this , in turn , implies that an extra term must be added to the result . in appendix [ appc ] such term is analytically evaluated starting from the definition in eq . ( [ terminant2 ] ) , and turns out to be @xmath84 } { ( -\beta)^{n-1 } }          \,e_{n}(-\beta\gamma ) . \end{array } \label{extraterm}\ ] ] + all subsequent numerical experiments will be done within the _ mathematica _ language . consider the evaluation of the airy function , defined as @xmath85\ , \mathrm{d}s , \label{airy.1}\ ] ] which is of the form given in eq . ( [ sd.1 ] ) with @xmath86 , @xmath87 , and @xmath88 . the detailed analysis of the saddle topology , as well as the expanding coefficients @xmath20 has been summarized in ref . @xcite , so that here it will not be given again . we only recall that the two saddles are @xmath89 and @xmath90 , and that @xmath91 , @xmath92 , @xmath93 . we focus our attention on the evaluation of the airy function across the stokes line@xcite , i.e. , for @xmath94 , in order to compare the performances of the 2nd - level h - wt with respect those displayed by the 1st - order h - wt in the same situation . more precisely , we write the argument of the airy function as @xmath95 , where @xmath96 is a real positive parameter , whose value coincides with the singulant @xmath97 . the study of the asymptotic evaluation of the airy function across its stokes line has played a pivotal role in the development of several asymptotic technique , mainly in light of the relative simplicity of the involved saddle topology . such a simplicity could help in grasping , whereas possible , some conceptual aspects related to the use of the h - wt . differently from what done in ref . @xcite , where the relative error values were displayed via the use of tables , in the present paper we are going to resort to graphical visualizations , due to the presence of the two `` free '' parameters @xmath8 and @xmath47 . in the first experiment , whose results are shown in figure  [ figerrorairyf16 ] , the airy function is evaluated for @xmath98 . note that the same experiment was carried out in ref . @xcite via the use of h. the values of the relative error , obtained through the 1st - level h - wt , are shown , as black dots , versus the values of @xmath8 , reported on the abscissa axis . for each value of @xmath8 , the values of the relative error obtained via the 2nd - level h - wt , with @xmath99 $ ] , are also plotted and , for the sake of clarity , are joined with lines of different color each of them , which corresponds to a different value of @xmath8 , departing from @xmath8 itself . this can be noted from the figure , where it is immediately seen how the higher the @xmath8 , the longer the corresponding coloured `` leg '' is . from a first look to the figure , it appears that the relative error , obtained with both the 1st- and the 2nd - level h - wt , is lower bounded . we shall find that all subsequent numerical experiments present the same characteristic . as a first remark , it should be noting that the improvement of the estimate accuracy induced by the 2nd - level h - wt with respect to that obtained via the 1st - level one of the same order appears , at least for values of @xmath8 not too large , not to adequately refund the unavoidable increase of the computational complexity required by the application of the 2nd - level transformation . for example , it is seen from fig . [ figerrorairyf16 ] that a relative error of the order of @xmath100 , achieved through the 2nd - level h - wt with @xmath101 and @xmath102 , would be reached via a 1st - order h - wt by letting @xmath103 , but with a considerable saving of computational effort . the above example clearly suggests the use of the 2nd - level h - wt only in those cases where the best accuracy attainable via the 1st - level h - wt turns out to be not adequate . this happens , for instance , when the integral is attempted to be evaluated beyond the asymptotic realm . to put into evidence this aspect , figure  [ figerrorairyf14to2 ] shows the same as in fig . [ figerrorairyf16 ] , but for a decreasing sequence of values of @xmath96 , namely 14 ( a ) , 10 ( b ) , 6 ( c ) , and 2 ( d ) . in particular , in fig . [ figerrorairyf14to2](d ) , where the airy function argument is located at a distance @xmath104 from the origin of the complex plane , the 1st - level h - wt provides a best error of the order of @xmath105 , achieved for @xmath106 . higher accuracies are not allowed because the information gained at the first h stage turns out to be no longer sufficient to generate wt - resummable sequences . the 2nd - level h - wt , on the other hand , provides a best error of the order of @xmath107 , which is attained for @xmath108 . some intuitive insights about the resummation process associated to the h - wt could be grasped by noting that a lower bounded error is an intrinsic imprint of superasymptotic and hyperasymptotic resummations , whereas it is not generally featured by the application of the wt to alternating factorial divergent series@xcite . accordingly , one should be inclined to think that such an error behavior could be ascribed to the presence of the  regularization " step operated by h on the raw input data . speaking within a more general context , this should be somewhat related to the possible presence on nonanalytic , nonperturbative correction terms which can not be grasped simply by resummation processes , but rather require the use of  generalized nonanalytic expansions " @xcite . + in a second experiment concerning the airy function , the asymptotic parameter @xmath96 is let running within the interval @xmath109 $ ] and , for each value of @xmath96 , an exhaustive search of the optimal values of the truncations @xmath8 and @xmath110 , which minimize the 1st- and 2nd - level relative errors , respectively , is operated . the results are shown in fig . [ figoptimalerrorairy ] , where the optimal relative errors obtained via the 1st- ( open circles ) and the 2nd - level ( dots ) h - wt are shown as functions of @xmath96 . the values of the optimal truncation @xmath8 for the 1st - level h - wt are also reported , versus @xmath96 , in fig . [ figoptimalnairy ] , while those of @xmath8 and @xmath47 , for the 2nd - level h - wt , in fig . [ figoptimalnmairy](a ) and ( b ) , respectively . we will come back later on the above results . + in concluding the present section , however , we want to provide a table of explicit values obtained through the use of the 2nd - level h - wt . we choose to evaluate the airy function for @xmath111 , for which the optimal setting of the truncation parameters turns out to be @xmath108 . the preliminary step is the evaluation , through a simple wt , of the contribution , to the airy integral , coming from the saddle @xmath112 . the result is shown in table [ table.1 ] . furthermore , the subsequent action of the 2nd - level h - wt on the saddle @xmath113 is shown in table [ table.2 ] , where the complete estimates of the airy function provided are reported together with the corresponding values of the truncation @xmath114 $ ] , with @xmath115 . the second numerical experiment concerns the evaluation of the instanton integral @xmath116\,\mathrm{d}s , \label{inst.1}\ ] ] with @xmath117 , already considered in ref . @xcite as a simplified prototype for the modeling of instanton tunneling between symmetric double wells . it was shown in ref . @xcite that the integral in eq . ( [ inst.1 ] ) can be written as@xcite @xmath118 where , by referring to eq . ( [ sd.1 ] ) , @xmath119 , @xmath120 , and where @xmath1 is the steepest descent path connecting the points @xmath121 and @xmath122 via the lines @xmath123 and @xmath124 . the complete saddle topology , as well as the expressions of the expanding coefficients associated to all saddles have been described in ref . in particular , there are three saddles , @xmath125 , @xmath126 , and @xmath127 , with @xmath128 , @xmath129 , and @xmath130 . the saddles involved in the evaluation of @xmath12 are @xmath112 and @xmath131 , but only the latter requires a h - wt treatment , since the associated singulant is @xmath132 , and the corresponding asymptotic series turns out to be nonalternating . furthermore , @xmath133 , @xmath134 , and @xmath135 , while we recall that the integral in eq . ( [ inst.1 ] ) can be expressed in closed form via @xmath136 , \label{inst.4}\ ] ] where @xmath137 denotes the @xmath22th - order modified bessel function of the first kind . the first experiment concerns the evaluation of @xmath138 via the 1st- and the 2nd - level h - wt . in fig . [ figinstantonek1ov2 ] it is seen how the 2nd - level relative error is bounded , with a minimum value of the order of @xmath139 , achieved for @xmath140 . on the opposite , the 1st - level h - wt turns out to be completely inadequate to provide a reasonably accurate estimate of the function , due to the very low value of @xmath5 . the searching for optimal values has also been carried out in the present case , but using , as the varying asymptotic parameter , @xmath141 $ ] , i.e. , where the 1st - order h - wt displays the worst results in terms of accuracy , as shown in fig . 2a of ref . the error values are shown in fig . [ figoptimalerrorinstanton ] , versus @xmath5 , while the optimal settings of @xmath8 and of @xmath110 are plotted , against @xmath5 , in fig . [ figoptimalninstanton ] and  [ figoptimalnminstanton ] , respectively . it is worth now comparing the results pertinent to the airy and the @xmath142 functions . what we are going to show can seem at first sight somewhat surprising , but gives a possible first hint toward the understanding of the h - wt mechanisms . for simplicity we shall refer to the 1st - level transformation , but the results will apply also to higher - order levels . in fig . [ figcomparisonairyinst ] , the values of the relative error obtained for the airy function ( dots ) and for the instanton function ( solid curve ) are plotted , versus @xmath8 , when the values of the parameter @xmath96 and the parameter @xmath5 are numerically equal . in particular , figure ( a ) corresponds to @xmath143 , ( b ) to 7 , ( c ) to 12 , and ( d ) to 20 . it is clearly seen that the behavior of the relative error follows basically the same law . to give a possible explanation of this , in fig . [ figcompairyinst ] a pictorial representation of the complete saddle network and the complex integration path involved in the evaluation of the airy ( a ) and of the instanton ( b ) functions is plotted . in both pictures , the black dot denotes the saddle for which the h - wt is required . although the two saddle distributions are clearly different , they present some common features that , together with eq . ( [ istage.complete ] ) , are enough to justify what happens in fig .  [ figcomparisonairyinst ] . each of the `` black '' saddles is adjacent to a single saddle . for the airy function @xmath144 is adjacent to @xmath145 , while for the instanton function @xmath146 is adjacent to @xmath145 . the values of the corresponding singulants are @xmath96 and 1 , respectively . the use of the resurgence relation in eq . ( [ resurgence.1 ] ) now gives , for the two `` black '' saddles , @xmath147 for the airy function and @xmath148 for the instanton function . from the above equation it is seen that the behavior of the expanding coefficients follows the same asymptotic law as soon as @xmath149 . at the same time , however , the above equality guarantees that also the asymptotic laws for the expanding coefficients corresponding to the adjacent saddles is identical . in fact , for the airy function the saddle adjacent to @xmath145 is @xmath144 itself , with a singulant value of @xmath150 . as far as the instanton function is concerned , the saddles adjacent to @xmath145 are @xmath144 and @xmath146 , but for both of them the singulant values are -1 . accordingly , the use of eq . ( [ resurgence.1 ] ) , together with the condition @xmath149 , provides again an equivalence between the asymptotic laws of @xmath151 for the airy function and @xmath151 for the instanton function . finally , on using eq . ( [ istage.complete ] ) it is not difficult to convince that the retrieving process is the same for the two functions at the 1st - level@xcite . leaving a deeper understanding of this phenomenon to future investigations , it is here worthwhile to point out how an immediate consequence of the above described  topological equivalence " could be the restriction of the study of the h - wt retrieving performances to a few classes of prototype test cases .      as a last numerical experiment we consider asymptotic evaluations of the so - called swallowtail diffraction catastrophe @xcite , which is defined through the following integral : @xmath152\ , \mathrm{d}s , \label{sw.1}\ ] ] which is of the form given in eq . ( [ sd.1 ] ) with @xmath119 , @xmath153 , and @xmath88 . the integration path @xmath1 can be thought of as the union of steepest descent paths approaching , for @xmath154 , the directions @xmath155 , with @xmath156 . although a systematic treatment of the swallowtail asymptotics , along the general classical rules recalled in sec . [ notations ] , can be derived by paralleling the analysis carried out , for the pearcey function , in ref . @xcite , up to our knowledge it is not present in the current literature . as shown in appendix [ swt ] , all expanding coefficients @xmath20 are given by @xmath157 where @xmath158 and where the polynomials @xmath159 are defined via the generating function formula @xmath160 it is not difficult to prove that eq . ( [ generatingc.3 ] ) allows the numerical evaluation of the polynomials @xmath161 to be efficiently performed via the use of the following recurrence rule , whose derivation is outlined in appendix [ recurrencerule ] : @xmath162 with the triggering values @xmath163 , @xmath164 , and @xmath165 . the numerical experiments we are going to illustrate concern asymptotic evaluations of @xmath166 at points belonging to the corresponding stokes set , which has been estensively studied in ref . @xcite ( see , in particular , fig . 3 of this reference ) . accordingly , the triplets @xmath167 have been chosen following the prescriptions given in ref . @xcite , in order to investigate points at the intersection between the stokes surface and the plane @xmath168 , along the branch corresponding to @xmath169 and @xmath170 . this leads to triplets of the form @xmath171 , with @xmath172 being a positive parameter@xcite . we start by considering the case @xmath173 . the saddle topology is constituted by four saddles , which are listed below together with the corresponding list of adjacent ones : @xmath174 with the orientation matrix being @xmath175 , \label{sw.1.2}\ ] ] where the values of @xmath176 which are not relevant for the present experiment have been replaced by dots . three of the four saddles , namely @xmath145 , @xmath146 and @xmath177 , do contribute to the integral . in particular , the integration path consists in the union of three steepest - descent arcs , the first connecting the points @xmath178 and @xmath179 passing through @xmath131 , the second connecting the point @xmath179 to the saddle @xmath112 , passing through @xmath180 , and the third connecting the saddle @xmath112 to the point @xmath181 . accordingly , the stokes phenomenon occurs via the so - called _ saddle connection _ between saddles @xmath180 and @xmath112@xcite . a pictorial representation of the topology above described is given , for the reader convenience , in fig . [ figswpaths ] .    before showing the numerical results about the performances of the 1st- and the 2nd - level h - wt , it is worth giving some details about the way the saddle topology of the swallowtail integral influences the retrieving capabilities of the wt . we refer , in particular , to the contribution , to the swallowtail integral , of @xmath146 and @xmath177 . as far as the former is concerned , since there is only one adjacent saddle , @xmath177 , due to the fact that the corresponding singulant @xmath182 is purely imaginary , it turns out that @xmath183 , thus allowing the wt to operate the resummation . similar considerations can be done for the contributive saddle @xmath145 . the situation is somewhat different for the saddle @xmath177 , which is connected to @xmath145 . in fact , from the above described topology , it turns out that the adjacent saddle @xmath146 is dominant , i.e. , presents the minimum value of @xmath184 , with @xmath185 . accordingly , one should conclude also for @xmath186 an asymptotic `` factorial divided by power '' law , similar to that corresponding to @xmath187 and , due to the fact that @xmath188 , one should expect the wt to be able in resumming the corresponding asymptotic series . but this , on the contrary , does not happen , because of the presence of the other two , nondominant , saddles @xmath144 and @xmath145 , symmetrically placed in the complex singulant space . this can be explained by expliciting eq . ( [ resurgence.1 ] ) as @xmath189 , \label{resurgence.1bis}\ ] ] where @xmath190^{*}$ ] . this example shows how the divergent asymptotic series generated on stokes sets do not necessarily display a strictly nonalternating sign pattern as , for example , happened for the airy function , but rather how the asymptotic behavior of their single terms can display more complex patterns , depending on the whole saddle topology . figure  [ figswk2 ] shows the relative errors obtained through the 1st- and the 2nd - level h - wts in the case of swallowtail evaluation across the stokes set defined above , when @xmath173 . as for figs . [ figerrorairyf16 ] ,  [ figerrorairyf14to2 ] , and   [ figinstantonek1ov2 ] , the errors are plotted versus the parameters @xmath8 and @xmath47 . in this case their optimal values turn out to be @xmath106 ( for the 1st - level ) and @xmath191 ( for the 2nd - level ) , with corresponding error values of @xmath192 and @xmath193 , respectively , evaluated with respect to the  exact " value of @xmath194 , obtained via the method recently proposed in ref . finally , an experiment about optimal resummation of swallowtail functions has been carried out , by using @xmath195 $ ] as the parameter representative of the  asymptoticity " features . the error values are shown , versus @xmath172 , in fig . [ figoptimalerrorsw ] , while the optimal settings of @xmath8 and of @xmath110 are plotted , against @xmath172 , in fig . [ figoptimalnsw ] and  [ figoptimalnmsw ] , respectively . the h - wt was introduced in ref . @xcite as a powerful and easily implementable refinement of h aimed at allowing the wt to successfully decode those divergent asymptotic series generated through the application of the steepest descent method to saddle point integrals evaluated across stokes sets , for which their single terms do not display a strictly alternating sign pattern . the scheme proposed in @xcite employed the wt only on the asymptotic series generated by the first - stage hyperasymptotic treatment of the corresponding diverging remainder . in the present sequel we reported on the possibility of combining h and the wt to higher orders in h. in particular , the full development of the 2nd - level h - wt has been detailed within the classical framework of h. the results obtained from the application of the 2nd - level h - wt , also in comparison to those obtained via the 1st - level one , on the different types of saddle - point integrals considered , showed how the increase of complexity and computational effort required by the new transformation be adequately refunded , in terms of accuracy of the estimate , particularly when the integrals are evaluated for values of their parameters which are beyond the asymptotic regime , whereas h turns out to be inapplicable and the 1st - order h - wt unavoidably lacks precision . at the same time , however , it should be noted how , for  ordinary " asymptotic evaluations , at least in the cases considered in the present work , the performances of the 1st- and the 2nd - level h - wts seem to be comparable in terms of the estimate accuracy , against a considerably difference in the computational efforts required by the two transformations . although the h - wt has been developed , here and in ref . @xcite , with reference to the evaluation of saddle - point integrals of the type in eq . ( [ sd.1 ] ) , we believe it could be useful also in dealing with problems of different nature like , for instance , the hyperasymptotic treatment of a wide class of linear and nonlinear ordinary differential equations , which has been recently considered  @xcite . the semi - analytical algorithms proposed for the numerical evaluation of the 2nd - level hyperterminants would reveal useful in this perspective .    before concluding , it is worth pointing out , as an important open problem , the need for an _ a priori _ algorithm for estimating the values of @xmath8 and @xmath47 that lead to optimal results . especially in cases where it is not convenient ( or possible ) to evaluate the original function , such an algorithm would certainly be of great help to a typical user . unfortunately , differently from h , for which the optimal settings of the hyperseries truncations are directly extracted from the singulant values@xcite , at present it does not seem possible to provide similar information for the h - wt . the difficulties in giving practical guidelines for the choice of @xmath8 and @xmath110 can also be appreciated from the results presented in sec . [ numericalexperiments ] and especially from figs . [ figoptimalnairy ] , [ figoptimalnmairy ] , [ figoptimalninstanton ] , [ figoptimalnminstanton ] , [ figoptimalnsw ] , and [ figoptimalnmsw ] , where it seems quite difficult to obtain general rules for the optimal settings of them . a possible hint , grasped from a quantitative comparison between the results obtained for the airy and the instanton functions , seems to be given by the strong connection between the h - wt retrieving performances and the saddle topology associated to the integral under consideration . what we found is that different saddle networks can share a sort of  topological equivalence " property , which is related to the set of the saddles adjacent to that under consideration and to the values of the relevant singulants . if two networks turn out to be equivalent at a certain hyperasymptotic level , this would result in the same computational effort , in terms of relative error , as far as the corresponding h - wt retrieved estimates are concerned . this , in particular , would imply that the study of the h - wt retrieving performances could be , in principle , carried out only on a restricted class of prototype functions . i would like to thank all anonymous reviewers for their constructive criticisms and suggestions . i am also grateful to turi maria spinozzi for his invaluable help during the preparation of the present work . the starting point is again the integral representation of the remainder @xmath32 given in eq . ( [ remainder.2 ] ) , where we let @xmath196 with @xmath47 being a _ new _ truncation order . substitution of eq . ( [ iistage.1 ] ) into eq . ( [ remainder.2 ] ) gives @xmath197 on using resurgence , i.e. , eq . ( [ resurgence.1 ] ) , on @xmath25 , namely @xmath198 and @xmath199 after straightforward algebra the second term in eq . ( [ iistage.2 ] ) becomes @xmath200 or , by interchanging the integral with the @xmath201-sum , @xmath202 eventually , on making the substitution @xmath203 , after some algebra the quantity in eq . ( [ iistage.3.2 ] ) can be written as @xmath204 where the 2nd - level hyperterminant @xmath52 , defined by @xmath205 has been introduced . accordingly , the complete expression of the 2nd - level hyperasymptotic expansion is @xmath206 which coincides with eq . ( [ iistage.complete ] ) . we first recast the double integral in eq . ( [ terminant2 ] ) as follows : @xmath208 which , if @xmath209 were not negative , would take the form @xmath210 ,    \end{array } \label{betanegativo.2}\ ] ] where the integrand turns out to be continuous even for real , positive values of @xmath209 , since @xmath211 if @xmath207 , eq . ( [ betanegativo.1 ] ) must be written as @xmath212 where @xmath213 denotes the principal part operator . on evaluating the integral in eq . ( [ betanegativo.4 ] ) we have @xmath214\\      \\- \mathrm{i}\pi\exp(\beta)\,\displaystyle\frac{(-\beta)^{m}}{\beta - v}.      \end{array } \label{betanegativo.5}\ ] ] finally , the direct substitution of eq . ( [ betanegativo.5 ] ) into eq . ( [ betanegativo.4 ] ) leads to the closed - form of the extra term as @xmath215 } { ( -\beta)^{n-1 } }          \,e_{n}(-\beta\gamma ) . \end{array } \label{betanegativo.6.1}\ ] ] the swallowtail function in eq . ( [ sd.1 ] ) corresponds to let @xmath216 , @xmath88 , and @xmath217 in eq . ( [ sw.1 ] ) . the ( four ) saddle points @xmath218 are solutions of the algebraic equation @xmath219 which can be exactly solved by using , for instance , cardano s formula . the evaluation of the expanding coefficients @xmath20 requires to evaluate the integral in eq . ( [ sdreview.5.0 ] ) . the first step is to expand @xmath220 around @xmath11 , i.e. , @xmath221=\\ \\ = \displaystyle\frac{(s - s_n)^2}{5\mathrm{i}}\,p(s ) , \end{array } \label{sw.3}\ ] ] where @xmath222 substitution into eq . ( [ sdreview.5.0 ] ) gives @xmath223^{r+1/2}}=\\ \\ = \displaystyle\frac{(r-1/2)!\,\,\,(5\,\mathrm{i})^{r+1/2}\,}{(2r)!}\ , \left\ { \displaystyle\frac{\mathrm{d}^{2r}}{\mathrm{d } s^{2r}}\ , \displaystyle\frac 1{[p(s)]^{r+1/2 } } \right\}_{s = s_n}. \end{array } \label{sw.5}\ ] ] equation  ( [ sw.5 ] ) can be further simplified by letting @xmath224 and then by changing @xmath225 into @xmath226 , which yields @xmath227_{t=0 } , \end{array } \label{generatingc.1.5.1}\ ] ] where @xmath228 note that eq . ( [ generatingc.1.5.1 ] ) can also be given the alternative following form : @xmath229 where the function @xmath230 is defined through @xmath231 which coincides with eq . ( [ generatingc.3 ] ) . the starting point is the definition of @xmath161 through the generating function in eq . ( [ generatingc.3.1 ] ) which , once derived with respect to @xmath232 , gives @xmath233 and , by taking eq . ( [ generatingc.3.1 ] ) into account once again , leads to @xmath234 by operating the products term by term , by rearranging the series indices , after long but straightforward algebra eq . ( [ generatingc.4 ] ) follows . 00    r. borghi ,  joint use of the weniger transformation and hyperasymptotics for accurate asymptotic evaluations of a class of saddle - point integrals , " phys . e * 78 , * 026703 ( 2008 ) . m. born and e. wolf , _ principles of optics _ ( cambridge university press , cambridge , 1999 ) .    m. v. berry ,  some quantum - to - classical asymptotics , " in les houches lecture series lii , eds . giannoni , a. voros , and j. zinn - justin , chaos and quantum physics , elsevier north - holland , amsterdam , 251 - 304 ( 1991 ) . f. j. wright ,  the stokes set of the cusp diffraction catastrophe , " j. phys . a * 13 , * 2913 - 2928 ( 1980 ) . m. v. berry and c. j. howls ,  stokes surfaces of diffraction catastrophes with codimension three , " nonlinearity * 3 , * 281 - 291 ( 1990 ) . e. j. weniger ,  nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series , " comput . phys . rep . * 10 , * 189 - 371 ( 1989 ) . e. j. weniger ,  mathematical properties of a new levin - type sequence transformation introduced by cizek , zamastil , and skala . i. algebraic theory " , j. math . * 45 , * 1209 - 1246 ( 2004 ) .    for an updated review about methods for decoding diverging series , see for instance e. caliceti , m. meyer - hermann , p. ribeca , a. surzhykov , and u. d. jentschura ,  from useful algorithms for slowly convergent series to physical predictions based on divergent perturbative expansions , " phys . * 446 , * 1 - 96 ( 2007 ) . arxiv:0707.1596v1 .    c. brezinski and m. redivo zaglia , _ extrapolation methods _ ( north - holland , amsterdam , 1991 ) . u. d. jentschura , p. j. mohr , g. soff , and e. j. weniger ,  convergence acceleration via combined nonlinear - condensation transformations , " comput . * 116 , * 28 - 54 ( 1999 ) . d. jentschura ,  resummation of nonalternating divergent perturbative expansions , " phys . d * 62 , * 076001 ( 2000 ) . s. v. aksenov , m. a. savageau , u. d. jentschura , j. becher , g. soff , and p. j. mohr ,  application of the combined nonlinear - condensation transformation to problems in statistical analysis and theoretical physics , " comput . * 150 , * 1 - 20 ( 2003 ) . m. v. berry and c. j. howls ,  hyperasymptotics , " proc . a * 430 , * 653 - 668 ( 1990 ) . m. v. berry and c. j. howls ,  hyperasymptotics for integrals with saddles , " proc . r. soc . a * 434 , * 657 - 675 ( 1991 ) .    m. v. berry and c. upstill ,  catastrophe optics : morphologies of caustics and their diffraction patterns , " prog * xviii , * 257 - 346 ( 1980 ) .    i. s. gradshteyn and i. m. ryzhik , _ table of integrals , series , and products , _ , a. jeffrey and d. zwillinger , eds . ( academic , 2000 ) . a. b. olde daalhuis ,  hyperterminants i , " j. comput . appl . math . * 76 , * 255 - 264 ( 1996 ) . a. b. olde daalhuis ,  hyperterminants ii , " j. comput . appl . math . * 89 , * 87 - 95 ( 1998 ) .    a. p. prudnikov , yu . a. brychkov , and o. i. marichev , _ integrals and series . vol . ( gordon breach publisher , new york , 1986 ) . j. p. boyd ,  the devil s invention : asymptotic , superasymptotic and hyperasymptotic series , " acta applicandae mathematicae * 56 , * 1 - 98 ( 1999 ) .    m. v. berry and c. j. howls ,  infinity interpreted , " phys . world * 6 , * 35 - 39 , ( 1993 ) . it could be proved that the above  equivalence " between the airy and the instanton functions keeps its validity also for the 2nd - level h - wt . j. n. l. connor , p. r. curtis , and d. farrelly , `` the uniform asymptotic swallowtail approximation : practical methods for oscillating integrals with four coalescing saddle points , '' j. phys . a : math * 17 , * 281 - 310 ( 1984 ) . j. f. nye , _ natural focusing and fine structures of light _ , ( institute of physics publishing , bristol , 1999 ) . j. f. nye ,  dislocation lines in the swallowtail diffraction catastrophe , " proc . a * 463 , * 343 - 355 ( 2007 ) . r. borghi , `` on the numerical evaluation of cuspoid diffraction catastrophes , '' j. opt . am . a * 25 , * 1682 - 1690 ( 2008 ) .    a. b. olde daalhuis ,  hyperasymptotics for nonlinear odes . i. a riccati equation , " proc . lond . a * 461 , * 2503 - 2520 ( 2005 ) . a. b. olde daalhuis ,  hyperasymptotics for nonlinear odes ii . the first painlev equation and a second - order riccati equation , " proc . a * 461 , * 3005 - 3021 ( 2005 ) . m. v. berry and a. shelankov ,  the aharonov - bohm wave and the cornu spiral , " j. phys . a * 32 , * l447-l455 ( 1999 ) . a. sidi ,  a new method for deriving pad approximants for some hypergeometric functions , " j. comput . appl . math . * 7 , * 37 - 40 ( 1981 ) . values , provided by the wt , of the contribution coming from the saddle @xmath112 in the evaluation of the airy function across its stokes set , when @xmath111 . [ cols="^,^",options="header " , ] fig . 1 - r. borghi
Abstract:

Answer: the use of hyperasymptotics and the weniger transformation has been proposed , in a joint fashion , for decoding the divergent asymptotic series generated by the steepest descent on a wide class of saddle - point integrals evaluated across stokes sets [ r. borghi , phys . rev . e * 78 , * 026703 ( 2008 ) ] . in the present sequel , the full development of the h - wt up to the second order in h is derived . numerical experiments , carried out on several classes of saddle - point integrals , including the swallowtail diffraction catastrophe , show the effectiveness of the 2nd - level h - wt , in particular when the integrals are evaluated beyond the asymptotic realm .