anchor.diff:    ! Maximum number of iterations
     integer, parameter :: MAXITER = 10000
     ! Maximum number of printing elements
  -  integer, parameter :: MAXNEL  = 3
  +  integer, parameter :: MAXNEL  = 20
     ! Minimum trust region radius
     real(8), parameter :: DELMIN = 1.0D-30
     ! Initial value of RHO
  
positive.code:        hxnorm = hynorm
         
         iter = iter + 1
  
         currfeas = max(epsfeas, min(hynorm, currfeas) * rho)
  
         ! Verify convergence conditions
  
  
negative.code: 
         xp(1:n) = x(1:n)
  
         ! Update current feasibility
         
         currfeas = max(epsfeas, min(currfeas, (1.0D0 - ALPHA) * hxnorm))
  
         if ( verbose ) write(*,905)
  