Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/dse35.ssp
There are 2 other files named dse35.ssp in the archive. Click here to see a list.
C                                                                       DSE3  10
C     ..................................................................DSE3  20
C                                                                       DSE3  30
C        SUBROUTINE DSE35                                               DSE3  40
C                                                                       DSE3  50
C        PURPOSE                                                        DSE3  60
C           TO COMPUTE A VECTOR OF SMOOTHED FUNCTION VALUES GIVEN A     DSE3  70
C           VECTOR OF FUNCTION VALUES WHOSE ENTRIES CORRESPOND TO       DSE3  80
C           EQUIDISTANTLY SPACED ARGUMENT VALUES.                       DSE3  90
C                                                                       DSE3 100
C        USAGE                                                          DSE3 110
C           CALL DSE35(Y,Z,NDIM,IER)                                    DSE3 120
C                                                                       DSE3 130
C        DESCRIPTION OF PARAMETERS                                      DSE3 140
C           Y     -  GIVEN VECTOR OF DOUBLE PRECISION FUNCTION VALUES   DSE3 150
C                    (DIMENSION NDIM)                                   DSE3 160
C           Z     -  RESULTING VECTOR OF DOUBLE PRECISION SMOOTHED      DSE3 170
C                    FUNCTION VALUES (DIMENSION NDIM)                   DSE3 180
C           NDIM  -  DIMENSION OF VECTORS Y AND Z                       DSE3 190
C           IER   -  RESULTING ERROR PARAMETER                          DSE3 200
C                    IER = -1  - NDIM IS LESS THAN 5                    DSE3 210
C                    IER =  0  - NO ERROR                               DSE3 220
C                                                                       DSE3 230
C        REMARKS                                                        DSE3 240
C           (1)  IF IER=-1 THERE HAS BEEN NO COMPUTATION.               DSE3 250
C           (2)   Z CAN HAVE THE SAME STORAGE ALLOCATION AS Y.  IF Y IS DSE3 260
C                 DISTINCT FROM Z, THEN IT IS NOT DESTROYED.            DSE3 270
C                                                                       DSE3 280
C        SUBROUTINE AND FUNCTION SUBPROGRAMS REQUIRED                   DSE3 290
C           NONE                                                        DSE3 300
C                                                                       DSE3 310
C        METHOD                                                         DSE3 320
C           IF X IS THE (SUPPRESSED) VECTOR OF ARGUMENT VALUES, THEN    DSE3 330
C           EXCEPT AT THE POINTS X(1),X(2),X(NDIM-1) AND X(NDIM), EACH  DSE3 340
C           SMOOTHED VALUE Z(I) IS OBTAINED BY EVALUATING AT X(I) THE   DSE3 350
C           LEAST-SQUARES POLYNOMIAL OF DEGREE 3 RELEVANT TO THE 5      DSE3 360
C           SUCCESSIVE POINTS (X(I+K),Y(I+K)) K = -2,-1,...,2.  (SEE    DSE3 370
C           HILDEBRAND, F.B., INTRODUCTION TO NUMERICAL ANALYSIS,       DSE3 380
C           MC GRAW-HILL, NEW YORK/TORONTO/LONDON, 1956, PP. 295-302.)  DSE3 390
C                                                                       DSE3 400
C     ..................................................................DSE3 410
C                                                                       DSE3 420
      SUBROUTINE DSE35(Y,Z,NDIM,IER)                                    DSE3 430
C                                                                       DSE3 440
C                                                                       DSE3 450
      DIMENSION Y(1),Z(1)                                               DSE3 460
      DOUBLE PRECISION Y,Z,A,B,C,D                                      DSE3 470
C                                                                       DSE3 480
C        TEST OF DIMENSION                                              DSE3 490
      IF(NDIM-5)4,1,1                                                   DSE3 500
C                                                                       DSE3 510
C        PREPARE LOOP                                                   DSE3 520
    1 B=Y(1)                                                            DSE3 530
      C=Y(2)                                                            DSE3 540
C                                                                       DSE3 550
C        START LOOP                                                     DSE3 560
      DO 3 I=5,NDIM                                                     DSE3 570
      A=B                                                               DSE3 580
      B=C                                                               DSE3 590
      C=Y(I-2)                                                          DSE3 600
C                                                                       DSE3 610
C        GENERATE FOURTH CENTRAL DIFFERENCE                             DSE3 620
      D=C-B-Y(I-1)                                                      DSE3 630
      D=D+D+C                                                           DSE3 640
      D=D+D+A+Y(I)                                                      DSE3 650
C                                                                       DSE3 660
C        CHECK FIRST TWO COMPONENTS                                     DSE3 670
      IF(I-5)2,2,3                                                      DSE3 680
    2 Z(1)=A-.014285714285714286D0*D                                    DSE3 690
      Z(2)=B+.057142857142857143D0*D                                    DSE3 700
    3 Z(I-2)=C-.08571428571428571D0*D                                   DSE3 710
C        END OF LOOP                                                    DSE3 720
C                                                                       DSE3 730
C        UPDATE LAST TWO COMPONENTS                                     DSE3 740
      Z(NDIM-1)=Y(NDIM-1)+.057142857142857143D0*D                       DSE3 750
      Z(NDIM)=Y(NDIM)-.014285714285714286D0*D                           DSE3 760
      IER=0                                                             DSE3 770
      RETURN                                                            DSE3 780
C                                                                       DSE3 790
C        ERROR EXIT IN CASE NDIM IS LESS THAN 5                         DSE3 800
    4 IER=-1                                                            DSE3 810
      RETURN                                                            DSE3 820
      END                                                               DSE3 830