tensorprod22.for

Go to the documentation of this file.

       SUBROUTINE tensorprod2(A,B,C,NDI)
c
       Implicit None
c
       INTEGER I,J,K,L,NDI
c
       DOUBLE PRECISION A(ndi,ndi),B(ndi,ndi),C(ndi,ndi,ndi,ndi)
c
      DO i=1,ndi
       DO j=1,ndi
         DO k=1,ndi
          DO l=1,ndi
          c(i,j,k,l)=a(i,j)*b(k,l)
          END DO
         END DO
       END DO
      END DO
c
      RETURN
c
      end SUBROUTINE tensorprod2