- rs_local = raysum_create (NULL, mpi_ct.local_work_units[mpi_ct.my_rank], mpi_ndet, TRUE);
- // rs_local->ndet = mpi_ndet;
- // rs_local->nview = mpi_nview;
+ if (opt_debug) {
+ fprintf(stdout, "Calc'd local work units process %d: nviews=%d, local_work_units=%d, start_work_units=%d\n",
+ mpi_ct.my_rank, mpi_nview, mpi_ct.local_work_units[mpi_ct.my_rank], mpi_ct.start_work_unit[mpi_ct.my_rank]);
+ MPI_Barrier(mpi_ct.comm);
+ }
+
+ rs_local = raysum_create (NULL, mpi_ct.local_work_units[mpi_ct.my_rank], mpi_ndet);
+ if (opt_debug) {
+ fprintf(stdout, "Created rs_local %lx for process %d: local views=%4d, local dets=%4d\n", (unsigned long int) rs_local, mpi_ct.my_rank, mpi_ct.local_work_units[mpi_ct.my_rank], mpi_ndet);
+ MPI_Barrier(mpi_ct.comm);
+ }
+
+ rs_local->ndet = mpi_ndet;
+ rs_local->nview = mpi_nview;