anchor.diff:        ierr = 10 * n
         go to 50
   c
  - 10    call  reduc(nm,n,a,b,fv2,ierr)
  + 10    call  util_reduc(nm,n,a,b,fv2,ierr)
         if (ierr .ne. 0) go to 50
         if (matz .ne. 0) go to 20
   c     .......... find eigenvalues only ..........
  -      call  tred1(nm,n,a,w,fv1,fv2)
  -      call  tqlrat(n,w,fv2,ierr)
  +      call  util_tred1(nm,n,a,w,fv1,fv2)
  +      call  util_tqlrat(n,w,fv2,ierr)
         go to 50
   c     .......... find both eigenvalues and eigenvectors ..........
  - 20    call  tred2(nm,n,a,w,fv1,z)
  -      call  tql2(nm,n,w,fv1,z,ierr)
  + 20    call  util_tred2(nm,n,a,w,fv1,z)
  +      call  util_tql2(nm,n,w,fv1,z,ierr)
         if (ierr .ne. 0) go to 50
  -      call  rebak(nm,n,b,fv2,n,z)
  +      call  util_rebak(nm,n,b,fv2,n,z)
    50    return
         end
  -      subroutine tql2(nm,n,d,e,z,ierr)
  +      subroutine util_tql2(nm,n,d,e,z,ierr)
   c
         integer i,j,k,l,m,n,ii,l1,l2,nm,mml,ierr
         double precision d(n),e(n),z(nm,n)
  
positive.code:       util_epslon = eps*dabs(x)
        return
        end
        subroutine tqlrat(n,d,e2,ierr)
  c
        integer i,j,l,m,n,ii,l1,mml,ierr
        double precision d(n),e2(n)
  
negative.code: c
   200                           return
        end
        subroutine util_reduc(nm,n,a,b,dl,ierr)
  c
        integer i,j,k,n,i1,j1,nm,nn,ierr
        double precision a(nm,n),b(nm,n),dl(n)
  