anchor.diff:            xp(i) = x(i)
          end do
   
  -       rho = max(10.0D0 * curropt, rho)
  -
          delta = max(DELMIN, rho, delta)
   
  +       rho = max(10.0D0 * epsopt, delta)
  +
          call qpsolver(n,x,l,u,me,mi,aevalf,aevalc,levalc,levaljac, &
  -            nf,ALPHA,ffilter,hfilter,currfeas,curropt,verbose,delta, &
  +            nf,ALPHA,ffilter,hfilter,currfeas,epsopt,verbose,delta, &
               fy,hynorm,rho,flag)
   
          dnorm = evalDist(n,xp,x)
  
positive.code: 
    ! PARAMETERS
    real(8), parameter :: BETA = 1.0D+04
    real(8), parameter :: ALPHA = 1.0D-1
    real(8), parameter :: MU = 1.0D-1
    real(8), parameter :: ETA = 2.5D-1
    ! Restoration reduction factor
  
negative.code:     ! LOCAL SCALARS
      logical :: isforb,isalph,isbeta
      integer :: i,j,k,iter,jcnnz,m,nf
      real(8) :: c,currfeas,delta,dnorm,dxznorm,fx,fy,fz,hxnorm,&
           hznorm,hynorm,rho
  
      iter = 1
  