anchor.diff:            
             integer :: c1, c2
             double precision :: hij
  -          c1 = ptr_futur_tmicrolist(p1)
  -          c2 = ptr_futur_tmicrolist(p2)
  +          c1 = ptr_tmicrolist(p1)
  +          c2 = ptr_tmicrolist(p2)
             do while(.true.)
  -            if(c1 >= ptr_tmicrolist(p1+1) .or. c2 >= ptr_tmicrolist(p2+1)) then
  -              if(ptr_tmicrolist(p1+1) /= c1) then
  -                call i_H_psi(exc_det,tmicrolist(1,1,c1),psi_coef_tmicrolist(c1, 1),N_int, ptr_tmicrolist(p1+1)-c1 ,psi_selectors_size*3,N_states,i_H_psi_value2)
  +            if(c1 >= ptr_futur_tmicrolist(p1) .or. c2 >= ptr_futur_tmicrolist(p2)) then
  +              if(ptr_futur_tmicrolist(p1) /= c1) then
  +                call i_H_psi(exc_det,tmicrolist(1,1,c1),psi_coef_tmicrolist(c1, 1),N_int, ptr_futur_tmicrolist(p1)-c1 ,psi_selectors_size*3,N_states,i_H_psi_value2)
                   i_H_psi_value = i_H_psi_value + i_H_psi_value2
   ! !                 !$OMP ATOMIC
   !                 d1 += dabs(i_H_psi_value2(1))
                 end if
                 
  -              if(ptr_tmicrolist(p2+1) /= c2) then
  -                call i_H_psi(exc_det,tmicrolist(1,1,c2),psi_coef_tmicrolist(c2, 1),N_int, ptr_tmicrolist(p2+1)-c2 ,psi_selectors_size*3,N_states,i_H_psi_value2)
  +              if(ptr_futur_tmicrolist(p2) /= c2) then
  +                call i_H_psi(exc_det,tmicrolist(1,1,c2),psi_coef_tmicrolist(c2, 1),N_int, ptr_futur_tmicrolist(p2)-c2 ,psi_selectors_size*3,N_states,i_H_psi_value2)
                   i_H_psi_value = i_H_psi_value + i_H_psi_value2
   !                 !$OMP ATOMIC
   !                 d1 += dabs(i_H_psi_value2(1))
  
positive.code:       
        
      if(nt <= 2) then
        if(i < i_cur) then
          N_microlist = 0 
          ptr_microlist = 1
          N_tmicrolist = 0 
  
negative.code:   do i=0,mo_tot_num*2
      ptr_futur_microlist(i) = ptr_microlist(i+1)
      do j=ptr_microlist(i), ptr_microlist(i+1) - 1
        if(idx_microlist(j) > i_generator) then
          ptr_futur_microlist(i) = j
          N_futur_microlist(i) = ptr_microlist(i+1) - j
          exit
  