********************************************************************************
*                                                                              *
*  2DHF version 1-2003                                                         *
*  Copyright (C) 1996  Jacek Kobus, Leif Laaksonen, Dage Sundholm              *
*                                                                              *
*  This software may be used and distributed according to the terms            *
*  of the GNU General Public License, see README and COPYING.                  *
*                                                                              *
********************************************************************************
c ### prod ###
c
c     a vector times a vector: dy(i)=dx(i)*dy(i), i=1..n
c
      subroutine prod(n,dx,dy)
      real*8	 dx(*),dy(*)
      integer*4 i,n

      do i = 1,n
         dy(i) = dy(i) * dx(i)
      enddo

      return
      end