 	if(brdf.y < 1e-5) brdf.y = 0.0;
 	if(brdf.z < 1e-5) brdf.z = 0.0;
 	
-	brdf =  brdf*1.0*gainF(_k1, _k2)*shadowF;
+	brdf =  brdf*0.1*gainF(_k1, _k2)*shadowF;
 	brdf.xyz = getBRDF_RGB_T_D65(M_Adobe_XRNew, brdf.xyz);
 	
 	



		// xyz value of color for current wavelength (regarding current wavenumber k).
//		vec2 coords = vec2((k*v/(Omega)) + bias/dimN, (k*u/(Omega)) + bias/dimN);
		
		P = taylorGaussWindow(coord22, k, w);
		

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

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