ELF>€@@ #line 2 "texture_noise_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. * ***************************************************************************/ //------------------------------------------------------------------------------ // Texture utility functions //------------------------------------------------------------------------------ // Perlin Noise Data #define NOISE_PERM_SIZE 256 __constant int NoisePerm[2 * NOISE_PERM_SIZE] = { 151, 160, 137, 91, 90, 15, 131, 13, 201, 95, 96, 53, 194, 233, 7, 225, 140, 36, 103, 30, 69, 142, // Rest of noise permutation table 8, 99, 37, 240, 21, 10, 23, 190, 6, 148, 247, 120, 234, 75, 0, 26, 197, 62, 94, 252, 219, 203, 117, 35, 11, 32, 57, 177, 33, 88, 237, 149, 56, 87, 174, 20, 125, 136, 171, 168, 68, 175, 74, 165, 71, 134, 139, 48, 27, 166, 77, 146, 158, 231, 83, 111, 229, 122, 60, 211, 133, 230, 220, 105, 92, 41, 55, 46, 245, 40, 244, 102, 143, 54, 65, 25, 63, 161, 1, 216, 80, 73, 209, 76, 132, 187, 208, 89, 18, 169, 200, 196, 135, 130, 116, 188, 159, 86, 164, 100, 109, 198, 173, 186, 3, 64, 52, 217, 226, 250, 124, 123, 5, 202, 38, 147, 118, 126, 255, 82, 85, 212, 207, 206, 59, 227, 47, 16, 58, 17, 182, 189, 28, 42, 223, 183, 170, 213, 119, 248, 152, 2, 44, 154, 163, 70, 221, 153, 101, 155, 167, 43, 172, 9, 129, 22, 39, 253, 19, 98, 108, 110, 79, 113, 224, 232, 178, 185, 112, 104, 218, 246, 97, 228, 251, 34, 242, 193, 238, 210, 144, 12, 191, 179, 162, 241, 81, 51, 145, 235, 249, 14, 239, 107, 49, 192, 214, 31, 181, 199, 106, 157, 184, 84, 204, 176, 115, 121, 50, 45, 127, 4, 150, 254, 138, 236, 205, 93, 222, 114, 67, 29, 24, 72, 243, 141, 128, 195, 78, 66, 215, 61, 156, 180, 151, 160, 137, 91, 90, 15, 131, 13, 201, 95, 96, 53, 194, 233, 7, 225, 140, 36, 103, 30, 69, 142, 8, 99, 37, 240, 21, 10, 23, 190, 6, 148, 247, 120, 234, 75, 0, 26, 197, 62, 94, 252, 219, 203, 117, 35, 11, 32, 57, 177, 33, 88, 237, 149, 56, 87, 174, 20, 125, 136, 171, 168, 68, 175, 74, 165, 71, 134, 139, 48, 27, 166, 77, 146, 158, 231, 83, 111, 229, 122, 60, 211, 133, 230, 220, 105, 92, 41, 55, 46, 245, 40, 244, 102, 143, 54, 65, 25, 63, 161, 1, 216, 80, 73, 209, 76, 132, 187, 208, 89, 18, 169, 200, 196, 135, 130, 116, 188, 159, 86, 164, 100, 109, 198, 173, 186, 3, 64, 52, 217, 226, 250, 124, 123, 5, 202, 38, 147, 118, 126, 255, 82, 85, 212, 207, 206, 59, 227, 47, 16, 58, 17, 182, 189, 28, 42, 223, 183, 170, 213, 119, 248, 152, 2, 44, 154, 163, 70, 221, 153, 101, 155, 167, 43, 172, 9, 129, 22, 39, 253, 19, 98, 108, 110, 79, 113, 224, 232, 178, 185, 112, 104, 218, 246, 97, 228, 251, 34, 242, 193, 238, 210, 144, 12, 191, 179, 162, 241, 81, 51, 145, 235, 249, 14, 239, 107, 49, 192, 214, 31, 181, 199, 106, 157, 184, 84, 204, 176, 115, 121, 50, 45, 127, 4, 150, 254, 138, 236, 205, 93, 222, 114, 67, 29, 24, 72, 243, 141, 128, 195, 78, 66, 215, 61, 156, 180 }; float Grad(int x, int y, int z, float dx, float dy, float dz) { const int h = NoisePerm[NoisePerm[NoisePerm[x] + y] + z] & 15; const float u = h < 8 || h == 12 || h == 13 ? dx : dy; const float v = h < 4 || h == 12 || h == 13 ? dy : dz; return ((h & 1) ? -u : u) + ((h & 2) ? -v : v); } float NoiseWeight(float t) { const float t3 = t * t * t; const float t4 = t3 * t; return 6.f * t4 * t - 15.f * t4 + 10.f * t3; } float Noise(float x, float y, float z) { // Compute noise cell coordinates and offsets int ix = Floor2Int(x); int iy = Floor2Int(y); int iz = Floor2Int(z); const float dx = x - ix, dy = y - iy, dz = z - iz; // Compute gradient weights ix &= (NOISE_PERM_SIZE - 1); iy &= (NOISE_PERM_SIZE - 1); iz &= (NOISE_PERM_SIZE - 1); const float w000 = Grad(ix, iy, iz, dx, dy, dz); const float w100 = Grad(ix + 1, iy, iz, dx - 1, dy, dz); const float w010 = Grad(ix, iy + 1, iz, dx, dy - 1, dz); const float w110 = Grad(ix + 1, iy + 1, iz, dx - 1, dy - 1, dz); const float w001 = Grad(ix, iy, iz + 1, dx, dy, dz - 1); const float w101 = Grad(ix + 1, iy, iz + 1, dx - 1, dy, dz - 1); const float w011 = Grad(ix, iy + 1, iz + 1, dx, dy - 1, dz - 1); const float w111 = Grad(ix + 1, iy + 1, iz + 1, dx - 1, dy - 1, dz - 1); // Compute trilinear interpolation of weights const float wx = NoiseWeight(dx); const float wy = NoiseWeight(dy); const float wz = NoiseWeight(dz); const float x00 = Lerp(wx, w000, w100); const float x10 = Lerp(wx, w010, w110); const float x01 = Lerp(wx, w001, w101); const float x11 = Lerp(wx, w011, w111); const float y0 = Lerp(wy, x00, x10); const float y1 = Lerp(wy, x01, x11); return Lerp(wz, y0, y1); } float Noise3(const float3 P) { return Noise(P.x, P.y, P.z); } float FBm(const float3 P, const float omega, const int maxOctaves) { // Compute number of octaves for anti-aliased FBm const float foctaves = (float)maxOctaves; const int octaves = Floor2Int(foctaves); // Compute sum of octaves of noise for FBm float sum = 0.f, lambda = 1.f, o = 1.f; for (int i = 0; i < octaves; ++i) { sum += o * Noise3(lambda * P); lambda *= 1.99f; o *= omega; } const float partialOctave = foctaves - (float)octaves; sum += o * SmoothStep(.3f, .7f, partialOctave) * Noise3(lambda * P); return sum; } float Turbulence(const float3 P, const float omega, const int maxOctaves) { // Compute number of octaves for anti-aliased FBm const float foctaves = (float)maxOctaves; const int octaves = Floor2Int(foctaves); // Compute sum of octaves of noise for turbulence float sum = 0.f, lambda = 1.f, o = 1.f; for (int i = 0; i < octaves; ++i) { sum += o * fabs(Noise3(lambda * P)); lambda *= 1.99f; o *= omega; } const float partialOctave = foctaves - (float)(octaves); sum += o * SmoothStep(.3f, .7f, partialOctave) * fabs(Noise3(lambda * P)); // finally, add in value to account for average value of fabsf(Noise()) // (~0.2) for the remaining octaves... sum += (maxOctaves - foctaves) * 0.2f; return sum; } float tex_sin(float a) { a = 0.5f + 0.5f * sin(a); return a; } float tex_saw(float a) { const float b = 2.f * M_PI_F; int n = (int) (a / b); a -= n*b; if (a < 0.f) a += b; return a / b; } float tex_tri(float a) { const float b = 2.f * M_PI_F; const float rmax = 1.f; a = rmax - 2.f * fabs(floor((a * (1.f / b)) + .5f) - (a * (1.f / b))); return a; } SH5HƒìH‹HT$H‰ßèH‹=H‰ÞHèHƒÄ[ÃGCC: (GNU) 4.7.1zRx <AƒK nA.symtab.strtab.shstrtab.text.data.bss.rodata.str1.8.rela.text.startup.rela.ctors.comment.note.GNU-stack.rela.eh_frame@!@'@,2@@`< ;%S N¨%Z0¨cºxÀ8sÀ% ø‚€"È  H$Ïñÿ <  N Si˜«µÂtexture_noise_funcs_kernel.cpp_GLOBAL__sub_I_texture_noise_funcs_kernel.cpp.LC0_GLOBAL_OFFSET_TABLE__ZN3slg3ocl32KernelSource_texture_noise_funcsE_ZNSsC1EPKcRKSaIcE_ZNSsD1Ev__dso_handle__cxa_atexit üÿÿÿÿÿÿÿ üÿÿÿÿÿÿÿüÿÿÿÿÿÿÿ# üÿÿÿÿÿÿÿ-üÿÿÿÿÿÿÿ2üÿÿÿÿÿÿÿ