 	const float normF = 1.0f;
 
 	// These are frequency increments
-	int anchorX = int(floor(bias + coords.x * (dimN - 0)));
-	int anchorY = int(floor(bias + coords.y * (dimN - 0)));
+	int anchorX = int(floor(orgIm + coords.x * (dimN - 0)));
+	int anchorY = int(floor(orgIm + coords.y * (dimN - 0)));
 	
 	vec3 fftMag = vec3(0.0f);
 	

	float u = V.x; float v = V.y; float w = V.z;
	

	float iterMax = 300.0;
	float lambdaStep = (lambda_max - lambda_min)/(iterMax-1.0);
	float F2 = fFByR0*fFByR0;
	
	
	float stepSize = 1.0;
	sigma_f_pix = ((2.0*dx) / (PI*dimX));
	float comp_sigma = sigma_f_pix;
	sigma_f_pix *= sigma_f_pix;

	vec2 precomputedFourier = vec2(0.0, 0.0);
	vec2 susum = vec2(0.0, 0.0);
	int lower = 0; int upper = int(approxSteps)+1;
	upper = 10;
	float reHeight = 0.0; float imHeight = 0.0;
	float real_part = 0.0; float imag_part = 0.0;
	float fourier_coefficients = 1.0;
