Thursday, 16 June 2011

algorithms - On Clenshaw's summation formula

When one has a finite sum of the form



S=sumnk=0ckFk(x)



where Fk(x) satisfies a two-term recurrence relation



Fk+1(x)=alphakFk(x)+betakFk1(x)



the standard algorithm used for computing S is Clenshaw summation:



yn+2=yn+1=0



yk=alphakyk+1+betak+1yk+2+ck;;;k=n;mathrmto;1;mathrmstep;1



S=y1F1(x)+(beta1y2+c0)F0(x)



I am also very familiar with the concept of "minimal" and "dominant" solutions of a two-term recurrence; to review, the dominant solution is the solution most stably propagated when the two-term recurrence is run in the forward direction, while the minimal solution is the solution that is stably computed when the recurrence is run backwards. As an example, for the recurrence relation of the Bessel functions, Jn(x), the function of the first kind, is the minimal solution while Yn(x), the other solution to the Bessel ODE, is the dominant solution to the two-term recurrence.



Mystifying to me, however, is this remark given in the book "Numerical Recipes" without justification of any sort:



"...Clenshaw's recurrence is always stable, independent of whether the recurrence for the functions... is stable in the upward or downward direction."



My question now, is why would the stability properties of Fk (i.e. its being minimal or dominant) not spoil the recurrence, if this is true? If not, I would be interested in seeing a counterexample.

No comments:

Post a Comment