anchor.diff:  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))
  +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))
  +           '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))
  +           3X,'||D|| =',44X,D16.8,/,3X,'Number of f evaluat. =',1X, &
  +           I44,/,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))
  +           I5,' elements):',/,(2X,4(1X,1PD16.8)))
   912 FORMAT(3X,'WARNING! Solver return FLAG ',I3,'!',/)
   
     end subroutine dfoirfalg
  
positive.code:        ! Verify convergence conditions
  
         if ( hynorm .le. epsfeas .and. &
              dnorm .le. epsopt ) then
            flag = 0
            exit
         end if
  
negative.code: 
      allocate(linrhs(m), linpos(m  + 1), linvar(m * n), linval(m * n))
  
      hxnorm = evalinfeas(n,x,l,u,me,mi,flag)
  
      currfeas = hxnorm
  
      dzynorm = 1.0D+20
  
      call aevalf(n,x,fx,flag)
  
  