anchor.diff:    logical, allocatable :: actually_computed(:)
     integer :: total_computed
     
  -  allocate(delta_det(N_states, N_det_non_ref, 0:comb_teeth, 2))
  +  allocate(delta_det(N_states, N_det_non_ref, 0:comb_teeth+1, 2))
     allocate(cp(N_states, N_det_non_ref, N_cp, 2), mrcc_detail(N_states, N_det_generators))
     allocate(delta_loc(N_states, N_det_non_ref, 2))
     mrcc_detail = 0d0
  
positive.code:         if(cps(ind(i), j) > 0d0) then
            if(tooth_of_det(ind(i)) < cp_first_tooth(j)) stop "coef on supposedely deterministic det"
            double precision :: fac
            fac = cps(ind(i), j) / cps_N(j) * mrcc_weight_inv(ind(i)) * comb_step
            do k=1,N_det_non_ref
            do i_state=1,N_states
  
negative.code:     integer :: iloc
      iloc = mrcc_find(stato, mrcc_cweight, N_det_generators, 1, iloc)
      first_det_of_teeth(i) = iloc
      fractage(i) = (mrcc_cweight(iloc) - stato) / mrcc_weight(iloc)
      stato -= comb_step
    end do
    first_det_of_teeth(comb_teeth+1) = N_det_generators + 1
  