anchor.diff:      uevaljac => evaljac_
   
       m = me + mi
  -
  -    curropt = sqrt(epsopt)
   
       rho = RHOINI
   
  
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
  