contraction42.for

Go to the documentation of this file.

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