{\\n let a = cro(ray.d,B.p0-2*B.p1+B.p2);\\n let b = cro(ray.d,-2*(B.p0-B.p1)); \\n let c = cro(ray.d,B.p0-ray.o);\\n if (abs(a) < 1e-8) {\\n if (abs(b) < 1e-8) { return vec4f(-1); }\\n let u = -c / b;\\n let t = select(-1,dot(bezier_quad(B,u)-ray.o,ray.d),0<=u&&u<=1);\\n return vec4f(t,t,u,u);\\n }\\n\\n let det = b * b - 4. * a * c;\\n if (det < 0.) { return vec4f(-1); }\\n\\n let t = (-b + sqrt(det)*vec2(1,-1))/(2*a);\\n let t0 = select(-1,dot(bezier_quad(B,t.x)-ray.o,ray.d),0<=t.x&&t.x<=1);\\n let t1 = select(-1,dot(bezier_quad(B,t.y)-ray.o,ray.d),0<=t.y&&t.y<=1);\\n\\n return vec4f(t0,t1,t);\\n}\\n\\n#define N_POINTS 5\\n#define N_CURVES 15\\n#define DURATION 1200\\n\\n#define N_CURVES_S 30\\nconst CURVES_S = array(\\nBezierQuad(vec2f(0.5247,0.5630),vec2f(0.3523,0.6821),vec2f(0.2406,0.6594)),\\nBezierQuad(vec2f(0.0447,0.6196),vec2f(-0.0987,0.5904),vec2f(-0.1375,0.4742)),\\nBezierQuad(vec2f(-0.1211,0.0095),vec2f(-0.1090,-0.0235),vec2f(0.0560,-0.0609)),\\nBezierQuad(vec2f(0.1846,-0.0930),vec2f(0.2794,-0.1187),vec2f(0.2986,-0.1474)),\\nBezierQuad(vec2f(0.3065,-0.3420),vec2f(0.2681,-0.3938),vec2f(0.2143,-0.3909)),\\nBezierQuad(vec2f(-0.1622,-0.5389),vec2f(-0.1255,-0.5785),vec2f(-0.0352,-0.6045)),\\nBezierQuad(vec2f(-0.1449,-0.3174),vec2f(-0.1692,-0.3508),vec2f(-0.2116,-0.4091)),\\nBezierQuad(vec2f(0.0434,-0.6272),vec2f(0.1430,-0.6559),vec2f(0.2274,-0.6411)),\\nBezierQuad(vec2f(0.3400,-0.6241),vec2f(0.3831,-0.6181),vec2f(0.4416,-0.5786)),\\nBezierQuad(vec2f(0.5007,-0.5387),vec2f(0.5632,-0.4964),vec2f(0.5824,-0.4415)),\\nBezierQuad(vec2f(0.6192,-0.3442),vec2f(0.6538,-0.2546),vec2f(0.6311,-0.1592)),\\nBezierQuad(vec2f(0.3569,0.1265),vec2f(0.3104,0.1507),vec2f(0.2402,0.1675)),\\nBezierQuad(vec2f(0.1308,0.1936),vec2f(0.0938,0.2024),vec2f(0.1011,0.2763)),\\nBezierQuad(vec2f(0.1125,0.3672),vec2f(0.1162,0.3946),vec2f(0.1827,0.3998)),\\nBezierQuad(vec2f(0.3208,0.4106),vec2f(0.3528,0.4131),vec2f(0.4243,0.3479)),\\nBezierQuad(vec2f(0.5053,0.3915),vec2f(0.5219,0.4186),vec2f(0.5564,0.4750)),\\nBezierQuad(vec2f(-0.3078,0.5687),vec2f(-0.3078,0.5831),vec2f(-0.3078,0.6112)),\\nBezierQuad(vec2f(-0.3436,0.6480),vec2f(-0.3698,0.6480),vec2f(-0.4015,0.6480)),\\nBezierQuad(vec2f(-0.4219,0.4404),vec2f(-0.4695,0.4875),vec2f(-0.5505,0.5676)),\\nBezierQuad(vec2f(-0.6318,0.5945),vec2f(-0.6318,0.5519),vec2f(-0.6318,0.4960)),\\nBezierQuad(vec2f(-0.5800,0.4041),vec2f(-0.4701,0.2955),vec2f(-0.3787,0.2051)),\\nBezierQuad(vec2f(-0.3078,0.1544),vec2f(-0.3078,0.2311),vec2f(-0.3078,0.3060)),\\nBezierQuad(vec2f(-0.3910,0.0228),vec2f(-0.4695,0.1005),vec2f(-0.5723,0.2022)),\\nBezierQuad(vec2f(-0.5074,-0.0546),vec2f(-0.4701,-0.0915),vec2f(-0.3420,-0.2182)),\\nBezierQuad(vec2f(-0.3078,-0.2255),vec2f(-0.3078,-0.1559),vec2f(-0.3078,-0.1124)),\\nBezierQuad(vec2f(-0.3429,-0.4116),vec2f(-0.4695,-0.2865),vec2f(-0.5669,-0.1902)),\\nBezierQuad(vec2f(-0.5114,-0.4377),vec2f(-0.4701,-0.4785),vec2f(-0.3906,-0.5571)),\\nBezierQuad(vec2f(-0.3078,-0.5907),vec2f(-0.3078,-0.5429),vec2f(-0.3078,-0.4839)),\\nBezierQuad(vec2f(-0.6318,-0.5337),vec2f(-0.6318,-0.5813),vec2f(-0.6318,-0.6193)),\\nBezierQuad(vec2f(-0.5524,-0.5960),vec2f(-0.5669,-0.5815),vec2f(-0.5933,-0.5551)),\\n);\\n\\n\\nstruct Curve {\\n pts: array,\\n t: array,\\n u: array,\\n A: array,\\n\\n p: BezierCubic,\\n n: u32,\\n err: f32,\\n}\\n#storage curves array\\n\\nfn get_t() -> f32 {\\n return quintic(saturate(min(f32(time.frame), DURATION)/DURATION));\\n}\\n\\n#workgroup_count init 1 1 1 \\n@compute @workgroup_size(N_CURVES)\\nfn init(@builtin(global_invocation_id) id: vec3u) {\\n if (curves[id.x].err > 1e-1 && curves[id.x].n < 60) { return; }\\n\\n curves[id.x].err = 0.;\\n curves[id.x].n = 0;\\n let seed = time.frame * N_CURVES + id.x;\\n\\n var pts: array;\\n for (var i = 0u; i < N_POINTS; i++) {\\n pts[i] = hash13(seed + i).xy*2-1;\\n }\\n curves[id.x].pts = pts;\\n\\n for (var i = 1;

@compute.toys