********************************************************************************
*                                                                              *
*  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 ### prodas ###
c
c     a vector times a vector: dr(i)=s*dx(i)*dy(i)+dr(i), i=1..n
c
      subroutine   prodas(n,s,dx,dy,dr)
      real*8	 dx(*),dy(*),dr(*),s
      integer*4 i,n

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

      return
      end