 
 
 float getBias(){
+	
 	float tmp_bias = 0.0;
-	if(int(N_2)%2 == 0){
-		tmp_bias = (N_2/2.0)/(N_2-1.0);
-	}else{
-		tmp_bias = 0.5; // old: 50.0/99.0;
-	}
+	// Set coordinates for the Origin
+	if (int(dimN) % 2 == 0)
+		tmp_bias = float(dimN )/ 2.0f  ; // -2 dur to rotational lochay
+	else
+		tmp_bias = float(dimN - 1.0) / 2.0f  ;
+		
+
+	
+	
+//	float tmp_bias = 0.0;
+//	if(int(N_2)%2 == 0){
+//		tmp_bias = (N_2/2.0)/(N_2-1.0);
+//	}else{
+//		tmp_bias = 0.5; // old: 50.0/99.0;
+//	}
 	return tmp_bias;
 }
 float bias = getBias();

	float reHeight = 0.0; float imHeight = 0.0;
	float real_part = 0.0; float imag_part = 0.0;
	float fourier_coefficients = 1.0;
	vec2 sum = vec2(0);
	
	// approximation till iteration 30 of fourier coefficient
	for(int n = lower; n <= upper; n++){

	distU = distU * distU / varX_InTxtUnits;
	distV = distV * distV / varY_InTxtUnits;
	
	return exp((-distU - distV)/2.0f)/sigma_f_pix;

}

vec2 taylorGaussWindow(vec2 coords, float k, float w){
