     for(int i=0; i<N; i++)
         printf("%lf ", numbering[i]);
     printf("\n");
+    printf("%d\n", numbering_copy.size());
 
     cudaFree(numbering_gpu);
     cudaFree(mask);

        flag = 0;
        delta /= 8;
        int oldsize = 0;
        print_array_global<<< 1, N >>>(numbering_gpu);
        cudaerr = cudaDeviceSynchronize();
        if (cudaerr != cudaSuccess)
            printf("PRINT kernel launch failed with error \"%s\".\n",

    cudaDeviceSynchronize();
    parallel_prefix(is_richer_neighbor, irn_sum, n);
    cudaDeviceSynchronize();

    if(irn_sum[n] == 0)
        stratify_none(numbering, is_class_component, indptr, indices, delta, n, icc_sum[n]);
    else{
