1 // This file has been automatically generated from resource file src/Shaders/RaytraceBase.fs
3 static const char Shaders_RaytraceBase_fs[] =
4 "#ifdef ADAPTIVE_SAMPLING\n"
5 " #extension GL_ARB_shader_image_load_store : require\n"
6 " #extension GL_NV_shader_atomic_float : require\n"
9 "#ifdef USE_TEXTURES\n"
10 " #extension GL_ARB_bindless_texture : require\n"
13 "//! Normalized pixel coordinates.\n"
16 "//! Sub-pixel offset in X direction for FSAA.\n"
17 "uniform float uOffsetX = 0.f;\n"
18 "//! Sub-pixel offset in Y direction for FSAA.\n"
19 "uniform float uOffsetY = 0.f;\n"
21 "//! Origin of viewing ray in left-top corner.\n"
22 "uniform vec3 uOriginLT;\n"
23 "//! Origin of viewing ray in left-bottom corner.\n"
24 "uniform vec3 uOriginLB;\n"
25 "//! Origin of viewing ray in right-top corner.\n"
26 "uniform vec3 uOriginRT;\n"
27 "//! Origin of viewing ray in right-bottom corner.\n"
28 "uniform vec3 uOriginRB;\n"
30 "//! Width of the rendering window.\n"
31 "uniform int uWinSizeX;\n"
32 "//! Height of the rendering window.\n"
33 "uniform int uWinSizeY;\n"
35 "//! Direction of viewing ray in left-top corner.\n"
36 "uniform vec3 uDirectLT;\n"
37 "//! Direction of viewing ray in left-bottom corner.\n"
38 "uniform vec3 uDirectLB;\n"
39 "//! Direction of viewing ray in right-top corner.\n"
40 "uniform vec3 uDirectRT;\n"
41 "//! Direction of viewing ray in right-bottom corner.\n"
42 "uniform vec3 uDirectRB;\n"
44 "//! Inverse model-view-projection matrix.\n"
45 "uniform mat4 uUnviewMat;\n"
47 "//! Model-view-projection matrix.\n"
48 "uniform mat4 uViewMat;\n"
50 "//! Texture buffer of data records of bottom-level BVH nodes.\n"
51 "uniform isamplerBuffer uSceneNodeInfoTexture;\n"
52 "//! Texture buffer of minimum points of bottom-level BVH nodes.\n"
53 "uniform samplerBuffer uSceneMinPointTexture;\n"
54 "//! Texture buffer of maximum points of bottom-level BVH nodes.\n"
55 "uniform samplerBuffer uSceneMaxPointTexture;\n"
56 "//! Texture buffer of transformations of high-level BVH nodes.\n"
57 "uniform samplerBuffer uSceneTransformTexture;\n"
59 "//! Texture buffer of vertex coords.\n"
60 "uniform samplerBuffer uGeometryVertexTexture;\n"
61 "//! Texture buffer of vertex normals.\n"
62 "uniform samplerBuffer uGeometryNormalTexture;\n"
63 "#ifdef USE_TEXTURES\n"
64 " //! Texture buffer of per-vertex UV-coordinates.\n"
65 " uniform samplerBuffer uGeometryTexCrdTexture;\n"
67 "//! Texture buffer of triangle indices.\n"
68 "uniform isamplerBuffer uGeometryTriangTexture;\n"
70 "//! Texture buffer of material properties.\n"
71 "uniform samplerBuffer uRaytraceMaterialTexture;\n"
72 "//! Texture buffer of light source properties.\n"
73 "uniform samplerBuffer uRaytraceLightSrcTexture;\n"
74 "//! Environment map texture.\n"
75 "uniform sampler2D uEnvironmentMapTexture;\n"
77 "//! Total number of light sources.\n"
78 "uniform int uLightCount;\n"
79 "//! Intensity of global ambient light.\n"
80 "uniform vec4 uGlobalAmbient;\n"
82 "//! Enables/disables hard shadows.\n"
83 "uniform int uShadowsEnabled;\n"
84 "//! Enables/disables specular reflections.\n"
85 "uniform int uReflectEnabled;\n"
86 "//! Enables/disables spherical environment map.\n"
87 "uniform int uSphereMapEnabled;\n"
88 "//! Enables/disables environment map background.\n"
89 "uniform int uSphereMapForBack;\n"
91 "//! Radius of bounding sphere of the scene.\n"
92 "uniform float uSceneRadius;\n"
93 "//! Scene epsilon to prevent self-intersections.\n"
94 "uniform float uSceneEpsilon;\n"
96 "#ifdef USE_TEXTURES\n"
97 " //! Unique 64-bit handles of OpenGL textures.\n"
98 " uniform uvec2 uTextureSamplers[MAX_TEX_NUMBER];\n"
101 "#ifdef ADAPTIVE_SAMPLING\n"
102 " //! OpenGL image used for accumulating rendering result.\n"
103 " volatile restrict layout(size1x32) uniform image2D uRenderImage;\n"
105 " //! OpenGL image storing offsets of sampled pixels blocks.\n"
106 " coherent restrict layout(size2x32) uniform iimage2D uOffsetImage;\n"
109 "//! Top color of gradient background.\n"
110 "uniform vec4 uBackColorTop = vec4 (0.0);\n"
111 "//! Bottom color of gradient background.\n"
112 "uniform vec4 uBackColorBot = vec4 (0.0);\n"
114 "/////////////////////////////////////////////////////////////////////////////////////////\n"
115 "// Specific data types\n"
117 "//! Stores ray parameters.\n"
125 "//! Stores intersection parameters.\n"
126 "struct SIntersect\n"
135 "/////////////////////////////////////////////////////////////////////////////////////////\n"
136 "// Some useful constants\n"
138 "#define MAXFLOAT 1e15f\n"
140 "#define SMALL vec3 (exp2 (-80.0f))\n"
142 "#define ZERO vec3 (0.0f, 0.0f, 0.0f)\n"
143 "#define UNIT vec3 (1.0f, 1.0f, 1.0f)\n"
145 "#define AXIS_X vec3 (1.0f, 0.0f, 0.0f)\n"
146 "#define AXIS_Y vec3 (0.0f, 1.0f, 0.0f)\n"
147 "#define AXIS_Z vec3 (0.0f, 0.0f, 1.0f)\n"
149 "#define M_PI 3.14159265f\n"
151 "#define LUMA vec3 (0.2126f, 0.7152f, 0.0722f)\n"
153 "// =======================================================================\n"
154 "// function : MatrixRowMultiplyDir\n"
155 "// purpose : Multiplies a vector by matrix\n"
156 "// =======================================================================\n"
157 "vec3 MatrixRowMultiplyDir (in vec3 v,\n"
162 " return vec3 (dot (m0.xyz, v),\n"
163 " dot (m1.xyz, v),\n"
164 " dot (m2.xyz, v));\n"
167 "//! 32-bit state of random number generator.\n"
170 "// =======================================================================\n"
171 "// function : SeedRand\n"
172 "// purpose : Applies hash function by Thomas Wang to randomize seeds\n"
173 "// (see http://www.burtleburtle.net/bob/hash/integer.html)\n"
174 "// =======================================================================\n"
175 "void SeedRand (in int theSeed, in int theSizeX, in int theRadius)\n"
177 " RandState = uint (int (gl_FragCoord.y) / theRadius * theSizeX + int (gl_FragCoord.x) / theRadius + theSeed);\n"
179 " RandState = (RandState + 0x479ab41du) + (RandState << 8);\n"
180 " RandState = (RandState ^ 0xe4aa10ceu) ^ (RandState >> 5);\n"
181 " RandState = (RandState + 0x9942f0a6u) - (RandState << 14);\n"
182 " RandState = (RandState ^ 0x5aedd67du) ^ (RandState >> 3);\n"
183 " RandState = (RandState + 0x17bea992u) + (RandState << 7);\n"
186 "// =======================================================================\n"
187 "// function : RandInt\n"
188 "// purpose : Generates integer using Xorshift algorithm by G. Marsaglia\n"
189 "// =======================================================================\n"
192 " RandState ^= (RandState << 13);\n"
193 " RandState ^= (RandState >> 17);\n"
194 " RandState ^= (RandState << 5);\n"
196 " return RandState;\n"
199 "// =======================================================================\n"
200 "// function : RandFloat\n"
201 "// purpose : Generates a random float in 0 <= x < 1 range\n"
202 "// =======================================================================\n"
203 "float RandFloat()\n"
205 " return float (RandInt()) * (1.f / 4294967296.f);\n"
208 "// =======================================================================\n"
209 "// function : MatrixColMultiplyPnt\n"
210 "// purpose : Multiplies a vector by matrix\n"
211 "// =======================================================================\n"
212 "vec3 MatrixColMultiplyPnt (in vec3 v,\n"
218 " return vec3 (m0.x * v.x + m1.x * v.y + m2.x * v.z + m3.x,\n"
219 " m0.y * v.x + m1.y * v.y + m2.y * v.z + m3.y,\n"
220 " m0.z * v.x + m1.z * v.y + m2.z * v.z + m3.z);\n"
223 "// =======================================================================\n"
224 "// function : MatrixColMultiplyDir\n"
225 "// purpose : Multiplies a vector by matrix\n"
226 "// =======================================================================\n"
227 "vec3 MatrixColMultiplyDir (in vec3 v,\n"
232 " return vec3 (m0.x * v.x + m1.x * v.y + m2.x * v.z,\n"
233 " m0.y * v.x + m1.y * v.y + m2.y * v.z,\n"
234 " m0.z * v.x + m1.z * v.y + m2.z * v.z);\n"
237 "//=======================================================================\n"
238 "// function : InverseDirection\n"
239 "// purpose : Returns safely inverted direction of the given one\n"
240 "//=======================================================================\n"
241 "vec3 InverseDirection (in vec3 theInput)\n"
243 " vec3 anInverse = 1.f / max (abs (theInput), SMALL);\n"
245 " return mix (-anInverse, anInverse, step (ZERO, theInput));\n"
248 "//=======================================================================\n"
249 "// function : BackgroundColor\n"
250 "// purpose : Returns color of gradient background\n"
251 "//=======================================================================\n"
252 "vec4 BackgroundColor()\n"
254 "#ifdef ADAPTIVE_SAMPLING\n"
256 " ivec2 aFragCoord = ivec2 (gl_FragCoord.xy);\n"
258 " ivec2 aTileXY = imageLoad (uOffsetImage, ivec2 (aFragCoord.x / BLOCK_SIZE,\n"
259 " aFragCoord.y / BLOCK_SIZE)).xy;\n"
261 " aTileXY.y += aFragCoord.y % min (uWinSizeY - aTileXY.y, BLOCK_SIZE);\n"
263 " return mix (uBackColorBot, uBackColorTop, float (aTileXY.y) / uWinSizeY);\n"
267 " return mix (uBackColorBot, uBackColorTop, vPixel.y);\n"
272 "/////////////////////////////////////////////////////////////////////////////////////////\n"
273 "// Functions for compute ray-object intersection\n"
275 "// =======================================================================\n"
276 "// function : GenerateRay\n"
278 "// =======================================================================\n"
279 "SRay GenerateRay (in vec2 thePixel)\n"
281 " vec3 aP0 = mix (uOriginLB, uOriginRB, thePixel.x);\n"
282 " vec3 aP1 = mix (uOriginLT, uOriginRT, thePixel.x);\n"
284 " vec3 aD0 = mix (uDirectLB, uDirectRB, thePixel.x);\n"
285 " vec3 aD1 = mix (uDirectLT, uDirectRT, thePixel.x);\n"
287 " vec3 aDirection = normalize (mix (aD0, aD1, thePixel.y));\n"
289 " return SRay (mix (aP0, aP1, thePixel.y), aDirection);\n"
292 "// =======================================================================\n"
293 "// function : IntersectSphere\n"
294 "// purpose : Computes ray-sphere intersection\n"
295 "// =======================================================================\n"
296 "float IntersectSphere (in SRay theRay, in float theRadius)\n"
298 " float aDdotD = dot (theRay.Direct, theRay.Direct);\n"
299 " float aDdotO = dot (theRay.Direct, theRay.Origin);\n"
300 " float aOdotO = dot (theRay.Origin, theRay.Origin);\n"
302 " float aD = aDdotO * aDdotO - aDdotD * (aOdotO - theRadius * theRadius);\n"
306 " float aTime = (sqrt (aD) - aDdotO) * (1.0f / aDdotD);\n"
308 " return aTime > 0.0f ? aTime : MAXFLOAT;\n"
311 " return MAXFLOAT;\n"
314 "// =======================================================================\n"
315 "// function : IntersectTriangle\n"
316 "// purpose : Computes ray-triangle intersection (branchless version)\n"
317 "// =======================================================================\n"
318 "void IntersectTriangle (in SRay theRay,\n"
319 " in vec3 thePnt0,\n"
320 " in vec3 thePnt1,\n"
321 " in vec3 thePnt2,\n"
322 " out vec3 theUVT,\n"
323 " out vec3 theNorm)\n"
325 " vec3 aToTrg = thePnt0 - theRay.Origin;\n"
327 " vec3 aEdge0 = thePnt1 - thePnt0;\n"
328 " vec3 aEdge1 = thePnt0 - thePnt2;\n"
330 " theNorm = cross (aEdge1, aEdge0);\n"
332 " vec3 theVect = cross (theRay.Direct, aToTrg);\n"
334 " theUVT = vec3 (dot (theNorm, aToTrg),\n"
335 " dot (theVect, aEdge1),\n"
336 " dot (theVect, aEdge0)) * (1.f / dot (theNorm, theRay.Direct));\n"
338 " theUVT.x = any (lessThan (theUVT, ZERO)) || (theUVT.y + theUVT.z) > 1.f ? MAXFLOAT : theUVT.x;\n"
341 "#define EMPTY_ROOT ivec4(0)\n"
343 "//! Utility structure containing information about\n"
344 "//! currently traversing sub-tree of scene's BVH.\n"
347 " //! Transformed ray.\n"
350 " //! Inversed ray direction.\n"
353 " //! Parameters of sub-root node.\n"
357 "#define MATERIAL_AMBN(index) (18 * index + 0)\n"
358 "#define MATERIAL_DIFF(index) (18 * index + 1)\n"
359 "#define MATERIAL_SPEC(index) (18 * index + 2)\n"
360 "#define MATERIAL_EMIS(index) (18 * index + 3)\n"
361 "#define MATERIAL_REFL(index) (18 * index + 4)\n"
362 "#define MATERIAL_REFR(index) (18 * index + 5)\n"
363 "#define MATERIAL_TRAN(index) (18 * index + 6)\n"
364 "#define MATERIAL_TRS1(index) (18 * index + 7)\n"
365 "#define MATERIAL_TRS2(index) (18 * index + 8)\n"
366 "#define MATERIAL_TRS3(index) (18 * index + 9)\n"
368 "#define TRS_OFFSET(treelet) treelet.SubData.x\n"
369 "#define BVH_OFFSET(treelet) treelet.SubData.y\n"
370 "#define VRT_OFFSET(treelet) treelet.SubData.z\n"
371 "#define TRG_OFFSET(treelet) treelet.SubData.w\n"
373 "//! Identifies the absence of intersection.\n"
374 "#define INALID_HIT ivec4 (-1)\n"
376 "//! Global stack shared between traversal functions.\n"
377 "int Stack[STACK_SIZE];\n"
379 "// =======================================================================\n"
380 "// function : pop\n"
382 "// =======================================================================\n"
383 "int pop (inout int theHead)\n"
385 " int aData = Stack[theHead];\n"
387 " int aMask = aData >> 26;\n"
388 " int aNode = aMask & 0x3;\n"
392 " if ((aMask & 0x3) == aNode)\n"
398 " aMask |= (aMask << 2) & 0x30;\n"
400 " Stack[theHead] = (aData & 0x03FFFFFF) | (aMask << 26);\n"
403 " return (aData & 0x03FFFFFF) + aNode;\n"
406 "// =======================================================================\n"
407 "// function : SceneNearestHit\n"
408 "// purpose : Finds intersection with nearest scene triangle\n"
409 "// =======================================================================\n"
410 "ivec4 SceneNearestHit (in SRay theRay, in vec3 theInverse, inout SIntersect theHit, out int theTrsfId)\n"
412 " ivec4 aTriIndex = INALID_HIT;\n"
414 " int aNode = 0; // node to traverse\n"
415 " int aHead = -1; // pointer of stack\n"
416 " int aStop = -1; // BVH level switch\n"
418 " SSubTree aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
420 " for (bool toContinue = true; toContinue; /* none */)\n"
422 " ivec4 aData = texelFetch (uSceneNodeInfoTexture, aNode);\n"
424 " if (aData.x == 0) // if inner node\n"
426 " aData.y += BVH_OFFSET (aSubTree);\n"
428 " vec4 aHitTimes = vec4 (MAXFLOAT,\n"
433 " vec3 aRayOriginInverse = -aSubTree.TrsfRay.Origin * aSubTree.Inverse;\n"
435 " vec3 aNodeMin0 = texelFetch (uSceneMinPointTexture, aData.y + 0).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
436 " vec3 aNodeMin1 = texelFetch (uSceneMinPointTexture, aData.y + 1).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
437 " vec3 aNodeMin2 = texelFetch (uSceneMinPointTexture, aData.y + min (2, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
438 " vec3 aNodeMin3 = texelFetch (uSceneMinPointTexture, aData.y + min (3, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
439 " vec3 aNodeMax0 = texelFetch (uSceneMaxPointTexture, aData.y + 0).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
440 " vec3 aNodeMax1 = texelFetch (uSceneMaxPointTexture, aData.y + 1).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
441 " vec3 aNodeMax2 = texelFetch (uSceneMaxPointTexture, aData.y + min (2, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
442 " vec3 aNodeMax3 = texelFetch (uSceneMaxPointTexture, aData.y + min (3, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
444 " vec3 aTimeMax = max (aNodeMin0, aNodeMax0);\n"
445 " vec3 aTimeMin = min (aNodeMin0, aNodeMax0);\n"
447 " float aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
448 " float aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
450 " aHitTimes.x = mix (MAXFLOAT, aTimeEnter,\n"
451 " aTimeEnter <= aTimeLeave && aTimeEnter <= theHit.Time && aTimeLeave >= 0.f);\n"
453 " aTimeMax = max (aNodeMin1, aNodeMax1);\n"
454 " aTimeMin = min (aNodeMin1, aNodeMax1);\n"
456 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
457 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
459 " aHitTimes.y = mix (MAXFLOAT, aTimeEnter,\n"
460 " aTimeEnter <= aTimeLeave && aTimeEnter <= theHit.Time && aTimeLeave >= 0.f);\n"
462 " aTimeMax = max (aNodeMin2, aNodeMax2);\n"
463 " aTimeMin = min (aNodeMin2, aNodeMax2);\n"
465 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
466 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
468 " aHitTimes.z = mix (MAXFLOAT, aTimeEnter,\n"
469 " aTimeEnter <= aTimeLeave && aTimeEnter <= theHit.Time && aTimeLeave >= 0.f && aData.z > 1);\n"
471 " aTimeMax = max (aNodeMin3, aNodeMax3);\n"
472 " aTimeMin = min (aNodeMin3, aNodeMax3);\n"
474 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
475 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
477 " aHitTimes.w = mix (MAXFLOAT, aTimeEnter,\n"
478 " aTimeEnter <= aTimeLeave && aTimeEnter <= theHit.Time && aTimeLeave >= 0.f && aData.z > 2);\n"
480 " ivec4 aChildren = ivec4 (0, 1, 2, 3);\n"
482 " aChildren.xy = aHitTimes.y < aHitTimes.x ? aChildren.yx : aChildren.xy;\n"
483 " aHitTimes.xy = aHitTimes.y < aHitTimes.x ? aHitTimes.yx : aHitTimes.xy;\n"
484 " aChildren.zw = aHitTimes.w < aHitTimes.z ? aChildren.wz : aChildren.zw;\n"
485 " aHitTimes.zw = aHitTimes.w < aHitTimes.z ? aHitTimes.wz : aHitTimes.zw;\n"
486 " aChildren.xz = aHitTimes.z < aHitTimes.x ? aChildren.zx : aChildren.xz;\n"
487 " aHitTimes.xz = aHitTimes.z < aHitTimes.x ? aHitTimes.zx : aHitTimes.xz;\n"
488 " aChildren.yw = aHitTimes.w < aHitTimes.y ? aChildren.wy : aChildren.yw;\n"
489 " aHitTimes.yw = aHitTimes.w < aHitTimes.y ? aHitTimes.wy : aHitTimes.yw;\n"
490 " aChildren.yz = aHitTimes.z < aHitTimes.y ? aChildren.zy : aChildren.yz;\n"
491 " aHitTimes.yz = aHitTimes.z < aHitTimes.y ? aHitTimes.zy : aHitTimes.yz;\n"
493 " if (aHitTimes.x != MAXFLOAT)\n"
495 " int aHitMask = (aHitTimes.w != MAXFLOAT ? aChildren.w : aChildren.z) << 2\n"
496 " | (aHitTimes.z != MAXFLOAT ? aChildren.z : aChildren.y);\n"
498 " if (aHitTimes.y != MAXFLOAT)\n"
499 " Stack[++aHead] = aData.y | (aHitMask << 2 | aChildren.y) << 26;\n"
501 " aNode = aData.y + aChildren.x;\n"
505 " toContinue = (aHead >= 0);\n"
507 " if (aHead == aStop) // go to top-level BVH\n"
509 " aStop = -1; aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
513 " aNode = pop (aHead);\n"
516 " else if (aData.x < 0) // leaf node (contains triangles)\n"
521 " for (int anIdx = aData.y; anIdx <= aData.z; ++anIdx)\n"
523 " ivec4 aTriangle = texelFetch (uGeometryTriangTexture, anIdx + TRG_OFFSET (aSubTree));\n"
525 " vec3 aPoint0 = texelFetch (uGeometryVertexTexture, aTriangle.x += VRT_OFFSET (aSubTree)).xyz;\n"
526 " vec3 aPoint1 = texelFetch (uGeometryVertexTexture, aTriangle.y += VRT_OFFSET (aSubTree)).xyz;\n"
527 " vec3 aPoint2 = texelFetch (uGeometryVertexTexture, aTriangle.z += VRT_OFFSET (aSubTree)).xyz;\n"
529 " IntersectTriangle (aSubTree.TrsfRay, aPoint0, aPoint1, aPoint2, aTimeUV, aNormal);\n"
531 " if (aTimeUV.x < theHit.Time)\n"
533 " aTriIndex = aTriangle;\n"
535 " theTrsfId = TRS_OFFSET (aSubTree);\n"
537 " theHit = SIntersect (aTimeUV.x, aTimeUV.yz, aNormal);\n"
541 " toContinue = (aHead >= 0);\n"
543 " if (aHead == aStop) // go to top-level BVH\n"
545 " aStop = -1; aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
549 " aNode = pop (aHead);\n"
551 " else if (aData.x > 0) // switch node\n"
553 " aSubTree.SubData = ivec4 (4 * aData.x - 4, aData.yzw); // store BVH sub-root\n"
555 " vec4 aInvTransf0 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 0);\n"
556 " vec4 aInvTransf1 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 1);\n"
557 " vec4 aInvTransf2 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 2);\n"
558 " vec4 aInvTransf3 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 3);\n"
560 " aSubTree.TrsfRay.Direct = MatrixColMultiplyDir (theRay.Direct,\n"
565 " aSubTree.Inverse = mix (-UNIT, UNIT, step (ZERO, aSubTree.TrsfRay.Direct)) /\n"
566 " max (abs (aSubTree.TrsfRay.Direct), SMALL);\n"
568 " aSubTree.TrsfRay.Origin = MatrixColMultiplyPnt (theRay.Origin,\n"
574 " aNode = BVH_OFFSET (aSubTree); // go to sub-root node\n"
576 " aStop = aHead; // store current stack pointer\n"
580 " return aTriIndex;\n"
583 "// =======================================================================\n"
584 "// function : SceneAnyHit\n"
585 "// purpose : Finds intersection with any scene triangle\n"
586 "// =======================================================================\n"
587 "float SceneAnyHit (in SRay theRay, in vec3 theInverse, in float theDistance)\n"
589 " float aFactor = 1.f;\n"
591 " int aNode = 0; // node to traverse\n"
592 " int aHead = -1; // pointer of stack\n"
593 " int aStop = -1; // BVH level switch\n"
595 " SSubTree aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
597 " for (bool toContinue = true; toContinue; /* none */)\n"
599 " ivec4 aData = texelFetch (uSceneNodeInfoTexture, aNode);\n"
601 " if (aData.x == 0) // if inner node\n"
603 " aData.y += BVH_OFFSET (aSubTree);\n"
605 " vec4 aHitTimes = vec4 (MAXFLOAT,\n"
610 " vec3 aRayOriginInverse = -aSubTree.TrsfRay.Origin * aSubTree.Inverse;\n"
612 " vec3 aNodeMin0 = texelFetch (uSceneMinPointTexture, aData.y + 0).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
613 " vec3 aNodeMin1 = texelFetch (uSceneMinPointTexture, aData.y + 1).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
614 " vec3 aNodeMin2 = texelFetch (uSceneMinPointTexture, aData.y + min (2, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
615 " vec3 aNodeMin3 = texelFetch (uSceneMinPointTexture, aData.y + min (3, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
616 " vec3 aNodeMax0 = texelFetch (uSceneMaxPointTexture, aData.y + 0).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
617 " vec3 aNodeMax1 = texelFetch (uSceneMaxPointTexture, aData.y + 1).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
618 " vec3 aNodeMax2 = texelFetch (uSceneMaxPointTexture, aData.y + min (2, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
619 " vec3 aNodeMax3 = texelFetch (uSceneMaxPointTexture, aData.y + min (3, aData.z)).xyz * aSubTree.Inverse + aRayOriginInverse;\n"
621 " vec3 aTimeMax = max (aNodeMin0, aNodeMax0);\n"
622 " vec3 aTimeMin = min (aNodeMin0, aNodeMax0);\n"
624 " float aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
625 " float aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
627 " aHitTimes.x = mix (MAXFLOAT, aTimeEnter,\n"
628 " aTimeEnter <= aTimeLeave && aTimeEnter <= theDistance && aTimeLeave >= 0.f);\n"
630 " aTimeMax = max (aNodeMin1, aNodeMax1);\n"
631 " aTimeMin = min (aNodeMin1, aNodeMax1);\n"
633 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
634 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
636 " aHitTimes.y = mix (MAXFLOAT, aTimeEnter,\n"
637 " aTimeEnter <= aTimeLeave && aTimeEnter <= theDistance && aTimeLeave >= 0.f);\n"
639 " aTimeMax = max (aNodeMin2, aNodeMax2);\n"
640 " aTimeMin = min (aNodeMin2, aNodeMax2);\n"
642 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
643 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
645 " aHitTimes.z = mix (MAXFLOAT, aTimeEnter,\n"
646 " aTimeEnter <= aTimeLeave && aTimeEnter <= theDistance && aTimeLeave >= 0.f && aData.z > 1);\n"
648 " aTimeMax = max (aNodeMin3, aNodeMax3);\n"
649 " aTimeMin = min (aNodeMin3, aNodeMax3);\n"
651 " aTimeLeave = min (aTimeMax.x, min (aTimeMax.y, aTimeMax.z));\n"
652 " aTimeEnter = max (aTimeMin.x, max (aTimeMin.y, aTimeMin.z));\n"
654 " aHitTimes.w = mix (MAXFLOAT, aTimeEnter,\n"
655 " aTimeEnter <= aTimeLeave && aTimeEnter <= theDistance && aTimeLeave >= 0.f && aData.z > 2);\n"
657 " ivec4 aChildren = ivec4 (0, 1, 2, 3);\n"
659 " aChildren.xy = aHitTimes.y < aHitTimes.x ? aChildren.yx : aChildren.xy;\n"
660 " aHitTimes.xy = aHitTimes.y < aHitTimes.x ? aHitTimes.yx : aHitTimes.xy;\n"
661 " aChildren.zw = aHitTimes.w < aHitTimes.z ? aChildren.wz : aChildren.zw;\n"
662 " aHitTimes.zw = aHitTimes.w < aHitTimes.z ? aHitTimes.wz : aHitTimes.zw;\n"
663 " aChildren.xz = aHitTimes.z < aHitTimes.x ? aChildren.zx : aChildren.xz;\n"
664 " aHitTimes.xz = aHitTimes.z < aHitTimes.x ? aHitTimes.zx : aHitTimes.xz;\n"
665 " aChildren.yw = aHitTimes.w < aHitTimes.y ? aChildren.wy : aChildren.yw;\n"
666 " aHitTimes.yw = aHitTimes.w < aHitTimes.y ? aHitTimes.wy : aHitTimes.yw;\n"
667 " aChildren.yz = aHitTimes.z < aHitTimes.y ? aChildren.zy : aChildren.yz;\n"
668 " aHitTimes.yz = aHitTimes.z < aHitTimes.y ? aHitTimes.zy : aHitTimes.yz;\n"
670 " if (aHitTimes.x != MAXFLOAT)\n"
672 " int aHitMask = (aHitTimes.w != MAXFLOAT ? aChildren.w : aChildren.z) << 2\n"
673 " | (aHitTimes.z != MAXFLOAT ? aChildren.z : aChildren.y);\n"
675 " if (aHitTimes.y != MAXFLOAT)\n"
676 " Stack[++aHead] = aData.y | (aHitMask << 2 | aChildren.y) << 26;\n"
678 " aNode = aData.y + aChildren.x;\n"
682 " toContinue = (aHead >= 0);\n"
684 " if (aHead == aStop) // go to top-level BVH\n"
686 " aStop = -1; aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
690 " aNode = pop (aHead);\n"
693 " else if (aData.x < 0) // leaf node\n"
698 " for (int anIdx = aData.y; anIdx <= aData.z; ++anIdx)\n"
700 " ivec4 aTriangle = texelFetch (uGeometryTriangTexture, anIdx + TRG_OFFSET (aSubTree));\n"
702 " vec3 aPoint0 = texelFetch (uGeometryVertexTexture, aTriangle.x += VRT_OFFSET (aSubTree)).xyz;\n"
703 " vec3 aPoint1 = texelFetch (uGeometryVertexTexture, aTriangle.y += VRT_OFFSET (aSubTree)).xyz;\n"
704 " vec3 aPoint2 = texelFetch (uGeometryVertexTexture, aTriangle.z += VRT_OFFSET (aSubTree)).xyz;\n"
706 " IntersectTriangle (aSubTree.TrsfRay, aPoint0, aPoint1, aPoint2, aTimeUV, aNormal);\n"
708 "#ifdef TRANSPARENT_SHADOWS\n"
709 " if (aTimeUV.x < theDistance)\n"
711 " aFactor *= 1.f - texelFetch (uRaytraceMaterialTexture, MATERIAL_TRAN (aTriangle.w)).x;\n"
714 " if (aTimeUV.x < theDistance)\n"
721 " toContinue = (aHead >= 0) && (aFactor > 0.1f);\n"
723 " if (aHead == aStop) // go to top-level BVH\n"
725 " aStop = -1; aSubTree = SSubTree (theRay, theInverse, EMPTY_ROOT);\n"
729 " aNode = pop (aHead);\n"
731 " else if (aData.x > 0) // switch node\n"
733 " aSubTree.SubData = ivec4 (4 * aData.x - 4, aData.yzw); // store BVH sub-root\n"
735 " vec4 aInvTransf0 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 0);\n"
736 " vec4 aInvTransf1 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 1);\n"
737 " vec4 aInvTransf2 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 2);\n"
738 " vec4 aInvTransf3 = texelFetch (uSceneTransformTexture, TRS_OFFSET (aSubTree) + 3);\n"
740 " aSubTree.TrsfRay.Direct = MatrixColMultiplyDir (theRay.Direct,\n"
745 " aSubTree.TrsfRay.Origin = MatrixColMultiplyPnt (theRay.Origin,\n"
751 " aSubTree.Inverse = mix (-UNIT, UNIT, step (ZERO, aSubTree.TrsfRay.Direct)) / max (abs (aSubTree.TrsfRay.Direct), SMALL);\n"
753 " aNode = BVH_OFFSET (aSubTree); // go to sub-root node\n"
755 " aStop = aHead; // store current stack pointer\n"
762 "#define PI 3.1415926f\n"
764 "// =======================================================================\n"
765 "// function : Latlong\n"
766 "// purpose : Converts world direction to environment texture coordinates\n"
767 "// =======================================================================\n"
768 "vec2 Latlong (in vec3 thePoint, in float theRadius)\n"
770 " float aPsi = acos (-thePoint.z / theRadius);\n"
772 " float aPhi = atan (thePoint.y, thePoint.x) + PI;\n"
774 " return vec2 (aPhi * 0.1591549f,\n"
775 " aPsi * 0.3183098f);\n"
778 "// =======================================================================\n"
779 "// function : SmoothNormal\n"
780 "// purpose : Interpolates normal across the triangle\n"
781 "// =======================================================================\n"
782 "vec3 SmoothNormal (in vec2 theUV, in ivec4 theTriangle)\n"
784 " vec3 aNormal0 = texelFetch (uGeometryNormalTexture, theTriangle.x).xyz;\n"
785 " vec3 aNormal1 = texelFetch (uGeometryNormalTexture, theTriangle.y).xyz;\n"
786 " vec3 aNormal2 = texelFetch (uGeometryNormalTexture, theTriangle.z).xyz;\n"
788 " return normalize (aNormal1 * theUV.x +\n"
789 " aNormal2 * theUV.y +\n"
790 " aNormal0 * (1.0f - theUV.x - theUV.y));\n"
793 "#define POLYGON_OFFSET_UNIT 0.f\n"
794 "#define POLYGON_OFFSET_FACTOR 1.f\n"
795 "#define POLYGON_OFFSET_SCALE 0.006f\n"
797 "// =======================================================================\n"
798 "// function : PolygonOffset\n"
799 "// purpose : Computes OpenGL polygon offset\n"
800 "// =======================================================================\n"
801 "float PolygonOffset (in vec3 theNormal, in vec3 thePoint)\n"
803 " vec4 aProjectedNorm = vec4 (theNormal, -dot (theNormal, thePoint)) * uUnviewMat;\n"
805 " float aPolygonOffset = POLYGON_OFFSET_UNIT;\n"
807 " if (aProjectedNorm.z * aProjectedNorm.z > 1e-20f)\n"
809 " aProjectedNorm.xy *= 1.f / aProjectedNorm.z;\n"
811 " aPolygonOffset += POLYGON_OFFSET_FACTOR * max (abs (aProjectedNorm.x),\n"
812 " abs (aProjectedNorm.y));\n"
815 " return aPolygonOffset;\n"
818 "// =======================================================================\n"
819 "// function : SmoothUV\n"
820 "// purpose : Interpolates UV coordinates across the triangle\n"
821 "// =======================================================================\n"
822 "#ifdef USE_TEXTURES\n"
823 "vec2 SmoothUV (in vec2 theUV, in ivec4 theTriangle)\n"
825 " vec2 aTexCrd0 = texelFetch (uGeometryTexCrdTexture, theTriangle.x).st;\n"
826 " vec2 aTexCrd1 = texelFetch (uGeometryTexCrdTexture, theTriangle.y).st;\n"
827 " vec2 aTexCrd2 = texelFetch (uGeometryTexCrdTexture, theTriangle.z).st;\n"
829 " return aTexCrd1 * theUV.x +\n"
830 " aTexCrd2 * theUV.y +\n"
831 " aTexCrd0 * (1.0f - theUV.x - theUV.y);\n"
835 "// =======================================================================\n"
836 "// function : FetchEnvironment\n"
838 "// =======================================================================\n"
839 "vec4 FetchEnvironment (in vec2 theTexCoord)\n"
841 " return mix (vec4 (0.0f, 0.0f, 0.0f, 1.0f),\n"
842 " textureLod (uEnvironmentMapTexture, theTexCoord, 0.0f), float (uSphereMapEnabled));\n"
845 "// =======================================================================\n"
846 "// function : Refract\n"
847 "// purpose : Computes refraction ray (also handles TIR)\n"
848 "// =======================================================================\n"
849 "#ifndef PATH_TRACING\n"
850 "vec3 Refract (in vec3 theInput,\n"
851 " in vec3 theNormal,\n"
852 " in float theRefractIndex,\n"
853 " in float theInvRefractIndex)\n"
855 " float aNdotI = dot (theInput, theNormal);\n"
857 " float anIndex = aNdotI < 0.0f\n"
858 " ? theInvRefractIndex\n"
859 " : theRefractIndex;\n"
861 " float aSquare = anIndex * anIndex * (1.0f - aNdotI * aNdotI);\n"
863 " if (aSquare > 1.0f)\n"
865 " return reflect (theInput, theNormal);\n"
868 " float aNdotT = sqrt (1.0f - aSquare);\n"
870 " return normalize (anIndex * theInput -\n"
871 " (anIndex * aNdotI + (aNdotI < 0.0f ? aNdotT : -aNdotT)) * theNormal);\n"
875 "#define MIN_SLOPE 0.0001f\n"
876 "#define EPS_SCALE 8.0000f\n"
878 "#define THRESHOLD vec3 (0.1f)\n"
880 "#define INVALID_BOUNCES 1000\n"
882 "#define LIGHT_POS(index) (2 * index + 1)\n"
883 "#define LIGHT_PWR(index) (2 * index + 0)\n"
885 "// =======================================================================\n"
886 "// function : Radiance\n"
887 "// purpose : Computes color along the given ray\n"
888 "// =======================================================================\n"
889 "#ifndef PATH_TRACING\n"
890 "vec4 Radiance (in SRay theRay, in vec3 theInverse)\n"
892 " vec3 aResult = vec3 (0.0f);\n"
893 " vec4 aWeight = vec4 (1.0f);\n"
897 " float aRaytraceDepth = MAXFLOAT;\n"
899 " for (int aDepth = 0; aDepth < NB_BOUNCES; ++aDepth)\n"
901 " SIntersect aHit = SIntersect (MAXFLOAT, vec2 (ZERO), ZERO);\n"
903 " ivec4 aTriIndex = SceneNearestHit (theRay, theInverse, aHit, aTrsfId);\n"
905 " if (aTriIndex.x == -1)\n"
907 " vec4 aColor = vec4 (0.0);\n"
909 " if (bool(uSphereMapForBack) || aWeight.w == 0.0f /* reflection */)\n"
911 " float aTime = IntersectSphere (theRay, uSceneRadius);\n"
913 " aColor = FetchEnvironment (Latlong (\n"
914 " theRay.Direct * aTime + theRay.Origin, uSceneRadius));\n"
918 " aColor = BackgroundColor();\n"
921 " aResult += aWeight.xyz * aColor.xyz; aWeight.w *= aColor.w;\n"
923 " break; // terminate path\n"
926 " vec3 aInvTransf0 = texelFetch (uSceneTransformTexture, aTrsfId + 0).xyz;\n"
927 " vec3 aInvTransf1 = texelFetch (uSceneTransformTexture, aTrsfId + 1).xyz;\n"
928 " vec3 aInvTransf2 = texelFetch (uSceneTransformTexture, aTrsfId + 2).xyz;\n"
930 " aHit.Normal = normalize (vec3 (dot (aInvTransf0, aHit.Normal),\n"
931 " dot (aInvTransf1, aHit.Normal),\n"
932 " dot (aInvTransf2, aHit.Normal)));\n"
934 " theRay.Origin += theRay.Direct * aHit.Time; // intersection point\n"
936 " // Evaluate depth on first hit\n"
937 " if (aDepth == 0)\n"
939 " vec4 aNDCPoint = uViewMat * vec4 (theRay.Origin, 1.f);\n"
941 " float aPolygonOffset = PolygonOffset (aHit.Normal, theRay.Origin);\n"
942 " aRaytraceDepth = (aNDCPoint.z / aNDCPoint.w + aPolygonOffset * POLYGON_OFFSET_SCALE) * 0.5f + 0.5f;\n"
945 " vec3 aNormal = SmoothNormal (aHit.UV, aTriIndex);\n"
947 " aNormal = normalize (vec3 (dot (aInvTransf0, aNormal),\n"
948 " dot (aInvTransf1, aNormal),\n"
949 " dot (aInvTransf2, aNormal)));\n"
951 " vec3 aAmbient = texelFetch (\n"
952 " uRaytraceMaterialTexture, MATERIAL_AMBN (aTriIndex.w)).rgb;\n"
953 " vec4 aDiffuse = texelFetch (\n"
954 " uRaytraceMaterialTexture, MATERIAL_DIFF (aTriIndex.w));\n"
955 " vec4 aSpecular = texelFetch (\n"
956 " uRaytraceMaterialTexture, MATERIAL_SPEC (aTriIndex.w));\n"
957 " vec4 aOpacity = texelFetch (\n"
958 " uRaytraceMaterialTexture, MATERIAL_TRAN (aTriIndex.w));\n"
960 "#ifdef USE_TEXTURES\n"
961 " if (aDiffuse.w >= 0.f)\n"
963 " vec4 aTexCoord = vec4 (SmoothUV (aHit.UV, aTriIndex), 0.f, 1.f);\n"
965 " vec4 aTrsfRow1 = texelFetch (\n"
966 " uRaytraceMaterialTexture, MATERIAL_TRS1 (aTriIndex.w));\n"
967 " vec4 aTrsfRow2 = texelFetch (\n"
968 " uRaytraceMaterialTexture, MATERIAL_TRS2 (aTriIndex.w));\n"
970 " aTexCoord.st = vec2 (dot (aTrsfRow1, aTexCoord),\n"
971 " dot (aTrsfRow2, aTexCoord));\n"
973 " vec3 aTexColor = textureLod (\n"
974 " sampler2D (uTextureSamplers[int(aDiffuse.w)]), aTexCoord.st, 0.f).rgb;\n"
976 " aDiffuse.rgb *= aTexColor;\n"
977 " aAmbient.rgb *= aTexColor;\n"
981 " vec3 aEmission = texelFetch (\n"
982 " uRaytraceMaterialTexture, MATERIAL_EMIS (aTriIndex.w)).rgb;\n"
984 " float aGeomFactor = dot (aNormal, theRay.Direct);\n"
986 " aResult.xyz += aWeight.xyz * aOpacity.x * (\n"
987 " uGlobalAmbient.xyz * aAmbient * max (abs (aGeomFactor), 0.5f) + aEmission);\n"
989 " vec3 aSidedNormal = mix (aNormal, -aNormal, step (0.0f, aGeomFactor));\n"
991 " for (int aLightIdx = 0; aLightIdx < uLightCount; ++aLightIdx)\n"
993 " vec4 aLight = texelFetch (\n"
994 " uRaytraceLightSrcTexture, LIGHT_POS (aLightIdx));\n"
996 " float aDistance = MAXFLOAT;\n"
998 " if (aLight.w != 0.0f) // point light source\n"
1000 " aDistance = length (aLight.xyz -= theRay.Origin);\n"
1002 " aLight.xyz *= 1.0f / aDistance;\n"
1005 " float aLdotN = dot (aLight.xyz, aSidedNormal);\n"
1007 " if (aLdotN > 0.0f) // first check if light source is important\n"
1009 " float aVisibility = 1.0f;\n"
1011 " if (bool(uShadowsEnabled))\n"
1013 " SRay aShadow = SRay (theRay.Origin, aLight.xyz);\n"
1015 " aShadow.Origin += uSceneEpsilon * (aLight.xyz +\n"
1016 " mix (-aHit.Normal, aHit.Normal, step (0.0f, dot (aHit.Normal, aLight.xyz))));\n"
1018 " vec3 aInverse = 1.0f / max (abs (aLight.xyz), SMALL);\n"
1020 " aVisibility = SceneAnyHit (\n"
1021 " aShadow, mix (-aInverse, aInverse, step (ZERO, aLight.xyz)), aDistance);\n"
1024 " if (aVisibility > 0.0f)\n"
1026 " vec3 aIntensity = vec3 (texelFetch (\n"
1027 " uRaytraceLightSrcTexture, LIGHT_PWR (aLightIdx)));\n"
1029 " float aRdotV = dot (reflect (aLight.xyz, aSidedNormal), theRay.Direct);\n"
1031 " aResult.xyz += aWeight.xyz * (aOpacity.x * aVisibility) * aIntensity *\n"
1032 " (aDiffuse.xyz * aLdotN + aSpecular.xyz * pow (max (0.f, aRdotV), aSpecular.w));\n"
1037 " if (aOpacity.x != 1.0f)\n"
1039 " aWeight *= aOpacity.y;\n"
1041 " if (aOpacity.z != 1.0f)\n"
1043 " theRay.Direct = Refract (theRay.Direct, aNormal, aOpacity.z, aOpacity.w);\n"
1048 " aWeight *= bool(uReflectEnabled) ?\n"
1049 " texelFetch (uRaytraceMaterialTexture, MATERIAL_REFL (aTriIndex.w)) : vec4 (0.0f);\n"
1051 " vec3 aReflect = reflect (theRay.Direct, aNormal);\n"
1053 " if (dot (aReflect, aHit.Normal) * dot (theRay.Direct, aHit.Normal) > 0.0f)\n"
1055 " aReflect = reflect (theRay.Direct, aHit.Normal);\n"
1058 " theRay.Direct = aReflect;\n"
1061 " if (all (lessThanEqual (aWeight.xyz, THRESHOLD)))\n"
1063 " aDepth = INVALID_BOUNCES;\n"
1065 " else if (aOpacity.x == 1.0f || aOpacity.z != 1.0f) // if no simple transparency\n"
1067 " theRay.Origin += aHit.Normal * mix (\n"
1068 " -uSceneEpsilon, uSceneEpsilon, step (0.0f, dot (aHit.Normal, theRay.Direct)));\n"
1070 " theInverse = 1.0f / max (abs (theRay.Direct), SMALL);\n"
1072 " theInverse = mix (-theInverse, theInverse, step (ZERO, theRay.Direct));\n"
1075 " theRay.Origin += theRay.Direct * uSceneEpsilon;\n"
1078 " gl_FragDepth = aRaytraceDepth;\n"
1080 " return vec4 (aResult.x,\n"