source/unres/src_MD/md-diff/mts/eigtransf.f

   1 c-------------------------------------------------------------
   2       subroutine eigtransf(n,ndim,a,b,d,c)
   3       implicit none
   4       integer n,ndim
   5       double precision a(ndim,n),b(ndim,n),c(ndim,n),d(ndim)
   6       integer i,j,k
   7       do i=1,n
   8         do j=1,n
   9           c(i,j)=0.0d0
  10           do k=1,n
  11             c(i,j)=c(i,j)+a(i,k)*b(k,j)*d(k)
  12           enddo
  13         enddo
  14       enddo
  15       return
  16       end