anchor.diff:  
     ! Call the solver - Second run
   
  -  verbose = .false.
  -
     epsfeas = 1.0D-04
   
  +  epsopt  = 1.0D-01
  +
  +  call set_epsfeas(epsfeas)
  +
  +  PEN = 1.0D+05
  +
  +  call fird(n,x,l,u,me,mi,evalf,evalc,evaljac,verbose,ftype, &
  +       epsfeas,epsopt,f,feas,fcnt,flag)
  +
  +  ftcnt = ftcnt + fcnt
  +
  +  ! Call the solver - Third run
  +
  +  epsfeas = 1.0D-04
  +
     epsopt  = 1.0D-04
   
     call set_epsfeas(epsfeas)
   
  -  PEN = 1.0D+10
  -
  -  ftype = 2
  +  PEN = 1.0D+09
   
     call fird(n,x,l,u,me,mi,evalf,evalc,evaljac,verbose,ftype, &
          epsfeas,epsopt,f,feas,fcnt,flag)
  
positive.code: 
    close(99)
  
    ! Call the solver - First run
  
    verbose = .false.
  
    epsfeas = 1.0D-04
  
    epsopt  = 1.0D-01
  
negative.code:   call set_epsfeas(epsfeas)
  
    PEN = 1.0D+03
  
    call fird(n,x,l,u,me,mi,evalf,evalc,evaljac,verbose,ftype, &
         epsfeas,epsopt,f,feas,fcnt,flag)
  