home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / SMOOTH.FOR < prev    next >
Encoding:
Text File  |  1990-02-20  |  9.0 KB  |  317 lines
  1.       SUBROUTINE DataSmoothSG (dataset, numdat, smoothnum, derivnum,
  2.      +   smoothdata)
  3.       INCLUDE 'STDHDR.FOR'
  4.       REAL dataset(0:maxv), smoothdata(0:maxv)
  5.       INTEGER numdat, smoothnum, derivnum
  6.  
  7.  
  8.       INTEGER normcoef(0: 2, 0: 4)
  9.       INTEGER smoothcoef[ALLOCATABLE](:,:,:)
  10.       INTEGER numcoef, smoothdeg, startconv
  11.       INTEGER stopconv, i, j, k, iErr
  12.  
  13.       ALLOCATE(smoothcoef(0 : 2, 0 : 4, 0 : 12), STAT=iErr)
  14.       CALL CheckMem(iErr)
  15.  
  16.       normcoef(0, 0) = 35
  17.       normcoef(0, 1) = 21
  18.       normcoef(0, 2) = 231
  19.       normcoef(0, 3) = 429
  20.       normcoef(0, 4) = 143
  21.       normcoef(1, 0) = 10
  22.       normcoef(1, 1) = 28
  23.       normcoef(1, 2) = 60
  24.       normcoef(1, 3) = 110
  25.       normcoef(1, 4) = 182
  26.       normcoef(2, 0) = 7
  27.       normcoef(2, 1) = 42
  28.       normcoef(2, 2) = 462
  29.       normcoef(2, 3) = 429
  30.       normcoef(2, 4) = 1001
  31.  
  32.  
  33.       smoothcoef(0, 0, 0) = -3
  34.       smoothcoef(0, 0, 1) = 12
  35.       smoothcoef(0, 0, 2) = 17
  36.       smoothcoef(0, 0, 3) = 12
  37.       smoothcoef(0, 0, 4) = -3
  38.       smoothcoef(0, 0, 5) = 0
  39.       smoothcoef(0, 0, 6) = 0
  40.       smoothcoef(0, 0, 7) = 0
  41.       smoothcoef(0, 0, 8) = 0
  42.       smoothcoef(0, 0, 9) = 0
  43.       smoothcoef(0, 0, 10) = 0
  44.       smoothcoef(0, 0, 11) = 0
  45.       smoothcoef(0, 0, 12) = 0
  46.  
  47.  
  48.       smoothcoef(0, 1, 0) = -2
  49.       smoothcoef(0, 1, 1) = 3
  50.       smoothcoef(0, 1, 2) = 6
  51.       smoothcoef(0, 1, 3) = 7
  52.       smoothcoef(0, 1, 4) = 6
  53.       smoothcoef(0, 1, 5) = 3
  54.       smoothcoef(0, 1, 6) = -2
  55.       smoothcoef(0, 1, 7) = 0
  56.       smoothcoef(0, 1, 8) = 0
  57.       smoothcoef(0, 1, 9) = 0
  58.       smoothcoef(0, 1, 11) = 0
  59.       smoothcoef(0, 1, 11) = 0
  60.       smoothcoef(0, 1, 12) = 0
  61.  
  62.       smoothcoef(0, 2, 0) = -21
  63.       smoothcoef(0, 2, 1) = 14
  64.       smoothcoef(0, 2, 2) = 39
  65.       smoothcoef(0, 2, 3) = 54
  66.       smoothcoef(0, 2, 4) = 59
  67.       smoothcoef(0, 2, 5) = 54
  68.       smoothcoef(0, 2, 6) = 39
  69.       smoothcoef(0, 2, 7) = 14
  70.       smoothcoef(0, 2, 8) = -21
  71.       smoothcoef(0, 2, 9) = 0
  72.       smoothcoef(0, 2, 10) = 0
  73.       smoothcoef(0, 2, 11) = 0
  74.       smoothcoef(0, 2, 12) = 0
  75.  
  76.       smoothcoef(0, 3, 0) = -36
  77.       smoothcoef(0, 3, 1) = 9
  78.       smoothcoef(0, 3, 2) = 44
  79.       smoothcoef(0, 3, 3) = 69
  80.       smoothcoef(0, 3, 4) = 84
  81.       smoothcoef(0, 3, 5) = 89
  82.       smoothcoef(0, 3, 6) = 84
  83.       smoothcoef(0, 3, 7) = 69
  84.       smoothcoef(0, 3, 8) = 44
  85.       smoothcoef(0, 3, 9) = 9
  86.       smoothcoef(0, 3, 10) = -36
  87.       smoothcoef(0, 3, 11) = 0
  88.       smoothcoef(0, 3, 12) = 0
  89.  
  90.  
  91.       smoothcoef(0, 4, 0) = -11
  92.       smoothcoef(0, 4, 1) = 0
  93.       smoothcoef(0, 4, 2) = 9
  94.       smoothcoef(0, 4, 3) = 16
  95.       smoothcoef(0, 4, 4) = 21
  96.       smoothcoef(0, 4, 5) = 24
  97.       smoothcoef(0, 4, 6) = 25
  98.       smoothcoef(0, 4, 7) = 24
  99.       smoothcoef(0, 4, 8) = 21
  100.       smoothcoef(0, 4, 9) = 16
  101.       smoothcoef(0, 4, 10) = 9
  102.       smoothcoef(0, 4, 11) = 0
  103.       smoothcoef(0, 4, 12) = -11
  104.  
  105.       smoothcoef(1, 0, 0) = -2
  106.       smoothcoef(1, 0, 1) = -1
  107.       smoothcoef(1, 0, 2) = 0
  108.       smoothcoef(1, 0, 3) = 1
  109.       smoothcoef(1, 0, 4) = 2
  110.       smoothcoef(1, 0, 5) = 0
  111.       smoothcoef(1, 0, 6) = 0
  112.       smoothcoef(1, 0, 7) = 0
  113.       smoothcoef(1, 0, 8) = 0
  114.       smoothcoef(1, 0, 9) = 0
  115.       smoothcoef(1, 0, 10) = 0
  116.       smoothcoef(1, 0, 11) = 0
  117.       smoothcoef(1, 0, 12) = 0
  118.  
  119.       smoothcoef(1, 1, 0) = -3
  120.       smoothcoef(1, 1, 1) = -2
  121.       smoothcoef(1, 1, 2) = -1
  122.       smoothcoef(1, 1, 3) = 0
  123.       smoothcoef(1, 1, 4) = 1
  124.       smoothcoef(1, 1, 5) = 2
  125.       smoothcoef(1, 1, 6) = 3
  126.       smoothcoef(1, 1, 7) = 0
  127.       smoothcoef(1, 1, 8) = 0
  128.       smoothcoef(1, 1, 9) = 0
  129.       smoothcoef(1, 1, 10) = 0
  130.       smoothcoef(1, 1, 11) = 0
  131.       smoothcoef(1, 1, 12) = 0
  132.  
  133.       smoothcoef(1, 2, 0) = -4
  134.       smoothcoef(1, 2, 1) = -3
  135.       smoothcoef(1, 2, 2) = -2
  136.       smoothcoef(1, 2, 3) = -1
  137.       smoothcoef(1, 2, 4) = 0
  138.       smoothcoef(1, 2, 5) = 1
  139.       smoothcoef(1, 2, 6) = 2
  140.       smoothcoef(1, 2, 7) = 3
  141.       smoothcoef(1, 2, 8) = 4
  142.       smoothcoef(1, 2, 0) = 0
  143.       smoothcoef(1, 2, 10) = 0
  144.       smoothcoef(1, 2, 11) = 0
  145.       smoothcoef(1, 2, 12) = 0
  146.  
  147.  
  148.       smoothcoef(1, 3, 0) = -5
  149.       smoothcoef(1, 3, 1) = -4
  150.       smoothcoef(1, 3, 2) = -3
  151.       smoothcoef(1, 3, 3) = -2
  152.       smoothcoef(1, 3, 4) = -1
  153.       smoothcoef(1, 3, 5) = 0
  154.       smoothcoef(1, 3, 6) = 1
  155.       smoothcoef(1, 3, 7) = 2
  156.       smoothcoef(1, 3, 8) = 3
  157.       smoothcoef(1, 3, 9) = 4
  158.       smoothcoef(1, 3, 10) = 5
  159.       smoothcoef(1, 3, 11) = 0
  160.       smoothcoef(1, 3, 12) = 0
  161.  
  162.  
  163.       smoothcoef(1, 4, 0) = -6
  164.       smoothcoef(1, 4, 1) = -5
  165.       smoothcoef(1, 4, 2) = -4
  166.       smoothcoef(1, 4, 3) = -3
  167.       smoothcoef(1, 4, 4) = -2
  168.       smoothcoef(1, 4, 5) = -1
  169.       smoothcoef(1, 4, 6) = 0
  170.       smoothcoef(1, 4, 7) = 1
  171.       smoothcoef(1, 4, 8) = 2
  172.       smoothcoef(1, 4, 9) = 3
  173.       smoothcoef(1, 4, 10) = 4
  174.       smoothcoef(1, 4, 11) = 5
  175.       smoothcoef(1, 4, 12) = 6
  176.  
  177.  
  178.       smoothcoef(2, 0, 0) = 2
  179.       smoothcoef(2, 0, 1) = -1
  180.       smoothcoef(2, 0, 2) = -2
  181.       smoothcoef(2, 0, 3) = -1
  182.       smoothcoef(2, 0, 4) = 2
  183.       smoothcoef(2, 0, 5) = 0
  184.       smoothcoef(2, 0, 6) = 0
  185.       smoothcoef(2, 0, 7) = 0
  186.       smoothcoef(2, 0, 8) = 0
  187.       smoothcoef(2, 0, 9) = 0
  188.       smoothcoef(2, 0, 10) = 0
  189.       smoothcoef(2, 0, 11) = 0
  190.       smoothcoef(2, 0, 12) = 0
  191.  
  192.  
  193.       smoothcoef(2, 1, 0) = 5
  194.       smoothcoef(2, 1, 1) = 0
  195.       smoothcoef(2, 1, 2) = -3
  196.       smoothcoef(2, 1, 3) = -4
  197.       smoothcoef(2, 1, 4) = -3
  198.       smoothcoef(2, 1, 5) = 0
  199.       smoothcoef(2, 1, 6) = 5
  200.       smoothcoef(2, 1, 7) = 0
  201.       smoothcoef(2, 1, 8) = 0
  202.       smoothcoef(2, 1, 9) = 0
  203.       smoothcoef(2, 1, 10) = 0
  204.       smoothcoef(2, 1, 11) = 0
  205.       smoothcoef(2, 1, 12) = 0
  206.  
  207.  
  208.       smoothcoef(2, 2, 0) = 28
  209.       smoothcoef(2, 2, 1) = 7
  210.       smoothcoef(2, 2, 2) = -8
  211.       smoothcoef(2, 2, 3) = -17
  212.       smoothcoef(2, 2, 4) = -20
  213.       smoothcoef(2, 2, 5) = -17
  214.       smoothcoef(2, 2, 6) = -8
  215.       smoothcoef(2, 2, 7) = 7
  216.       smoothcoef(2, 2, 8) = 28
  217.       smoothcoef(2, 2, 9) = 0
  218.       smoothcoef(2, 2, 10) = 0
  219.       smoothcoef(2, 2, 11) = 0
  220.       smoothcoef(2, 2, 12) = 0
  221.  
  222.  
  223.       smoothcoef(2, 3, 0) = 15
  224.       smoothcoef(2, 3, 1) = 6
  225.       smoothcoef(2, 3, 2) = -1
  226.       smoothcoef(2, 3, 3) = -6
  227.       smoothcoef(2, 3, 4) = -9
  228.       smoothcoef(2, 3, 5) = -10
  229.       smoothcoef(2, 3, 6) = -9
  230.       smoothcoef(2, 3, 7) = -6
  231.       smoothcoef(2, 3, 8) = -1
  232.       smoothcoef(2, 3, 9) = 6
  233.       smoothcoef(2, 3, 10) = 15
  234.       smoothcoef(2, 3, 11) = 0
  235.       smoothcoef(2, 3, 12) = 0
  236.  
  237.  
  238.       smoothcoef(2, 4, 0) = 22
  239.       smoothcoef(2, 4, 1) = 11
  240.       smoothcoef(2, 4, 2) = 2
  241.       smoothcoef(2, 4, 3) = -5
  242.       smoothcoef(2, 4, 4) = -10
  243.       smoothcoef(2, 4, 5) = -13
  244.       smoothcoef(2, 4, 6) = -14
  245.       smoothcoef(2, 4, 7) = -13
  246.       smoothcoef(2, 4, 8) = -10
  247.       smoothcoef(2, 4, 9) = -5
  248.       smoothcoef(2, 4, 10) = 2
  249.       smoothcoef(2, 4, 11) = 11
  250.       smoothcoef(2, 4, 12) = 22
  251.  
  252.  
  253.       IF (smoothnum .ge. 1 .AND. smoothnum .LE. 5
  254.      +   .AND. derivnum .ge. 0 .AND. derivnum .LE. 2) THEN
  255.         numcoef = 2 * smoothnum + 3
  256.         smoothdeg = (numcoef - 1) / 2
  257.         startconv = smoothdeg
  258.         stopconv = numdat - smoothdeg - 1
  259.         DO i = 0, numdat - 1
  260.           smoothdata(i) = 0.0
  261.         END DO
  262.         DO i = startconv, stopconv
  263.           DO j = 0, numcoef - 1
  264.            smoothdata(i) = smoothdata(i) + dataset(i-smoothdeg+j)
  265.      +               * smoothcoef(derivnum, smoothnum - 1, j)
  266.           END DO
  267.           smoothdata(i) =
  268.      +          smoothdata(i)/REAL(normcoef(derivnum, smoothnum-1))
  269.         END DO
  270.         j = smoothdeg - 1
  271.         k = numdat - smoothdeg
  272.         DO i = 0, smoothdeg - 1
  273.           smoothdata(j) = smoothdata(j + 1)
  274.           smoothdata(k) = smoothdata(k - 1)
  275.           j = j - 1
  276.           k = k + 1
  277.         END DO
  278.       END IF
  279.  
  280.       DEALLOCATE(smoothcoef, STAT=iErr)
  281.       CALL CheckDealloc(iErr)
  282.       END !SUB
  283.  
  284.       SUBROUTINE DataSmoothWeights(dataset, numdat, smoothnum, weights,
  285.      +  wdivisor, smoothdata)
  286.       INCLUDE 'STDHDR.FOR'
  287.       REAL dataset(0:maxv), weights(0:maxc)
  288.       REAL smoothdata(0:maxv), wdivisor
  289.       INTEGER numdat, smoothnum
  290.       INTEGER smoothdeg, startconv, stopconv,i, j, k
  291.  
  292.       IF (smoothnum .ge. 1 .AND. MOD(smoothnum, 2) .EQ. 1) THEN
  293.         smoothdeg = smoothnum / 2
  294.         startconv = smoothdeg
  295.         stopconv = numdat - smoothdeg - 1
  296.         DO i = 0, numdat - 1
  297.           smoothdata(i) = 0.0
  298.         END DO
  299.         DO i = startconv, stopconv
  300.           DO j = 0, smoothnum - 1
  301.             smoothdata(i) = smoothdata(i) +
  302.      +             dataset(i - smoothdeg + j) * weights(j)
  303.           END DO
  304.           smoothdata(i) = smoothdata(i) / wdivisor
  305.         END DO
  306.         j = smoothdeg - 1
  307.         k = numdat - smoothdeg
  308.         DO i = 1, smoothdeg
  309.           smoothdata(j) = smoothdata(j + 1)
  310.           smoothdata(k) = smoothdata(k - 1)
  311.           j = j - 1
  312.           k = k + 1
  313.         END DO
  314.       END IF
  315.       END !SUB
  316.  
  317.