Trailing-Edge
-
PDP-10 Archives
-
decuslib10-02
-
43,50145/dteas.doc
There are 2 other files named dteas.doc in the archive. Click here to see a list.
SUBROUTINE DTEAS
PURPOSE
CALCULATE THE LIMIT OF A GIVEN SEQUENCE BY MEANS OF THE
EPSILON-ALGORITHM.
USAGE
CALL DTEAS(X,N,FIN,EPS,IER)
DESCRIPTION OF PARAMETERS
X - DOUBLE PRECISION VECTOR WHOSE COMPONENTS ARE TERMS
OF THE GIVEN SEQUENCE. ON RETURN THE COMPONENTS OF
VECTOR X ARE DESTROYED.
N - DIMENSION OF INPUT VECTOR X.
FIN - RESULTANT SCALAR IN DOUBLE PRECISION CONTAINING ON
RETURN THE LIMIT OF THE GIVEN SEQUENCE.
EPS - SINGLE PRECISION INPUT VALUE, WHICH SPECIFIES THE
UPPER BOUND OF THE RELATIVE (ABSOLUTE) ERROR IF THE
COMPONENTS OF X ARE ABSOLUTELY GREATER (LESS) THAN
ONE.
CALCULATION IS TERMINATED AS SOON AS THREE TIMES IN
SUCCESSION THE RELATIVE (ABSOLUTE) DIFFERENCE
BETWEEN NEIGHBOURING TERMS IS NOT GREATER THAN EPS.
IER - RESULTANT ERROR PARAMETER CODED IN THE FOLLOWING
FORM
IER=0 - NO ERROR
IER=1 - REQUIRED ACCURACY NOT REACHED WITH
MAXIMAL NUMBER OF ITERATIONS
IER=-1 - INTEGER N IS LESS THAN TEN.
REMARKS
NO ACTION BESIDES ERROR MESSAGE IN CASE N LESS THAN TEN.
THE CHARACTER OF THE GIVEN INFINITE SEQUENCE MUST BE
RECOGNIZABLE BY THOSE N COMPONENTS OF THE INPUT VECTOR X.
SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
NONE
METHOD
THE CONVERGENCE OF THE GIVEN SEQUENCE IS ACCELERATED BY
MEANS OF THE E(2)-TRANSFORMATION, USED IN AN ITERATIVE WAY.
FOR REFERENCE, SEE
ALGORITHM 215,SHANKS, CACM 1963, NO. 11, PP. 662. AND
P. WYNN, SINGULAR RULES FOR CERTAIN NON-LINEAR ALGORITHMS
BIT VOL. 3, 1963, PP. 175-195.