ELF          >                               @     @   #line 2 "materialdefs_funcs_velvet.cl"

/***************************************************************************
 * Copyright 1998-2013 by authors (see AUTHORS.txt)                        *
 *                                                                         *
 *   This file is part of LuxRender.                                       *
 *                                                                         *
 * Licensed under the Apache License, Version 2.0 (the "License");         *
 * you may not use this file except in compliance with the License.        *
 * You may obtain a copy of the License at                                 *
 *                                                                         *
 *     http://www.apache.org/licenses/LICENSE-2.0                          *
 *                                                                         *
 * Unless required by applicable law or agreed to in writing, software     *
 * distributed under the License is distributed on an "AS IS" BASIS,       *
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.*
 * See the License for the specific language governing permissions and     *
 * limitations under the License.                                          *
 ***************************************************************************/

//------------------------------------------------------------------------------
// Velvet material
//------------------------------------------------------------------------------

#if defined (PARAM_ENABLE_MAT_VELVET)

float3 VelvetMaterial_Evaluate(__global Material *material,
		__global HitPoint *hitPoint, const float3 lightDir, const float3 eyeDir,
		BSDFEvent *event, float *directPdfW
		TEXTURES_PARAM_DECL) {
		
	if (directPdfW)
		*directPdfW = fabs(lightDir.z * M_1_PI_F);

	*event = DIFFUSE | REFLECT;

	const float3 kd = Spectrum_Clamp(Texture_GetSpectrumValue(&texs[material->velvet.kdTexIndex], hitPoint
			TEXTURES_PARAM));

	const float A1 = Texture_GetFloatValue(&texs[material->velvet.p1TexIndex], hitPoint
			TEXTURES_PARAM);

	const float A2 = Texture_GetFloatValue(&texs[material->velvet.p2TexIndex], hitPoint
			TEXTURES_PARAM);

	const float A3 = Texture_GetFloatValue(&texs[material->velvet.p3TexIndex], hitPoint
			TEXTURES_PARAM);

	const float delta = Texture_GetFloatValue(&texs[material->velvet.thicknessTexIndex], hitPoint
			TEXTURES_PARAM);

	const float cosv = -dot(lightDir, eyeDir);

	// Compute phase function

	const float B = 3.0f * cosv;

	float p = 1.0f + A1 * cosv + A2 * 0.5f * (B * cosv - 1.0f) + A3 * 0.5 * (5.0f * cosv * cosv * cosv - B);
	p = p / (4.0f * M_PI_F);
 
	p = (p * delta) / fabs(eyeDir.z);

	// Clamp the BRDF (page 7)
	if (p > 1.0f)
		p = 1.0f;
	else if (p < 0.0f)
		p = 0.0f;

	return p * kd;
}

float3 VelvetMaterial_Sample(__global Material *material,
		__global HitPoint *hitPoint, const float3 fixedDir, float3 *sampledDir,
		const float u0, const float u1, 
		float *pdfW, float *cosSampledDir, BSDFEvent *event,
		const BSDFEvent requestedEvent
		TEXTURES_PARAM_DECL) {
	if (!(requestedEvent & (DIFFUSE | REFLECT)) ||
			(fabs(fixedDir.z) < DEFAULT_COS_EPSILON_STATIC))
		return BLACK;

	*sampledDir = (signbit(fixedDir.z) ? -1.f : 1.f) * CosineSampleHemisphereWithPdf(u0, u1, pdfW);

	*cosSampledDir = fabs((*sampledDir).z);
	if (*cosSampledDir < DEFAULT_COS_EPSILON_STATIC)
		return BLACK;

	*event = DIFFUSE | REFLECT;

	const float3 kd = Spectrum_Clamp(Texture_GetSpectrumValue(&texs[material->velvet.kdTexIndex], hitPoint
			TEXTURES_PARAM));

	const float A1 = Texture_GetFloatValue(&texs[material->velvet.p1TexIndex], hitPoint
			TEXTURES_PARAM);

	const float A2 = Texture_GetFloatValue(&texs[material->velvet.p2TexIndex], hitPoint
			TEXTURES_PARAM);

	const float A3 = Texture_GetFloatValue(&texs[material->velvet.p3TexIndex], hitPoint
			TEXTURES_PARAM);

	const float delta = Texture_GetFloatValue(&texs[material->velvet.thicknessTexIndex], hitPoint
			TEXTURES_PARAM);

	const float cosv = dot(-fixedDir, *sampledDir);;

	// Compute phase function

	const float B = 3.0f * cosv;

	float p = 1.0f + A1 * cosv + A2 * 0.5f * (B * cosv - 1.0f) + A3 * 0.5 * (5.0f * cosv * cosv * cosv - B);
	p = p / (4.0f * M_PI_F);
 
	p = (p * delta) / fabs(fixedDir.z);

	// Clamp the BRDF (page 7)
	if (p > 1.0f)
		p = 1.0f;
	else if (p < 0.0f)
		p = 0.0f;

	return kd * (p / *pdfW);
}
#endif
           SH5    HƒìH‹    HT$H‰ßè    H‹=    H‰ÞH    è    HƒÄ[Ã             GCC: (GNU) 4.7.1              zR x            <    AƒK nA  .symtab .strtab .shstrtab .text .data .bss .rodata.str1.8 .rela.text.startup .rela.ctors .comment .note.GNU-stack .rela.eh_frame                                                                                            @                                      !                     @                                      '                     @                                     ,      2               @       6                            @                     €      <                               ;                      P                                 S                     À                                    N                      à                                 Z      0               È                                   c                      Ú                                     x                     à      8                              s                      ø                                                             ‚                                                           È                          	                      h      á                                                          ñÿ                                                                                                                              &             <                                  
                                       Z                           	                 _                      u                    ª                      ½                      Ç                     Ô                       materialdefs_funcs_velvet_kernel.cpp _GLOBAL__sub_I_materialdefs_funcs_velvet_kernel.cpp .LC0 _GLOBAL_OFFSET_TABLE_ _ZN3slg3ocl38KernelSource_materialdefs_funcs_velvetE _ZNSsC1EPKcRKSaIcE _ZNSsD1Ev __dso_handle __cxa_atexit                     üÿÿÿÿÿÿÿ       	      üÿÿÿÿÿÿÿ             üÿÿÿÿÿÿÿ#       	      üÿÿÿÿÿÿÿ-             üÿÿÿÿÿÿÿ2             üÿÿÿÿÿÿÿ                                            