anchor.diff:         hxnorm = hynorm
          
          iter = iter + 1
  -
  -       currfeas = max(epsfeas, min(hynorm, currfeas) * rho)
   
          ! Verify convergence conditions
   
  
positive.code:   implicit none
  
    ! PARAMETERS
    real(8), parameter :: BETA = 1.0D+04
    real(8), parameter :: ALPHA = 1.0D-4
    real(8), parameter :: MU = 1.0D-1
    real(8), parameter :: ETA = 2.5D-1
    ! Restoration reduction factor
  
negative.code: 
         xp(1:n) = x(1:n)
  
         ! Update current feasibility
         
         currfeas = max(epsfeas, min(currfeas, (1.0D0 - ALPHA) * hxnorm))
  
         if ( verbose ) write(*,905)
  