princ_u.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	princ_u (nel, es1, es2, es3, es4, es5, es6, ev, dirprv)

Function/Subroutine Documentation

princ_u()

subroutine princ_u	(	integer	nel,
			es1,
			es2,
			es3,
			es4,
			es5,
			es6,
			ev,
			dirprv )

Definition at line 30 of file princ_u.F.

C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include "implicit_f.inc"
C-----------------------------------------------
C   G L O B A L   P A R A M E T E R S
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C O M M O N 
C-----------------------------------------------
C----------------------------------------------------------------
C  I N P U T   A R G U M E N T S
C----------------------------------------------------------------
      INTEGER       NEL     
      my_real
     .   es1(mvsiz), es2(mvsiz), es3(mvsiz), 
     .   es4(mvsiz), es5(mvsiz), es6(mvsiz),
     .   ev(mvsiz,3),dirprv(mvsiz,3,3)
C----------------------------------------------------------------
C  L O C A L  V A R I B L E S
C----------------------------------------------------------------
      INTEGER    I,J,K
      my_real
     .            av(mvsiz,6),evv(mvsiz,3)
C----------------------------------------------------------------
C--------- [B]=[F][F]^t -1 strain-----
      DO i=1,nel
       av(i,1) = es1(i)
       av(i,2) = es2(i)
       av(i,3) = es3(i)
       av(i,4) = es4(i)
       av(i,5) = es5(i)
       av(i,6) = es6(i)
      ENDDO
C
C         for a simple precision executing
c      IF (IRESP==1) THEN
c        CALL VALPVECDP(AV,EVV,DIRPRV,NEL)
c      ELSE
c        CALL VALPVEC(AV,EVV,DIRPRV,NEL)
c      ENDIF
        CALL valpvecop_v(av,evv,dirprv,nel)
 
      DO i=1,nel
        ev(i,1)=half*log(evv(i,1)+ one)
        ev(i,2)=half*log(evv(i,2)+ one)
        ev(i,3)=half*log(evv(i,3)+ one)
      ENDDO
C----------------
C      
      RETURN