diff options
Diffstat (limited to 'cuda/3d/cone_bp.cu')
| -rw-r--r-- | cuda/3d/cone_bp.cu | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/cuda/3d/cone_bp.cu b/cuda/3d/cone_bp.cu index df1e19f..3525eb4 100644 --- a/cuda/3d/cone_bp.cu +++ b/cuda/3d/cone_bp.cu @@ -213,15 +213,15 @@ __global__ void dev_cone_BP_SS(void* D_volData, unsigned int volPitch, int start float fZs = fZ; for (int iSubZ = 0; iSubZ < iRaysPerVoxelDim; ++iSubZ) { - const float fUNum = fCu.w + fX * fCu.x + fY * fCu.y + fZ * fCu.z; - const float fVNum = fCv.w + fX * fCv.x + fY * fCv.y + fZ * fCv.z; - const float fDen = fCd.w + fX * fCd.x + fY * fCd.y + fZ * fCd.z; + const float fUNum = fCu.w + fXs * fCu.x + fYs * fCu.y + fZs * fCu.z; + const float fVNum = fCv.w + fXs * fCv.x + fYs * fCv.y + fZs * fCv.z; + const float fDen = fCd.w + fXs * fCd.x + fYs * fCd.y + fZs * fCd.z; const float fr = __fdividef(1.0f, fDen); const float fU = fUNum * fr; const float fV = fVNum * fr; - fVal += tex3D(gT_coneProjTexture, fU, fV, fAngle) * fr; + fVal += tex3D(gT_coneProjTexture, fU, fAngle, fV) * fr * fr; fZs += fSubStep; } |