pythiaeRHIC/pythia/pybksb.f

0001
0002 C*********************************************************************
0003
0004 C...PYBKSB
0005 C...Auxiliary to PYSIGH, for technicolor corrections to QCD 2 -> 2
0006 C...processes.
0007
0008       SUBROUTINE PYBKSB(A,N,NP,INDX,B)
0009       IMPLICIT NONE
0010       INTEGER N,NP,INDX(N)
0011       COMPLEX*16 A(NP,NP),B(N)
0012       INTEGER I,II,J,LL
0013       COMPLEX*16 SUM
0014
0015       II=0
0016       DO 110 I=1,N
0017         LL=INDX(I)
0018         SUM=B(LL)
0019         B(LL)=B(I)
0020         IF (II.NE.0)THEN
0021           DO 100 J=II,I-1
0022             SUM=SUM-A(I,J)*B(J)
0023   100     CONTINUE
0024         ELSE IF (ABS(SUM).NE.0D0) THEN
0025           II=I
0026         ENDIF
0027         B(I)=SUM
0028   110 CONTINUE
0029       DO 130 I=N,1,-1
0030         SUM=B(I)
0031         DO 120 J=I+1,N
0032           SUM=SUM-A(I,J)*B(J)
0033   120   CONTINUE
0034         B(I)=SUM/A(I,I)
0035   130 CONTINUE
0036       RETURN
0037       END