anchor.diff:  !
   ! Start of main SCF loop
   !
  -  do while((max_error_DIIS > threshold_DIIS_nonzero) .and. (iteration_SCF < n_it_SCF_max) .and. dabs(Delta_energy_SCF) > thresh_SCF)
  +  do while(( (max_error_DIIS > threshold_DIIS_nonzero).or.(dabs(Delta_energy_SCF) > thresh_SCF) ) .and. (iteration_SCF < n_it_SCF_max))
   
   ! Increment cycle number
   
  
positive.code: 
      energy_SCF = HF_energy
      Delta_Energy_SCF = energy_SCF - energy_SCF_previous
      energy_SCF_previous = energy_SCF
  
  !   Print results at the end of each iteration
  
negative.code:           )
  
        Fock_matrix_AO_alpha = Fock_matrix_AO*0.5d0
        Fock_matrix_AO_beta  = Fock_matrix_AO*0.5d0
        TOUCH Fock_matrix_AO_alpha Fock_matrix_AO_beta
  
      endif
  
      MO_coef = eigenvectors_Fock_matrix_MO
  
      TOUCH MO_coef
  
  !   Calculate error vectors
  
  