|
@@ -71,7 +71,7 @@ lamino_bp_varied (global float *proj,
|
|
|
const float newx = (float) vX - param->vol_ox;
|
|
|
|
|
|
/* Vary alpha depending on the z height */
|
|
|
- const float alpha = -3 * PI_2 + param->theta + newz / param->vol_sz * param->proj_ox_variation;
|
|
|
+ const float alpha = -3 * PI_2 + param->theta;
|
|
|
|
|
|
const float cf = cos(param->phi);
|
|
|
const float ct = cos(alpha);
|
|
@@ -87,7 +87,7 @@ lamino_bp_varied (global float *proj,
|
|
|
const float mat_4 = -sg * sf + cg * st * cf;
|
|
|
const float mat_5 = cg * ct;
|
|
|
|
|
|
- const float oldx = newx * mat_0 + newy * mat_1 + newz * mat_2 + param->proj_ox;
|
|
|
+ const float oldx = newx * mat_0 + newy * mat_1 + newz * mat_2 + param->proj_ox + newz / param->vol_sz * param->proj_ox_variation;
|
|
|
const float oldy = newx * mat_3 + newy * mat_4 + newz * mat_5 + param->proj_oy;
|
|
|
|
|
|
volume[idx] += interpolate (proj, param, oldx, oldy);
|