- if (iDetPos < 0 || iDetPos >= nDet - 1)
- ; // errorIndexOutsideDetector (ix, iy, beta, r[ix][iy], phi[ix][iy], dDetPos, iDetPos);
- else
- pImCol[iy] += (filteredProj[iDetPos] + frac * (filteredProj[iDetPos+1] - filteredProj[iDetPos])) / (dU * dU);
+ if (iDetPos >= 0 && iDetPos < nDet - 1)
+ pImCol[iy] += (filteredProj[iDetPos] + frac * (filteredProj[iDetPos+1] - filteredProj[iDetPos]))
+ / (dU * dU);
+ } else if (interpType == Backprojector::INTERP_CUBIC) {
+ double d = iDetCenter + dPos; // position along detector
+ if (d >= 0 && d < nDet)
+ pImCol[iy] += pCubicInterp->interpolate (d) / (dU * dU);