anchor.diff:      ! Compute Hij for all alpha singles
       ! ----------------------------------
   
  -if (.False.) then
       tmp_det2(1:$N_int,2) = psi_det_beta_unique (1:$N_int, kcol)
       do i=1,n_singles_a
         l_a = singles_a(i)
  
positive.code:     
      ! Loop inside the alpha row to gather all the connected betas
      l_b = k_b+1
      nmax = min(N_det_beta_unique, N_det - l_b+1)
      do i=1,nmax
        lrow = psi_bilinear_matrix_transp_rows(l_b)
        if (lrow /= krow) exit
        lcol = psi_bilinear_matrix_transp_columns(l_b)
        buffer(1:$N_int,i) = psi_det_beta_unique(1:$N_int, lcol)
        idx(i) = l_b
        l_b = l_b+1
      enddo
      i = i-1
    
  
negative.code:   integer, allocatable           :: singles_a(:)
    integer, allocatable           :: singles_b(:)
    integer, allocatable           :: idx(:), idx0(:)
    integer                        :: maxab, n_singles_a, n_singles_b, kcol_prev
    integer*8                      :: k8
    double precision, allocatable  :: v_t(:,:), s_t(:,:)
    !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: v_t, s_t
  