NUM(IR)=NUM(K)
G(IR)=G(K)
CALL GG1(IR,IP+1,IP,N1)
C..
IF(ISIMIN.EQ.0) GO TO 1111
1000
CONTINUE
RETURN
END
C------------------------------------------------------------------C
C------------------------------------------------------------------C
SUBROUTINE GG1(ICOL,IP1,IPS,N1)
DOUBLE PRECISION A,G,X
DIMENSION A(10,40),G(40),NUM(40)
COMMON /MASS1/ A,G,NUM,X(40)
DO 1 I=IP1,N1
A(ICOL,I)=A(ICOL,I)/A(ICOL,ICOL)
1
CONTINUE
A(ICOL,ICOL)=1.
DO 2 J=1,IPS
IF(J.EQ.ICOL) GO TO 2
DO 3 I=IP1,N1
A(J,I)=A(J,I)-A(ICOL,I)*A(J,ICOL)
3
CONTINUE
A(J,ICOL)=0.
2
CONTINUE
RETURN
END
C------------------------------------------------------------------C
C------------------------------------------------------------------C
SUBROUTINE GG(INI,NP,NP1)
DOUBLE PRECISION P1,AM,AT
COMMON/COMGG/ P1(20,21)
DO 201 I=INI,NP
AM=P1(I,I)
I1=I
IJ=I+1
DO 202 K=IJ,NP
IF(ABS(AM).GE.DABS(P1(K,I))) GO TO 202
AM=P1(K,I)
I1=K
202
CONTINUE
28
Previous Page
