 	float comp_sigma = sigma_f_pix;
 	sigma_f_pix *= sigma_f_pix;
 	sigma_f_pix *= 2.0;
-	
+	t1 = u;
+	t2 = v;
 	
 	for(float iter = 0; iter < iterMax; iter = iter + 1.0){
 		

		
		float uv_N_n_hat = (kk*dx*t2);
		float uv_N_n = floor(uv_N_n_hat);
		float uu_N_n = (kk*dx*t1);
			
		float uu_N_base = uu_N_n - neighborRadius;
		float uv_N_base = uv_N_n - neighborRadius;
		

				if(coords.x < 0.0 || coords.x > 1.0 || coords.y < 0.0 || coords.y > 1.0) continue;
				
				// complex valued frequency contribution of current pixel.
				P = taylorApproximation(coords, kk, w);
				
				// compute gaussion weight of current pixel
				float w_ij = getGaussianWeight(dist2, sigma_f_pix);
