anchor.diff:         ! Verify convergence conditions
   
          if ( hynorm .le. epsfeas .and. &
  -            dnorm .le. epsopt ) then
  +            dzynorm .le. epsopt ) then
             flag = 0
             exit
          end if
  
positive.code: 901 FORMAT('H-iteration: the pair (',1PD17.8,',',1PD17.8,') was added.',/)
  902 FORMAT('F-iteration.',/)
  903 FORMAT(/,'Filter update',/,13('-'))
  904 FORMAT('Current point (first',1X,I5,' elements):',/,2X, &
             4(1X,1PD16.8))
  905 FORMAT(/,'Restoration Phase',/,17('-'),/)
  906 FORMAT(3X,'Forbidden?',56X,L,/,3X,'Alpha?',60X,L,/,3X,'Beta?',61X,L,/)
  907 FORMAT(3X,'Updated feasibility requirement to',27X,E9.1,/)
  908 FORMAT(3X,'F(Z) =',45X,1PD16.8,/,3X,'H(Z) =',45X,D16.8,/,3X, &
             'Restored point (first',1X,I5,' elements):',/,3X,16X, &
             3(1X,1PD16.8))
  909 FORMAT(/,'Optimization Phase',/,18('-'),/)
  910 FORMAT(3X,'F(Z+D) =',43X,1PD16.8,/,3X,'H(Z+D) =',43X,D16.8,/,   &
             3X,'DELTA =',44X,1PD16.8,/,3X,'RHO =',46X,1PD16.8,/,     &
             3X,'||D|| =',44X,D16.8,/,3X,'Optimized point (first',1X, &
             I5,' elements):',/,3X,16X,3(1X,1PD16.8))
  911 FORMAT(/,70('-'),/,'Final Iteration',/,70('-'),/,'F(X) =',48X, &
             1PD16.8,/,'H(X) =',48X,D16.8,/,   &
             'Function Evaluations =',28X,I20,/,'Solution (first',1X, &
             I5,' elements):',/,2X,4(1X,1PD16.8))
  912 FORMAT(3X,'WARNING! Solver return FLAG ',I3,'!',/)
  
    end subroutine dfoirfalg
  
negative.code:     ! LOCAL SCALARS
      logical :: isforb,isalph,isbeta
      integer :: i,j,k,iter,jcnnz,m,nf
      real(8) :: c,currfeas,delta,dzynorm,dxznorm,fx,fy,fz,hxnorm,&
           hznorm,hynorm,rho
  
      iter = 1
  