-pow((l - vec3f(590, 550, 500))/vec3f(50,70,50), vec3f(2)))\\n // + exp(-pow((l - vec3f(460, 0, 0))/vec3f(40,1,1), vec3f(2)))*.2\\n // ,vec3f(2.2)\\n // );\\n return spectral_zucconi(l);\\n}\\n\\nfn sd_sunlight(l: f32) -> float { // spectral distribution, blackbody radiation formula. Might be accurate?\\n let f = 170.;\\n // magic numbers trying to fit the curve\\n return pow(l/f, 3.) / (exp(l/f) - 1.);\\n}\\n\\nfn refr(\\n i: vec3f, // Incoming vector\\n n: vec3f, // Normal vector, pointing from high to low ior\\n mu: f32, // Index of refraction ratio high/low\\n tir: ptr // Total internal reflection\\n) -> vec3f {\\n\\n *tir = false;\\n\\n var t: vec3f;\\n if (dot(i, n) > 0.){\\n \\n if (dot(i, n) > sqrt(mu*mu - 1.)/mu){\\n t = refract(i, -n, mu);\\n } else {\\n t = reflect(i, n);\\n *tir = true;\\n }\\n\\n } else {\\n t = refract(i, n, 1./mu);\\n }\\n\\n return t;\\n}\\n\\nconst k = 1103515245u;\\nfn hash(y: uint3 ) -> vec3f {\\n var x = y;\\n x = (uint3(x.x>>8u, x.y>>8u, x.z>>8u)^x.yzx)*k;\\n x = (uint3(x.x>>8u, x.y>>8u, x.z>>8u)^x.yzx)*k;\\n x = (uint3(x.x>>8u, x.y>>8u, x.z>>8u)^x.yzx)*k;\\n \\n return vec3f(x)*(1.0/float(0xffffffffu));\\n} //hash by IQ https://www.shadertoy.com/view/XlXcW4. Ported to wgsl by ahs3n.\\n\\nfn df(p: vec3f) -> float {\\n let n = i32(custom.geometry*5);\\n switch (n) {\\n\\n case 0: {\\n return length(p-vec3f(0,0,1.))-.6 + sin(p.x*30.)*.01;\\n }\\n\\n case 1: {\\n return max(\\n max(-p.z + .2, abs(p.x)-.9), \\n max(abs(p.y)-.9, p.z-.5 - sin((p.x + p.y)*20.)*.02)\\n );\\n }\\n\\n case 2: {\\n return max(abs(p.x)-.5, max(abs(p.y)-.3, abs(p.z-.7) -.2- p.x*.3));\\n }\\n\\n case 3: {\\n return max(abs(p.x), max(abs(p.y), abs(p.z-.7)))-.3;\\n }\\n\\n case 4: {\\n return max(length(p.xy)-.2, -p.z-.1)// + sin(p.z*50.)*.002\\n ;\\n }\\n\\n default: {\\n return length(p - vec3f(0,0,1)) - .4;\\n }\\n }\\n}\\n\\nfn baseDF(p: vec3f) -> float {\\n return max(p.z, max(abs(p.x), abs(p.y))-1);\\n}\\n\\nfn dfSum(p: vec3f) -> float {\\n return min(df(p), baseDF(p));\\n}\\n\\nfn normal(p: vec3f) -> vec3f {\\n let d = dfSum(p);\\n return normalize(\\n vec3f(\\n d - dfSum(p - vec3f(eps,0,0)),\\n d - dfSum(p - vec3f(0,eps,0)),\\n d - dfSum(p - vec3f(0,0,eps))\\n )\\n );\\n}\\n\\nfn ray(o: vec3f, dir: vec3f, negative: bool) -> vec4f { \\n /*\\n vec4f(\\n depth,\\n hit?1:0,\\n ground?1:0,\\n [empty]\\n )\\n */\\n var p = o; \\n var t = 0.; \\n var d = 0.;\\n var ground = false;\\n var hit = false;\\n\\n for (var i = 0; i < 1<<7; i++){\\n p = o+dir*t;\\n d = min(\\n df(p) * select(1.,-1.,negative),\\n baseDF(p)\\n );\\n\\n t += d;\\n }\\n\\n if (d < eps){\\n hit = true;\\n }\\n\\n if (baseDF(p) < df(p) * select(1.,-1.,negative) && p.z < eps && hit){\\n ground = true;\\n }\\n\\n return vec4f(\\n t,\\n select(0.,1.,hit),\\n select(0.,1.,ground),\\n 0\\n );\\n}\\n\\nfn path(\\n o: vec3f, \\n dir: vec3f, \\n mu: f32,\\n rayDirOut: ptr, \\n rayOutPos: ptr \\n) -> vec4f {\\n\\n var hit = ray(o, dir, false);\\n\\n /*\\n vec4f(\\n depth,\\n hit?1:0,\\n ground?1:0,\\n [empty]\\n )\\n */\\n\\n var p = o + hit.x * dir;\\n var rd = dir;\\n\\n if (hit.y > .5 && hit.z < .5){ // hit and not ground\\n var norm = normal(p);\\n p -= norm*eps*3;\\n var tir: bool;\\n rd = refr(dir, norm, mu, &tir);\\n var tHit: vec4f;\\n var inside = true;\\n\\n for (var i = 0; i < 8; i++){\\n\\n

@compute.toys