anchor.diff:        pt2_detail(:, index(i)) += pt2_mwen(:,i)
         parts_to_get(index(i)) -= 1
         if(parts_to_get(index(i)) < 0) then 
  +        print *, "PARTS ??"
           stop "PARTS ??"
         end if
         if(parts_to_get(index(i)) == 0) actually_computed(index(i)) = .true.
  
positive.code:   
    comb_step = 1d0/dfloat(comb_teeth)
    do i=1,N_det_generators
      if(pt2_weight(i)/norm_left < comb_step/2d0) then
        first_det_of_comb = i
        exit
      end if
  
negative.code:   end do
    first_det_of_teeth(comb_teeth+1) = N_det_generators + 1
    first_det_of_teeth(1) = first_det_of_comb
    if(first_det_of_teeth(1) /= first_det_of_comb) then
       print *, 'Error in ', irp_here
       stop "comb provider"
    endif
    
  END_PROVIDER
  
  