anchor.diff:  
          ! Update current feasibility
          
  -       currfeas = max(epsfeas, min(currfeas, (1.0D0 - ALPHA) * hxnorm))
  +       currfeas = max(epsfeas, min((1.0D0 - ALPHA) * hxnorm, hxnorm * dzynorm))
   
          if ( verbose ) write(*,905)
   
  
positive.code: 
            call restore(n,x,l,u,me,mi,aevalc,aevaljac,currfeas,verbose, &
                 hznorm,flag)
  
            call aevalf(n,x,fz,flag)
  
  
negative.code:           end if
  
         end do
  
         if ( flag .ne. 0 ) exit
  
         if ( verbose ) WRITE(*,908) fz,hznorm,min(n,MAXNEL), &
                        (x(i), i = 1,min(MAXNEL,n))
  