anchor.diff:  
       IF ( FLAG .NE. 0 ) GOTO 31
   
  -    IF ((F-FOPT) .GT. (0.1D0*VQUAD)) THEN
  -       DELTA= 0.5D0*DELTA                 
  -    ELSE IF ((F-FOPT) .GE. (0.7D0*VQUAD)) THEN
  -       DELTA=DELTA  
  +    IF ( (F - FOPT) .GT. (0.1D0 * VQUAD) ) THEN
  +       DELTA = 0.5D0 * DELTA                 
  +    ELSE IF ( (F - FOPT) .GE. (0.7D0 * VQUAD) ) THEN
  +       DELTA = DELTA  
       ELSE
  -       DELTA =   DELTA + DELTA  
  +       DELTA = 2.0D0 * DELTA
       END IF
   
       ! CHOOSE WHO LEAVE Y CALCULATING THE VALUE OF SIGMA. THE VARIABLE
  -    ! IT' IS CHOOSEN FOR DEFINE WHO LEAVE.
  -
  +    ! 'IT' IS CHOOSEN FOR DEFINE WHO LEAVE.
  +
  +    ! TODO: Take a close look at IT when there is no update
       CALL SIGMA(H,N,NPT,Y,X,VETOR1,SIGM,ALFA,BETA,TAU,IT,DELTA)         
   
       ! IF ANY REDUCTION IN F, PUT X IN INTERPOLATION SET.
       IF (F .LE. FOPT) THEN  
          IF ( OUTPUT ) WRITE(*,1005) IT
  -       DO I=1, N            
  +       DO I = 1, N            
             Y(IT,I) = X(I) 
          END DO
  +       FF(IT) = F
       ELSE
          GO TO 23
       END IF
   
       ! UPDATE H              
  -    CALL INVERSAH(H, N, NPT,VETOR1, SIGM, IT, ALFA, BETA,TAU)
  -
  -    CALL ATUALIZAQ(H, N, NPT, Q, DELTA, Y, X, F, IT) 
  -
  -23  IF (F  .LE.  FOPT + 0.1D0*VQUAD) THEN                 
  +
  +    CALL INVERSAH(H, N, NPT, VETOR1, SIGM, IT, ALFA, BETA,TAU)
  +
  +    ! Update SVR Model
  +
  +    call qsvm(Y, FF, n, npt, HQ, g, b)
  +
  +    ! Update TRDF's data
  +
  +    call svmToTRDF(n, HQ, g, b, XBASE_A, HQ_A, GOPT_A, MOD_A, Q)
  +
  +23  IF ( F .LE. FOPT + 0.1D0 * VQUAD ) THEN                 
          FOPT = F
  -       DO I=1, N
  +       DO I = 1, N
             XNOVO(I) = X(I) 
          END DO
          GO TO 11  
       END IF
  -    IF (sigm .le. 0d0) go to 4
  -    IF ((RHO .LE. RHOEND) .OR. (IC == MAXIC)) GO TO 31
  +
  +    IF ( sigm .le. 0.0d0 ) go to 4
  +
  +    IF ( (RHO .LE. RHOEND) .OR. (IC == MAXIC) ) GO TO 31
  +
       KN=0
  -    DO   K=1,NPT
  -       SUM=0D0
  -       DO   J=1,N
  -          SUM=SUM+(Y(K,J)-X(J))**2
  -       END DO
  -       IF (SUM .GT. DISTSQ) THEN
  -          KN=K
  +    DO K = 1,NPT
  +       SUM = 0.0D0
  +       DO J = 1,N
  +          SUM = SUM + (Y(K,J) - X(J))**2.0D0
  +       END DO
  +
  +       IF ( SUM .GT. DISTSQ ) THEN
  +          KN = K
             exit
          END IF
       END DO
   
  -    IF (KN .GT. 0)  GOTO 4       
  +    IF ( KN .GT. 0 )  GOTO 4       
   
       DELTA = RHO  
  -    RHO = GAMA * RHO      
  +    RHO   = GAMA * RHO      
   
       GO TO 11    
   
  
positive.code:     !
      ! FCNT - number of function evaluations
  
      IMPLICIT NONE
  
      ! SCALAR ARGUMENTS
  
negative.code:     evaljac => evaljac_
      evalhc  => evalhc_
  
      F      = 1.0D300
      IC     = 0
      MAXIC  = maxfcnt ! MAXIMUM NUMBER OF FUNCTION EVALUATIONS
      RHOBEG = RBEG
      RHO    = RHOBEG
      DELTA  = RHO
      RHOEND = REND
      GAMA   = 0.1D0
  
      !     ---------------------------
      !     Feasibility phase - Phase 0
  