anchor.diff:                 ISEED( I ) = 1
    122        CONTINUE
   
  -            CALL ODRNV(2, ISEED, 2*IN-1, WORK(1))
  +            CALL DLARNV(2, ISEED, 2*IN-1, WORK(1))
               DO 125 I = 1,IN-1
                  D(IBEGIN+I-1) = D(IBEGIN+I-1)*(ONE+EPS*PERT*WORK(I))
                  E(IBEGIN+I-1) = E(IBEGIN+I-1)*(ONE+EPS*PERT*WORK(IN+I))
  
positive.code: *        Compute the negcount at the 1/4 and 3/4 points
           IF(MB.GT.1) THEN
              CALL ODRRC( 'T', IN, S1, S2, D(IBEGIN),
       $                   E(IBEGIN), PIVMIN, CNT, CNT1, CNT2, IINFO)
           ENDIF
  
           IF(MB.EQ.1) THEN
  
negative.code:          IF ( .NOT.USEDQD ) THEN
  *           If ODRRD has been used, shift the eigenvalue approximations
  *           according to their representation. This is necessary for
  *           a uniform DLARRV since dqds computes eigenvalues of the
  *           shifted representation. In DLARRV, W will always hold the
  *           UNshifted eigenvalue approximation.
              DO 134 J=WBEGIN,WEND
                 W(J) = W(J) - SIGMA
  