rotation.for

Go to the documentation of this file.

      SUBROUTINE rotation(F,R,U,NDI)
C>    COMPUTES ROTATION TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI
      DOUBLE PRECISION F(ndi,ndi),R(ndi,ndi),U(ndi,ndi),UINV(ndi,ndi)
C
      CALL matinv3d(u,uinv,ndi)
C
      r = matmul(f,uinv)
      RETURN
      END SUBROUTINE rotation