anchor.diff:  
          delta = max(DELMIN, rho, delta)
   
  -       rho = max(10.0D0 * epsopt, delta)
  +!       rho = max(10.0D0 * epsopt, delta)
  +       rho = max(10.0D0 * epsopt, min(delta, dzynorm))
   
          call qpsolver(n,x,l,u,me,mi,aevalf,aevalc,levalc,levaljac, &
               nf,ALPHA,ffilter,hfilter,currfeas,epsopt,verbose,delta, &
               fy,hynorm,rho,flag)
   
  -       dnorm = evalDist(n,xp,x)
  +       dzynorm = evalDist(n,xp,x)
   
          if ( flag .ne. 0 ) write(*,912) flag
   
  -       if ( verbose ) write(*,910) fy,hynorm,delta,rho,dnorm, &
  -            min(n,MAXNEL),(x(i), i=1,min(n,MAXNEL))
  +       if ( verbose ) write(*,910) fy,hynorm,delta,rho,dzynorm, &
  +            nfev,min(n,MAXNEL),(x(i), i=1,min(n,MAXNEL))
   
          ! ------------- !
          ! Filter Update !
  
positive.code: 
            ! Just for inequalities
            if ( i .gt. me ) then
               call aevalc(n,x,i,linrhs(i),flag)
               linrhs(i) = linrhs(i) - max(0.0D0, linrhs(i))
            end if
  
            call aevaljac(n,x,i,linvar(k),linval(k),jcnnz,flag)
            linpos(i) = k
  
            do j = k,k + jcnnz - 1
               linrhs(i) = linrhs(i) - linval(j) * x(linvar(j))
            end do
  
            k = k + jcnnz
  
negative.code:     ! LOCAL SCALARS
      logical :: isforb,isalph,isbeta
      integer :: i,j,k,iter,jcnnz,m,nf
      real(8) :: c,currfeas,delta,dzynorm,dxznorm,fx,fy,fz,hxnorm,&
           hznorm,hynorm,rho
  
      iter = 1
  