(\\n cos(t), sin(t), \\n - sin(t), cos(t), \\n );\\n}\\n\\nfn RotXY(x: float3, t: float) -> float3\\n{\\n return float3(Rotate(t)*x.xy, x.z);\\n} \\n\\nfn GetCameraMatrix(ang: float2) -> float3x3\\n{\\n let x_dir = float3(cos(ang.x)*sin(ang.y), cos(ang.y), sin(ang.x)*sin(ang.y));\\n let y_dir = normalize(cross(x_dir, float3(0.0,1.0,0.0)));\\n let z_dir = normalize(cross(y_dir, x_dir));\\n return float3x3(-x_dir, y_dir, z_dir);\\n}\\n\\n\\n//project to clip space\\nfn Project(cam: Camera, p: float3) -> float3\\n{\\n let td = distance(cam.pos, p);\\n let dir = (p - cam.pos)/td;\\n let screen = dir*cam.cam;\\n return float3(screen.xy*cam.size.y/(cam.fov*screen.z) + 0.5*cam.size,screen.z*td);\\n}\\n\\nconst max_iterations = 256;\\nconst color_thresholds = float4(255.0, 130.0, 80.0, 255.0);\\n\\n\\nfn AdditiveBlend(color: float3, projected: float2)\\n{\\n let screen_size = int2(camera.size);\\n let screen_coord = int2(projected.xy);\\n let index = screen_coord.x + screen_size.x * screen_coord.y;\\n let scaledColor = int3(floor(ATOMIC_SCALE*color + rand4().xyz));\\n\\n if(scaledColor.x>0)\\n {\\n atomicAdd(&atomic_storage[index*4+0], scaledColor.x);\\n }\\n \\n if(scaledColor.y>0)\\n {\\n atomicAdd(&atomic_storage[index*4+1], scaledColor.y);\\n }\\n\\n if(scaledColor.z>0)\\n {\\n atomicAdd(&atomic_storage[index*4+2], scaledColor.z);\\n }\\n}\\n\\nfn RasterizePoint(pos: float3, color: float3)\\n{\\n let screen_size = int2(camera.size);\\n let projectedPos = Project(camera, pos);\\n let screenCoord = int2(projectedPos.xy);\\n\\n //outside of our view\\n if(screenCoord.x < 0 || screenCoord.x >= screen_size.x || \\n screenCoord.y < 0 || screenCoord.y >= screen_size.y || projectedPos.z < 0.0)\\n {\\n return;\\n }\\n\\n let depth = projectedPos.z;\\n let new_color = color * exp(-0.0*depth)/(depth*depth+1e-8);\\n\\n var rad = clamp(0.1/depth, 1.0, 8.0);\\n var radi = round(rad);\\n rad = max(rad*0.15, 0.33);\\n if(depth < 0.2) {radi += 1.0;}\\n\\n let screenCoordf = float2(screenCoord);\\n \\n for(var x = -radi; x <= radi; x+=1.0)\\n {\\n for(var y = -radi; y <= radi; y+=1.0)\\n {\\n let dx = vec2(x,y);\\n let projPos = screenCoordf + dx;\\n let dist = distance(projPos, projectedPos.xy)/rad;\\n let w = exp(-dist)/(rad*rad);\\n AdditiveBlend(new_color*w, projPos);\\n }\\n }\\n\\n}\\n\\n\\nfn saturate(x: f32) -> f32 {\\n return min(1.0, max(0.0, x));\\n}\\n\\nfn saturate_vec3(x: vec3) -> vec3 {\\n return min(vec3(1.0, 1.0, 1.0), max(vec3(0.0, 0.0, 0.0), x));\\n}\\n\\nfn bump3y(x: vec3, yoffset: vec3) -> vec3 {\\n var y: vec3 = vec3(1.0, 1.0, 1.0) - x * x;\\n y = saturate_vec3(y - yoffset);\\n return y;\\n}\\n\\nfn spectral_zucconi(w: f32) -> vec3 {\\n let x: f32 = saturate((w - 400.0) / 300.0);\\n\\n let cs: vec3 = vec3(3.54541723, 2.86670055, 2.29421995);\\n let xs: vec3 = vec3(0.69548916, 0.49416934, 0.28269708);\\n let ys: vec3 = vec3(0.02320775, 0.15936245, 0.53520021);\\n\\n return bump3y(cs * (x - xs), ys);\\n}\\n\\nfn hue(v: float) -> float3 {\\n return .6 + .6 * cos(6.3 * v + float3(0.,23.,21.));\\n}\\n\\nfn buddhabrot(iters: i32) -> vec4 {\\n var z: vec3 = vec3(0.01, 0.01, 0.01);\\n var c: vec3 = (rand4().xyz - 0.5) * 1.2;\\n\\n var i: i32 = 0;\\n let BULB_POWER = custom.BulbPower * BULB_MAX_POWER + 5.0*(0.5*sin(0.25*custom.AnimationSpeed*time.elapsed)+0.5);\\n let BULB_POWER_DELTA = (custom.PowerDelta - 0.5) * BULB_POWER_DELTA_MAX;\\n loop {\\n if (i >= iters) { break; }\\n\\n let r: f32 = length(z);\\n let b: f32 = (BULB_POWER + BULB_POWER_DELTA) * acos(z.y / r);\\n let a: f32 =(BULB_POWER + BULB_POWER_DELTA) * atan2(z.x, z.z);\\n z = c + pow(r, (BULB_POWER + BULB_POWER_DELTA)) * vec3(sin(b

@compute.toys