pull4.for

Go to the documentation of this file.

      SUBROUTINE pull4(MAT,SPATIAL,FINV,DET,NDI)
C>        PULL-BACK TIMES DET OF 4TH ORDER TENSOR
       IMPLICIT NONE
       include 'PARAM_UMAT.INC'
C
       INTEGER I1,J1,K1,L1,II1,JJ1,KK1,LL1,NDI
       DOUBLE PRECISION MAT(ndi,ndi,ndi,ndi),FINV(ndi,ndi)
       DOUBLE PRECISION SPATIAL(ndi,ndi,ndi,ndi)
       DOUBLE PRECISION AUX,DET
C
       DO i1=1,ndi
        DO j1=1,ndi
         DO k1=1,ndi
          DO l1=1,ndi
           aux=zero
           DO ii1=1,ndi
            DO jj1=1,ndi
             DO kk1=1,ndi
              DO ll1=1,ndi
              aux=aux+det*
     +        finv(i1,ii1)*finv(j1,jj1)*
     +        finv(k1,kk1)*finv(l1,ll1)*spatial(ii1,jj1,kk1,ll1)
              END DO
             END DO
            END DO
           END DO
           mat(i1,j1,k1,l1)=aux
          END DO
         END DO
        END DO
       END DO
C
       RETURN
      END SUBROUTINE pull4