contraction44.for

Go to the documentation of this file.

      SUBROUTINE contraction44(S,LT,RT,NDI)
C>       DOUBLE CONTRACTION BETWEEN 4TH ORDER TENSORS
C>      INPUT:
C>       LT - RIGHT 4TH ORDER TENSOR
C>       RT - LEFT  4TH ORDER TENSOR
C>      OUTPUT:
C>       S - DOUBLE CONTRACTED TENSOR (4TH ORDER)
       IMPLICIT NONE
       include 'PARAM_UMAT.INC'
C
       INTEGER I1,J1,K1,L1,M1,N1,NDI
C
       DOUBLE PRECISION LT(ndi,ndi,ndi,ndi),RT(ndi,ndi,ndi,ndi)
       DOUBLE PRECISION S(ndi,ndi,ndi,ndi)
       DOUBLE PRECISION AUX
C                     
       DO i1=1,ndi
        DO j1=1,ndi
         DO k1=1,ndi
          DO l1=1,ndi
           aux=zero
           DO m1=1,ndi
            DO n1=1,ndi
                aux=aux+lt(i1,j1,m1,n1)*rt(m1,n1,k1,l1)
            END DO
           END DO
           s(i1,j1,k1,l1)=aux
          END DO
         END DO
        END DO
       END DO
C
       RETURN
      END SUBROUTINE contraction44