+++ /dev/null
-c-------------------------------------------------------------
- subroutine eigtransf1(n,ndim,ab,d,c)
- implicit none
- integer n,ndim
- double precision ab(ndim,ndim,n),c(ndim,n),d(ndim)
- integer i,j,k
- do i=1,n
- do j=1,n
- c(i,j)=0.0d0
- do k=1,n
- c(i,j)=c(i,j)+ab(k,j,i)*d(k)
- enddo
- enddo
- enddo
- return
- end