X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=libctsim%2Fbackprojectors.cpp;h=695582aa86d493c6c8ba85615f9c4989bf3acc7f;hb=0219ede69e1c3afc6e160b8f276bfd4617acbc08;hp=2b9eba45cdc4ec7ce4b300030c19a4fafaf5af1e;hpb=352a0691b9bd67f6c93ea822b353be3c101f4adb;p=ctsim.git

diff --git a/libctsim/backprojectors.cpp b/libctsim/backprojectors.cpp
index 2b9eba4..695582a 100644
--- a/libctsim/backprojectors.cpp
+++ b/libctsim/backprojectors.cpp
@@ -330,6 +330,9 @@ Backproject::PostProcessing()
 void
 Backproject::ScaleImageByRotIncrement ()
 {
+#if HAVE_OPENMP
+  #pragma omp parallel for
+#endif
   for (int ix = 0; ix < nx; ix++)
     for (int iy = 0; iy < ny; iy++)
       v[ix][iy] *= rotScale;
@@ -548,6 +551,9 @@ BackprojectDiff::BackprojectView (const double* const filteredProj, const double
   if (interpType == Backprojector::INTERP_LINEAR) {
     // precalculate scaled difference for linear interpolation
     deltaFilteredProj = new double [nDet];
+#if HAVE_OPENMP
+  #pragma omp parallel for
+#endif
     for (int i = 0; i < nDet - 1; i++)
       deltaFilteredProj[i] = filteredProj[i+1] - filteredProj[i];
     deltaFilteredProj[nDet - 1] = 0;  // last detector
@@ -624,6 +630,9 @@ BackprojectIntDiff::BackprojectView (const double* const filteredProj, const dou
   if (interpType == Backprojector::INTERP_LINEAR) {
     // precalculate scaled difference for linear interpolation
     deltaFilteredProj = new double [nDet];
+#if HAVE_OPENMP
+    #pragma omp parallel for
+#endif
     for (int i = 0; i < nDet - 1; i++)
       deltaFilteredProj[i] = (filteredProj[i+1] - filteredProj[i]) * dInvScale;
     deltaFilteredProj[nDet - 1] = 0;  // last detector