CTITLESAPK12 - SUBROUTINE USED BY SDPICK 00010000 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA 00020000 CA DESIGNER P. C. LUH 00030000 CA AUTHOR P. C. LUH 00040000 CA LANGUAGE FORTRAN 00050000 CA SYSTEM IBM ONLY 00060000 CA DATE MM/DD/YY 00070000 C REVISED 03/06/86 LBL OLD NAME = MORDR 00080000 C 00090000 C 00100000 C SUBROUTINE MORDR(N,B,A,D1,D2) 00110000 C 00120000 SUBROUTINE SAPK12(N,B,A,D1,D2) 00130000 C 00140000 DIMENSION B(1),A(1) 00150000 C 00160000 D1=B(1) 00170000 IF(N-2) 10,20,20 00180000 10 D2=B(1) 00190000 RETURN 00200000 C 00210000 20 M2=N/2 00220000 NP=(N+1)/2 00230000 IS=MOD(N,2)+1 00240000 NNP=N 00250000 C 00260000 DO 40 I=1,M2 00270000 X=B(I) 00280000 YK=B(I) 00290000 KX=I 00300000 KN=I 00310000 K=I+1 00320000 C 00330000 DO 30 L=K,NNP 00340000 IF(B(L).LT.X) GO TO 25 00350000 X=B(L) 00360000 KX=L 00370000 25 IF(YK.LE.B(L)) GO TO 30 00380000 YK=B(L) 00390000 KN=L 00400000 30 CONTINUE 00410000 C 00420000 Q=B(I) 00430000 B(I)=YK 00440000 B(KN)=Q 00450000 Q=A(I) 00460000 A(I)=A(KN) 00470000 A(KN)=Q 00480000 IF(KX.EQ.I) KX=KN 00490000 Q=B(NNP) 00500000 B(NNP)=X 00510000 B(KX)=Q 00520000 Q=A(NNP) 00530000 A(NNP)=A(KX) 00540000 A(KX)=Q 00550000 40 NNP=NNP-1 00560000 C 00570000 D1=B(NP) 00580000 GO TO (50,60),IS 00590000 50 D2=B(NP+1) 00600000 RETURN 00610000 C 00620000 60 D2=B(NP) 00630000 C 00640000 RETURN 00650000 C 00660000 END 00670000