anchor.diff:  
       time = omp_get_wtime()
     
  -    if(time - timeLast > 30.0 .or. more /= 1) then
  +    if(time - timeLast > 10.0 .or. more /= 1) then
         timeLast = time
         do i=1, first_det_of_teeth(1)-1
           if(not(actually_computed(i))) then
  
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
  
  