diff options
author | Willem Jan Palenstijn <wjp@usecode.org> | 2019-03-29 21:21:29 +0100 |
---|---|---|
committer | Willem Jan Palenstijn <Willem.Jan.Palenstijn@cwi.nl> | 2019-09-25 14:10:08 +0200 |
commit | 3cf63d335ebe392a8c77f0c90395c18150647aeb (patch) | |
tree | 60a96bfe45c62aa9e57c1c5751e72d3794934845 /cuda/2d/algo.cu | |
parent | 87715885f3b4a80693493e37aa8293899a6b987e (diff) | |
download | astra-3cf63d335ebe392a8c77f0c90395c18150647aeb.tar.gz astra-3cf63d335ebe392a8c77f0c90395c18150647aeb.tar.bz2 astra-3cf63d335ebe392a8c77f0c90395c18150647aeb.tar.xz astra-3cf63d335ebe392a8c77f0c90395c18150647aeb.zip |
Adjust adjoint to line integral scaling
Diffstat (limited to 'cuda/2d/algo.cu')
-rw-r--r-- | cuda/2d/algo.cu | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/cuda/2d/algo.cu b/cuda/2d/algo.cu index b4c2864..11422ff 100644 --- a/cuda/2d/algo.cu +++ b/cuda/2d/algo.cu @@ -258,10 +258,10 @@ bool ReconAlgo::copyDataToGPU(const float* pfSinogram, unsigned int iSinogramPit if (!ok) return false; - // rescale sinogram to adjust for pixel size - processSino<opMul>(D_sinoData, fSinogramScale, - //1.0f/(fPixelSize*fPixelSize), - sinoPitch, dims); + // rescale sinogram + if (fSinogramScale != 1.0f) + processSino<opMul>(D_sinoData, fSinogramScale, + sinoPitch, dims); ok = copyVolumeToDevice(pfReconstruction, iReconstructionPitch, dims, @@ -331,11 +331,11 @@ bool ReconAlgo::callBP(float* D_volumeData, unsigned int volumePitch, if (parProjs) { assert(!fanProjs); return BP(D_volumeData, volumePitch, D_projData, projPitch, - dims, parProjs, outputScale); + dims, parProjs, fOutputScale * outputScale); } else { assert(fanProjs); return FanBP(D_volumeData, volumePitch, D_projData, projPitch, - dims, fanProjs, outputScale); + dims, fanProjs, fOutputScale * outputScale); } } |