SUBROUTINE DSVRGP(N,RA,RB,IPERM) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DIMENSION RA(N),RB(N),IPERM(N) DO 10 I=1,N RB(I)=RA(I) 10 CONTINUE CALL DPSORT(RB,N,IPERM,2,IER) RETURN END