pull2.for

Go to the documentation of this file.

      SUBROUTINE pull2(PK,SIG,FINV,DET,NDI)
C>       PULL-BACK TIMES DET OF A 2ND ORDER TENSOR
       IMPLICIT NONE
       include 'PARAM_UMAT.INC'
C
       INTEGER I1,J1,II1,JJ1,NDI
       DOUBLE PRECISION PK(ndi,ndi),FINV(ndi,ndi)
       DOUBLE PRECISION SIG(ndi,ndi)
       DOUBLE PRECISION AUX,DET
C
       DO i1=1,ndi
        DO j1=1,ndi
          aux=zero
         DO ii1=1,ndi
          DO jj1=1,ndi
            aux=aux+det*finv(i1,ii1)*finv(j1,jj1)*sig(ii1,jj1)
         END DO
        END DO
           pk(i1,j1)=aux
        END DO
       END DO
C
       RETURN
      END SUBROUTINE pull2