anchor.diff:    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) stop "comb provider"
  +  if(first_det_of_teeth(1) /= first_det_of_comb) then
  +     print *, 'Error in ', irp_here
  +     stop "comb provider"
  +  endif
     
   END_PROVIDER
   
  
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: 
      time = omp_get_wtime()
    
      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
  