home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libconv / TRY / simple_dcorm1d.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  984 b   |  37 lines
  1.       subroutine simple_dcorm1d(   n_fx, ny, n_gx, f, g, h)
  2.       double precision f(0:(n_fx-1),ny)
  3.       double precision g(0:(n_gx-1))
  4.       double precision h(0:(n_fx+n_gx-2),ny)
  5.       integer n_fx, n_gx, ny
  6. c-----------------------------------------------------------------------------
  7. c
  8. c   This Fortran Subroutine written by
  9. c    Jean-Pierre Panziera
  10. c    Silicon Graphics Inc
  11. c    September 12, 1991
  12. c
  13. c-----------------------------------------------------------------------------
  14.     integer j, ihx, ifx, igx
  15.     double precision zero
  16.     parameter ( zero = 0.)
  17. c-----------------------------------------------------------------------------
  18. c
  19. c   Compute the N 1D correlations
  20. c
  21. c-----------------------------------------------------------------------------
  22.     do j = 1, ny
  23.         do ihx = 0, n_fx+n_gx-2
  24.         h(ihx,j) = zero
  25.         end do
  26.         do ifx = 0, n_fx-1
  27.         do  igx =  0, n_gx-1
  28.             h(ifx+igx,j) = h(ifx+igx,j) + f(ifx,j) * g(n_gx-1-igx)
  29.         end do
  30.         end do
  31.     end do
  32.  
  33.       return
  34.       end
  35.  
  36.  
  37.