黄永刚单晶塑性有限元umat子程序

DUM=A(IMAX,K) A(IMAX,K)=A(J,K) A(J,K)=DUM END DO

D=-D

VV(IMAX)=VV(J) END IF

INDX(J)=IMAX

IF (A(J,J).EQ.0.) A(J,J)=TINY IF (J.NE.N) THEN DUM=1./A(J,J) DO I=J+1,N

A(I,J)=A(I,J)*DUM END DO END IF

END DO

RETURN END

C----------------------------------------------------------------------

SUBROUTINE LUBKSB (A, N, NP, INDX, B)

C----- Linear equation solver based on LU decomposition

C----- Use single precision on cray C

IMPLICIT REAL*8 (A-H,O-Z)

DIMENSION A(NP,NP), INDX(N), B(N)

II=0

DO I=1,N

LL=INDX(I) SUM=B(LL) B(LL)=B(I)

IF (II.NE.0) THEN DO J=II,I-1

SUM=SUM-A(I,J)*B(J) END DO

ELSE IF (SUM.NE.0.) THEN II=I END IF

B(I)=SUM END DO

DO I=N,1,-1 SUM=B(I)

IF (I.LT.N) THEN DO J=I+1,N

SUM=SUM-A(I,J)*B(J) END DO END IF

B(I)=SUM/A(I,I) END DO

RETURN END