anchor.diff:         
          iter = iter + 1
   
  -       currfeas = max(epsfeas, min(hynorm, currfeas) * curropt)
  -
  -       curropt = max(epsopt, min(rho, dnorm) / (1.1D0 ** iter))
  +       currfeas = max(epsfeas, min(hynorm, currfeas) * rho)
   
          ! Verify convergence conditions
   
  
positive.code:     ! LOCAL SCALARS
      logical :: isforb,isalph,isbeta
      integer :: i,j,k,iter,jcnnz,m,nf
      real(8) :: c,currfeas,curropt,delta,dnorm,dxznorm,fx,fy,fz,hxnorm,&
           hznorm,hynorm,rho
  
      iter = 1
  
negative.code:     uevaljac => evaljac_
  
      m = me + mi
  
      rho = RHOINI
  
  