proj_eulerian.for

Go to the documentation of this file.

      SUBROUTINE projeul(A,AA,PE,NDI)
C>    EULERIAN PROJECTION TENSOR
C      INPUTS:
C          IDENTITY TENSORS - A, AA
C      OUTPUTS:
C          4TH ORDER SYMMETRIC EULERIAN PROJECTION TENSOR - PE
C
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER I,J,K,L,NDI
C
      DOUBLE PRECISION A(ndi,ndi),AA(ndi,ndi,ndi,ndi),
     1                 pe(ndi,ndi,ndi,ndi)
C
      DO i=1,ndi
         DO j=1,ndi
          DO k=1,ndi
             DO l=1,ndi
              pe(i,j,k,l)=aa(i,j,k,l)-(one/three)*(a(i,j)*a(k,l))
           END DO
          END DO
        END DO
      END DO
C
      RETURN
      END SUBROUTINE projeul