home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / itrpapsm / sg13.for < prev   
Encoding:
Text File  |  1985-11-29  |  2.4 KB  |  79 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE SG13
  5. C
  6. C        PURPOSE
  7. C           TO COMPUTE A VECTOR OF SMOOTHED FUNCTION VALUES GIVEN
  8. C           VECTORS OF ARGUMENT VALUES AND CORRESPONDING FUNCTION
  9. C           VALUES.
  10. C
  11. C        USAGE
  12. C           CALL SG13(X,Y,Z,NDIM,IER)
  13. C
  14. C        DESCRIPTION OF PARAMETERS
  15. C           X     -  GIVEN VECTOR OF ARGUMENT VALUES (DIMENSION NDIM)
  16. C           Y     -  GIVEN VECTOR OF FUNCTION VALUES CORRESPONDING TO X
  17. C                    (DIMENSION NDIM)
  18. C           Z     -  RESULTING VECTOR OF SMOOTHED FUNCTION VALUES
  19. C                    (DIMENSION NDIM)
  20. C           NDIM  -  DIMENSION OF VECTORS X,Y,AND Z
  21. C           IER   -  RESULTING ERROR PARAMETER
  22. C                    IER = -1  - NDIM IS LESS THAN 3
  23. C                    IER =  0  - NO ERROR
  24. C
  25. C        REMARKS
  26. C           (1)  IF IER=-1 THERE HAS BEEN NO COMPUTATION.
  27. C           (2)   Z CAN HAVE THE SAME STORAGE ALLOCATION AS X OR Y.  IF
  28. C                 X OR Y IS DISTINCT FROM Z, THEN IT IS NOT DESTROYED.
  29. C
  30. C        SUBROUTINES AND SUBPROGRAMS REQUIRED
  31. C           NONE
  32. C
  33. C        METHOD
  34. C           EXCEPT AT THE ENDPOINTS X(1) AND X(NDIM), EACH SMOOTHED
  35. C           VALUE Z(I) IS OBTAINED BY EVALUATING AT X(I) THE LEAST-
  36. C           SQUARES POLYNOMIAL OF DEGREE 1 RELEVANT TO THE 3 SUCCESSIVE
  37. C           POINTS (X(I+K),Y(I+K)) K = -1,0,1.(SEE HILDEBRAND, F.B.,
  38. C           INTRODUCTION  TO NUMERICAL ANALYSIS, MC GRAW-HILL, NEW YORK/
  39. C           TORONTO/LONDON, 1956, PP.258-311.)
  40. C
  41. C     ..................................................................
  42. C
  43.       SUBROUTINE SG13(X,Y,Z,NDIM,IER)
  44. C
  45. C
  46.       DIMENSION X(1),Y(1),Z(1)
  47. C
  48. C        TEST OF DIMENSION
  49.       IF(NDIM-3)7,1,1
  50. C
  51. C        START LOOP
  52.     1 DO 6 I=3,NDIM
  53.       XM=.3333333*(X(I-2)+X(I-1)+X(I))
  54.       YM=.3333333*(Y(I-2)+Y(I-1)+Y(I))
  55.       T1=X(I-2)-XM
  56.       T2=X(I-1)-XM
  57.       T3=X(I)-XM
  58.       XM=T1*T1+T2*T2+T3*T3
  59.       IF(XM)3,3,2
  60.     2 XM=(T1*(Y(I-2)-YM)+T2*(Y(I-1)-YM)+T3*(Y(I)-YM))/XM
  61. C
  62. C        CHECK FIRST POINT
  63.     3 IF(I-3)4,4,5
  64.     4 H=XM*T1+YM
  65.     5 Z(I-2)=H
  66.     6 H=XM*T2+YM
  67. C        END OF LOOP
  68. C
  69. C        UPDATE LAST TWO COMPONENTS
  70.       Z(NDIM-1)=H
  71.       Z(NDIM)=XM*T3+YM
  72.       IER=0
  73.       RETURN
  74. C
  75. C        ERROR EXIT IN CASE NDIM IS LESS THAN 3
  76.     7 IER=-1
  77.       RETURN
  78.       END
  79.