program cache_blocking
c
c     Este es un programa para probar como funciona el CACHE
c     mediante el producto de dos vectores a y b cuyo resultado se guarda
c     en prod. Cada arreglo es de dmiension DIM <= 30 000 000.
c
      integer DIM, TOT
      parameter (DIM= 34 000,TOT=2*DIM)
      double precision a(DIM), b(DIM), prod
      data a,b/TOT*0.d0/
c
      prod = 0.d0
      do i=1,dim
c         write(*,*) 'a=', a(i),' b=', b(i)
         prod = prod + a(i)*b(i)
      enddo
      write(*,*) prod
c
      end