x_ite 8.11.5 → 8.11.6

package/dist/assets/components/ParticleSystems.js

@@ -1,11 +1,11 @@
-/* X_ITE v8.11.5 */(() => { // webpackBootstrap
+/* X_ITE v8.11.6 */(() => { // webpackBootstrap
 /******/ 	"use strict";
 /******/ 	var __webpack_modules__ = ({
 
 /***/ 355:
 /***/ ((module) => {
 
-module.exports = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("lib/jquery");
+module.exports = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("lib/jquery");
 
 /***/ })
 
@@ -73,22 +73,22 @@ var __webpack_exports__ = {};
 // UNUSED EXPORTS: default
 
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Components\")"
-const Components_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Components");
+const Components_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Components");
 var Components_default = /*#__PURE__*/__webpack_require__.n(Components_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Fields\")"
-const Fields_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Fields");
+const Fields_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Fields");
 var Fields_default = /*#__PURE__*/__webpack_require__.n(Fields_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Base/X3DFieldDefinition\")"
-const X3DFieldDefinition_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Base/X3DFieldDefinition");
+const X3DFieldDefinition_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Base/X3DFieldDefinition");
 var X3DFieldDefinition_default = /*#__PURE__*/__webpack_require__.n(X3DFieldDefinition_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Base/FieldDefinitionArray\")"
-const FieldDefinitionArray_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Base/FieldDefinitionArray");
+const FieldDefinitionArray_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Base/FieldDefinitionArray");
 var FieldDefinitionArray_default = /*#__PURE__*/__webpack_require__.n(FieldDefinitionArray_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Components/Core/X3DNode\")"
-const X3DNode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Components/Core/X3DNode");
+const X3DNode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Components/Core/X3DNode");
 var X3DNode_default = /*#__PURE__*/__webpack_require__.n(X3DNode_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Namespace\")"
-const Namespace_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Namespace");
+const Namespace_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Namespace");
 var Namespace_default = /*#__PURE__*/__webpack_require__.n(Namespace_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Browser/ParticleSystems/GeometryTypes.js
 /*******************************************************************************
@@ -155,477 +155,35 @@ const __default__ = GeometryTypes;
 Namespace_default().set ("x_ite/Browser/ParticleSystems/GeometryTypes", __default__);
 /* harmony default export */ const ParticleSystems_GeometryTypes = (__default__);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Base/X3DConstants\")"
-const X3DConstants_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Base/X3DConstants");
+const X3DConstants_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Base/X3DConstants");
 var X3DConstants_default = /*#__PURE__*/__webpack_require__.n(X3DConstants_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Browser/ParticleSystems/Line3.glsl.js
-const Line3_glsl_default_ = /* glsl*/ `
-struct Line3 {
-   vec3 point;
-   vec3 direction;
-};
-
-// Line3
-// line3 (const in vec3 point1, const in vec3 point2)
-// {
-//    return Line3 (point1, normalize (point2 - point1));
-// }
-
-/* Line intersect triangle */
-
-bool
-intersects (const in Line3 line, const in vec3 a, const in vec3 b, const in vec3 c, out vec3 r)
-{
-   // find vectors for two edges sharing vert0
-   vec3 edge1 = b - a;
-   vec3 edge2 = c - a;
-
-   // begin calculating determinant - also used to calculate U parameter
-   vec3 pvec = cross (line .direction, edge2);
-
-   // if determinant is near zero, ray lies in plane of triangle
-   float det = dot (edge1, pvec);
-
-   // Non culling intersection
-
-   if (det == 0.0)
-      return false;
-
-   float inv_det = 1.0 / det;
-
-   // calculate distance from vert0 to ray point
-   vec3 tvec = line .point - a;
-
-   // calculate U parameter and test bounds
-   float u = dot (tvec, pvec) * inv_det;
-
-   if (u < 0.0 || u > 1.0)
-      return false;
-
-   // prepare to test V parameter
-   vec3 qvec = cross (tvec, edge1);
-
-   // calculate V parameter and test bounds
-   float v = dot (line .direction, qvec) * inv_det;
-
-   if (v < 0.0 || u + v > 1.0)
-      return false;
-
-   r = vec3 (u, v, 1.0 - u - v);
-
-   return true;
-}
-`;
+const Line3_glsl_default_ = /* glsl */ `struct Line3{vec3 point;vec3 direction;};bool intersects(const in Line3 line,const in vec3 a,const in vec3 b,const in vec3 c,out vec3 r){vec3 edge1=b-a;vec3 edge2=c-a;vec3 pvec=cross(line.direction,edge2);float det=dot(edge1,pvec);if(det==0.0)return false;float inv_det=1.0/det;vec3 tvec=line.point-a;float u=dot(tvec,pvec)*inv_det;if(u<0.0||u>1.0)return false;vec3 qvec=cross(tvec,edge1);float v=dot(line.direction,qvec)*inv_det;if(v<0.0||u+v>1.0)return false;r=vec3(u,v,1.0-u-v);return true;}
+`
 ;
 
 Namespace_default().set ("x_ite/Browser/ParticleSystems/Line3.glsl", Line3_glsl_default_);
 /* harmony default export */ const Line3_glsl = (Line3_glsl_default_);
 ;// CONCATENATED MODULE: ./src/x_ite/Browser/ParticleSystems/Plane3.glsl.js
-const Plane3_glsl_default_ = /* glsl*/ `
-struct Plane3
-{
-   vec3  normal;
-   float distanceFromOrigin;
-};
-
-Plane3
-plane3 (const in vec3 point, const in vec3 normal)
-{
-   return Plane3 (normal, dot (normal, point));
-}
-
-float
-plane_distance (const in Plane3 plane, const in vec3 point)
-{
-   return dot (point, plane .normal) - plane .distanceFromOrigin;
-}
-
-/* Plane intersect line */
-bool
-intersects (const in Plane3 plane, const in Line3 line, out vec3 point)
-{
-   // Check if the line is parallel to the plane.
-   float theta = dot (line .direction, plane .normal);
-
-   // Plane and line are parallel.
-   if (theta == 0.0)
-      return false;
-
-   // Plane and line are not parallel. The intersection point can be calculated now.
-   float t = (plane .distanceFromOrigin - dot (plane .normal, line .point)) / theta;
-
-   point = line .point + line .direction * t;
-
-   return true;
-}
-
-/* Find find the first point that is farther to the plane than value. */
-// int
-// upper_bound (const in vec4 points [ARRAY_SIZE], in int count, const in float value, const in Plane3 plane)
-// {
-//    int first = 0;
-//    int step  = 0;
-
-//    while (count > 0)
-//    {
-//       int index = first;
-
-//       step = count >> 1;
-
-//       index += step;
-
-//       if (value < plane_distance (plane, points [index] .xyz))
-//       {
-//          count = step;
-//       }
-//       else
-//       {
-//          first  = ++ index;
-//          count -= step + 1;
-//       }
-//    }
-
-//    return first;
-// }
-
-/* CombSort: sort points in distance to a plane. */
-void
-sort (inout vec4 points [ARRAY_SIZE], const in int count, const in Plane3 plane)
-{
-   const float shrink = 1.0 / 1.3;
-
-   int  gap       = count;
-   bool exchanged = true;
-
-   while (exchanged)
-   {
-      gap = int (float (gap) * shrink);
-
-      if (gap <= 1)
-      {
-         exchanged = false;
-         gap       = 1;
-      }
-
-      for (int i = 0, l = count - gap; i < l; ++ i)
-      {
-         int j = gap + i;
-
-         if (plane_distance (plane, points [i] .xyz) > plane_distance (plane, points [j] .xyz))
-         {
-            vec4 tmp1 = points [i];
-            points [i] = points [j];
-            points [j] = tmp1;
-
-            exchanged = true;
-         }
-      }
-   }
-}
-
-
-// /* CombSort: sort points and normals in distance to a plane. */
-// void
-// sort (inout vec4 points [ARRAY_SIZE], inout vec3 normals [ARRAY_SIZE], const in int count, const in Plane3 plane)
-// {
-//    const float shrink = 1.0 / 1.3;
-
-//    int  gap       = count;
-//    bool exchanged = true;
-
-//    while (exchanged)
-//    {
-//       gap = int (float (gap) * shrink);
-
-//       if (gap <= 1)
-//       {
-//          exchanged = false;
-//          gap       = 1;
-//       }
-
-//       for (int i = 0, l = count - gap; i < l; ++ i)
-//       {
-//          int j = gap + i;
-
-//          if (plane_distance (plane, points [i] .xyz) > plane_distance (plane, points [j] .xyz))
-//          {
-//             vec4 tmp1 = points [i];
-//             points [i] = points [j];
-//             points [j] = tmp1;
-
-//             vec3 tmp2   = normals [i];
-//             normals [i] = normals [j];
-//             normals [j] = tmp2;
-
-//             exchanged = true;
-//          }
-//       }
-//    }
-// }
-
-int
-min_index (const in vec4 points [ARRAY_SIZE], const in int count, const in float value, const in Plane3 plane)
-{
-   int   index = -1;
-   float dist  = 1000000.0;
-
-   for (int i = 0; i < count; ++ i)
-   {
-      float d = plane_distance (plane, points [i] .xyz);
-
-      if (d >= value && d < dist)
-      {
-         dist  = d;
-         index = i;
-      }
-   }
-
-   return index;
-}
-`;
+const Plane3_glsl_default_ = /* glsl */ `struct Plane3{vec3 normal;float distanceFromOrigin;};Plane3 plane3(const in vec3 point,const in vec3 normal){return Plane3(normal,dot(normal,point));}float plane_distance(const in Plane3 plane,const in vec3 point){return dot(point,plane.normal)-plane.distanceFromOrigin;}bool intersects(const in Plane3 plane,const in Line3 line,out vec3 point){float theta=dot(line.direction,plane.normal);if(theta==0.0)return false;float t=(plane.distanceFromOrigin-dot(plane.normal,line.point))/theta;point=line.point+line.direction*t;return true;}void sort(inout vec4 points[ARRAY_SIZE],const in int count,const in Plane3 plane){const float shrink=1.0/1.3;int gap=count;bool exchanged=true;while(exchanged){gap=int(float(gap)*shrink);if(gap<=1){exchanged=false;gap=1;}for(int i=0,l=count-gap;i<l;++i){int j=gap+i;if(plane_distance(plane,points[i].xyz)>plane_distance(plane,points[j].xyz)){vec4 tmp1=points[i];points[i]=points[j];points[j]=tmp1;exchanged=true;}}}}int min_index(const in vec4 points[ARRAY_SIZE],const in int count,const in float value,const in Plane3 plane){int index=-1;float dist=1000000.0;for(int i=0;i<count;++i){float d=plane_distance(plane,points[i].xyz);if(d>=value&&d<dist){dist=d;index=i;}}return index;}
+`
 ;
 
 Namespace_default().set ("x_ite/Browser/ParticleSystems/Plane3.glsl", Plane3_glsl_default_);
 /* harmony default export */ const Plane3_glsl = (Plane3_glsl_default_);
 ;// CONCATENATED MODULE: ./src/x_ite/Browser/ParticleSystems/Box3.glsl.js
-const Box3_glsl_default_ = /* glsl*/ `
-bool
-intersects (const in vec3 min, const in vec3 max, const in Line3 line)
-{
-   vec3 intersection;
-
-   // front
-
-   if (intersects (plane3 (max, vec3 (0.0, 0.0, 1.0)), line, intersection))
-   {
-      if (all (greaterThanEqual (vec4 (intersection .xy, max .xy), vec4 (min .xy, intersection .xy))))
-         return true;
-   }
-
-   // back
-
-   if (intersects (plane3 (min, vec3 (0.0, 0.0, -1.0)), line, intersection))
-   {
-      if (all (greaterThanEqual (vec4 (intersection .xy, max .xy), vec4 (min .xy, intersection .xy))))
-         return true;
-   }
-
-   // top
-
-   if (intersects (plane3 (max, vec3 (0.0, 1.0, 0.0)), line, intersection))
-   {
-      if (all (greaterThanEqual (vec4 (intersection .xz, max .xz), vec4 (min .xz, intersection .xz))))
-         return true;
-   }
-
-   // bottom
-
-   if (intersects (plane3 (min, vec3 (0.0, -1.0, 0.0)), line, intersection))
-   {
-      if (all (greaterThanEqual (vec4 (intersection .xz, max .xz), vec4 (min .xz, intersection .xz))))
-         return true;
-   }
-
-   // right
-
-   if (intersects (plane3 (max, vec3 (1.0, 0.0, 0.0)), line, intersection))
-   {
-      if (all (greaterThanEqual (vec4 (intersection .yz, max .yz), vec4 (min .yz, intersection .yz))))
-         return true;
-   }
-
-   return false;
-}
-`;
+const Box3_glsl_default_ = /* glsl */ `bool intersects(const in vec3 min,const in vec3 max,const in Line3 line){vec3 intersection;if(intersects(plane3(max,vec3(0.0,0.0,1.0)),line,intersection)){if(all(greaterThanEqual(vec4(intersection.xy,max.xy),vec4(min.xy,intersection.xy))))return true;}if(intersects(plane3(min,vec3(0.0,0.0,-1.0)),line,intersection)){if(all(greaterThanEqual(vec4(intersection.xy,max.xy),vec4(min.xy,intersection.xy))))return true;}if(intersects(plane3(max,vec3(0.0,1.0,0.0)),line,intersection)){if(all(greaterThanEqual(vec4(intersection.xz,max.xz),vec4(min.xz,intersection.xz))))return true;}if(intersects(plane3(min,vec3(0.0,-1.0,0.0)),line,intersection)){if(all(greaterThanEqual(vec4(intersection.xz,max.xz),vec4(min.xz,intersection.xz))))return true;}if(intersects(plane3(max,vec3(1.0,0.0,0.0)),line,intersection)){if(all(greaterThanEqual(vec4(intersection.yz,max.yz),vec4(min.yz,intersection.yz))))return true;}return false;}
+`
 ;
 
 Namespace_default().set ("x_ite/Browser/ParticleSystems/Box3.glsl", Box3_glsl_default_);
 /* harmony default export */ const Box3_glsl = (Box3_glsl_default_);
 ;// CONCATENATED MODULE: ./src/x_ite/Browser/ParticleSystems/BVH.glsl.js
-const BVH_glsl_default_ = /* glsl*/ `
-#define BVH_NODE        0
-#define BVH_TRIANGLE    1
-#define BVH_STACK_SIZE  32
-
-int bvhNodeIndex = 0;
-
-void
-setBVHIndex (const in int index)
-{
-   bvhNodeIndex = index;
-}
-
-int
-getBVHRoot (const in sampler2D volume, const in int hierarchyIndex, const in int rootIndex)
-{
-   return int (texelFetch (volume, rootIndex, 0) .x) + hierarchyIndex;
-}
-
-int
-getBVHType (const in sampler2D volume)
-{
-   return int (texelFetch (volume, bvhNodeIndex, 0) .x);
-}
-
-vec3
-getBVHMin (const in sampler2D volume)
-{
-   return texelFetch (volume, bvhNodeIndex + 1, 0) .xyz;
-}
-
-vec3
-getBVHMax (const in sampler2D volume)
-{
-   return texelFetch (volume, bvhNodeIndex + 2, 0) .xyz;
-}
-
-int
-getBVHLeft (const in sampler2D volume, const in int hierarchyIndex)
-{
-   return int (texelFetch (volume, bvhNodeIndex, 0) .y) + hierarchyIndex;
-}
-
-int
-getBVHRight (const in sampler2D volume, const in int hierarchyIndex)
-{
-   return int (texelFetch (volume, bvhNodeIndex, 0) .z) + hierarchyIndex;
-}
-
-int
-getBVHTriangle (const in sampler2D volume)
-{
-   return int (texelFetch (volume, bvhNodeIndex, 0) .y);
-}
-
-/* Ray triangle intersection test */
-
-int
-getIntersections (const in sampler2D volume, const in int verticesIndex, const in int hierarchyIndex, const in int rootIndex, const in Line3 line, out vec4 points [ARRAY_SIZE])
-{
-   int current    = getBVHRoot (volume, hierarchyIndex, rootIndex);
-   int count      = 0;
-   int stackIndex = -1;
-   int stack [BVH_STACK_SIZE];
-
-   while (stackIndex >= 0 || current >= 0)
-   {
-      if (current >= 0)
-      {
-         setBVHIndex (current);
-
-         if (getBVHType (volume) == BVH_NODE)
-         {
-            // Node
-
-            if (intersects (getBVHMin (volume), getBVHMax (volume), line))
-            {
-               stack [++ stackIndex] = current;
-
-               current = getBVHLeft (volume, hierarchyIndex);
-            }
-            else
-            {
-               current = -1;
-            }
-         }
-         else
-         {
-            // Triangle
-
-            int  t = getBVHTriangle (volume);
-            int  v = verticesIndex + t;
-            vec3 r = vec3 (0.0);
-
-            vec3 a = texelFetch (volume, v,     0) .xyz;
-            vec3 b = texelFetch (volume, v + 1, 0) .xyz;
-            vec3 c = texelFetch (volume, v + 2, 0) .xyz;
-
-            if (intersects (line, a, b, c, r))
-               points [count ++] = vec4 (r .z * a + r .x * b + r .y * c, 1.0);
-
-            current = -1;
-         }
-      }
-      else
-      {
-         setBVHIndex (stack [stackIndex --]);
-
-         current = getBVHRight (volume, hierarchyIndex);
-      }
-   }
-
-   return count;
-}
-
-int
-getIntersections (const in sampler2D volume, const in int verticesIndex, const in int normalsIndex, const in int hierarchyIndex, const in int rootIndex, const in Line3 line, out vec4 points [ARRAY_SIZE], out vec3 normals [ARRAY_SIZE])
-{
-   int current    = getBVHRoot (volume, hierarchyIndex, rootIndex);
-   int count      = 0;
-   int stackIndex = -1;
-   int stack [BVH_STACK_SIZE];
-
-   while (stackIndex >= 0 || current >= 0)
-   {
-      if (current >= 0)
-      {
-         setBVHIndex (current);
-
-         if (getBVHType (volume) == BVH_NODE)
-         {
-            // Node
-
-            if (intersects (getBVHMin (volume), getBVHMax (volume), line))
-            {
-               stack [++ stackIndex] = current;
-
-               current = getBVHLeft (volume, hierarchyIndex);
-            }
-            else
-            {
-               current = -1;
-            }
-         }
-         else
-         {
-            // Triangle
-
-            int  t = getBVHTriangle (volume);
-            int  v = verticesIndex + t;
-            vec3 r = vec3 (0.0);
-
-            vec3 a = texelFetch (volume, v,     0) .xyz;
-            vec3 b = texelFetch (volume, v + 1, 0) .xyz;
-            vec3 c = texelFetch (volume, v + 2, 0) .xyz;
-
-            if (intersects (line, a, b, c, r))
-            {
-               points [count] = vec4 (r .z * a + r .x * b + r .y * c, 1.0);
-
-               int n = normalsIndex + t;
-
-               vec3 n0 = texelFetch (volume, n,     0) .xyz;
-               vec3 n1 = texelFetch (volume, n + 1, 0) .xyz;
-               vec3 n2 = texelFetch (volume, n + 2, 0) .xyz;
-
-               normals [count] = save_normalize (r .z * n0 + r .x * n1 + r .y * n2);
-
-               ++ count;
-            }
-
-            current = -1;
-         }
-      }
-      else
-      {
-         setBVHIndex (stack [stackIndex --]);
-
-         current = getBVHRight (volume, hierarchyIndex);
-      }
-   }
-
-   return count;
-}
-`;
+const BVH_glsl_default_ = /* glsl */ `#define BVH_NODE 0
+#define BVH_TRIANGLE 1
+#define BVH_STACK_SIZE 32
+int bvhNodeIndex=0;void setBVHIndex(const in int index){bvhNodeIndex=index;}int getBVHRoot(const in sampler2D volume,const in int hierarchyIndex,const in int rootIndex){return int(texelFetch(volume,rootIndex,0).x)+hierarchyIndex;}int getBVHType(const in sampler2D volume){return int(texelFetch(volume,bvhNodeIndex,0).x);}vec3 getBVHMin(const in sampler2D volume){return texelFetch(volume,bvhNodeIndex+1,0).xyz;}vec3 getBVHMax(const in sampler2D volume){return texelFetch(volume,bvhNodeIndex+2,0).xyz;}int getBVHLeft(const in sampler2D volume,const in int hierarchyIndex){return int(texelFetch(volume,bvhNodeIndex,0).y)+hierarchyIndex;}int getBVHRight(const in sampler2D volume,const in int hierarchyIndex){return int(texelFetch(volume,bvhNodeIndex,0).z)+hierarchyIndex;}int getBVHTriangle(const in sampler2D volume){return int(texelFetch(volume,bvhNodeIndex,0).y);}int getIntersections(const in sampler2D volume,const in int verticesIndex,const in int hierarchyIndex,const in int rootIndex,const in Line3 line,out vec4 points[ARRAY_SIZE]){int current=getBVHRoot(volume,hierarchyIndex,rootIndex);int count=0;int stackIndex=-1;int stack[BVH_STACK_SIZE];while(stackIndex>=0||current>=0){if(current>=0){setBVHIndex(current);if(getBVHType(volume)==BVH_NODE){if(intersects(getBVHMin(volume),getBVHMax(volume),line)){stack[++stackIndex]=current;current=getBVHLeft(volume,hierarchyIndex);}else{current=-1;}}else{int t=getBVHTriangle(volume);int v=verticesIndex+t;vec3 r=vec3(0.0);vec3 a=texelFetch(volume,v,0).xyz;vec3 b=texelFetch(volume,v+1,0).xyz;vec3 c=texelFetch(volume,v+2,0).xyz;if(intersects(line,a,b,c,r))points[count++]=vec4(r.z*a+r.x*b+r.y*c,1.0);current=-1;}}else{setBVHIndex(stack[stackIndex--]);current=getBVHRight(volume,hierarchyIndex);}}return count;}int getIntersections(const in sampler2D volume,const in int verticesIndex,const in int normalsIndex,const in int hierarchyIndex,const in int rootIndex,const in Line3 line,out vec4 points[ARRAY_SIZE],out vec3 normals[ARRAY_SIZE]){int current=getBVHRoot(volume,hierarchyIndex,rootIndex);int count=0;int stackIndex=-1;int stack[BVH_STACK_SIZE];while(stackIndex>=0||current>=0){if(current>=0){setBVHIndex(current);if(getBVHType(volume)==BVH_NODE){if(intersects(getBVHMin(volume),getBVHMax(volume),line)){stack[++stackIndex]=current;current=getBVHLeft(volume,hierarchyIndex);}else{current=-1;}}else{int t=getBVHTriangle(volume);int v=verticesIndex+t;vec3 r=vec3(0.0);vec3 a=texelFetch(volume,v,0).xyz;vec3 b=texelFetch(volume,v+1,0).xyz;vec3 c=texelFetch(volume,v+2,0).xyz;if(intersects(line,a,b,c,r)){points[count]=vec4(r.z*a+r.x*b+r.y*c,1.0);int n=normalsIndex+t;vec3 n0=texelFetch(volume,n,0).xyz;vec3 n1=texelFetch(volume,n+1,0).xyz;vec3 n2=texelFetch(volume,n+2,0).xyz;normals[count]=save_normalize(r.z*n0+r.x*n1+r.y*n2);++count;}current=-1;}}else{setBVHIndex(stack[stackIndex--]);current=getBVHRight(volume,hierarchyIndex);}}return count;}
+`
 ;
 
 Namespace_default().set ("x_ite/Browser/ParticleSystems/BVH.glsl", BVH_glsl_default_);
@@ -920,593 +478,18 @@ Object .assign (Object .setPrototypeOf (X3DParticleEmitterNode .prototype, (X3DN
          gl      = browser .getContext ();
 
       const vertexShaderSource = /* glsl */ `#version 300 es
-
-      precision highp float;
-      precision highp int;
-      precision highp sampler2D;
-
-      uniform int   randomSeed;
-      uniform int   geometryType;
-      uniform bool  createParticles;
-      uniform float particleLifetime;
-      uniform float lifetimeVariation;
-      uniform float deltaTime;
-      uniform vec2  particleSize;
-
-      uniform int       numForces;
-      uniform sampler2D forces;
-
-      uniform int       boundedVerticesIndex;
-      uniform int       boundedNormalsIndex;
-      uniform int       boundedHierarchyIndex;
-      uniform int       boundedHierarchyRoot;
-      uniform sampler2D boundedVolume;
-
-      uniform int       numColors;
-      uniform sampler2D colorRamp;
-
-      uniform int       texCoordCount;
-      uniform int       numTexCoords;
-      uniform sampler2D texCoordRamp;
-
-      ${Object .values (this .uniforms) .join ("\n")}
-
-      in vec4 input0; // (life, lifetime, elapsedTime, texCoordIndex0)
-      in vec4 input2; // (velocity, 0.0)
-      in vec4 input6; // position
-
-      out vec4 output0; // (life, lifetime, elapsedTime, texCoordIndex0)
-      out vec4 output1; // (color)
-      out vec4 output2; // (velocity, 0.0)
-
-      out vec4 output3; // scale rotation matrix
-      out vec4 output4; // scale rotation matrix
-      out vec4 output5; // scale rotation matrix
-      out vec4 output6; // position
-
-      // Constants
-
-      ${Object .entries (ParticleSystems_GeometryTypes) .map (([k, v]) => `#define ${k} ${v}`) .join ("\n")}
-
-      const int   ARRAY_SIZE = 32;
-      const float M_PI       = 3.14159265359;
-
-      uniform float NaN;
-
-      // Texture
-
-      vec4
-      texelFetch (const in sampler2D sampler, const in int index, const in int lod)
-      {
-         int   x = textureSize (sampler, lod) .x;
-         ivec2 p = ivec2 (index % x, index / x);
-         vec4  t = texelFetch (sampler, p, lod);
-
-         return t;
-      }
-
-      // Math
-
-      // Save normalize, that will not divide by zero.
-      vec3
-      save_normalize (const in vec3 vector)
-      {
-         float l = length (vector);
-
-         if (l == 0.0)
-            return vec3 (0.0);
-
-         return vector / l;
-      }
-
-      // Quaternion
-
-      vec4
-      Quaternion (const in vec3 fromVector, const in vec3 toVector)
-      {
-         vec3 from = save_normalize (fromVector);
-         vec3 to   = save_normalize (toVector);
-
-         float cos_angle = dot (from, to);
-         vec3  cross_vec = cross (from, to);
-         float cross_len = length (cross_vec);
-
-         if (cross_len == 0.0)
-         {
-            if (cos_angle > 0.0)
-            {
-               return vec4 (0.0, 0.0, 0.0, 1.0);
-            }
-            else
-            {
-               vec3 t = cross (from, vec3 (1.0, 0.0, 0.0));
-
-               if (dot (t, t) == 0.0)
-                  t = cross (from, vec3 (0.0, 1.0, 0.0));
-
-               t = save_normalize (t);
-
-               return vec4 (t, 0.0);
-            }
-         }
-         else
-         {
-            float s = sqrt (abs (1.0 - cos_angle) * 0.5);
-
-            cross_vec = save_normalize (cross_vec);
-
-            return vec4 (cross_vec * s, sqrt (abs (1.0 + cos_angle) * 0.5));
-         }
-      }
-
-      vec3
-      multVecQuat (const in vec3 v, const in vec4 q)
-      {
-         float a = q .w * q .w - q .x * q .x - q .y * q .y - q .z * q .z;
-         float b = 2.0 * (v .x * q .x + v .y * q .y + v .z * q .z);
-         float c = 2.0 * q .w;
-         vec3  r = a * v .xyz + b * q .xyz + c * (q .yzx * v .zxy - q .zxy * v .yzx);
-
-         return r;
-      }
-
-      mat3
-      Matrix3 (const in vec4 quaternion)
-      {
-         float x = quaternion .x;
-         float y = quaternion .y;
-         float z = quaternion .z;
-         float w = quaternion .w;
-         float A = y * y;
-         float B = z * z;
-         float C = x * y;
-         float D = z * w;
-         float E = z * x;
-         float F = y * w;
-         float G = x * x;
-         float H = y * z;
-         float I = x * w;
-
-         return mat3 (1.0 - 2.0 * (A + B),
-                      2.0 * (C + D),
-                      2.0 * (E - F),
-                      2.0 * (C - D),
-                      1.0 - 2.0 * (B + G),
-                      2.0 * (H + I),
-                      2.0 * (E + F),
-                      2.0 * (H - I),
-                      1.0 - 2.0 * (A + G));
-      }
-
-      /* Random number generation */
-
-      uint seed = 1u;
-
-      void
-      srand (const in int value)
-      {
-         seed = uint (value);
-      }
-
-      // Return a uniform distributed random floating point number in the interval [0, 1].
-      float
-      random ()
-      {
-         seed = seed * 1103515245u + 12345u;
-
-         return float (seed) / 4294967295.0;
-      }
-
-      float
-      getRandomValue (const in float min, const in float max)
-      {
-         return min + random () * (max - min);
-      }
-
-      float
-      getRandomLifetime ()
-      {
-         float v    = particleLifetime * lifetimeVariation;
-         float min_ = max (0.0, particleLifetime - v);
-         float max_ = particleLifetime + v;
-
-         return getRandomValue (min_, max_);
-      }
-
-      float
-      getRandomSpeed ()
-      {
-         float v    = speed * variation;
-         float min_ = max (0.0, speed - v);
-         float max_ = speed + v;
-
-         return getRandomValue (min_, max_);
-      }
-
-      vec3
-      getRandomNormal ()
-      {
-         float theta = getRandomValue (-M_PI, M_PI);
-         float cphi  = getRandomValue (-1.0, 1.0);
-         float r     = sqrt (1.0 - cphi * cphi); // sin (acos (cphi));
-
-         return vec3 (sin (theta) * r, cos (theta) * r, cphi);
-      }
-
-      vec3
-      getRandomNormalWithAngle (const in float angle)
-      {
-         float theta = getRandomValue (-M_PI, M_PI);
-         float cphi  = getRandomValue (cos (angle), 1.0);
-         float r     = sqrt (1.0 - cphi * cphi); // sin (acos (cphi));
-
-         return vec3 (sin (theta) * r, cos (theta) * r, cphi);
-      }
-
-      vec3
-      getRandomNormalWithDirectionAndAngle (const in vec3 direction, const in float angle)
-      {
-         vec4 rotation = Quaternion (vec3 (0.0, 0.0, 1.0), direction);
-         vec3 normal   = getRandomNormalWithAngle (angle);
-
-         return multVecQuat (normal, rotation);
-      }
-
-      vec3
-      getRandomSurfaceNormal (const in vec3 direction)
-      {
-         float theta    = getRandomValue (-M_PI, M_PI);
-         float cphi     = pow (random (), 1.0 / 3.0);
-         float r        = sqrt (1.0 - cphi * cphi); // sin (acos (cphi));
-         vec3  normal   = vec3 (sin (theta) * r, cos (theta) * r, cphi);
-         vec4  rotation = Quaternion (vec3 (0.0, 0.0, 1.0), direction);
-
-         return multVecQuat (normal, rotation);
-      }
-
-      vec3
-      getRandomSphericalVelocity ()
-      {
-         vec3  normal = getRandomNormal ();
-         float speed  = getRandomSpeed ();
-
-         return normal * speed;
-      }
-
-      // Algorithms
-
-      int
-      upperBound (const in sampler2D sampler, in int count, const in float value)
-      {
-         int first = 0;
-         int step  = 0;
-
-         while (count > 0)
-         {
-            int index = first;
-
-            step = count >> 1;
-
-            index += step;
-
-            if (value < texelFetch (sampler, index, 0) .x)
-            {
-               count = step;
-            }
-            else
-            {
-               first  = ++ index;
-               count -= step + 1;
-            }
-         }
-
-         return first;
-      }
-
-      void
-      interpolate (const in sampler2D sampler, const in int count, const in float fraction, out int index0, out int index1, out float weight)
-      {
-         // Determine index0, index1 and weight.
-
-         if (count == 1 || fraction <= texelFetch (sampler, 0, 0) .x)
-         {
-            index0 = 0;
-            index1 = 0;
-            weight = 0.0;
-         }
-         else if (fraction >= texelFetch (sampler, count - 1, 0) .x)
-         {
-            index0 = count - 2;
-            index1 = count - 1;
-            weight = 1.0;
-         }
-         else
-         {
-            int index = upperBound (sampler, count, fraction);
-
-            if (index < count)
-            {
-               index1 = index;
-               index0 = index - 1;
-
-               float key0 = texelFetch (sampler, index0, 0) .x;
-               float key1 = texelFetch (sampler, index1, 0) .x;
-
-               weight = clamp ((fraction - key0) / (key1 - key0), 0.0, 1.0);
-            }
-            else
-            {
-               index0 = 0;
-               index1 = 0;
-               weight = 0.0;
-            }
-         }
-      }
-
-      void
-      interpolate (const in sampler2D sampler, const in int count, const in float fraction, out int index0)
-      {
-         // Determine index0.
-
-         if (count == 1 || fraction <= texelFetch (sampler, 0, 0) .x)
-         {
-            index0 = 0;
-         }
-         else if (fraction >= texelFetch (sampler, count - 1, 0) .x)
-         {
-            index0 = count - 2;
-         }
-         else
-         {
-            int index = upperBound (sampler, count, fraction);
-
-            if (index < count)
-               index0 = index - 1;
-            else
-               index0 = 0;
-         }
-      }
-
-      vec3
-      getRandomBarycentricCoord ()
-      {
-         // Random barycentric coordinates.
-
-         float u = random ();
-         float v = random ();
-
-         if (u + v > 1.0)
-         {
-            u = 1.0 - u;
-            v = 1.0 - v;
-         }
-
-         float t = 1.0 - u - v;
-
-         return vec3 (t, u, v);
-      }
-
-      void
-      getRandomPointOnSurface (const in sampler2D surface, const in int verticesIndex, const in int normalsIndex, out vec4 position, out vec3 normal)
-      {
-         // Determine index0, index1 and weight.
-
-         float lastAreaSoFar = texelFetch (surface, verticesIndex - 1, 0) .x;
-         float fraction      = random () * lastAreaSoFar;
-
-         int   index0;
-         int   index1;
-         int   index2;
-         float weight;
-
-         interpolate (surface, verticesIndex, fraction, index0, index1, weight);
-
-         // Interpolate and return position.
-
-         index0 *= 3;
-         index1  = index0 + 1;
-         index2  = index0 + 2;
-
-         vec4 vertex0 = texelFetch (surface, verticesIndex + index0, 0);
-         vec4 vertex1 = texelFetch (surface, verticesIndex + index1, 0);
-         vec4 vertex2 = texelFetch (surface, verticesIndex + index2, 0);
-
-         vec3 normal0 = texelFetch (surface, normalsIndex + index0, 0) .xyz;
-         vec3 normal1 = texelFetch (surface, normalsIndex + index1, 0) .xyz;
-         vec3 normal2 = texelFetch (surface, normalsIndex + index2, 0) .xyz;
-
-         // Random barycentric coordinates.
-
-         vec3 r = getRandomBarycentricCoord ();
-
-         // Calculate position and direction.
-
-         position = r .z * vertex0 + r .x * vertex1 + r .y * vertex2;
-         normal   = save_normalize (r .z * normal0 + r .x * normal1 + r .y * normal2);
-      }
-
-      // Functions
-
-      ${this .functions .join ("\n")}
-
-      // Current values
-
-      vec4
-      getColor (const in float lifetime, const in float elapsedTime)
-      {
-         if (numColors > 0)
-         {
-            // Determine index0, index1 and weight.
-
-            float fraction = elapsedTime / lifetime;
-
-            int   index0;
-            int   index1;
-            float weight;
-
-            interpolate (colorRamp, numColors, fraction, index0, index1, weight);
-
-            // Interpolate and return color.
-
-            vec4 color0 = texelFetch (colorRamp, numColors + index0, 0);
-            vec4 color1 = texelFetch (colorRamp, numColors + index1, 0);
-
-            return mix (color0, color1, weight);
-         }
-         else
-         {
-            return vec4 (1.0);
-         }
-      }
-
-      void
-      bounce (const in vec4 fromPosition, inout vec4 toPosition, inout vec3 velocity)
-      {
-         if (boundedHierarchyRoot < 0)
-            return;
-
-         Line3 line = Line3 (fromPosition .xyz, save_normalize (velocity));
-
-         vec4 points  [ARRAY_SIZE];
-         vec3 normals [ARRAY_SIZE];
-
-         int numIntersections = getIntersections (boundedVolume, boundedVerticesIndex, boundedNormalsIndex, boundedHierarchyIndex, boundedHierarchyRoot, line, points, normals);
-
-         if (numIntersections == 0)
-            return;
-
-         Plane3 plane1 = plane3 (line .point, line .direction);
-
-         int index = min_index (points, numIntersections, 0.0, plane1);
-
-         if (index == -1)
-            return;
-
-         Plane3 plane2 = plane3 (points [index] .xyz, normals [index]);
-
-         if (sign (plane_distance (plane2, fromPosition .xyz)) == sign (plane_distance (plane2, toPosition .xyz)))
-            return;
-
-         velocity   = reflect (velocity, normals [index]);
-         toPosition = vec4 (points [index] .xyz + reflect (points [index] .xyz - fromPosition .xyz, normals [index]), 1.0);
-      }
-
-      int
-      getTexCoordIndex0 (const in float lifetime, const in float elapsedTime)
-      {
-         if (numTexCoords == 0)
-         {
-            return -1;
-         }
-         else
-         {
-            float fraction = elapsedTime / lifetime;
-            int   index0   = 0;
-
-            interpolate (texCoordRamp, numTexCoords, fraction, index0);
-
-            return numTexCoords + index0 * texCoordCount;
-         }
-      }
-
-      void
-      main ()
-      {
-         int   life        = int (input0 [0]);
-         float lifetime    = input0 [1];
-         float elapsedTime = input0 [2] + deltaTime;
-
-         srand ((gl_VertexID + randomSeed) * randomSeed);
-
-         if (elapsedTime > lifetime)
-         {
-            // Create new particle or hide particle.
-
-            lifetime    = getRandomLifetime ();
-            elapsedTime = 0.0;
-
-            output0 = vec4 (max (life + 1, 1), lifetime, elapsedTime, getTexCoordIndex0 (lifetime, elapsedTime));
-
-            if (createParticles)
-            {
-               output1 = getColor (lifetime, elapsedTime);
-               output2 = vec4 (getRandomVelocity (), 0.0);
-               output6 = getRandomPosition ();
-            }
-            else
-            {
-               output1 = vec4 (0.0);
-               output2 = vec4 (0.0);
-               output6 = vec4 (NaN);
-            }
-         }
-         else
-         {
-            // Animate particle.
-
-            vec3 velocity = input2 .xyz;
-            vec4 position = input6;
-
-            for (int i = 0; i < numForces; ++ i)
-            {
-               vec4  force      = texelFetch (forces, i, 0);
-               float turbulence = force .w;
-               vec3  normal     = getRandomNormalWithDirectionAndAngle (force .xyz, turbulence);
-               float speed      = length (force .xyz);
-
-               velocity += normal * speed;
-            }
-
-            position .xyz += velocity * deltaTime;
-
-            bounce (input6, position, velocity);
-
-            output0 = vec4 (life, lifetime, elapsedTime, getTexCoordIndex0 (lifetime, elapsedTime));
-            output1 = getColor (lifetime, elapsedTime);
-            output2 = vec4 (velocity, 0.0);
-            output6 = position;
-         }
-
-         switch (geometryType)
-         {
-            case POINT:
-            case SPRITE:
-            case GEOMETRY:
-            {
-               output3 = vec4 (1.0, 0.0, 0.0, 0.0);
-               output4 = vec4 (0.0, 1.0, 0.0, 0.0);
-               output5 = vec4 (0.0, 0.0, 1.0, 0.0);
-               break;
-            }
-            case LINE:
-            {
-               mat3 r = Matrix3 (Quaternion (vec3 (0.0, 0.0, 1.0), output2 .xyz));
-               mat3 s = mat3 (1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, particleSize .y);
-               mat3 m = r * s;
-
-               output3 = vec4 (m [0], 0.0);
-               output4 = vec4 (m [1], 0.0);
-               output5 = vec4 (m [2], 0.0);
-               break;
-            }
-            default: // QUAD, TRIANGLE
-            {
-               output3 = vec4 (particleSize .x, 0.0, 0.0, 0.0);
-               output4 = vec4 (0.0, particleSize .y, 0.0, 0.0);
-               output5 = vec4 (0.0, 0.0, 1.0, 0.0);
-               break;
-            }
-         }
-      }
-      `;
+precision highp float;precision highp int;precision highp sampler2D;uniform int randomSeed;uniform int geometryType;uniform bool createParticles;uniform float particleLifetime;uniform float lifetimeVariation;uniform float deltaTime;uniform vec2 particleSize;uniform int numForces;uniform sampler2D forces;uniform int boundedVerticesIndex;uniform int boundedNormalsIndex;uniform int boundedHierarchyIndex;uniform int boundedHierarchyRoot;uniform sampler2D boundedVolume;uniform int numColors;uniform sampler2D colorRamp;uniform int texCoordCount;uniform int numTexCoords;uniform sampler2D texCoordRamp;
+${Object .values (this .uniforms) .join ("\n")}
+ in vec4 input0;in vec4 input2;in vec4 input6;out vec4 output0;out vec4 output1;out vec4 output2;out vec4 output3;out vec4 output4;out vec4 output5;out vec4 output6;
+${Object .entries (ParticleSystems_GeometryTypes) .map (([k, v]) => `#define ${k} ${v}`) .join ("\n")}
+ const int ARRAY_SIZE=32;const float M_PI=3.14159265359;uniform float NaN;vec4 texelFetch(const in sampler2D sampler,const in int index,const in int lod){int x=textureSize(sampler,lod).x;ivec2 p=ivec2(index % x,index/x);vec4 t=texelFetch(sampler,p,lod);return t;}vec3 save_normalize(const in vec3 vector){float l=length(vector);if(l==0.0)return vec3(0.0);return vector/l;}vec4 Quaternion(const in vec3 fromVector,const in vec3 toVector){vec3 from=save_normalize(fromVector);vec3 to=save_normalize(toVector);float cos_angle=dot(from,to);vec3 cross_vec=cross(from,to);float cross_len=length(cross_vec);if(cross_len==0.0){if(cos_angle>0.0){return vec4(0.0,0.0,0.0,1.0);}else{vec3 t=cross(from,vec3(1.0,0.0,0.0));if(dot(t,t)==0.0)t=cross(from,vec3(0.0,1.0,0.0));t=save_normalize(t);return vec4(t,0.0);}}else{float s=sqrt(abs(1.0-cos_angle)*0.5);cross_vec=save_normalize(cross_vec);return vec4(cross_vec*s,sqrt(abs(1.0+cos_angle)*0.5));}}vec3 multVecQuat(const in vec3 v,const in vec4 q){float a=q.w*q.w-q.x*q.x-q.y*q.y-q.z*q.z;float b=2.0*(v.x*q.x+v.y*q.y+v.z*q.z);float c=2.0*q.w;vec3 r=a*v.xyz+b*q.xyz+c*(q.yzx*v.zxy-q.zxy*v.yzx);return r;}mat3 Matrix3(const in vec4 quaternion){float x=quaternion.x;float y=quaternion.y;float z=quaternion.z;float w=quaternion.w;float A=y*y;float B=z*z;float C=x*y;float D=z*w;float E=z*x;float F=y*w;float G=x*x;float H=y*z;float I=x*w;return mat3(1.0-2.0*(A+B),2.0*(C+D),2.0*(E-F),2.0*(C-D),1.0-2.0*(B+G),2.0*(H+I),2.0*(E+F),2.0*(H-I),1.0-2.0*(A+G));}uint seed=1u;void srand(const in int value){seed=uint(value);}float random(){seed=seed*1103515245u+12345u;return float(seed)/4294967295.0;}float getRandomValue(const in float min,const in float max){return min+random()*(max-min);}float getRandomLifetime(){float v=particleLifetime*lifetimeVariation;float min_=max(0.0,particleLifetime-v);float max_=particleLifetime+v;return getRandomValue(min_,max_);}float getRandomSpeed(){float v=speed*variation;float min_=max(0.0,speed-v);float max_=speed+v;return getRandomValue(min_,max_);}vec3 getRandomNormal(){float theta=getRandomValue(-M_PI,M_PI);float cphi=getRandomValue(-1.0,1.0);float r=sqrt(1.0-cphi*cphi);return vec3(sin(theta)*r,cos(theta)*r,cphi);}vec3 getRandomNormalWithAngle(const in float angle){float theta=getRandomValue(-M_PI,M_PI);float cphi=getRandomValue(cos(angle),1.0);float r=sqrt(1.0-cphi*cphi);return vec3(sin(theta)*r,cos(theta)*r,cphi);}vec3 getRandomNormalWithDirectionAndAngle(const in vec3 direction,const in float angle){vec4 rotation=Quaternion(vec3(0.0,0.0,1.0),direction);vec3 normal=getRandomNormalWithAngle(angle);return multVecQuat(normal,rotation);}vec3 getRandomSurfaceNormal(const in vec3 direction){float theta=getRandomValue(-M_PI,M_PI);float cphi=pow(random(),1.0/3.0);float r=sqrt(1.0-cphi*cphi);vec3 normal=vec3(sin(theta)*r,cos(theta)*r,cphi);vec4 rotation=Quaternion(vec3(0.0,0.0,1.0),direction);return multVecQuat(normal,rotation);}vec3 getRandomSphericalVelocity(){vec3 normal=getRandomNormal();float speed=getRandomSpeed();return normal*speed;}int upperBound(const in sampler2D sampler,in int count,const in float value){int first=0;int step=0;while(count>0){int index=first;step=count>>1;index+=step;if(value<texelFetch(sampler,index,0).x){count=step;}else{first=++index;count-=step+1;}}return first;}void interpolate(const in sampler2D sampler,const in int count,const in float fraction,out int index0,out int index1,out float weight){if(count==1||fraction<=texelFetch(sampler,0,0).x){index0=0;index1=0;weight=0.0;}else if(fraction>=texelFetch(sampler,count-1,0).x){index0=count-2;index1=count-1;weight=1.0;}else{int index=upperBound(sampler,count,fraction);if(index<count){index1=index;index0=index-1;float key0=texelFetch(sampler,index0,0).x;float key1=texelFetch(sampler,index1,0).x;weight=clamp((fraction-key0)/(key1-key0),0.0,1.0);}else{index0=0;index1=0;weight=0.0;}}}void interpolate(const in sampler2D sampler,const in int count,const in float fraction,out int index0){if(count==1||fraction<=texelFetch(sampler,0,0).x){index0=0;}else if(fraction>=texelFetch(sampler,count-1,0).x){index0=count-2;}else{int index=upperBound(sampler,count,fraction);if(index<count)index0=index-1;else index0=0;}}vec3 getRandomBarycentricCoord(){float u=random();float v=random();if(u+v>1.0){u=1.0-u;v=1.0-v;}float t=1.0-u-v;return vec3(t,u,v);}void getRandomPointOnSurface(const in sampler2D surface,const in int verticesIndex,const in int normalsIndex,out vec4 position,out vec3 normal){float lastAreaSoFar=texelFetch(surface,verticesIndex-1,0).x;float fraction=random()*lastAreaSoFar;int index0;int index1;int index2;float weight;interpolate(surface,verticesIndex,fraction,index0,index1,weight);index0*=3;index1=index0+1;index2=index0+2;vec4 vertex0=texelFetch(surface,verticesIndex+index0,0);vec4 vertex1=texelFetch(surface,verticesIndex+index1,0);vec4 vertex2=texelFetch(surface,verticesIndex+index2,0);vec3 normal0=texelFetch(surface,normalsIndex+index0,0).xyz;vec3 normal1=texelFetch(surface,normalsIndex+index1,0).xyz;vec3 normal2=texelFetch(surface,normalsIndex+index2,0).xyz;vec3 r=getRandomBarycentricCoord();position=r.z*vertex0+r.x*vertex1+r.y*vertex2;normal=save_normalize(r.z*normal0+r.x*normal1+r.y*normal2);}
+${this .functions .join ("\n")}
+ vec4 getColor(const in float lifetime,const in float elapsedTime){if(numColors>0){float fraction=elapsedTime/lifetime;int index0;int index1;float weight;interpolate(colorRamp,numColors,fraction,index0,index1,weight);vec4 color0=texelFetch(colorRamp,numColors+index0,0);vec4 color1=texelFetch(colorRamp,numColors+index1,0);return mix(color0,color1,weight);}else{return vec4(1.0);}}void bounce(const in vec4 fromPosition,inout vec4 toPosition,inout vec3 velocity){if(boundedHierarchyRoot<0)return;Line3 line=Line3(fromPosition.xyz,save_normalize(velocity));vec4 points[ARRAY_SIZE];vec3 normals[ARRAY_SIZE];int numIntersections=getIntersections(boundedVolume,boundedVerticesIndex,boundedNormalsIndex,boundedHierarchyIndex,boundedHierarchyRoot,line,points,normals);if(numIntersections==0)return;Plane3 plane1=plane3(line.point,line.direction);int index=min_index(points,numIntersections,0.0,plane1);if(index==-1)return;Plane3 plane2=plane3(points[index].xyz,normals[index]);if(sign(plane_distance(plane2,fromPosition.xyz))==sign(plane_distance(plane2,toPosition.xyz)))return;velocity=reflect(velocity,normals[index]);toPosition=vec4(points[index].xyz+reflect(points[index].xyz-fromPosition.xyz,normals[index]),1.0);}int getTexCoordIndex0(const in float lifetime,const in float elapsedTime){if(numTexCoords==0){return-1;}else{float fraction=elapsedTime/lifetime;int index0=0;interpolate(texCoordRamp,numTexCoords,fraction,index0);return numTexCoords+index0*texCoordCount;}}void main(){int life=int(input0[0]);float lifetime=input0[1];float elapsedTime=input0[2]+deltaTime;srand((gl_VertexID+randomSeed)*randomSeed);if(elapsedTime>lifetime){lifetime=getRandomLifetime();elapsedTime=0.0;output0=vec4(max(life+1,1),lifetime,elapsedTime,getTexCoordIndex0(lifetime,elapsedTime));if(createParticles){output1=getColor(lifetime,elapsedTime);output2=vec4(getRandomVelocity(),0.0);output6=getRandomPosition();}else{output1=vec4(0.0);output2=vec4(0.0);output6=vec4(NaN);}}else{vec3 velocity=input2.xyz;vec4 position=input6;for(int i=0;i<numForces;++i){vec4 force=texelFetch(forces,i,0);float turbulence=force.w;vec3 normal=getRandomNormalWithDirectionAndAngle(force.xyz,turbulence);float speed=length(force.xyz);velocity+=normal*speed;}position.xyz+=velocity*deltaTime;bounce(input6,position,velocity);output0=vec4(life,lifetime,elapsedTime,getTexCoordIndex0(lifetime,elapsedTime));output1=getColor(lifetime,elapsedTime);output2=vec4(velocity,0.0);output6=position;}switch(geometryType){case POINT:case SPRITE:case GEOMETRY:{output3=vec4(1.0,0.0,0.0,0.0);output4=vec4(0.0,1.0,0.0,0.0);output5=vec4(0.0,0.0,1.0,0.0);break;}case LINE:{mat3 r=Matrix3(Quaternion(vec3(0.0,0.0,1.0),output2.xyz));mat3 s=mat3(1.0,0.0,0.0,0.0,1.0,0.0,0.0,0.0,particleSize.y);mat3 m=r*s;output3=vec4(m[0],0.0);output4=vec4(m[1],0.0);output5=vec4(m[2],0.0);break;}default:{output3=vec4(particleSize.x,0.0,0.0,0.0);output4=vec4(0.0,particleSize.y,0.0,0.0);output5=vec4(0.0,0.0,1.0,0.0);break;}}}
+`
 
       const fragmentShaderSource = /* glsl */ `#version 300 es
-
-      precision highp float;
-
-      void
-      main () { }
-      `;
+precision highp float;void main(){}
+`
 
       // Vertex shader
 
@@ -1515,6 +498,9 @@ Object .assign (Object .setPrototypeOf (X3DParticleEmitterNode .prototype, (X3DN
       gl .shaderSource (vertexShader, vertexShaderSource);
       gl .compileShader (vertexShader);
 
+      if (!gl .getShaderParameter (vertexShader, gl .COMPILE_STATUS))
+         console .error (gl .getShaderInfoLog (vertexShader));
+
       // Fragment shader
 
       const fragmentShader = gl .createShader (gl .FRAGMENT_SHADER);
@@ -1522,6 +508,9 @@ Object .assign (Object .setPrototypeOf (X3DParticleEmitterNode .prototype, (X3DN
       gl .shaderSource (fragmentShader, fragmentShaderSource);
       gl .compileShader (fragmentShader);
 
+      if (!gl .getShaderParameter (fragmentShader, gl .COMPILE_STATUS))
+         console .error (gl .getShaderInfoLog (fragmentShader));
+
       // Program
 
       const program = gl .createProgram ();
@@ -1634,7 +623,7 @@ const X3DParticleEmitterNode_default_ = X3DParticleEmitterNode;
 Namespace_default().set ("x_ite/Components/ParticleSystems/X3DParticleEmitterNode", X3DParticleEmitterNode_default_);
 /* harmony default export */ const ParticleSystems_X3DParticleEmitterNode = (X3DParticleEmitterNode_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Numbers/Vector3\")"
-const Vector3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Numbers/Vector3");
+const Vector3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Numbers/Vector3");
 var Vector3_default = /*#__PURE__*/__webpack_require__.n(Vector3_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/PointEmitter.js
 /*******************************************************************************
@@ -1955,7 +944,7 @@ const X3DParticlePhysicsModelNode_default_ = X3DParticlePhysicsModelNode;
 Namespace_default().set ("x_ite/Components/ParticleSystems/X3DParticlePhysicsModelNode", X3DParticlePhysicsModelNode_default_);
 /* harmony default export */ const ParticleSystems_X3DParticlePhysicsModelNode = (X3DParticlePhysicsModelNode_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Base/X3DCast\")"
-const X3DCast_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Base/X3DCast");
+const X3DCast_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Base/X3DCast");
 var X3DCast_default = /*#__PURE__*/__webpack_require__.n(X3DCast_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/BoundedPhysicsModel.js
 /*******************************************************************************
@@ -2529,28 +1518,28 @@ const ForcePhysicsModel_default_ = ForcePhysicsModel;
 Namespace_default().set ("x_ite/Components/ParticleSystems/ForcePhysicsModel", ForcePhysicsModel_default_);
 /* harmony default export */ const ParticleSystems_ForcePhysicsModel = (ForcePhysicsModel_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Components/Shape/X3DShapeNode\")"
-const X3DShapeNode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Components/Shape/X3DShapeNode");
+const X3DShapeNode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Components/Shape/X3DShapeNode");
 var X3DShapeNode_default = /*#__PURE__*/__webpack_require__.n(X3DShapeNode_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Browser/Rendering/GeometryContext\")"
-const GeometryContext_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Browser/Rendering/GeometryContext");
+const GeometryContext_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Browser/Rendering/GeometryContext");
 var GeometryContext_default = /*#__PURE__*/__webpack_require__.n(GeometryContext_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Rendering/VertexArray\")"
-const VertexArray_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Rendering/VertexArray");
+const VertexArray_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Rendering/VertexArray");
 var VertexArray_default = /*#__PURE__*/__webpack_require__.n(VertexArray_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Rendering/TraverseType\")"
-const TraverseType_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Rendering/TraverseType");
+const TraverseType_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Rendering/TraverseType");
 var TraverseType_default = /*#__PURE__*/__webpack_require__.n(TraverseType_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Browser/Shape/AlphaMode\")"
-const AlphaMode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Browser/Shape/AlphaMode");
+const AlphaMode_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Browser/Shape/AlphaMode");
 var AlphaMode_default = /*#__PURE__*/__webpack_require__.n(AlphaMode_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Numbers/Matrix4\")"
-const Matrix4_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Numbers/Matrix4");
+const Matrix4_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Numbers/Matrix4");
 var Matrix4_default = /*#__PURE__*/__webpack_require__.n(Matrix4_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Numbers/Matrix3\")"
-const Matrix3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Numbers/Matrix3");
+const Matrix3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Numbers/Matrix3");
 var Matrix3_default = /*#__PURE__*/__webpack_require__.n(Matrix3_namespaceObject);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Geometry/Plane3\")"
-const Plane3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Geometry/Plane3");
+const Plane3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Geometry/Plane3");
 var Plane3_default = /*#__PURE__*/__webpack_require__.n(Plane3_namespaceObject);
 ;// CONCATENATED MODULE: ./src/standard/Math/Algorithms/QuickSort.js
 /*******************************************************************************
@@ -4134,7 +3123,7 @@ const ParticleSystem_default_ = ParticleSystem;
 Namespace_default().set ("x_ite/Components/ParticleSystems/ParticleSystem", ParticleSystem_default_);
 /* harmony default export */ const ParticleSystems_ParticleSystem = (ParticleSystem_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Components/Rendering/IndexedLineSet\")"
-const IndexedLineSet_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Components/Rendering/IndexedLineSet");
+const IndexedLineSet_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Components/Rendering/IndexedLineSet");
 var IndexedLineSet_default = /*#__PURE__*/__webpack_require__.n(IndexedLineSet_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/PolylineEmitter.js
 /*******************************************************************************
@@ -4386,7 +3375,7 @@ const PolylineEmitter_default_ = PolylineEmitter;
 Namespace_default().set ("x_ite/Components/ParticleSystems/PolylineEmitter", PolylineEmitter_default_);
 /* harmony default export */ const ParticleSystems_PolylineEmitter = (PolylineEmitter_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Geometry/Triangle3\")"
-const Triangle3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Geometry/Triangle3");
+const Triangle3_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Geometry/Triangle3");
 var Triangle3_default = /*#__PURE__*/__webpack_require__.n(Triangle3_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/SurfaceEmitter.js
 /*******************************************************************************
@@ -4646,7 +3635,7 @@ const SurfaceEmitter_default_ = SurfaceEmitter;
 Namespace_default().set ("x_ite/Components/ParticleSystems/SurfaceEmitter", SurfaceEmitter_default_);
 /* harmony default export */ const ParticleSystems_SurfaceEmitter = (SurfaceEmitter_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"x_ite/Components/Geometry3D/IndexedFaceSet\")"
-const IndexedFaceSet_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("x_ite/Components/Geometry3D/IndexedFaceSet");
+const IndexedFaceSet_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("x_ite/Components/Geometry3D/IndexedFaceSet");
 var IndexedFaceSet_default = /*#__PURE__*/__webpack_require__.n(IndexedFaceSet_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/VolumeEmitter.js
 /*******************************************************************************
@@ -4943,7 +3932,7 @@ const VolumeEmitter_default_ = VolumeEmitter;
 Namespace_default().set ("x_ite/Components/ParticleSystems/VolumeEmitter", VolumeEmitter_default_);
 /* harmony default export */ const ParticleSystems_VolumeEmitter = (VolumeEmitter_default_);
 ;// CONCATENATED MODULE: external "window [Symbol .for (\"X_ITE.X3D\")] .require (\"standard/Math/Algorithm\")"
-const Algorithm_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.5")] .require ("standard/Math/Algorithm");
+const Algorithm_namespaceObject = window [Symbol .for ("X_ITE.X3D-8.11.6")] .require ("standard/Math/Algorithm");
 var Algorithm_default = /*#__PURE__*/__webpack_require__.n(Algorithm_namespaceObject);
 ;// CONCATENATED MODULE: ./src/x_ite/Components/ParticleSystems/WindPhysicsModel.js
 /*******************************************************************************