 		float fourier_re = fourier_coefficients*precomputedFourier.x;
 		float fourier_im = fourier_coefficients*precomputedFourier.y;
 		
-		if(n % 4 == 0){
-			real_part += fourier_re;
-			imag_part += fourier_im;
-			
-		}else if(n % 4 == 1){
-			real_part -= fourier_im;
-			imag_part += fourier_re;
-			
-		}else if(n % 4 == 2){
-			real_part -= fourier_re;
-			imag_part -= fourier_im;
-			
-		}else{
-			real_part += fourier_im;
-			imag_part -= fourier_re;
-		}
+		
+		real_part += fourier_re;
+		imag_part += fourier_im;
+		
+//		if(n % 4 == 0){
+//			real_part += fourier_re;
+//			imag_part += fourier_im;
+//			
+//		}else if(n % 4 == 1){
+//			real_part -= fourier_im;
+//			imag_part += fourier_re;
+//			
+//		}else if(n % 4 == 2){
+//			real_part -= fourier_re;
+//			imag_part -= fourier_im;
+//			
+//		}else{
+//			real_part += fourier_im;
+//			imag_part -= fourier_re;
+//		}
 	}
 	
 	return vec2(real_part, imag_part);

	
	float uv_sqr = pow(u*u+v*v, 0.5);

	float iterMax = 200.0;
	float lambdaStep = (lambda_max - lambda_min)/(iterMax-1.0);
	float F2 = fFByR0*fFByR0;
	

float T_2 = t_0 * N_1;
float periods = periodCount-1.0; // 26 // number of patch periods along surface
float Omega = ((N_1/N_2)*2.0*PI)/t_0; // (N_1/N_2)*2*PI/t_0, before 8.0*PI*pow(10.0,7.0);
//float bias = (N_2/2.0)/(N_2-1.0); // old: 50.0/99.0;
float bias = (1.0/(N_1-1.0))*((N_1-1.0)/2.0); // old: 50.0/99.0;
//float bias = 0.5; // old: 50.0/99.0;


float neighborRadius = (neigh_rad < 5 && neigh_rad > -1) ? float(neigh_rad) : 0.0;

//transformation constant
