ELF>`"@@ #line 2 "mc_funcs.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. * ***************************************************************************/ void ConcentricSampleDisk(const float u0, const float u1, float *dx, float *dy) { float r, theta; // Map uniform random numbers to $[-1,1]^2$ float sx = 2.f * u0 - 1.f; float sy = 2.f * u1 - 1.f; // Map square to $(r,\theta)$ // Handle degeneracy at the origin if (sx == 0.f && sy == 0.f) { *dx = 0.f; *dy = 0.f; return; } if (sx >= -sy) { if (sx > sy) { // Handle first region of disk r = sx; if (sy > 0.f) theta = sy / r; else theta = 8.f + sy / r; } else { // Handle second region of disk r = sy; theta = 2.f - sx / r; } } else { if (sx <= sy) { // Handle third region of disk r = -sx; theta = 4.f - sy / r; } else { // Handle fourth region of disk r = -sy; theta = 6.f + sx / r; } } theta *= M_PI_F / 4.f; *dx = r * cos(theta); *dy = r * sin(theta); } float3 CosineSampleHemisphere(const float u0, const float u1) { float x, y; ConcentricSampleDisk(u0, u1, &x, &y); const float z = sqrt(fmax(0.f, 1.f - x * x - y * y)); return (float3)(x, y, z); } float3 CosineSampleHemisphereWithPdf(const float u0, const float u1, float *pdfW) { float x, y; ConcentricSampleDisk(u0, u1, &x, &y); const float z = sqrt(fmax(0.f, 1.f - x * x - y * y)); *pdfW = z * M_1_PI_F; return (float3)(x, y, z); } float3 UniformSampleSphere(const float u1, const float u2) { float z = 1.f - 2.f * u1; float r = sqrt(max(0.f, 1.f - z * z)); float phi = 2.f * M_PI_F * u2; float x = r * cos(phi); float y = r * sin(phi); return (float3)(x, y, z); } float3 UniformSampleCone(const float u0, const float u1, const float costhetamax, const float3 x, const float3 y, const float3 z) { const float costheta = mix(costhetamax, 1.f, u0); const float sintheta = sqrt(1.f - costheta * costheta); const float phi = u1 * 2.f * M_PI_F; const float kx = cos(phi) * sintheta; const float ky = sin(phi) * sintheta; const float kz = costheta; return (float3)(kx * x.x + ky * y.x + kz * z.x, kx * x.y + ky * y.y + kz * z.y, kx * x.z + ky * y.z + kz * z.z); } float UniformConePdf(const float costhetamax) { return 1.f / (2.f * M_PI_F * (1.f - costhetamax)); } float PowerHeuristic(const float fPdf, const float gPdf) { const float f = fPdf; const float g = gPdf; return (f * f) / (f * f + g * g); } float PdfWtoA(const float pdfW, const float dist, const float cosThere) { return pdfW * fabs(cosThere) / (dist * dist); } float PdfAtoW(const float pdfA, const float dist, const float cosThere) { return pdfA * dist * dist / fabs(cosThere); } //------------------------------------------------------------------------------ // Distribution1D //------------------------------------------------------------------------------ // Implementation of std::upper_bound() __global float *std_upper_bound(__global float *first, __global float *last, const float val) { __global float *it; uint count = last - first; uint step; while (count > 0) { it = first; step = count / 2; it += step; if (!(val < *it)) { first = ++it; count -= step + 1; } else count = step; } return first; } //__global float *std_upper_bound(__global float *first, __global float *last, const float val) { // __global float *it = first; // // while ((it <= last) && (*it <= val)) { // it++; // } // // return it; //} float Distribution1D_SampleContinuous(__global float *distribution1D, const float u, float *pdf, uint *off) { const uint count = as_uint(distribution1D[0]); __global float *func = &distribution1D[1]; __global float *cdf = &distribution1D[1 + count]; // Find surrounding CDF segments and _offset_ if (u <= cdf[0]) { *pdf = func[0]; if (off) *off = 0; return 0.f; } if (u >= cdf[count]) { *pdf = func[count - 1]; if (off) *off = count - 1; return 1.f; } __global float *ptr = std_upper_bound(cdf, cdf + count + 1, u); const uint offset = ptr - cdf - 1; // Compute offset along CDF segment const float du = (u - cdf[offset]) / (cdf[offset + 1] - cdf[offset]); // Compute PDF for sampled offset *pdf = func[offset]; // Save offset if (off) *off = offset; return (offset + du) / count; } uint Distribution1D_SampleDiscrete(__global float *distribution1D, const float u, float *pdf) { const uint count = as_uint(distribution1D[0]); __global float *func = &distribution1D[1]; __global float *cdf = &distribution1D[1 + count]; // Find surrounding CDF segments and _offset_ if (u <= cdf[0]) { *pdf = func[0] / count; return 0; } if (u >= cdf[count]) { *pdf = func[count - 1] / count; return count - 1; } __global float *ptr = std_upper_bound(cdf, cdf + count + 1, u); const uint offset = ptr - cdf - 1; // Compute PDF for sampled offset *pdf = func[offset] / count; return offset; } uint Distribution1D_Offset(__global float *distribution1D, const float u) { const uint count = as_uint(distribution1D[0]); return min(count - 1, Floor2UInt(u * count)); } float Distribution1D_Pdf_UINT(__global float *distribution1D, const uint offset) { const uint count = as_uint(distribution1D[0]); __global float *func = &distribution1D[1]; return func[offset] / count; } float Distribution1D_Pdf_FLOAT(__global float *distribution1D, const float u) { const uint count = as_uint(distribution1D[0]); __global float *func = &distribution1D[1]; return func[Distribution1D_Offset(distribution1D, u)] / count; } //------------------------------------------------------------------------------ // Distribution2D //------------------------------------------------------------------------------ void Distribution2D_SampleContinuous(__global float *distribution2D, const float u0, const float u1, float2 *uv, float *pdf) { const uint width = as_uint(distribution2D[0]); const uint height = as_uint(distribution2D[1]); __global float *marginal = &distribution2D[2]; // size of Distribution1Dsize of count + size of func + size of cdf const uint marginalSize = 1 + height + height + 1; // size of Distribution1Dsize of count + size of func + size of cdf const uint conditionalSize = 1 + width + width + 1; float pdf1; uint index; (*uv).s1 = Distribution1D_SampleContinuous(marginal, u1, &pdf1, &index); float pdf0; __global float *conditional = &distribution2D[2 + marginalSize + index * conditionalSize]; (*uv).s0 = Distribution1D_SampleContinuous(conditional, u0, &pdf0, NULL); *pdf = pdf0 * pdf1; } float Distribution2D_Pdf(__global float *distribution2D, const float u0, const float u1) { const uint width = as_uint(distribution2D[0]); const uint height = as_uint(distribution2D[1]); __global float *marginal = &distribution2D[2]; // size of Distribution1Dsize of count + size of func + size of cdf const uint marginalSize = 1 + height + height + 1; // size of Distribution1Dsize of count + size of func + size of cdf const uint conditionalSize = 1 + width + width + 1; const uint index = Distribution1D_Offset(marginal, u1); __global float *conditional = &distribution2D[2 + marginalSize + index * conditionalSize]; return Distribution1D_Pdf_FLOAT(conditional, u0) * Distribution1D_Pdf_FLOAT(marginal, u1); } SH5HHHT$HH=HHH[GCC: (GNU) 4.7.1zRx <AK nA.symtab.strtab.shstrtab.text.data.bss.rodata.str1.8.rela.text.startup.rela.ctors.comment.note.GNU-stack.rela.eh_frame@!@'@,2@ @@!< ;(S!Nh)Z0!c!x!8s) !`&  ((<  8 =Swmc_funcs_kernel.cpp_GLOBAL__sub_I_mc_funcs_kernel.cpp.LC0_GLOBAL_OFFSET_TABLE__ZN3slg3ocl21KernelSource_mc_funcsE_ZNSsC1EPKcRKSaIcE_ZNSsD1Ev__dso_handle__cxa_atexit  # -2