anchor.diff:        tmp_det2(1:N_int,1) = psi_det_alpha_unique(1:N_int, lrow)
         call i_H_j_mono_spin( tmp_det, tmp_det2, N_int, 1, hij)
         do l=1,N_st
  -        v_t(l,l_a) += hij * u_t(l,k_a)
  -        v_t(l,k_a) += hij * u_t(l,l_a)
  +        v_t(l,l_a) = v_t(l,l_a) + hij * u_t(l,k_a)
  +        v_t(l,k_a) = v_t(l,k_a) + hij * u_t(l,l_a)
           ! single => sij = 0 
         enddo
       enddo
  
positive.code:     tmp_det(1:N_int,1) = psi_det_alpha_unique(1:N_int, krow)
      tmp_det(1:N_int,2) = psi_det_beta_unique (1:N_int, kcol)
      
      is_single_a = .False.
      do k=1,singles_a(0,krow)
        is_single_a( singles_a(k,krow) ) = .True.
      enddo
  
negative.code:       enddo
        call i_H_j_double_spin( tmp_det(1,1), psi_det_alpha_unique(1, doubles(i)), N_int, hij)
        do l=1,N_st
          v_t(l,l_a) = v_t(l,l_a) + hij * u_t(l,k_a)
          v_t(l,k_a) = v_t(l,k_a) + hij * u_t(l,l_a)
          ! same spin => sij = 0 
        enddo
      enddo
  