2(p.x) + u32(m.y)*rez + u32(m.z)*rez*rez;\\n var w = u32(p.x);\\n workgroupBarrier();\\n C2[w] = C[r];\\n workgroupBarrier();\\n}\\n#workgroup_count cnv 1 rez rez\\n@compute @workgroup_size(rez,1,1)\\nfn cnv(@builtin(global_invocation_id) id3: vec3u)\\n{\\n var id = c31(id3);\\n var d = D[id]; //voxel of simulationCube\\n var p = vec3f(id3)+.5f;\\n var v = 0f;\\n var et = 0f;\\n for(var z=-rZ; z<=rZ; z+=1f){\\n for(var y=-rZ; y<=rZ; y+=1f){ loadLineD(p, vec3f(0,y,z));\\n for(var x=-rZ; x<=rZ; x+=1f){\\n var xyz = vec3f(x,y,z);\\n var l2 = dot(xyz,xyz);\\n if(l2<.1f || l2>rZ*rZ+.1f){continue;}\\n var i2 = u32(mod1(p.x+x , f32(rez)));\\n var v0 = D2[i2]; //neighbor vector\\n var l = sqrt(l2);\\n var xyzl = xyz/l;\\n v0 = v0-d;\\n var e = wave(l); et += abs(e);\\n v += e*length(v0-xyzl*dot(v0,xyzl));//rotations\\n }}}\\n C[id] = v/et;\\n}\\n#workgroup_count fun 1 rez rez\\n@compute @workgroup_size(rez,1,1)\\nfn fun(@builtin(global_invocation_id) id3: vec3u)\\n{\\n var id = c31(id3);\\n var d = D[id]; //voxel of simulationCube\\n var p = vec3f(id3)+.5f;\\n var v = vec3f(0);\\n var et = 0f;\\n for(var z=-rZ; z<=rZ; z+=1f){\\n for(var y=-rZ; y<=rZ; y+=1f){ loadLineC(p, vec3f(0,y,z));\\n for(var x=-rZ; x<=rZ; x+=1f){\\n var xyz = vec3f(x,y,z);\\n var l2 = dot(xyz,xyz);\\n if(l2<.1f || l2>rZ*rZ+.1f){continue;}\\n var i2 = u32(mod1(p.x+x , f32(rez)));\\n var c = C2[i2];\\n var l = sqrt(l2);\\n var xyzl = xyz/l;\\n var e = wave(l); et += abs(e);\\n v += e*xyzl*c;\\n }}}\\n var u = (p/f32(rez)-.5f)*f32(time.frame)*.01f;\\n v = v/et*custom.c*33f +\\n u/exp(dot(u,u))*3f;\\n\\n D[id] += (v-D[id])*.4f;\\n}\\n@compute @workgroup_size(8,8,1)\\nfn main_image(@builtin(global_invocation_id) id: vec3u)\\n{\\n var screen_size = textureDimensions(screen);\\n if(id.x >= screen_size.x || id.y >= screen_size.y){ return; }\\n var fragCoord = float2(id.xy) + .5f;\\n var iResolution = float2(screen_size);\\n var clk = f32(mouse.click);\\n var iTime = time.elapsed;\\n\\n var u = (2f*fragCoord -iResolution)/iResolution.y;\\n var m = (2f*float2(mouse.pos)-iResolution)/iResolution.y*2f;\\n var camPos = (1f-clk)*vec3f(cos(sin(iTime*.01f)*3.f),\\n cos(iTime*.23f)*.3f,\\n sin(cos(iTime*.07f)*3.f))\\n + clk*vec3f(cos(m.x),m.y,sin(m.x));\\n camPos = normalize(camPos)*3f;\\n var camDir = -normalize(camPos);\\n \\n var mtx0 = normalize(vec3f(camDir.z,0.f,-camDir.x));\\n var mtx = mat3x3(mtx0, cross(camDir,mtx0), camDir);\\n var ray = mtx*normalize(vec3f(u,2.f)); //direction of ray from camera\\n var ray2= 1.f/ray;\\n var ray3= step(vec3f(0),ray)*2.f-1.f;\\n \\n var x3 = f32(rez);\\n var x3d= 2.f/x3;\\n var x4 = rez*2u; //max voxels the ray will transverse inside simulationCube\\n \\n var inc = step(camPos,vec3f( 1))* //inside cube\\n step(vec3f(-1),camPos);\\n var tMin = (vec3f(-1)-camPos)*ray2;\\n var tMax = (vec3f( 1)-camPos)*ray2; ray2 = abs(ray2);\\n var t1 = min(tMin, tMax);\\n var t2 = max(tMin, tMax);\\n var tN = max(max(t1.x, t1.y), t1.z); //length of ray between camera and simulationCube\\n var tF = min(min(t2.x, t2.y), t2.z);\\n tF = f32(tF>tN);\\n \\n var p = camPos+ray*tN*(1.f-inc.x*inc.y*inc.z)*1.001f; //collision position of ray on simulationCube \\n p = (p*.5f+.5f)*x3; //transform to voxels coordinates, range 0 to x3\\n var lig = vec3f(1); //light created by camera going to voxels\\n var rif = vec3f(0); //reflected light by voxels going to camera\\n for(var i=0u; i=1f||o.y>=1f

@compute.toys