anchor.diff:      if(time - timeLast > 1d0 .or. more /= 1) then
         timeLast = time
         cur_cp = N_cp
  -      if(.not. actually_computed(1)) cycle pullLoop
  +      if(.not. actually_computed(mrcc_jobs(1))) cycle pullLoop
   
         do i=2,N_det_generators
  -        if(.not. actually_computed(i)) then
  +        if(.not. actually_computed(mrcc_jobs(i))) then
  +          print *, "first not comp", i
             cur_cp = done_cp_at(i-1)
             exit
           end if
  
positive.code:   logical, allocatable         :: computed(:)
    integer                        :: i, j, last_full, dets(comb_teeth)
    integer                        :: k, l, cur_cp, under_det(comb_teeth+1)
    
    allocate(computed(N_det_generators))
    cps = 0d0
    cur_cp = 1
    done_cp_at = 0
  
negative.code:         E0 = E0 + mrcc_detail(1, first_det_of_teeth(cp_first_tooth(cur_cp))) * (1d0-fractage(cp_first_tooth(cur_cp)))
        end if
        call wall_time(time)
        !if ((dabs(eqt) < relative_error*0d0 .and. cps_N(cur_cp) >= 30)  .or. total_computed == N_det_generators) then
        if(cur_cp > 10) then 
          ! Termination
          !print '(G10.3, 2X, F16.7, 2X, G16.3, 2X, F16.4, A20)', Nabove(tooth), avg+E, eqt, time-time0, ''
          print *, "GREPME", cur_cp, E+E0+avg, eqt, time-time0, total_computed
  