anchor.diff:  !                 endif
   !               end do
   !             else
  -              call i_H_psi(exc_det,microlist(1,1,ptr_futur_microlist(sporb)),psi_coef_microlist(ptr_futur_microlist(sporb), 1),N_int,N_futur_microlist(sporb),psi_selectors_size*4,N_states,i_H_psi_value2)
  +              call i_H_psi(exc_det,microlist(1,1,ptr_microlist(sporb)),psi_coef_microlist(ptr_microlist(sporb), 1),N_int,N_microlist(sporb)-N_futur_microlist(sporb),psi_selectors_size*4,N_states,i_H_psi_value2)
                 i_H_psi_value = i_H_psi_value + i_H_psi_value2
   !               !$OMP ATOMIC
   !               d2 += dabs(i_H_psi_value2(1))
  
positive.code: !               
            nok = .false.
            sporb = i_particle + (ispin - 1) * mo_tot_num                                              
  ! ! !           
            if(N_microlist(sporb) > 0) call check_past(exc_det, microlist(1,1,ptr_microlist(sporb)), idx_microlist(ptr_microlist(sporb)), N_microlist(sporb), i_generator, nok, N_int)
            if(nok) cycle
  !         
            if(N_microlist(0) > 0) call i_H_psi(exc_det,microlist,psi_coef_microlist,N_int,N_microlist(0),psi_selectors_size*3,N_states,i_H_psi_value)
            if(N_microlist(sporb) > 0) call i_H_psi(exc_det,microlist(1,1,ptr_microlist(sporb)),psi_coef_microlist(ptr_microlist(sporb), 1),N_int,N_microlist(sporb),psi_selectors_size*3,N_states,i_H_psi_value2)
            i_H_psi_value(:) = i_H_psi_value(:) + i_H_psi_value2(:)
            double precision :: Hii, diag_H_mat_elem_fock
            Hii = diag_H_mat_elem_fock(psi_det_generators(1,1,i_generator),exc_det,fock_diag_tmp,N_int)
  
negative.code:         
          
          p2 = i_particle2 + (ispin2 - 1) * mo_tot_num
          if(N_tmicrolist(p2) > 0 .and. idx_tmicrolist(ptr_tmicrolist(p2+1)-1) > i_generator) cycle
          if(isinwf(p1, p2)) cycle
          exc_det = ion_det
          
  