anchor.diff:  
       do i=1,n_singles
         lcol = singles(i)
  -      ! TODO cycle if lcol <= kcol
  +      if (lcol <= kcol) cycle
  +
         tmp_det2(1:N_int,2) = psi_det_beta_unique(1:N_int, lcol)
   
         l_a = psi_bilinear_matrix_columns_loc(lcol)
  -      ! TODO loop
  +      do while (l_a <= k_a)
  +        l_a += 1
  +      enddo
   
         do while ( l_a < psi_bilinear_matrix_columns_loc(lcol+1) )
           lrow = psi_bilinear_matrix_rows(l_a)
  
positive.code: 
            call i_H_j_double_alpha_beta(tmp_det,tmp_det2,N_int,hij)
            v_0(k_a, 1:N_st) += hij * psi_bilinear_matrix_values(l_a,1:N_st)
          endif
          l_a += 1
  
  
negative.code: 
    enddo
  
  end
  
  
  