 		
 	if(isnan(brdf.x) ||isnan(brdf.y) ||isnan(brdf.z)) o_col = vec4(1.0, 0.0, 0.0, 1.0);
 	else if(isinf(brdf.x) ||isinf(brdf.y) ||isinf(brdf.z)) o_col = vec4(0.0, 0.0, 1.0, 1.0);
-//	else o_col = brdf+vec4(ambient,ambient,ambient,0.0);
-	else o_col = vec4(ambient,ambient,ambient,0.0);
+	else o_col = brdf+vec4(ambient,ambient,ambient,0.0);
+//	else o_col = vec4(ambient,ambient,ambient,0.0);
 	frag_shaded	= o_col;
 }

				t1 = v;
				t2 = u;
			}
				
			for(float iter = lower; iter <= upper; iter++){
				if(iter == 0.0) continue;
				lambda_iter = (dx*t1)/iter;
				k = 2.0*PI / lambda_iter;

					
						P = taylorApproximation(coords, k, w);

						

														
						float norm_fact = sigma_f_pix*PI;
							
