in a previous paper @xcite , a family of fixed - point algorithms for the numerical approximation of nonlinear systems of the form @xmath0 was introduced . in ( [ mm1 ] ) , @xmath1 is a nonsingular @xmath2 real matrix and @xmath3 is an homogeneous function with degree @xmath4 ( this means that @xmath5 ) . among other applications , systems of this form
are very typical in the numerical generation of traveling waves and ground states in water wave problems and nonlinear optics . for the numerical approximation to solutions of ( [ mm1 ] ) ,
the use of the classical fixed - point algorithm is not suitable .
this is due to the fact that if @xmath6 is a solution and @xmath7 stands for the iteration matrix at @xmath6 , then , since @xmath8 is homogeneous of degree @xmath9 then @xmath10 and therefore @xmath11 that is , @xmath9 is an eigenvalue of @xmath12 with magnitude above one .
this makes the iteration not convergent in general .    as an alternative and
based on the petviashvili method , @xcite , the following fixed - point algorithms were considered in @xcite : @xmath13 from @xmath14 and where @xmath15 is a @xmath16 function satisfying the following properties :    * a set of fixed points of the iteration operator @xmath17 coincides with a set of fixed points of ( [ mm1 ] ) .
this means that : ( a ) if @xmath6 is a solution of ( [ mm1 ] ) then @xmath18 ; ( b ) inversely , if the sequence @xmath19 , generated by ( [ mm2 ] ) , converges to some @xmath20 , then @xmath21 ( and , consequently , @xmath20 is a solution of ( [ mm1 ] ) ) .
* @xmath22 is homogeneous with degree @xmath23 such that @xmath24 .
the function @xmath22 is called the stabilizing factor of the method , inheriting the nomenclature of the petviashvili method .
actually , formula ( [ mm2 ] ) generalizes the petviashvili scheme , which corresponds to the choice @xmath25 the first part of the work , carried out in @xcite ( see also @xcite ) , analyzed the convergence of ( [ mm2 ] ) .
the main conclusion was that , compared to the classical fixed - point algorithm , the stabilizing factor acts like a filter of the spectrum of the matrix @xmath12 in the sense that :    * the eigenvalue @xmath26 of @xmath12 is transformed to the eigenvalue @xmath27 of the iteration matrix @xmath28 of ( [ iterop ] ) .
* the rest of the spectrum of @xmath28 is included into the spectrum of @xmath12 .
thus , the convergence of ( [ mm2 ] ) depends of the spectrum of @xmath12 different from @xmath9 . from these conclusions ,
several results of convergence can be derived , see @xcite for details .
the motivation of this paper is two - fold .
first , several numerical experiments in the literature show that the petviashvili type algorithms are in sometimes computationally slower than other alternatives . in order to continue to benefit from the easy implementation ( one of the advantages of the methods )
the algorithms should improve their performance with the inclusion of some acceleration technique .
a further motivation comes from the known mechanism of some extrapolation methods , @xcite , to transform divergent into convergent cases .
the application of this property to these petviashvili type methods may extend their use to compute traveling waves under more demanding conditions , for example in two dimensions or / and in case of highly oscillatory waves .
the literature on acceleration techniques is very rich with many different families and strategies , @xcite .
this paper will be focused on two types of procedures : the vector extrapolation methods , @xcite and the anderson mixing , @xcite .
we think that the first one is the most widely studied group ; in particular , known convergence results for some of these methods will serve us to justify several examples of transformation from divergence to convergence when generating traveling waves iteratively .
the second family accelerates the convergence by introducing the strategy of minimization of the residual in some norms at each step .
they have been revealed efficient in , for example , electronic structure computations , @xcite ( see also @xcite and references therein ) and , to our knowledge , this is the first time they are applied to the numerical generation of traveling waves .
the main purpose of this paper is then exploring by numerical means the application of these two acceleration methods to the generation of traveling waves through several problems of interest and from the petviashvili type methods ( [ mm2 ] ) and their extended versions derived in @xcite .
with the case studies presented here we have tried to cover different situations of hard computation of the waves as an attempt to give some guidelines of application . in this sense
, the paper provides several conclusions to be emphasized :    * the use of acceleration techniques is highly recommended here since it improves the performance in general and allows to extend the application of the methods to computationally harder situations , with especial emphasis on two - dimensional simulations and highly oscillatory wave generation . * by comparing the two families of acceleration techniques considered in this study , the vector extrapolation methods are in general more competitive for these problems compared to the anderson acceleration methods . among the vector extrapolation methods ,
the polynomial methods provide a better performance in general ( some exceptions can be seen in the experiments below ) . *
the main drawback of the anderson acceleration methods concerns the numerical treatment of the associated minimization problem , since most of the difficulties come from ill - conditioning .
this might be improved by including suitable preconditioning techniques ( here the methods were implemented in a standard way , @xcite ) . however , it is remarkable that when the anderson acceleration methods work , their performance is in general comparable to that of some vector extrapolation methods .
the structure of the paper is as follows .
section [ se2 ] is devoted to a description of the two families of acceleration techniques considered in this study .
this also includes some comments on the implementation and convergence results .
the application of both techniques to the methods ( [ mm2 ] ) is studied in section [ se3 ] through a plethora of numerical experiments involving the computation of different types of traveling waves : ground states , classical and generalized solitary waves as well as periodic traveling waves .
the numerical study will be focused on the two main motivations of the paper : the improvement of the efficiency and the extension of application of the methods to computationally harder problems and where the iteration is initially not convergent .
finally , section [ se4 ] completes the computational study with some illustrations of the application of the acceleration to the extended versions of the algorithms ( [ mm2 ] ) , treated in @xcite and suitable when the nonlinearity in ( [ mm1 ] ) contains several homogeneous terms with different degree .
some concluding remarks are in section [ se5 ] .
besides the local character of the convergence , in some cases and compared to other alternatives , fixed point algorithms has the additional disadvantage of a slow performance . in
what follows , several techniques of acceleration will be considered and applied to the methods ( [ mm2 ] ) , with the aim of improving their efficiency .
furthermore , as in the case of the classical algorithm , @xcite , some cases of divergence will be transformed to convergent iterations .
this section introduces two families of acceleration techniques : the vector extrapolation methods ( vem from now on ) and the anderson acceleration methods ( aam ) .
we will include a description of the schemes ( including some convergence results ) and some comments on implementation .
the first group of acceleration techniques consists of vector extrapolation methods ( vem ) . for a more detailed analysis and implementation of the methods
see @xcite and references therein . here
we will describe the general features of the procedures and their application to ( [ mm2 ] ) .
two families of vem are typically emphasized in the literature .
the first one covers the so - called polynomial methods ; they include , as the most widely cited , the minimal polynomial extrapolation ( mpe ) , the reduced rank extrapolation ( rre ) and the modified minimal polynomial extrapolation ( mmpe ) methods , @xcite .
the second family consists of the so - called @xmath29-algorithms ; typical examples are the scalar and vector @xmath29-algorithms and the topological @xmath29-algorithm , @xcite .
all the methods share of course the idea of introducing the extrapolation as a procedure to transform the original sequence @xmath30 of the involved iterative process by some strategy .
the polynomial methods are usually described in terms of the transformation ( @xmath31 ) @xmath32 where    * @xmath33 . * @xmath34 ( @xmath35 ) denotes the @xmath36 matrix of columns @xmath37 . *
@xmath38 stands for the @xmath36 matrix of some columns @xmath39 with @xmath40 as the adjoint matrix of @xmath38 ( conjugate transpose ) .    in ( [ fsec211 ] ) , @xmath41 stands for the moore - penrose generalized inverse of @xmath42 , defined as @xmath43 , @xcite .
different choices of the vectors @xmath44 lead to the most widely used polynomial methods :    1 .
minimal polynomial extrapolation ( mpe ) : @xmath45 .
2 .   reduced rank extrapolation ( rre ) : @xmath46 .
3 .   modified minimal polynomial extrapolation ( mmpe ) : @xmath47 , for arbitrary , fixed , linearly independent vectors @xmath48 .
the formulation of the vem may follow an alternative approach , @xcite .
the transformation ( [ fsec211 ] ) can be computed in the form @xmath49 where the coefficients @xmath50 are obtained from the resolution ( in some sense ) of overdetermined , inconsistent systems @xmath51 for some vectors @xmath52 .
different methods emerge by combining different choices of the norm where the residual vector @xmath53 is minimized with suitable vectors @xmath52 .
thus , for example , assuming @xmath54 , we have :    * rre is obtained by writing ( [ fsec214 ] ) in the form @xmath55 where the @xmath56 solve ( [ fsec215 ] ) with @xmath57 and the euclidean norm with equal weights is used .
* mpe is obtained by using ( [ fsec214 ] ) with @xmath58 where @xmath59 and the @xmath60 solve ( [ fsec215 ] ) with @xmath61 in the sense of minimization with the euclidean norm with equal weights .
* mmpe is obtained from ( [ fsec214 ] ) but where instead of ( [ fsec215 ] ) a system of the form @xmath62 is used . in ( [ fsec216 ] )
@xmath63 , being @xmath64
.    these formulations can be unified by a representation with determinants , @xcite .
this writes the extrapolation steps @xmath65 in the form @xmath66 with @xmath67 where the @xmath68 are scalars that depend on the extrapolation method and where the expansion of ( [ fsec218 ] ) is in the sense @xmath69 with @xmath70 the cofactor of @xmath71 in the first row .
thus in the case of the numerator in ( [ fsec217 ] ) , formula ( [ fsec219 ] ) is a vector , while in the case of the denominator in ( [ fsec217 ] ) , formula ( [ fsec219 ] ) is a scalar .
( see @xcite for the interpretation in terms of the schur complement of a matrix . )
the three previously mentioned polynomial methods correspond to the following choices of @xmath68 :    * mpe : @xmath72 .
* rre : @xmath73 . * mmpe : @xmath74 , where @xmath75 are linearly independent vectors in @xmath76 .    a second family of vem
is called the @xmath29-algorithms .
a description of them may start from the scalar @xmath29-algorithm of wynn , @xcite .
this scalar extrapolation method can be derived from the representation ( [ fsec217 ] ) , ( [ fsec218 ] ) ( in the scalar case ) with @xmath77 ( which are scalars in the scalar case ) .
the corresponding ratio of determinants @xmath78 is called the classical @xmath79- ( or shanks schmidt @xmath80 ) transform , @xcite .
this ratio can be evaluated recursively for increasing @xmath81 and @xmath82 without the computation of determinants or schur complements .
the corresponding formulation is @xmath83 where @xmath84 and works along diagonals on @xmath85 constant .
thus , from ( [ fsec2111 ] ) , formulas ( [ fsec2112 ] ) , ( [ fsec2113 ] ) compute each entry of a triangular array in terms of the previous entries .
the extension of the scalar @xmath29-algorithm to the vectorial case was carried out by brezinski , @xcite , and wynn , @xcite , by using different definitions of inverse  of a vector , see @xcite .
wynn suggests to consider the transpose of the moore - penrose generalized inverse of a vector , @xmath86 leading to the vector @xmath29-algorithm ( vea ) , whose formulas are of the form ( [ fsec2112 ] ) , ( [ fsec2113 ] ) where the scalars are substituted by vectors and ( [ fsec2113 ] ) makes use of ( [ fsec2114 ] ) .
this implies that the @xmath79-transform ( [ fsec2111 ] ) is understood in the above described vectorial sense .
this is called the generalized shanks schmidt ( gss ) transform , @xcite . on the other hand , brezinski defines the inverse of pair of vectors @xmath87 such that @xmath88 as the pair of vectors @xmath89 where @xmath90 thus , @xmath91 is called the inverse of @xmath92 with respect to @xmath93 and viceversa .
this definition leads to the so - called topological @xmath29-algorithm ( tea ) , when an arbitrary vector @xmath20 is fixed and the inverses of @xmath94 and @xmath95 are considered with respect to @xmath20 , that is @xmath96 the recursive formulas are @xmath97 brezinski proved , @xcite , the connection with the gss - transform , showing that @xmath98 for an efficient implementation of ( [ fsec2115 ] ) see @xcite .
thus ( tea ) corresponds to take @xmath99 in ( [ fsec217 ] ) , ( [ fsec218 ] ) .
the mechanism of working of the vem can be described as follows , see @xcite for details .
one starts from assuming an asymptotic expression for the sequence @xmath100 of the form @xmath101 where @xmath102 and @xmath103 , ordered such that @xmath104 if @xmath105 with only a finite number of @xmath106 having the same modulus .
the expansion ( [ fsec212 ] ) can be generalized by considering , instead of constant vectors @xmath107 , polynomials @xmath108 in @xmath82 with vector coefficients of the form @xmath109 with @xmath110 linearly independent in @xmath76 and where if @xmath111 then @xmath112 , @xcite . for simplicity
, the description below will make use of ( [ fsec212 ] ) .    the asymptotic expansion ( [ fsec212 ] ) is considered in a general vector space ( finite or infinite dimensional ) where the iteration is defined .
it is understood in the sense that any truncation differs from @xmath100 in less than some power of the next @xmath113 .
this means that for any positive integer @xmath8 there are @xmath114 and a positive integer @xmath115 that only depend on @xmath8 such that for every @xmath116 @xmath117 allows to identify @xmath118 as limit or anti - limit of the sequence @xmath100 , according to the size of @xmath119 .
( that is , if @xmath120 then @xmath121 exists and equals @xmath118 .
if @xmath122 then @xmath121 does not exist and @xmath118 is called the anti - limit of the sequence @xmath100 . ) under these conditions , several results of convergence for mpe , rre , mmpe and tea are obtained in the literature , @xcite and references therein . for these methods , one can find an extrapolation step @xmath123 such that @xmath124 ) may explain the convergent behaviour of the extrapolation in some cases . if the @xmath113 s are identified as the eigenvalues of the linearization operator of the iteration at the limit ( or anti - limit ) @xmath118 , then the extrapolation has the effect of translating the behaviour of the iteration to an eigenvalue @xmath125 that may be into the unit disk , even if the previous ones are out of it .
hence , @xmath118 may be anti - limit for the original iteration and the extrapolation converges to it .
these results are extended to the defective linear case with more general polynomials ( [ fsec212b ] ) in @xcite .
the integer @xmath123 is related to the concept of minimal polynomial @xmath126 of a matrix @xmath42 with respect to a vector @xmath92 , @xcite ; this is the unique polynomial of least degree such that @xmath127 thus in the case of linear iteration with matrix @xmath42 , @xmath123 is taken to be the degree of the minimal polynomial of @xmath42 with respect to the first iteration @xmath128 . in the case of a nonlinear system written in fixed point form
@xmath129 then @xmath123 is theoretically defined as the degree of the minimal polynomial of @xmath130 with respect to @xmath128 , where @xmath6 is a solution of ( [ fsec2110 ] ) .
in contrast with the linear case , there is no way to determine @xmath123 in advance for the nonlinear case .
this forces to consider several strategies for the choice and the corresponding implementation , see the discussion in @xcite and the comments here below .
we also mention that in the linear case , the extrapolation methods mpe and rre are mathematically equivalent to the method of arnoldi , @xcite , and the gmres , @xcite , respectively , see @xcite , while the mmpe is mathematically equivalent to the hessenberg method , @xcite and tea to the method of lanczos , @xcite , see @xcite .
efficient and stable implementation of the rre , mpe and mmpe methods by using qr and lu factorizations can be seen in @xcite .
for the case of tea , see @xcite and @xcite for vea .
the implementation is usually carried out in a cycling mode .
a cycle of the iteration is performed by the following steps : consider a method ( [ mm2 ] ) with a stabilizing factor @xmath22 satisfying ( p1 ) , ( p2 ) .
given @xmath14 and a width of extrapolation @xmath131 , for @xmath132 , the advance @xmath133 is :    * set @xmath134 and compute @xmath135 steps of the fixed - point algorithm : @xmath136 * compute the extrapolation steps ( [ fsec217 ] ) with any of the methods described above and @xmath137 .
* set @xmath138 and go to step ( a ) .    the cycle ( a)-(b)-(c ) is repeated until the error ( residual or between two consecutive iterations ) is below a prefixed tolerance , a maximum number of iterations is attained or the discrepancy between the stabilizing factor at the iterations and one is below a prefixed tolerance .    the width of extrapolation @xmath135 depends on the choice of the technique : @xmath139 for mpe , rre or mmpe and @xmath140 for vea or tea , @xcite .
since @xmath123 is generally unknown , some strategy for the implementation must be adopted . in practice , as discussed in @xcite , the methods are implemented with some ( small ) values of @xmath135 and take that with the best performance
. it may be also different for each cycle , although quadratic convergence is not expected if @xmath123 is too small .
this choice of @xmath123 will be computationally studied in some examples in section [ se3 ] .
the hypotheses for the expansion ( [ fsec212 ] ) include the conditions @xmath141 .
in many problems for traveling wave generation , @xmath142 appears as eigenvalue of the iterative technique ( of fixed point type ) although under especial circumstances that allow to extend the convergence results in some sense .
this especial situation is related to the presence of symmetries in the equations for traveling waves . in order to extend the results of convergence to this case
, one has to consider the orbits by the symmetry group of the equations and interpret the convergence in the orbital sense , that is a convergence for the orbits . the description of this orbital convergence can be seen in @xcite .
finally , local quadratic convergence is proved in @xcite ( see also @xcite ) for the four methods and vea for a general nonlinear system ( [ fsec2110 ] ) and under the following hypotheses on @xmath143 :    * the jacobian matrix @xmath144 does not have @xmath142 as eigenvalue .
* @xmath81 is the degree of the minimal polynomial of @xmath144 with respect to @xmath145 . *
the algorithm is implemented in the cycling mode where @xmath123 is chosen in the @xmath146-th cycle as the degree of the minimal polynomial of @xmath144 with respect to @xmath147 .
this result can be extended to the case where @xmath148 admits a @xmath149-parameter ( @xmath150 ) group of symmetries and , consequently , @xmath142 is an eigenvalue of @xmath144 , by using the reduced system for the orbits of the group and in the orbital sense , @xcite .
this would lead to quadratic orbital convergence but linear local convergence .
a second family of acceleration techniques considered here is the so - called anderson family or anderson mixing , @xcite .
it is widely used in electronic structure computations and only recently it has been analyzed in a more general context , @xcite .
( to our knowledge , this is the first time that aam are applied to accelerate traveling wave computations . )
the main goal of the approach consists of combining the iteration with a minimization problem for the residual at each step . for linear problems ,
this technique is essentially equivalent to the gmres method , @xcite .
the stages for an iteration step are as follows : given @xmath151 , set @xmath152 . for @xmath153
* set @xmath154 * set @xmath155 where @xmath156 * determine @xmath157 that solves @xmath158 * set @xmath159    the resolution of the optimization problem ( [ minimize ] ) is the source of the additional computational work of the acceleration .
one way to reduce this extra effort is the so - called multisecant updating @xcite .
( this also clarifies the connection with quasi - newton methods . ) this technique consists of writing the problem in an equivalent form @xmath160 but with a more direct resolution , and determining the acceleration from it .
the general step becomes :    * set @xmath154 * determine @xmath161 by solving ( [ minip ] ) .
* set @xmath162 * set @xmath163    as mentioned in @xcite , if @xmath164 is full - rank , the solution of the minimization problem can be written as @xmath165 and the anderson acceleration has the alternative form @xmath166 ( note that @xmath167 can be viewed as an approximate inverse of the jacobian of @xmath168 ) . the formulation ( [ aaalternative ] ) motivates the generalization of the anderson mixing , @xcite .
this is performed replacing @xmath169 in ( [ aaalternative ] ) by some @xmath170 satisfying @xmath171 and ( [ aaalternative ] ) becomes @xmath172 the resulting methods are collected in the so - called anderson s family , @xcite .
two particular members are emphasized : the type - i method ( denoted by aa - i from now on ) , which corresponds to @xmath173 in ( [ aaalternative2 ] ) and type - ii method ( or aa - ii ) , which is the original anderson mixing ( [ aaalternative ] ) .    to our knowledge
, some convergence results can be seen in @xcite . in @xcite
the authors identify some anderson methods for linear problems and in some sense with the gmres method and the arnoldi ( fom ) method ; some convergence results can be derived from this identification . in @xcite the equivalence with gmres for linear problems is completely characterized .
finally , @xcite gives some proofs of convergence of the anderson acceleration when applied to contractive mappings : @xmath23-linear convergence of the residual for linear problems under certain conditions when @xmath174 and local @xmath175-linear convergence in the nonlinear case .
( these types of convergence are defined in the paper . )    on the other hand , as observed in @xcite , the implementation of aam should be carried out by attending to three main points : a convenient formulation of the minimization problem ( [ minimize ] ) , a numerical method for its efficient resolution and , finally , the parameter @xmath176 , which plays a similar role to that of the extrapolation width @xmath135 in the vem . in our computations below , we have followed the treatment described in @xcite .
this is based on the use of the unconstrained form ( [ minip ] ) and its numerical resolution with @xmath177 decomposition . for other alternatives in both problems , see the discussion in @xcite and the references cited there .
( according to our results below , the use of alternative preconditioning techniques might be recommendable in some cases . )
as far as the choice of @xmath176 is concerned , a similar strategy to that of @xmath135 will be used , since our experiments , @xcite , suggest that @xmath176 ( as @xmath135 ) strongly depends on the problem under study and large values are not recommended .
finally the codes are implemented by retaining the definition of @xmath178 since other alternatives , @xcite , did not improve the results in a relevant way .
presented here is a comparative study on the use of vem and aam as acceleration techniques from the petviashvili type methods ( [ mm2 ] ) in traveling wave computations .
the comparison is organized according to two main points : the type of traveling wave to be generated and the elements of each family of methods to be used for the generation .
the first point includes the following case studies :    1 .
classical solitary waves , generalized solitary waves and periodic traveling waves of the four - parameter boussinesq system , @xcite .
localized ground state solutions of nls type equations , @xcite .
highly oscillatory solitary waves of the one- and two - dimensional benjamin equation , @xcite .    this plethora of waves attempts to discuss and overcome different computational difficulties and with the aim of establishing as more general conclusions as possible . as for the second point
, each family of techniques has been represented by the following methods :    * mpe and rre standing for polynomial extrapolation methods . * vea and tea standing for @xmath29 algorithms . * aa - i and aa - ii standing for the aam .    for simplicity
, the acceleration will be applied to the petviashvili method ( [ mm2 ] ) , ( [ mm3c ] ) with @xmath179 . due to the similar behaviour of the methods of the family ( [ mm2 ] ) , illustrated in @xcite , the conclusions from the corresponding results
can reasonable serve when the petviashvili method is substituted by any of ( [ mm2 ] ) .    in all the cases
considered , the traveling wave profiles appeared as solutions of initial value problems of ode s .
their discretization to generate approximations to the profiles , was carried out in a common and standard way .
the corresponding initial periodic boundary value problem ( on a sufficiently long interval ) was discretized by using fourier collocation techniques , @xcite ; the discretization leads to a nonlinear system of algebraic equations for the approximate values of the profile at the collocation points or for the discrete fourier coefficients of the approximation .
this system is iteratively solved with the classical petviashvili method ( [ mm2 ] ) , ( [ mm3c ] ) along with the selected acceleration technique .
this will be described in each equation considered below .
several stopping criteria for the iterations are implemented :    * a maximum number of iterations . *
the iteration stops when one of the following quantities are below a prefixed , small tolerance @xmath180 : * * the difference in euclidean norm between two consecutive iterations @xmath181 * * the residual error ( also in euclidean norm ) @xmath182 * * the discrepancy between the stabilizing factor and ( in case of convergence ) its limit one @xmath183    the numerical experiments that form the comparative study are of different type :    * for several values of @xmath123 ( in the case of vem ) and @xmath176 ( in the case of aam ) we have computed the number of iterations required by each method to achieve a residual error below the prefixed tolerance .
this allows to compare some performance of the methods between the two families , among different techniques within a same family and indeed with the petviashvili method without acceleration . *
some eigenvalues of the iteration matrices for the classical fixed point algorithm and the petviashvili method have been computed ( with the corresponding standard matlab function ) in order to explain the behaviour of the second one , @xcite and how the acceleration eventually changes it . *
the form of the approximate profiles and some experiments to check their accuracy are also displayed .      in this first example
we study the numerical generation of traveling wave solutions of the four - parameter family of boussinesq system @xmath184 where @xmath185 and the four parameters @xmath186 satisfy @xmath187 with some constant @xmath188 $ ] , @xcite . system ( [ fsec311 ] ) appears as one of the alternatives to model the bidirectional propagation of the irrotational free surface flow of an incompressible , inviscid fluid in a uniform horizontal channel under the effects of gravity when the surface tension and cross - channel variations of the fluid are assumed to be negligible .
if @xmath189 denotes the undisturbed water depth , then @xmath190 stands for the deviation of the free surface from @xmath189 at the point @xmath191 and time @xmath192 , while @xmath193 is the horizontal velocity of the fluid at @xmath191 and at the height @xmath194 ( where @xmath195 corresponds to the channel bottom ) at time @xmath192 . for the derivation of ( [ fsec311 ] ) from the two - dimensional euler equations and the mathematical theory see @xcite . for the modification of ( [ fsec312 ] ) to include the influence of surface tension see @xcite .
the boussinesq system ( [ fsec311 ] ) admits different types of traveling wave solutions .
first , being an approximation to the corresponding two - dimensional euler equations in the theory of surface waves , it is expected to have classical solitary wave solutions .
they are solutions of the initial value problem of ( [ fsec311 ] ) of smooth traveling wave form @xmath196 with some speed @xmath197 and decaying to zero as @xmath198 .
substitution into ( [ fsec311 ] ) and after one integration the profiles @xmath199 must satisfy the ode system @xmath200 the problems of existence , asymptotic decay and stability of solutions of ( [ fsec313 ] ) have been analyzed in many references and for particular values of @xmath186 , see @xcite and references therein .
furthermore , in the same reference , linearly well - posed systems ( [ fsec313 ] ) may be studied as a first order ode system and , based on normal form theory , a discussion on the values of the parameters leading to boussinesq systems admitting solitary wave solutions with speed @xmath201 is established . according to it , two classes of systems can be distinguished .
the first one admits classical ( in the sense above defined ) solitary wave solutions .
this group contains the bona - smith system ( @xmath202 ) , @xcite , or the bbm - bbm system ( @xmath203 ) , @xcite .
( the classical boussinesq system , which corresponds to @xmath204 , is also in this group , although it is out of the general discussion of @xcite and has been studied separately , @xcite . )
a second class of boussinesq systems admits generalized solitary wave solutions , that is traveling wave profiles which are not homoclinic to zero at infinity but to small amplitude periodic waves , @xcite .
the kdv - kdv system ( @xmath205 ) , @xcite , is an example of this second group .
finally , the existence of periodic traveling wave solutions of ( [ fsec311 ] ) is studied in @xcite by applying topological degree theory for positive operator to the corresponding periodic initial value problem posed on an interval @xmath206 and some cnoidal wave solutions of the bbm - bbm system are computed .
a smooth periodic traveling wave solution @xmath207 with some speed @xmath197 must satisfy the ode system @xmath208 for some real constants @xmath209 .
these are related to the period parameter @xmath210 .
the resolution of ( [ fsec313b ] ) involves a modified system for which these constants of integration are set to zero , @xcite .
this is briefly described as follows .
one first searches for constant solutions @xmath211 of ( [ fsec313b ] ) .
this leads to the system @xmath212 which can be solved as a cubic equation for @xmath213 : @xmath214 once @xmath215 is obtained ( there may be more than one solution indeed ) the differences @xmath216 must satisfy @xmath217 this strategy will be considered in the numerical generation of the profiles in ( [ fsec313b ] ) : the system ( [ fsec313d ] ) wil be discretized to compute approximations to the variables @xmath218 and to the variables @xmath219 from them .    in order to generate numerically classical and generalized solitary wave solutions of ( [ fsec311 ] ) the corresponding periodic value problem of ( [ fsec313 ] ) on a long enough interval
@xmath206 is discretized with a fourier collocation method leading to a discrete system of the form @xmath220 where @xmath221 are approximations @xmath222 to the values of a solution of ( [ fsec313b ] ) at the grid points @xmath223 , @xmath224 is the pseudospectral differentiation matrix , @xcite , @xmath225 is the @xmath2 identity matrix and the nonlinear term @xmath8 , which is homogeneous of degree @xmath226 , involves hadamard products . in the case of periodic traveling waves and as was mentioned above , system ( [ fsec314 ] ) will be substituted in the implementation by @xmath227 for the approximations @xmath228 to the @xmath229 variables at the grid points and where @xmath230 are previously known from the resolution of ( [ fsec313c ] ) , ( [ fsec313c1 ] ) .
then @xmath231 .    the methods ( [ mm2 ] ) along with the corresponding acceleration technique
are then applied to the discrete systems ( [ fsec314 ] ) and ( [ fsec314b ] ) .
the implementation is performed in the fourier space ; for example ( [ fsec314 ] ) becomes @xmath232 thus the @xmath233 system ( [ fsec314 ] ) is divided into @xmath234 blocks of @xmath235 systems for the corresponding @xmath9-th discrete fourier coefficients @xmath236 ( for simplicity , we assume that @xmath237 for some @xmath238 . ) alternatively , ( [ fsec314 ] ) can be written in the form @xmath239 where @xmath240 denotes periodic convolution and if @xmath241 , the vectors @xmath242 have discrete fourier coefficients @xmath243 in order to explain the behaviour of the iteration , the size of the eigenvalues of the iteration matrix will be relevant in the numerical study . in this case , the corresponding iteration matrix of the classical fixed point iteration at a solution @xmath244 has the form    @xmath245    ( where @xmath246 stands for the diagonal matrix with diagonal entries given by the components of @xmath247 ) . some information on the spectrum of @xmath12 is known .
we already have the eigenvalue @xmath248 , corresponding to the degree of homogeneity of the nonlinear part , with @xmath244 as an eigenvector .
also , the application of @xmath224 to ( [ fsec314 ] ) leads to @xmath249 which means that @xmath142 is an eigenvalue of @xmath12 and @xmath250 is an associated eigenvector .
this corresponds to the translational  invariance of ( [ fsec311 ] ) .
three particular systems of ( [ fsec311 ] ) will be taken to illustrate the numerical generation of traveling waves .
the first one is the classical boussinesq system ( @xmath204 ) , @xcite @xmath251 which is known to have classical solitary wave solutions , @xcite .
the second one is the so - called kdv - kdv system ( @xmath252 ) @xmath253 that admits generalized solitary wave solutions , @xcite . finally , in order to illustrate the numerical generation of periodic traveling waves , the bbm - bbm system ( @xmath203 )
, @xmath254 will be taken , @xcite .      in the case of system ( [ fsec317 ] ) a first experiment of comparison of the acceleration techniques has been made by taking @xmath255 and a hyperbolic secant profile as initial iteration with @xmath256 and @xmath257 collocation points . the petviashvili method ( [ mm2 ] ) , ( [ mm3c ] ) with @xmath258 was first run , generating approximate @xmath259 and @xmath118 profiles as shown in figures [ figuresw0](a ) and ( b ) while figures [ figuresw0](c ) and ( d ) stand for the corresponding phase portraits of the approximate profiles in ( a ) and ( b ) .
( they show the classical character of the solitary waves , represented as homoclinic to zero orbits with exponential decay , @xcite . )
the accuracy of the iteration is checked in figure [ figuresw1 ] .
figure [ figuresw1](a ) illustrates the convergence of the sequence @xmath260 of stabilizing factors , computed with the corresponding formula ( [ mm3c ] ) and the optimal value @xmath258 . the discrepancy ( [ fsec33 ] )
is below the tolerance @xmath261 in @xmath262 iterations , while the first residual error below @xmath180 is @xmath263 at @xmath264 .
( this also happens in the rest of the experiments : when the procedure is convergent , the error @xmath265 achieves the tolerance before the residual error ; therefore , the control on this last one is a harder test and will be adopted as the main one to stop the iteration . )    convergence is also confirmed by table [ tav1b ] .
this shows the six largest magnitude eigenvalues of the iteration matrix ( [ fsec316 ] ) ( first column ) and of the iteration matrix ( at the same iterate ) of the petviashvili procedure ( second column ) both at the last computed iterate @xmath266 .
the first column reveals the dominant eigenvalues @xmath267 , both simple , while the rest is below one .
the filtering effect of the petviashvili method , @xcite , is observed in the second column ; the dominant eigenvalue is filtered to zero ( recall that @xmath258 ) and the rest is preserved . since @xmath268 corresponds to the translational symmetry of ( [ fsec311 ] ) , this guarantees the local convergence of the method ( also in the orbital sense mentioned above ) .
.classical solitary wave generation of ( [ fsec317 ] ) .
six largest magnitude eigenvalues of the approximated iteration matrix @xmath269 ( first column ) and of the iteration matrix @xmath270 , generated by the petviashvili method ( [ mm2 ] ) , ( [ mm3c ] ) with @xmath258 , both evaluated at the last computed iterate @xmath266 . [ cols="^,^",options="header " , ]
in this paper we have studied numerically the use of acceleration techniques applied to fixed point algorithms of petviashvili type to generate numerically traveling waves in nonlinear dispersive wave equations .
the comparison has been established between vector extrapolation methods and anderson acceleration methods for different types of traveling waves . from the plethora of numerical experiments ,
our main conclusions are :    * the use of acceleration techniques improves the performance of the petviashvili type methods in all the cases .
this improvement is observed in two main points : first , when the petviashvili type method is convergent , the acceleration reduces the number of iterations in a relevant way .
( in some cases , this is really important : in some one - dimensional problems , the reduction is at least of @xmath271 and attains up to @xmath272 . ) on the other hand , the mechanism of acceleration , especially in the case of vem , allows to transform initially divergent sequences into convergent processes .
this is particularly relevant in traveling waves with high oscillations .
furthermore , acceleration has been shown to be more efficient than other alternatives for some cases , like numerical continuation . * in general , vem provide better results and among them , polynomial methods such as mpe and rre are more efficient than @xmath29-algorithms like vea and tea , although in some convergent cases the acceleration in terms of the number of iterations is very similar among all the methods while in computational time the @xmath29-algorithms work better . * the aam are competitive in some cases but they are mostly affected by ill - conditioning and a more computational effort due to their longer implementation . the best results of these methods are obtained when generating numerically periodic traveling waves in some nonlinear dispersive systems and ground state profiles in nls type equations , while their performance is poor when computing highly oscillatory traveling waves .
their application to the e - petviashvili type methods suggests an uncertain behaviour with respect to relevant quantities of the problem through the iteration .
the main question to that this comparative study has not been able to answer is in the authors opinion finding a deeper understanding of the way how acceleration techniques ( especially vem ) work on these problems . in particular , we miss some conclusions about the width of extrapolation ( which is related to the extrapolation step for convergence ) to be used a priori ( if it is possible to do that ) . we have observed that this looks to be strongly dependent on the problem under study .
this might be though a good starting point for a future research .
this research has been supported by project mtm2014 - 54710-p .
ablowitzm l. a. abramyan , y. a. stepanyants , the structure of two - dimensional solitons in media with anomalous small dispersion , siv .
jetp , 61(5)(1985 ) 963 - 966 .
albert , j.l .
bona , j.m .
restrepo , solitary - wave solutions of the benjamin equation , siam j. appl .
, 59 ( 1999 ) 21392161 .
j. alvarez , a. duran , petviashvili type methods for traveling wave computations : i. analysis of convergence , j. comp .
, 266 ( 2014 ) 39 - 51 .
j. alvarez , a. duran , an extended petviashvili method for the numerical generation of traveling and localized waves , commun .
nonlinear sci .
simulat . , 19(2014 ) 2272 - 2283 .
j. alvarez , a. duran , corrigendum to `` petviashvili type methods for traveling wave computations : i. analysis of convergence '' [ j. comput .
( 2014 ) 3951 ] , j. comp .
math . , 277(2015 ) 215 - 216 .
d. g. anderson , iterative procedures for nonlinear integral equations , j. assoc .
, 12 ( 1965 ) 547 - 560 .
t. b. benjamin , internal waves of permanent form in fluids of great depth , j. fluid mech .
29 ( 1967 ) 559 - 592 .
t. b. benjamin . a new kind of solitary wave , j. fluid mech . , 245 ( 1992 ) 401 - 411
. t. b. benjamin , solitary and periodic waves of a new kind , phil .
a , 354 ( 1996 ) 1775 - 1806 .
j. l. bona and m. chen , a boussinesq system for two - way propagation of nonlinear dispersive waves , physica d , 116 ( 1998 ) , 191224 .
j. l. bona , m. chen , j .- c .
saut , boussinesq equations and other systems for small - amplitude long waves in nonlinear dispersive media : i. derivation and linear theory , j. nonlin .
12 ( 2002 ) , 283 - 318 . j. l. bona , m. chen , j .- c .
saut , boussinesq equations and other systems for small - amplitude long waves in nonlinear dispersive media : ii . the nonlinear theory , nonlinearity 17 ( 2004 ) , 925 - 952 . j. l. bona , v. a. dougalis , d. e. mitsotakis , numerical solution of kdv - kdv systems of boussinesq equations : i. the numerical scheme and generalized solitary waves , math .
simul . , 74(2007 ) 214 - 228 .
j. l. bona , v. a. dougalis , d. e. mitsotakis , numerical solution of kdv - kdv systems of boussinesq equations : ii .
generation and evolution of radiating solitary waves , nonlinearity 21 ( 2008 ) 2825 - 2848 .
j. l. bona and r. smith , a model for the two - way propagation of water waves in a channel , math .
79(1976 ) , 167182
. j. v. boussinesq , thorie des ondes et des remous qui se propagent le long dun canal rectangulaire horizontal , en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond , j. math .
pures appl .
17 ( 1872 ) 55 - 108 . , chebyshev and fourier spectral methods , 2nd ed .
dover publications , new york , 2000 .
c. brezinski , gnralization de la transformation de shanks , de la table de pad et de lepsilon - algorithm , calcolo , 12 ( 1975 ) 317 - 360 . c. brezinski , a general extrapolation algorithm , numer .
, 35(1980 ) 175 - 187 .
c. brezinski , convergence acceleration during the 20th century , j. comput .
, 122 ( 2000 ) 1 - 21 . c. brezinski , a. c. rieu , the solution of systems of equations usiung the vector @xmath29-algorithm and an application to boundary value problems , math .
comp . , 28(1974 )
731 - 741 .
c. brezinski , m. redivo zaglia , extrapolation methods , theory and practice , north - holland , amsterdam , 1991 . c. brezinski , m. redivo zaglia , a schur complement approach to a general extrapolation algorithm , linear algebra and appl . , 368 ( 2003 ) 279 - 301 .
, spectral methods in fluid dynamics .
springer - verlag , new york - heidelberg - berlin , 1988 .
s. cabay , l. w. jackson , a polynomial extrapolation method for finding limits and antilimits for vector sequences , siam j. numer .
, 13(1976 ) 734 - 752
. m. chen , m. chen , n. v. nguyen , cnoidal wave solutions to boussinesq systems , nonlinearity , 20 ( 2007 ) , 1443 - 1461 . m. chen , n. v. nguyen , s .-
sun , solitary - wave solutions to boussinesq systems with large surface tension , disc .
systems , ser .
s , 2(1)(2009 ) 37 - 53 .
m. chen , n. v. nguyen , s .-
sun , existence of traveling - wave solutions to boussinesq systems , diff .
eq . , 24(2011 ) 895 - 908 .
p. daripa , r. k. dash , a class of model equations for bi - directional propagation of capillary - gravity waves .
sci . , 41(2003 ) 201 - 218 .
j. demmel , applied numerical linear algebra , siam philadelphia , 1997 .
v. a. dougalis , a. duran , d. e. mitsotakis , numerical approximation of solitary waves of the benjamin equation , accepted in math .
simul . , 2012 . v. a. dougalis , a. duran , d. e. mitsotakis , numerical solution of the benjamin equation , wave motion 52(2015 ) 194 - 215 . v. a. dougalis , d. e. mitsotakis , theory and numerical analysis of boussinesq systems : a review , in : effective computational methods in wave propagation , n. a. kamparis , v. a. dougalis , j. a. ekaterinaris ( eds . ) crc press 2008 , 63 - 110 .
h. fang , y. saad , two classes of multisecant methods for nonlinear acceleration , numer .
linear algebra appl . , 16 ( 2009 ) 197 - 221 .
r. p. eddy , extrapolation to the limit of a vector sequence , in : p.c.c .
wang ( ed . ) , information linkage between applied mathematics and industry , academic press , new york , 1979 , 387 - 396 .
e. eyert , a comparative study on methods for convergece acceleration of iterative vector sequences , j. comput .
, 124(1996 ) 271 - 285 .
e. gekeler , on the solution of systems of equations by the epsilon algorithm of wynn , math .
, 26(1972 ) 427 - 437
. g. h. golub , ch .
f. van loan , matrix computations , j. h. u. press , baltimore , 1996 .
k. jbilou , h. sadok , some results about vector extrapolation methods and related fixed point iterations , j. comput .
appl . math . , 36(1991 ) 385 - 398 .
k. jbilou , h. sadok , lu - implementation of the modified minimal polynomial extrapolation method , i m a j. numer .
, 19(1999 ) 549 - 561 .
k. jbilou , h. sadok , vector extrapolation methods
. applications and numerical comparisons , j. comput .
math . , 122 ( 2000 ) 149 - 165 .
b. b. kadomtsev , v. i. petviashvili , on the stability of solitary waves in weakly dispersive media , sov .
phys . dokl . 15
( 1970 ) 539 - 541 .
b. kim , three - dimensional solitary waves in dispersive wave systems , doctoral dissertation , department of mathematics , mit , 2006 .
b. kim , t. r. akylas , on gravity - capillary lumps , j. fluid mech .
540 ( 2005 ) 337 - 351 .
b. kim , t. r. akylas , on gravity - capillary lumps .
two - dimensional benjamin equation , j. fluid mech . 557 ( 2006 ) 3 + 237 - 256 .
, a generalized petviashvili method for scalar and vector hamiltonian equations with arbitrary form of nonlinearity , j.  comput .
( 2007 ) 1668 - 1692 .
lakoba , j. yang , a mode elimination technique to improve convergence of iteration methods for finding solitary waves , j. comp .
phys . 226 ( 2007 ) 1693 - 1709 . c. lanczos , solutions of systems of linear equations by minimized iterations , j. res .
, 49(1952 ) 33 - 53
. h. le ferrand , the quadratic convergence of the topological @xmath29-algorithm for systems of nonlinear equations , numer .
algorithms , 3 ( 1992 ) 273 - 284 .
e. lombardi , oscillatory integrals and phenomena beyond all algebraic orders , with applications to homoclinic orbits in reversible systems , lecture notes in mathematics , springer - verlag , berlin , 2000 .
m. mesina , convergence acceleration for the iterative solution of the equations @xmath273 , comput .
methods appl .
engrg . , 10(1977 ) 165 - 173 .
c. d. meyer , matrix analysis and applied linear algebra , siam philadelphia , 2000 .
p. ni , anderson acceleration of fixed - point iteration with applications to electronic structure computations , ph .
d. thesis , worcester polytechnic institute , worcester , ma , 2009 .
r. l. pego , m. i. weinstein , convective linear stability of solitary waves for boussinesq equations , l stud .
appl . math .
99(1997 ) 311 - 375 .
, convergence of petviashvili s iteration method for numerical approximation of stationary solutions of nonlinear wave equations , siam j.  numer .
anal . 42 ( 2004 )
1110 - 1127 . ,
soviet j.  plasma phys . 2 ( 1976 ) 257 - 258 .
f. potra , h. engler , a characterization of the behaviour of the anderson acceleration on linear problems , linear algebra and its applications , 438(3)(2013 ) 1002 - 1011 .
p. pulay , improved scf convergence , j. comput .
, 3 ( 1982 ) 556 - 560 .
y. saad , krylov subspace methods for solving large unsymmetric linear systems , math .
, 37(1981 ) 105 - 126 . y. saad , m. h. schultz , gmres : a generalized minimal residual algorithm for solving nonsymmetric linear systems , siam j. sci .
comput . , 7(1986 ) 856 - 869 .
a. sidi , convergence and stability of minimal polynomial and reduced rank extrapolation algorothms , siam j. numer .
, 23 ( 1986 ) 197 - 209 .
a. sidi , extrapolation vs. projection methods for linear systems of equations , j. comp .
math . , 22(1)(1988 ) 71 - 88 .
a. sidi , efficient implementation of minimal polynomial and reduced rank extrapolation methods , j. comput .
math . , 36(1991 ) 305 - 337 .
a. sidi , practical extrapolation methods , theory and applications , cambridge university press , new york , 2003 .
a. sidi , j. bridger , convergence and stability analyses for some vector extrapolation methods in the presence of defective iteration matrices , j. comp .
math . , 22(1988 ) 35 - 61 .
a. sidi , w. f. ford , d. a. smith , acceleration of convergence of vector sequences , siam j. numer .
, 23 ( 1986 ) 178 - 196 .
j. r. schmidt , on the numerical solution of linear simultaneous equations by an iterative method , phil .
, ser . 7 , 32(1941 ) 369 - 383 .
d. shanks , non - linear transformations of divergent and slowly convergent sequences , j. math . and phys .
, 34(1955 ) 1 - 42 .
s. skelboe , computation of the periodic steady - state response of nonlinear networks by extrapolation methods , ieee trans .
circuits and systems , 27(1980 ) 161 - 175 .
d. a. smith , w. f. ford , a. sidi , extrapolation methods for vector sequences , siam rev .
, 29 ( 1987 ) 199 - 233 .
r. c. e. tan , implementation of the topologicasl @xmath29.algorithm , siam j. sci .
9 ( 1988 ) 839 - 848 .
a. toth , c. t. kelley , convergence analysis for anderson acceleration , siam j. numer .
, 53(2015 ) 805 - 819
. j. van iseghem , convergence of vectorial sequences .
applications , numer . math . ,
68(1994 ) 549 - 562 . v. v. voronovich , v. i. shira , y. a. stepanyants , two - dimensional models for nonlinear voriticy waves in shear flows , stud . appl .
math . , 100(1998 ) 1 - 32 .
h. f. walker , p. ni , anderson acceleration for fixed - point iterations , siam j. numer .
, 49 ( 2011 ) 1715 - 1735 .
p. wynn , on a device for computing the @xmath274 transformation , mathematical tables and other aids to computation , 10(1956 ) 91 - 96 .
p. wynn , acceleration techniques in numerical analysis , with particular reference to problems in one independent variable , proceedings of the ifip congress 1962 , north - holland , 149 - 156 .
p. wynn , general purpose vector epsilon algorithm algol procedures , numer .
math . , 6(1964 ) 22 - 36 .
p. wynn , on the convergence and stability of the epsilon algorithm , siam j. numer .
, 3(1966 ) 91 - 122 .
j. yang , _ nonlinear waves in integrable and nonintegrable systems _ , siam , philadelphia , 2010 .
j. yang , classification of solitary wave bifurcations in generalized nonlinear schrdinger equation , stud .
129(2012 ) 133162 .
c. yang , j. c. meza , b. lee , l. w. wang , kssolv- a matlab toolbox for solving the kohn - sham equations , acm trans .
software , 36 ( 2009 ) 1 - 35 .