anchor.diff:      
       ! 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
  +    do i=1,N_det_beta_unique
         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
  +      if (l_b > N_det) exit
       enddo
       i = i-1
     
  
positive.code:     
      ! Loop inside the beta column to gather all the connected alphas
      l_a = k_a+1
      nmax = min(N_det_alpha_unique, N_det - l_a+1)
      do i=1,nmax
        lcol = psi_bilinear_matrix_columns(l_a)
        if (lcol /= kcol) exit
        lrow = psi_bilinear_matrix_rows(l_a)
        buffer(1:$N_int,i) = psi_det_alpha_unique(1:$N_int, lrow)
        idx(i) = l_a
        l_a = l_a+1
      enddo
      i = i-1
      
  
negative.code:       !$OMP          istart, iend, istep,        &
        !$OMP          ishift, idx0, u_t, maxab, v_0, s_0)             &
        !$OMP   PRIVATE(krow, kcol, tmp_det, spindet, k_a, k_b, i,     &
        !$OMP          lcol, lrow, l_a, l_b,                           &
        !$OMP          buffer, doubles, n_doubles, &
        !$OMP          tmp_det2, hij, sij, idx, l, kcol_prev, v_t,     &
        !$OMP          singles_a, n_singles_a, singles_b,              &
  