- im_global = new ImageFile (im_filename, nx, ny);
- im_global->fileCreate();
-#endif
-
-#ifdef HAVE_MPI
- mpi_t1 = MPI::Wtime();
- proj_reconst (*im_local, rs_local, opt_filter, opt_filter_param,
- opt_interp, opt_interp_param, opt_backproj, opt_trace);
-
- mpi_t2 = MPI::Wtime();
- mpi_t = mpi_t2 - mpi_t1;
- mpiWorld.getComm().Reduce(&mpi_t, &mpi_t_g, 1, MPI::DOUBLE, MPI::MAX, 0);
- if (mpiWorld.getRank() == 0 && opt_verbose)
- printf("Time to reconstruct = %f, Max time = %f\n", mpi_t, mpi_t_g);
-#else
- proj_reconst (*im_global, rs_global, opt_filter, opt_filter_param,
- opt_interp, opt_interp_param, opt_backproj, opt_trace);
+ imGlobal = new ImageFile (im_filename, nx, ny);
+ imGlobal->fileCreate();