CD Spectral Inversion

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.9%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22%2C%20app_title%3D%22CD%20Spectral%20Inversion%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20import%20json%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%0A%20%20%20%20return%20json%2C%20mo%2C%20np%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Circular%20dichroism%3A%20reading%20protein%20structure%20with%20a%20physics-constrained%20network%0A%0A%20%20%20%20Far-UV%20circular%20dichroism%20reads%20out%20a%20protein's%20secondary%20structure%2C%20but%20only%20forwards%3A%20a%0A%20%20%20%20spectrum%20is%20the%20weighted%20sum%20of%20fixed%20basis%20curves%2C%20%24y(%5Clambda)%3Df_H%20B_H%2Bf_S%20B_S%2Bf_C%20B_C%24.%0A%20%20%20%20Going%20backwards%2C%20recovering%20the%20fractions%20from%20a%20noisy%20spectrum%2C%20is%20ill-posed%20because%20the%0A%20%20%20%20basis%20curves%20overlap.%20Here%20a%20small%20MLP%20learns%20the%20inverse%20map%2C%20with%20a%20**softmax%20output**%20that%0A%20%20%20%20forces%20the%20prediction%20onto%20the%20simplex%20(non-negative%2C%20summing%20to%20one)%20and%20a%20**physics%20loss**%0A%20%20%20%20that%20penalises%20predictions%20whose%20reconstructed%20spectrum%20does%20not%20match%20the%20input.%20Watch%20the%0A%20%20%20%20reconstruction%20converge%20onto%20the%20measurement%20as%20the%20network%20trains.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(np)%3A%0A%20%20%20%20WL%20%3D%20np.arange(190.0%2C%20251.0%2C%201.0)%0A%20%20%20%20N_WL%20%3D%20len(WL)%0A%0A%20%20%20%20def%20_basis()%3A%0A%20%20%20%20%20%20%20%20helix%20%3D%20(-30%20*%20np.exp(-0.5%20*%20((WL%20-%20208)%20%2F%205)%20**%202)%20-%2025%20*%20np.exp(-0.5%20*%20((WL%20-%20222)%20%2F%204)%20**%202)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%2050%20*%20np.exp(-0.5%20*%20((WL%20-%20193)%20%2F%208)%20**%202))%0A%20%20%20%20%20%20%20%20sheet%20%3D%20(-15%20*%20np.exp(-0.5%20*%20((WL%20-%20216)%20%2F%207)%20**%202)%20%2B%2020%20*%20np.exp(-0.5%20*%20((WL%20-%20195)%20%2F%206)%20**%202))%0A%20%20%20%20%20%20%20%20coil%20%3D%20(-8%20*%20np.exp(-0.5%20*%20((WL%20-%20200)%20%2F%206)%20**%202)%20%2B%205%20*%20np.exp(-0.5%20*%20((WL%20-%20218)%20%2F%208)%20**%202))%0A%20%20%20%20%20%20%20%20return%20np.column_stack(%5Bhelix%2C%20sheet%2C%20coil%5D)%0A%0A%20%20%20%20B%20%3D%20_basis()%0A%20%20%20%20SCALE%20%3D%2030.0%20%20%20%23%20input%20spectra%20are%20O(30)%3B%20scale%20to%20O(1)%20for%20training%20stability%0A%0A%20%20%20%20def%20spectrum(frac)%3A%0A%20%20%20%20%20%20%20%20return%20B%20%40%20np.asarray(frac%2C%20float)%0A%0A%20%20%20%20return%20B%2C%20N_WL%2C%20SCALE%2C%20spectrum%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(B%2C%20N_WL%2C%20SCALE%2C%20np)%3A%0A%20%20%20%20%23%20%E2%94%80%E2%94%80%20physics-informed%20MLP%3A%2061%20-%3E%20H%20-%3E%20H%20-%3E%203%20logits%20-%3E%20softmax%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%0A%20%20%20%20def%20init_mlp(H%2C%20seed%3D42)%3A%0A%20%20%20%20%20%20%20%20rng%20%3D%20np.random.default_rng(seed)%0A%20%20%20%20%20%20%20%20k1%20%3D%20np.sqrt(2.0%20%2F%20N_WL)%3B%20k2%20%3D%20np.sqrt(2.0%20%2F%20H)%0A%20%20%20%20%20%20%20%20return%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22W1%22%3A%20rng.standard_normal((H%2C%20N_WL))%20*%20k1%2C%20%22b1%22%3A%20np.zeros(H)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22W2%22%3A%20rng.standard_normal((H%2C%20H))%20*%20k2%2C%20%20%20%20%22b2%22%3A%20np.zeros(H)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22W3%22%3A%20rng.standard_normal((3%2C%20H))%20*%20k2%2C%20%20%20%20%22b3%22%3A%20np.zeros(3)%2C%0A%20%20%20%20%20%20%20%20%7D%0A%0A%20%20%20%20def%20_softmax(Z)%3A%0A%20%20%20%20%20%20%20%20e%20%3D%20np.exp(Z%20-%20Z.max(1%2C%20keepdims%3DTrue))%0A%20%20%20%20%20%20%20%20return%20e%20%2F%20e.sum(1%2C%20keepdims%3DTrue)%0A%0A%20%20%20%20def%20forward(p%2C%20X)%3A%0A%20%20%20%20%20%20%20%20z1%20%3D%20X%20%40%20p%5B%22W1%22%5D.T%20%2B%20p%5B%22b1%22%5D%3B%20a1%20%3D%20np.maximum(0%2C%20z1)%0A%20%20%20%20%20%20%20%20z2%20%3D%20a1%20%40%20p%5B%22W2%22%5D.T%20%2B%20p%5B%22b2%22%5D%3B%20a2%20%3D%20np.maximum(0%2C%20z2)%0A%20%20%20%20%20%20%20%20logit%20%3D%20a2%20%40%20p%5B%22W3%22%5D.T%20%2B%20p%5B%22b3%22%5D%0A%20%20%20%20%20%20%20%20return%20_softmax(logit)%2C%20(z1%2C%20a1%2C%20z2%2C%20a2)%0A%0A%20%20%20%20def%20train(p%2C%20X%2C%20Y%2C%20x_query%2C%20y_meas%2C%20epochs%2C%20lr%2C%20lam%2C%20B_arr)%3A%0A%20%20%20%20%20%20%20%20loss_hist%2C%20pred_hist%2C%20recon_hist%20%3D%20%5B%5D%2C%20%5B%5D%2C%20%5B%5D%0A%20%20%20%20%20%20%20%20_scale%20%3D%20SCALE%0A%20%20%20%20%20%20%20%20_var%20%3D%20max(float(np.var(y_meas))%2C%201e-9)%0A%20%20%20%20%20%20%20%20for%20_%20in%20range(epochs)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pred%2C%20(z1%2C%20a1%2C%20z2%2C%20a2)%20%3D%20forward(p%2C%20X)%0A%20%20%20%20%20%20%20%20%20%20%20%20_Bn%20%3D%20X.shape%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20MSE%20gradient%20through%20softmax%20(analytic%2C%20vectorised)%0A%20%20%20%20%20%20%20%20%20%20%20%20d%20%3D%20(2.0%20%2F%203.0)%20*%20(pred%20-%20Y)%0A%20%20%20%20%20%20%20%20%20%20%20%20dlogit%20%3D%20pred%20*%20(d%20-%20(pred%20*%20d).sum(1%2C%20keepdims%3DTrue))%0A%20%20%20%20%20%20%20%20%20%20%20%20dW3%20%3D%20dlogit.T%20%40%20a2%20%2F%20_Bn%3B%20db3%20%3D%20dlogit.mean(0)%0A%20%20%20%20%20%20%20%20%20%20%20%20da2%20%3D%20(dlogit%20%40%20p%5B%22W3%22%5D)%20*%20(z2%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20dW2%20%3D%20da2.T%20%40%20a1%20%2F%20_Bn%3B%20db2%20%3D%20da2.mean(0)%0A%20%20%20%20%20%20%20%20%20%20%20%20da1%20%3D%20(da2%20%40%20p%5B%22W2%22%5D)%20*%20(z1%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20dW1%20%3D%20da1.T%20%40%20X%20%2F%20_Bn%3B%20db1%20%3D%20da1.mean(0)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_k%2C%20_g%20in%20((%22W1%22%2C%20dW1)%2C%20(%22b1%22%2C%20db1)%2C%20(%22W2%22%2C%20dW2)%2C%20(%22b2%22%2C%20db2)%2C%20(%22W3%22%2C%20dW3)%2C%20(%22b3%22%2C%20db3))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20p%5B_k%5D%20-%3D%20lr%20*%20_g%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20track%3A%20supervised%20MSE%20%2B%20physics%20consistency%20on%20the%20query%20spectrum%0A%20%20%20%20%20%20%20%20%20%20%20%20_q%2C%20_%20%3D%20forward(p%2C%20x_query%5BNone%2C%20%3A%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_qc%20%3D%20_q%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_recon%20%3D%20B_arr%20%40%20_qc%0A%20%20%20%20%20%20%20%20%20%20%20%20_mse%20%3D%20float(np.mean((pred%20-%20Y)%20**%202))%0A%20%20%20%20%20%20%20%20%20%20%20%20_phys%20%3D%20float(np.mean((_recon%20-%20y_meas)%20**%202)%20%2F%20_var)%0A%20%20%20%20%20%20%20%20%20%20%20%20loss_hist.append(_mse%20%2B%20lam%20*%20_phys)%0A%20%20%20%20%20%20%20%20%20%20%20%20pred_hist.append(_qc.tolist())%0A%20%20%20%20%20%20%20%20%20%20%20%20recon_hist.append(np.round(_recon%2C%203).tolist())%0A%20%20%20%20%20%20%20%20return%20loss_hist%2C%20pred_hist%2C%20recon_hist%0A%0A%20%20%20%20def%20gen_dataset(n%2C%20sigma%2C%20seed)%3A%0A%20%20%20%20%20%20%20%20rng%20%3D%20np.random.default_rng(seed)%0A%20%20%20%20%20%20%20%20u%20%3D%20np.sort(rng.uniform(0%2C%201%2C%20(n%2C%202))%2C%20axis%3D1)%0A%20%20%20%20%20%20%20%20comps%20%3D%20np.column_stack(%5Bu%5B%3A%2C%200%5D%2C%20u%5B%3A%2C%201%5D%20-%20u%5B%3A%2C%200%5D%2C%201%20-%20u%5B%3A%2C%201%5D%5D)%0A%20%20%20%20%20%20%20%20specs%20%3D%20comps%20%40%20B.T%20%2B%20sigma%20*%20rng.standard_normal((n%2C%20N_WL))%0A%20%20%20%20%20%20%20%20return%20specs%20%2F%20SCALE%2C%20comps%0A%0A%20%20%20%20return%20gen_dataset%2C%20init_mlp%2C%20train%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20fh%20%3D%20mo.ui.slider(0.0%2C%201.0%2C%20value%3D0.5%2C%20step%3D0.05%2C%20label%3D%22true%20helix%22%2C%20show_value%3DTrue)%0A%20%20%20%20fs%20%3D%20mo.ui.slider(0.0%2C%201.0%2C%20value%3D0.2%2C%20step%3D0.05%2C%20label%3D%22true%20sheet%22%2C%20show_value%3DTrue)%0A%20%20%20%20sig%20%3D%20mo.ui.slider(0.0%2C%203.0%2C%20value%3D0.8%2C%20step%3D0.1%2C%20label%3D%22measurement%20noise%20%CF%83%22%2C%20show_value%3DTrue)%0A%20%20%20%20lam%20%3D%20mo.ui.slider(0.0%2C%202.0%2C%20value%3D0.5%2C%20step%3D0.1%2C%20label%3D%22physics%20weight%20%CE%BB%22%2C%20show_value%3DTrue)%0A%20%20%20%20ep%20%3D%20mo.ui.slider(20%2C%2080%2C%20value%3D45%2C%20step%3D5%2C%20label%3D%22training%20epochs%22%2C%20show_value%3DTrue)%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20**Set%20a%20true%20structure%2C%20measure%20it%20noisily%2C%20and%20let%20the%20network%20learn%20the%20inverse%20map.**%0A%20%20%20%20The%20physics%20weight%20%CE%BB%20controls%20how%20much%20the%20loss%20rewards%20a%20self-consistent%20reconstruction.%0A%20%20%20%20%7Bmo.hstack(%5Bfh%2C%20fs%2C%20sig%2C%20lam%2C%20ep%5D%2C%20justify%3D%22start%22%2C%20gap%3D2)%7D%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%20ep%2C%20fh%2C%20fs%2C%20lam%2C%20sig%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(B%2C%20ep%2C%20fh%2C%20fs%2C%20gen_dataset%2C%20init_mlp%2C%20lam%2C%20np%2C%20sig%2C%20spectrum%2C%20train)%3A%0A%20%20%20%20_h%20%3D%20fh.value%3B%20_s%20%3D%20min(fs.value%2C%201%20-%20_h)%3B%20_c%20%3D%20max(0.0%2C%201%20-%20_h%20-%20_s)%0A%20%20%20%20true_frac%20%3D%20np.array(%5B_h%2C%20_s%2C%20_c%5D)%0A%20%20%20%20_rng%20%3D%20np.random.default_rng(0)%0A%20%20%20%20y_meas%20%3D%20spectrum(true_frac)%20%2B%20sig.value%20*%20_rng.standard_normal(len(B))%0A%20%20%20%20_x_query%20%3D%20y_meas%20%2F%2030.0%0A%0A%20%20%20%20_X%2C%20_Y%20%3D%20gen_dataset(450%2C%20sig.value%2C%201)%0A%20%20%20%20_p%20%3D%20init_mlp(32%2C%20seed%3D42)%0A%20%20%20%20loss_hist%2C%20pred_hist%2C%20recon_hist%20%3D%20train(%0A%20%20%20%20%20%20%20%20_p%2C%20_X%2C%20_Y%2C%20_x_query%2C%20y_meas%2C%20int(ep.value)%2C%20lr%3D0.15%2C%20lam%3Dlam.value%2C%20B_arr%3DB)%0A%0A%20%20%20%20rec_frac%20%3D%20np.array(pred_hist%5B-1%5D)%0A%20%20%20%20%23%20spectral%20R%5E2%20of%20the%20final%20reconstruction%0A%20%20%20%20_recon%20%3D%20np.array(recon_hist%5B-1%5D)%0A%20%20%20%20_ss_res%20%3D%20float(np.sum((_recon%20-%20y_meas)%20**%202))%0A%20%20%20%20_ss_tot%20%3D%20float(np.sum((y_meas%20-%20y_meas.mean())%20**%202))%20%2B%201e-12%0A%20%20%20%20cd_stats%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22true%22%3A%20true_frac.tolist()%2C%20%22rec%22%3A%20rec_frac.tolist()%2C%0A%20%20%20%20%20%20%20%20%22r2%22%3A%201%20-%20_ss_res%20%2F%20_ss_tot%2C%20%22final_loss%22%3A%20float(loss_hist%5B-1%5D)%2C%0A%20%20%20%20%7D%0A%20%20%20%20cd_anim%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22frames%22%3A%20len(loss_hist)%2C%0A%20%20%20%20%20%20%20%20%22loss%22%3A%20%5Bfloat(v)%20for%20v%20in%20loss_hist%5D%2C%0A%20%20%20%20%20%20%20%20%22recon%22%3A%20recon_hist%2C%0A%20%20%20%20%20%20%20%20%22pred%22%3A%20pred_hist%2C%0A%20%20%20%20%20%20%20%20%22meas%22%3A%20np.round(y_meas%2C%203).tolist()%2C%0A%20%20%20%20%20%20%20%20%22helix_basis%22%3A%20np.round(B%5B%3A%2C%200%5D%2C%203).tolist()%2C%0A%20%20%20%20%20%20%20%20%22sheet_basis%22%3A%20np.round(B%5B%3A%2C%201%5D%2C%203).tolist()%2C%0A%20%20%20%20%20%20%20%20%22true%22%3A%20true_frac.tolist()%2C%0A%20%20%20%20%20%20%20%20%22wl0%22%3A%20float(190.0)%2C%20%22wl1%22%3A%20float(250.0)%2C%0A%20%20%20%20%7D%0A%20%20%20%20return%20cd_anim%2C%20cd_stats%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(cd_anim%2C%20cd_stats%2C%20json%2C%20mo)%3A%0A%20%20%20%20_d%20%3D%20cd_anim%0A%20%20%20%20_JS%20%3D%20r%22%22%22%3C!doctype%20html%3E%3Chtml%3E%3Chead%3E%3Cmeta%20charset%3D%22utf-8%22%3E%3Cstyle%3E%0A%20%20%20%20html%2Cbody%7Bmargin%3A0%3Bbackground%3A%230f1117%3Boverflow%3Ahidden%7Dcanvas%7Bdisplay%3Ablock%3Bwidth%3A100vw%3Bheight%3A100vh%7D%0A%20%20%20%20%3C%2Fstyle%3E%3C%2Fhead%3E%3Cbody%3E%3Ccanvas%20id%3D%22c%22%3E%3C%2Fcanvas%3E%3Cscript%3E%0A%20%20%20%20%2F*%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20CANVAS%20ENGINE%20v3%20%E2%80%94%20shared%20template%20for%20all%20marimo%20notebooks%0A%20%20%20%20%20%20%20Injected%20by%20replacing%20__SPEC__%20with%20json.dumps(spec).%0A%0A%20%20%20%20%20%20%20Spec%20root%20fields%3A%0A%20%20%20%20%20frames%2C%20dt%2C%20panels%5B%5D%2C%20layout%5B%5D%2C%20bg%20(optional%20%22%23rrggbb%22)%0A%0A%20%20%20%20%20%20%20Per-panel%20fields%20(all%20optional%20where%20sensible)%3A%0A%20%20%20%20%20kind%20%20%20%20%20%20%20%E2%80%93%20%22timeseries%22%7C%22scatter%22%7C%22heatmap%22%7C%22particles%22%7C...%0A%20%20%20%20%20title%20%20%20%20%20%20%E2%80%93%20panel%20title%20shown%20at%20top%0A%20%20%20%20%20xlabel%20%20%20%20%20%E2%80%93%20x-axis%20label%20(may%20include%20units%20e.g.%20%22position%20x%20(nm)%22)%0A%20%20%20%20%20ylabel%20%20%20%20%20%E2%80%93%20y-axis%20label%0A%20%20%20%20%20xmin%2Cxmax%20%20%E2%80%93%20explicit%20data%20range%20(auto-computed%20if%20absent)%0A%20%20%20%20%20ymin%2Cymax%20%20%E2%80%93%20explicit%20data%20range%20(auto-computed%20if%20absent)%0A%20%20%20%20%20logx%2Clogy%20%20%E2%80%93%20true%20for%20log-scale%20axis%0A%20%20%20%20%20%20%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%20*%2F%0A%0A%20%20%20%20const%20S%3D__SPEC__%2Ccv%3Ddocument.getElementById(%22c%22)%2Cctx%3Dcv.getContext(%222d%22)%3B%0A%20%20%20%20const%20DPR%3DMath.min(window.devicePixelRatio%7C%7C1%2C2)%3B%0A%20%20%20%20const%20BG%3DS.bg%7C%7C%22%230f1117%22%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Canvas%20sizing%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20fit()%7Bconst%20w%3Dcv.offsetWidth%7C%7Cwindow.innerWidth%7C%7C900%2Ch%3Dcv.offsetHeight%7C%7Cwindow.innerHeight%7C%7C500%3Bif(w%3E0%26%26h%3E0)%7Bcv.width%3Dw*DPR%3Bcv.height%3Dh*DPR%3B%7D%7D%0A%20%20%20%20addEventListener(%22resize%22%2Cfit)%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Animation%20state%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20let%20frame%3D0%2Clast%3D0%2Cplaying%3Dtrue%2Cspeed%3D1.0%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Per-panel%20zoom%2Fpan%20state%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20const%20VP%3DS.panels.map(()%3D%3E(%7Bzoom%3A1%2Cpx%3A0%2Cpy%3A0%7D))%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Mouse%2Ftouch%20state%20for%20zoom%2Bpan%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20let%20drag%3Dfalse%2CdragPi%3D-1%2CdragX%3D0%2CdragY%3D0%3B%0A%20%20%20%20function%20panelAt(cx%2Ccy)%7B%0A%20%20%20%20%20%20const%20W%3Dcv.width%2CH%3Dcv.height%3B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3CS.layout.length%3Bi%2B%2B)%7B%0A%20%20%20%20const%5Blx%2Cly%2Clw%2Clh%5D%3DS.layout%5Bi%5D%3B%0A%20%20%20%20if(cx%3E%3Dlx*W%26%26cx%3C%3D(lx%2Blw)*W%26%26cy%3E%3Dly*H%26%26cy%3C%3D(ly%2Blh)*H)return%20i%3B%0A%20%20%20%20%20%20%7Dreturn%20-1%3B%0A%20%20%20%20%7D%0A%20%20%20%20cv.addEventListener(%22mousedown%22%2Ce%3D%3E%7Bdrag%3Dtrue%3Bconst%20r%3Dcv.getBoundingClientRect()%3BdragX%3De.clientX-r.left%3BdragY%3De.clientY-r.top%3BdragPi%3DpanelAt(dragX*DPR%2CdragY*DPR)%3B%7D)%3B%0A%20%20%20%20cv.addEventListener(%22mousemove%22%2Ce%3D%3E%7B%0A%20%20%20%20%20%20if(!drag%7C%7CdragPi%3C0)return%3B%0A%20%20%20%20%20%20const%20r%3Dcv.getBoundingClientRect()%2Cnx%3D(e.clientX-r.left)%2Cny%3D(e.clientY-r.top)%3B%0A%20%20%20%20%20%20const%20dx%3D(nx-dragX)%2Fcv.clientWidth%2Cdy%3D(ny-dragY)%2Fcv.clientHeight%3B%0A%20%20%20%20%20%20VP%5BdragPi%5D.px%2B%3Ddx%3BVP%5BdragPi%5D.py%2B%3Ddy%3BdragX%3Dnx%3BdragY%3Dny%3B%0A%20%20%20%20%7D)%3B%0A%20%20%20%20cv.addEventListener(%22mouseup%22%2C()%3D%3Edrag%3Dfalse)%3B%0A%20%20%20%20cv.addEventListener(%22mouseleave%22%2C()%3D%3Edrag%3Dfalse)%3B%0A%20%20%20%20cv.addEventListener(%22wheel%22%2Ce%3D%3E%7B%0A%20%20%20%20%20%20const%20r%3Dcv.getBoundingClientRect()%2Ccx%3D(e.clientX-r.left)*DPR%2Ccy%3D(e.clientY-r.top)*DPR%3B%0A%20%20%20%20%20%20const%20pi%3DpanelAt(cx%2Ccy)%3Bif(pi%3C0)return%3B%0A%20%20%20%20%20%20const%20factor%3DMath.exp(-e.deltaY*0.001)%3B%0A%20%20%20%20%20%20VP%5Bpi%5D.zoom%3DMath.max(0.15%2CMath.min(20%2CVP%5Bpi%5D.zoom*factor))%3B%0A%20%20%20%20%20%20e.preventDefault()%3B%0A%20%20%20%20%7D%2C%7Bpassive%3Afalse%7D)%3B%0A%20%20%20%20cv.addEventListener(%22dblclick%22%2Ce%3D%3E%7B%0A%20%20%20%20%20%20const%20r%3Dcv.getBoundingClientRect()%2Ccx%3D(e.clientX-r.left)*DPR%2Ccy%3D(e.clientY-r.top)*DPR%3B%0A%20%20%20%20%20%20const%20pi%3DpanelAt(cx%2Ccy)%3Bif(pi%3E%3D0)%7BVP%5Bpi%5D.zoom%3D1%3BVP%5Bpi%5D.px%3D0%3BVP%5Bpi%5D.py%3D0%3B%7D%0A%20%20%20%20%7D)%3B%0A%20%20%20%20cv.addEventListener(%22touchstart%22%2Ce%3D%3E%7Bconst%20t%3De.touches%5B0%5D%3Bconst%20r%3Dcv.getBoundingClientRect()%3BdragX%3Dt.clientX-r.left%3BdragY%3Dt.clientY-r.top%3BdragPi%3DpanelAt(dragX*DPR%2CdragY*DPR)%3Bdrag%3Dtrue%3Be.preventDefault()%3B%7D%2C%7Bpassive%3Afalse%7D)%3B%0A%20%20%20%20cv.addEventListener(%22touchmove%22%2Ce%3D%3E%7Bif(!drag%7C%7CdragPi%3C0)return%3Bconst%20t%3De.touches%5B0%5D%3Bconst%20r%3Dcv.getBoundingClientRect()%3Bconst%20nx%3Dt.clientX-r.left%2Cny%3Dt.clientY-r.top%3BVP%5BdragPi%5D.px%2B%3D(nx-dragX)%2Fcv.clientWidth%3BVP%5BdragPi%5D.py%2B%3D(ny-dragY)%2Fcv.clientHeight%3BdragX%3Dnx%3BdragY%3Dny%3Be.preventDefault()%3B%7D%2C%7Bpassive%3Afalse%7D)%3B%0A%20%20%20%20cv.addEventListener(%22touchend%22%2C()%3D%3Edrag%3Dfalse)%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Tick%20generation%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20niceNum(x%2Cround)%7Bif(x%3D%3D%3D0)return%200%3Bconst%20s%3Dx%3C0%3F-1%3A1%3Bx%3DMath.abs(x)%3Bconst%20e%3DMath.floor(Math.log10(x))%2Cf%3Dx%2FMath.pow(10%2Ce)%3Blet%20nf%3Bif(round)nf%3Df%3C1.5%3F1%3Af%3C3%3F2%3Af%3C7%3F5%3A10%3Belse%20nf%3Df%3C%3D1%3F1%3Af%3C%3D2%3F2%3Af%3C%3D5%3F5%3A10%3Breturn%20s*nf*Math.pow(10%2Ce)%3B%7D%0A%20%20%20%20function%20genTicks(vmin%2Cvmax%2Cn%3D5)%7B%0A%20%20%20%20%20%20if(vmax%3C%3Dvmin)return%5Bvmin%2Cvmax%5D%3B%0A%20%20%20%20%20%20const%20d%3DniceNum((vmax-vmin)%2F(n-1)%2Ctrue)%3B%0A%20%20%20%20%20%20const%20lo%3DMath.ceil(vmin%2Fd)*d%2Chi%3DMath.floor(vmax%2Fd)*d%3B%0A%20%20%20%20%20%20const%20t%3D%5B%5D%3Bfor(let%20v%3Dlo%3Bv%3C%3Dhi%2Bd*0.5%3Bv%2B%3Dd)t.push(parseFloat(v.toPrecision(10)))%3B%0A%20%20%20%20%20%20return%20t.length%3Ft%3A%5Bvmin%2Cvmax%5D%3B%0A%20%20%20%20%7D%0A%20%20%20%20function%20logTicks(vmin%2Cvmax)%7B%0A%20%20%20%20%20%20const%20lo%3DMath.floor(Math.log10(Math.max(vmin%2C1e-30)))%2Chi%3DMath.ceil(Math.log10(Math.max(vmax%2C1e-30)))%3B%0A%20%20%20%20%20%20const%20t%3D%5B%5D%3Bfor(let%20e%3Dlo%3Be%3C%3Dhi%3Be%2B%2B)%7Bt.push(Math.pow(10%2Ce))%3Bif(hi-lo%3C3)%7Bt.push(2*Math.pow(10%2Ce))%3Bt.push(5*Math.pow(10%2Ce))%3B%7D%7D%0A%20%20%20%20%20%20return%20t.filter(v%3D%3Ev%3E%3Dvmin*0.99%26%26v%3C%3Dvmax*1.01)%3B%0A%20%20%20%20%7D%0A%20%20%20%20function%20fmtN(v)%7B%0A%20%20%20%20%20%20if(v%3D%3D%3D0)return%220%22%3Bconst%20a%3DMath.abs(v)%3B%0A%20%20%20%20%20%20if(a%3E%3D1e4%7C%7Ca%3C1e-3)return%20v.toExponential(1).replace(%22e%2B%22%2C%22e%22).replace(%22e0%22%2C%22e%22).replace(%22e-0%22%2C%22e-%22)%3B%0A%20%20%20%20%20%20if(a%3E%3D100)return%20v.toFixed(0)%3Bif(a%3E%3D10)return%20v.toFixed(1)%3Bif(a%3E%3D1)return%20v.toFixed(2)%3B%0A%20%20%20%20%20%20return%20v.toFixed(3).replace(%2F0%2B%24%2F%2C%22%22).replace(%2F%5C.%24%2F%2C%22%22)%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Core%202D%20axis%20renderer%20with%20grid%20%2B%20ticks%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20axes2D(p%2Clx%2Cly%2Clw%2Clh%2Cxd%2Cyd%2Cvp)%7B%0A%20%20%20%20%20%20const%20W%3Dcv.width%2CH%3Dcv.height%3B%0A%20%20%20%20%20%20const%20px%3DMath.round(lx*W)%2Cpy%3DMath.round(ly*H)%2Cpw%3DMath.round(lw*W)%2Cph%3DMath.round(lh*H)%3B%0A%20%20%20%20%20%20const%20hasTitle%3Dp.title%26%26p.title.length%3E0%3B%0A%20%20%20%20%20%20const%20ML%3D46*DPR%2CMR%3D8*DPR%2CMT%3D(hasTitle%3F22%3A8)*DPR%2CMB%3D28*DPR%3B%0A%20%20%20%20%20%20const%20mx%3Dpx%2BML%2Cmy%3Dpy%2BMT%2Cmw%3Dpw-ML-MR%2Cmh%3Dph-MT-MB%3B%0A%0A%20%20%20%20%20%20%2F*%20data%E2%86%92screen%20transforms%20with%20zoom%2Bpan%20*%2F%0A%20%20%20%20%20%20const%20z%3Dvp.zoom%3B%0A%20%20%20%20%20%20const%20toX%3Dv%3D%3Emx%2Bmw*((v-xd.min)%2F(xd.max-xd.min)*z%2Bvp.px)%3B%0A%20%20%20%20%20%20const%20toY%3Dv%3D%3Emy%2Bmh*(1-(v-yd.min)%2F(yd.max-yd.min)*z%2Bvp.py)%3B%0A%20%20%20%20%20%20%2F*%20screen%E2%86%92data%20*%2F%0A%20%20%20%20%20%20const%20fromX%3Dsx%3D%3E(sx-mx)%2Fmw%2Fz*(xd.max-xd.min)%2Bxd.min-vp.px*(xd.max-xd.min)%2Fz%3B%0A%0A%20%20%20%20%20%20%2F*%20clip%20to%20plot%20area%20*%2F%0A%20%20%20%20%20%20ctx.save()%3Bctx.beginPath()%3Bctx.rect(mx%2Cmy%2Cmw%2Cmh)%3Bctx.clip()%3B%0A%0A%20%20%20%20%20%20%2F*%20grid%20*%2F%0A%20%20%20%20%20%20const%20xticks%3Dxd.log%3FlogTicks(xd.min%2Cxd.max)%3AgenTicks(xd.min%2Cxd.max%2C6)%3B%0A%20%20%20%20%20%20const%20yticks%3Dyd.log%3FlogTicks(yd.min%2Cyd.max)%3AgenTicks(yd.min%2Cyd.max%2C5)%3B%0A%20%20%20%20%20%20ctx.strokeStyle%3D%22rgba(255%2C255%2C255%2C0.055)%22%3Bctx.lineWidth%3D1%3B%0A%20%20%20%20%20%20for(const%20v%20of%20xticks)%7Bconst%20x%3Dxd.log%3FtoX(Math.log10(v))%3AtoX(v)%3Bif(x%3Cmx%7C%7Cx%3Emx%2Bmw)continue%3Bctx.beginPath()%3Bctx.moveTo(x%2Cmy)%3Bctx.lineTo(x%2Cmy%2Bmh)%3Bctx.stroke()%3B%7D%0A%20%20%20%20%20%20for(const%20v%20of%20yticks)%7Bconst%20y%3Dyd.log%3FtoY(Math.log10(v))%3AtoY(v)%3Bif(y%3Cmy%7C%7Cy%3Emy%2Bmh)continue%3Bctx.beginPath()%3Bctx.moveTo(mx%2Cy)%3Bctx.lineTo(mx%2Bmw%2Cy)%3Bctx.stroke()%3B%7D%0A%20%20%20%20%20%20ctx.restore()%3B%0A%0A%20%20%20%20%20%20%2F*%20axis%20border%20*%2F%0A%20%20%20%20%20%20ctx.strokeStyle%3D%22%23243040%22%3Bctx.lineWidth%3D1.2%3B%0A%20%20%20%20%20%20ctx.strokeRect(mx%2Cmy%2Cmw%2Cmh)%3B%0A%0A%20%20%20%20%20%20%2F*%20tick%20marks%20%2B%20labels%20*%2F%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.font%3D(9.5*DPR)%2B%22px%20monospace%22%3B%0A%20%20%20%20%20%20ctx.textAlign%3D%22center%22%3B%0A%20%20%20%20%20%20for(const%20v%20of%20xticks)%7B%0A%20%20%20%20const%20x%3Dxd.log%3FtoX(Math.log10(v))%3AtoX(v)%3Bif(x%3Cmx-1%7C%7Cx%3Emx%2Bmw%2B1)continue%3B%0A%20%20%20%20ctx.fillStyle%3D%22%231e2535%22%3Bctx.fillRect(x-0.5%2Cmy%2Bmh%2C1%2C4*DPR)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.fillText(xd.log%3F(%2210%5E%22%2BMath.round(Math.log10(v)))%3AfmtN(v)%2Cx%2Cpy%2Bph-2*DPR)%3B%0A%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20ctx.textAlign%3D%22right%22%3B%0A%20%20%20%20%20%20for(const%20v%20of%20yticks)%7B%0A%20%20%20%20const%20y%3Dyd.log%3FtoY(Math.log10(v))%3AtoY(v)%3Bif(y%3Cmy-1%7C%7Cy%3Emy%2Bmh%2B1)continue%3B%0A%20%20%20%20ctx.fillStyle%3D%22%231e2535%22%3Bctx.fillRect(mx-4*DPR%2Cy-0.5%2C4*DPR%2C1)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.fillText(yd.log%3F(%2210%5E%22%2BMath.round(Math.log10(v)))%3AfmtN(v)%2Cmx-5*DPR%2Cy%2B3.5*DPR)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20axis%20labels%20%2B%20title%20*%2F%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%239aa3b0%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3B%0A%20%20%20%20%20%20if(p.xlabel)%7Bctx.textAlign%3D%22center%22%3Bctx.fillText(p.xlabel%2Cmx%2Bmw%2F2%2Cpy%2Bph)%3B%7D%0A%20%20%20%20%20%20if(p.ylabel)%7Bctx.save()%3Bctx.translate(px%2B11*DPR%2Cmy%2Bmh%2F2)%3Bctx.rotate(-Math.PI%2F2)%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(p.ylabel%2C0%2C0)%3Bctx.restore()%3B%7D%0A%20%20%20%20%20%20if(p.title)%7Bctx.fillStyle%3D%22%23b3bac8%22%3Bctx.font%3D(10.5*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(p.title%2Cmx%2Bmw%2F2%2Cpy%2B13*DPR)%3B%7D%0A%0A%20%20%20%20%20%20return%7Bmx%2Cmy%2Cmw%2Cmh%2CtoX%2CtoY%2CfromX%7D%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Legend%20helper%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20legend(px%2Cpy%2Cpw%2Cph%2Centries)%7B%0A%20%20%20%20%20%20const%20lx%3Dpx%2Bpw-8*DPR%2Cly%3Dpy%2B14*DPR%3B%0A%20%20%20%20%20%20ctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22right%22%3B%0A%20%20%20%20%20%20entries.forEach((%5Blabel%2Ccol%5D%2Ci)%3D%3E%7B%0A%20%20%20%20const%20y%3Dly%2Bi*13*DPR%3B%0A%20%20%20%20ctx.fillStyle%3Dcol%3Bctx.fillRect(lx-28*DPR%2Cy-6*DPR%2C22*DPR%2C2.5*DPR)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.fillText(label%2Clx%2Cy)%3B%0A%20%20%20%20%20%20%7D)%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20HSV%20helper%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20hsv(h%2Cs%2Cv)%7Blet%20r%2Cg%2Cb%2Ci%3DMath.floor(h*6)%2Cf%3Dh*6-i%2Cp%3Dv*(1-s)%2Cq%3Dv*(1-f*s)%2Ct%3Dv*(1-(1-f)*s)%3Bswitch(i%256)%7Bcase%200%3Ar%3Dv%2Cg%3Dt%2Cb%3Dp%3Bbreak%3Bcase%201%3Ar%3Dq%2Cg%3Dv%2Cb%3Dp%3Bbreak%3Bcase%202%3Ar%3Dp%2Cg%3Dv%2Cb%3Dt%3Bbreak%3Bcase%203%3Ar%3Dp%2Cg%3Dq%2Cb%3Dv%3Bbreak%3Bcase%204%3Ar%3Dt%2Cg%3Dp%2Cb%3Dv%3Bbreak%3Bcase%205%3Ar%3Dv%2Cg%3Dp%2Cb%3Dq%3Bbreak%3B%7Dreturn%60rgb(%24%7BMath.round(r*255)%7D%2C%24%7BMath.round(g*255)%7D%2C%24%7BMath.round(b*255)%7D)%60%3B%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%203D%20orbit%20camera%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20let%20az3%3D0.52%2Cel3%3D0.38%2Czoom3%3D1.0%2Cdrag3%3Dfalse%2Clx3%3D0%2Cly3%3D0%3B%0A%20%20%20%20function%20P3(wx%2Cwy%2Cwz%2Ccxp%2Ccyp%2Csc)%7Bconst%20ca%3DMath.cos(az3)%2Csa%3DMath.sin(az3)%3Bconst%20x1%3Dwx*ca%2Bwy*sa%2Cy1%3D-wx*sa%2Bwy*ca%2Cz1%3Dwz%3Bconst%20ce%3DMath.cos(el3)%2Cse%3DMath.sin(el3)%3Breturn%5Bcxp%2Bx1*sc*zoom3%2Ccyp-(-y1*se%2Bz1*ce)*sc*zoom3%5D%3B%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Panel%20rendering%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20drawPanel(p%2Cpi%2Clx%2Cly%2Clw%2Clh)%7B%0A%20%20%20%20%20%20const%20W%3Dcv.width%2CH%3Dcv.height%3B%0A%20%20%20%20%20%20const%20px%3DMath.round(lx*W)%2Cpy%3DMath.round(ly*H)%2Cpw%3DMath.round(lw*W)%2Cph%3DMath.round(lh*H)%3B%0A%20%20%20%20%20%20ctx.save()%3Bctx.beginPath()%3Bctx.rect(px%2Cpy%2Cpw%2Cph)%3Bctx.clip()%3B%0A%20%20%20%20%20%20const%20f%3DMath.min(frame%2CS.frames-1)%3B%0A%20%20%20%20%20%20const%20vp%3DVP%5Bpi%5D%3B%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20PARTICLES%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20if(p.kind%3D%3D%3D%22particles%22)%7B%0A%20%20%20%20const%20xs%3Dp.x%5Bf%5D%2Cys%3Dp.y%5Bf%5D%2Cths%3Dp.theta%3Fp.theta%5Bf%5D%3Anull%2CL%3Dp.L%3B%0A%20%20%20%20const%20sc%3DMath.min(pw%2Cph)*0.9%2FL%2Cox%3Dpx%2B(pw-L*sc)%2F2%2Coy%3Dpy%2B(ph-L*sc)%2F2%2Cr%3DMath.max(1.5%2C(p.r%7C%7C0.3)*sc)%3B%0A%20%20%20%20ctx.strokeStyle%3D%22%23243040%22%3Bctx.lineWidth%3D1%3Bctx.strokeRect(ox%2Coy%2CL*sc%2CL*sc)%3B%0A%20%20%20%20for(let%20i%3D0%3Bi%3Cxs.length%3Bi%2B%2B)%7B%0A%20%20%20%20%20%20const%20clr%3Dths%3Fhsv(((ths%5Bi%5D%25(2*Math.PI))%2B2*Math.PI)%25(2*Math.PI)%2F(2*Math.PI)%2C0.85%2C0.95)%3A%22%237dd3fc%22%3B%0A%20%20%20%20%20%20ctx.fillStyle%3Dclr%3Bctx.beginPath()%3Bctx.arc(ox%2Bxs%5Bi%5D*sc%2Coy%2Bys%5Bi%5D*sc%2Cr%2C0%2C7)%3Bctx.fill()%3B%0A%20%20%20%20%7D%0A%20%20%20%20%2F*%20scale%20bar%20*%2F%0A%20%20%20%20const%20bar%3DMath.round(L%2F4)*1%2CbarPx%3Dbar*sc%3B%0A%20%20%20%20ctx.strokeStyle%3D%22%235a6880%22%3Bctx.lineWidth%3D1.5*DPR%3B%0A%20%20%20%20ctx.beginPath()%3Bctx.moveTo(ox%2B4*DPR%2Coy%2BL*sc-8*DPR)%3Bctx.lineTo(ox%2B4*DPR%2BbarPx%2Coy%2BL*sc-8*DPR)%3Bctx.stroke()%3B%0A%20%20%20%20ctx.fillStyle%3D%22%235a6880%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(fmtN(bar)%2Cox%2B4*DPR%2Coy%2BL*sc-12*DPR)%3B%0A%20%20%20%20if(p.title)%7Bctx.fillStyle%3D%22%239aa3b0%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(p.title%2Cpx%2B4*DPR%2Cpy%2B12*DPR)%3B%7D%0A%20%20%20%20if(p.eta_val)%7Bctx.fillStyle%3D%22%23626b7c%22%3Bctx.fillText(%22%CE%B7%3D%22%2Bp.eta_val%5Bf%5D.toFixed(2)%2B%22%20%20N%3D%22%2Bxs.length%2Cpx%2B4*DPR%2Cpy%2B22*DPR)%3B%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20TIMESERIES%20(single%20or%20multi-series%20growing%20line)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22timeseries%22%7C%7Cp.kind%3D%3D%3D%22phi_t%22%7C%7Cp.kind%3D%3D%3D%22r_t%22%7C%7Cp.kind%3D%3D%3D%22grow_curve%22)%7B%0A%20%20%20%20const%20ys_arr%3Dp.ys%7C%7Cp.ys_arr%7C%7C(p.phi%3F%5Bp.phi%5D%3A(p.r%3F%5Bp.r%5D%3A(p.y%3F%5Bp.y%5D%3A%5B%5D)))%3B%0A%20%20%20%20const%20n%3Dys_arr%5B0%5D%3Fys_arr%5B0%5D.length%3A0%3B%0A%20%20%20%20const%20ts%3Dp.t%7C%7C(Array.from(%7Blength%3An%7D%2C(_%2Ci)%3D%3Ei))%3B%0A%20%20%20%20let%20ymin%3Dp.ymin!%3D%3Dundefined%3Fp.ymin%3AInfinity%2Cymax%3Dp.ymax!%3D%3Dundefined%3Fp.ymax%3A-Infinity%3B%0A%20%20%20%20if(ymin%3D%3D%3DInfinity)for(const%20ya%20of%20ys_arr)for(const%20v%20of%20ya)%7Bif(v%3Cymin)ymin%3Dv%3Bif(v%3Eymax)ymax%3Dv%3B%7D%0A%20%20%20%20ymin%3Dp.ymin!%3D%3Dundefined%3Fp.ymin%3Aymin-Math.abs(ymax-ymin)*0.05%3B%0A%20%20%20%20ymax%3Dp.ymax!%3D%3Dundefined%3Fp.ymax%3Aymax%2BMath.abs(ymax-ymin)*0.05%3B%0A%20%20%20%20const%20xd%3D%7Bmin%3Ats%5B0%5D%2Cmax%3Ats%5Bn-1%5D%7C%7C1%2Clog%3Afalse%7D%2Cyd%3D%7Bmin%3Aymin%2Cmax%3Aymax%7C%7C1%2Clog%3A!!p.logy%7D%3B%0A%20%20%20%20const%7Bmx%2Cmy%2Cmw%2Cmh%2CtoX%2CtoY%7D%3Daxes2D(p%2Clx%2Cly%2Clw%2Clh%2Cxd%2Cyd%2Cvp)%3B%0A%20%20%20%20ctx.save()%3Bctx.beginPath()%3Bctx.rect(mx%2Cmy%2Cmw%2Cmh)%3Bctx.clip()%3B%0A%20%20%20%20const%20colors%3Dp.colors%7C%7C%5B%22%237dd3fc%22%2C%22%23f97316%22%2C%22%23a594f2%22%2C%22%235ed28c%22%2C%22%239c3b2c%22%5D%3B%0A%20%20%20%20const%20upTo%3Dp.grow%3D%3D%3Dfalse%3Fn%3AMath.min(f%2B1%2Cn)%3B%0A%20%20%20%20for(let%20si%3D0%3Bsi%3Cys_arr.length%3Bsi%2B%2B)%7B%0A%20%20%20%20%20%20const%20ya%3Dys_arr%5Bsi%5D%3Bctx.strokeStyle%3Dcolors%5Bsi%25colors.length%5D%3Bctx.lineWidth%3D2*DPR%3Bctx.beginPath()%3B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3CupTo%3Bi%2B%2B)%7Bconst%20x%3DtoX(ts%5Bi%5D)%2Cy%3Dyd.log%3FtoY(Math.log10(Math.max(ya%5Bi%5D%2C1e-30)))%3AtoY(ya%5Bi%5D)%3Bi%3D%3D%3D0%3Fctx.moveTo(x%2Cy)%3Actx.lineTo(x%2Cy)%3B%7D%0A%20%20%20%20%20%20ctx.stroke()%3B%0A%20%20%20%20%7D%0A%20%20%20%20%2F*%20moving%20dot%20on%20first%20series%20*%2F%0A%20%20%20%20const%20di%3Dp.dot!%3D%3Dundefined%3Fp.dot%3A0%3Bif(di%3Cys_arr.length%26%26upTo%3E0)%7B%0A%20%20%20%20%20%20const%20iv%3DMath.min(upTo-1%2Cn-1)%3Bconst%20x%3DtoX(ts%5Biv%5D)%2Cy%3Dyd.log%3FtoY(Math.log10(Math.max(ys_arr%5Bdi%5D%5Biv%5D%2C1e-30)))%3AtoY(ys_arr%5Bdi%5D%5Biv%5D)%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%23f97316%22%3Bctx.beginPath()%3Bctx.arc(x%2Cy%2C4.5*DPR%2C0%2C7)%3Bctx.fill()%3B%0A%20%20%20%20%7D%0A%20%20%20%20ctx.restore()%3B%0A%20%20%20%20%2F*%20vlines%20*%2F%0A%20%20%20%20if(p.vlines)for(const%20vl%20of%20p.vlines)%7Bconst%20x%3DtoX(vl.x)%3Bif(x%3Cmx%7C%7Cx%3Emx%2Bmw)continue%3Bctx.strokeStyle%3Dvl.color%7C%7C%22%233a4255%22%3Bctx.lineWidth%3D1%3Bctx.setLineDash(%5B4%2C3%5D)%3Bctx.beginPath()%3Bctx.moveTo(x%2Cmy)%3Bctx.lineTo(x%2Cmy%2Bmh)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3Bif(vl.label)%7Bctx.fillStyle%3Dvl.color%7C%7C%22%235a6880%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(vl.label%2Cx%2B3*DPR%2Cmy%2B10*DPR)%3B%7D%7D%0A%20%20%20%20if(p.labels%7C%7Cp.colors)legend(px%2Cpy%2Cpw%2Cph%2C(p.labels%7C%7C%5B%5D).map((l%2Ci)%3D%3E%5Bl%2C(p.colors%7C%7Ccolors)%5Bi%25colors.length%5D%5D))%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20SCATTER%20%2F%20PHASE_DIAG%20%2F%20BIFURCATION%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22scatter%22%7C%7Cp.kind%3D%3D%3D%22scatter_vline%22%7C%7Cp.kind%3D%3D%3D%22phase_diag%22%7C%7Cp.kind%3D%3D%3D%22bifurcation%22)%7B%0A%20%20%20%20const%20xs%3Dp.x%7C%7Cp.eta_sweep%7C%7Cp.K_sweep%7C%7C%5B%5D%3B%0A%20%20%20%20const%20ys%3Dp.y%7C%7Cp.phi_sweep%7C%7Cp.r_sweep%7C%7C%5B%5D%3B%0A%20%20%20%20let%20xmin%3Dp.xmin!%3D%3Dundefined%3Fp.xmin%3AMath.min(...xs)%3B%0A%20%20%20%20let%20xmax%3Dp.xmax!%3D%3Dundefined%3Fp.xmax%3AMath.max(...xs)%3B%0A%20%20%20%20let%20ymin%3Dp.ymin!%3D%3Dundefined%3Fp.ymin%3AMath.min(...ys)*0.95%3B%0A%20%20%20%20let%20ymax%3Dp.ymax!%3D%3Dundefined%3Fp.ymax%3AMath.max(...ys)*1.05%7C%7C1.05%3B%0A%20%20%20%20const%20xd%3D%7Bmin%3Axmin%2Cmax%3Axmax%2Clog%3Afalse%7D%2Cyd%3D%7Bmin%3Aymin%2Cmax%3Aymax%2Clog%3Afalse%7D%3B%0A%20%20%20%20const%7Bmx%2Cmy%2Cmw%2Cmh%2CtoX%2CtoY%7D%3Daxes2D(p%2Clx%2Cly%2Clw%2Clh%2Cxd%2Cyd%2Cvp)%3B%0A%20%20%20%20ctx.save()%3Bctx.beginPath()%3Bctx.rect(mx%2Cmy%2Cmw%2Cmh)%3Bctx.clip()%3B%0A%20%20%20%20%2F*%20sweep%20curve%20*%2F%0A%20%20%20%20if(xs.length%3E1)%7Bctx.strokeStyle%3Dp.color%7C%7C%22%239c3b2c%22%3Bctx.lineWidth%3D2.5*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3Cxs.length%3Bi%2B%2B)%7Bconst%20x%3DtoX(xs%5Bi%5D)%2Cy%3DtoY(ys%5Bi%5D)%3Bi%3D%3D%3D0%3Fctx.moveTo(x%2Cy)%3Actx.lineTo(x%2Cy)%3B%7Dctx.stroke()%3B%7D%0A%20%20%20%20%2F*%20static%20vlines%20(theory%20markers)%20*%2F%0A%20%20%20%20if(p.Kc!%3D%3Dundefined)%7Bconst%20x%3DtoX(p.Kc)%3Bctx.strokeStyle%3D%22%234a5568%22%3Bctx.lineWidth%3D1%3Bctx.setLineDash(%5B3%2C3%5D)%3Bctx.beginPath()%3Bctx.moveTo(x%2Cmy)%3Bctx.lineTo(x%2Cmy%2Bmh)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3Bctx.fillStyle%3D%22%234a5568%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(%22Kc%3D%22%2BfmtN(p.Kc)%2Cx%2Cmy%2B10*DPR)%3B%7D%0A%20%20%20%20if(p.vlines)for(const%20vl%20of%20p.vlines)%7Bconst%20x%3DtoX(vl.x)%3Bif(x%3Cmx%7C%7Cx%3Emx%2Bmw)continue%3Bctx.strokeStyle%3Dvl.color%7C%7C%22%234a5568%22%3Bctx.lineWidth%3D1%3Bctx.setLineDash(%5B3%2C3%5D)%3Bctx.beginPath()%3Bctx.moveTo(x%2Cmy)%3Bctx.lineTo(x%2Cmy%2Bmh)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3Bif(vl.label)%7Bctx.fillStyle%3Dvl.color%7C%7C%22%235a6880%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(vl.label%2Cx%2Cmy%2B10*DPR)%3B%7D%7D%0A%20%20%20%20%2F*%20current-frame%20dot%20%2B%20vline%20*%2F%0A%20%20%20%20const%20dotx%3Dp.dot_x%3Fp.dot_x%5Bf%5D%3A(p.eta_now%3Fp.eta_now%5Bf%5D%3A(p.K_now%3Fp.K_now%5Bf%5D%3Anull))%3B%0A%20%20%20%20const%20doty%3Dp.dot_y%3Fp.dot_y%5Bf%5D%3A(p.phi_now%3Fp.phi_now%5Bf%5D%3A(p.r_now%3Fp.r_now%5Bf%5D%3Anull))%3B%0A%20%20%20%20if(dotx!%3D%3Dnull%26%26dotx!%3D%3Dundefined)%7Bconst%20x%3DtoX(dotx)%3Bctx.strokeStyle%3D%22%23f97316%22%3Bctx.lineWidth%3D1%3Bctx.setLineDash(%5B2%2C2%5D)%3Bctx.beginPath()%3Bctx.moveTo(x%2Cmy)%3Bctx.lineTo(x%2Cmy%2Bmh)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3Bif(doty!%3D%3Dnull)%7Bconst%20y%3DtoY(doty)%3Bctx.fillStyle%3D%22%23f97316%22%3Bctx.beginPath()%3Bctx.arc(x%2Cy%2C5*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%7D%0A%20%20%20%20ctx.restore()%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20PHASES%20(Kuramoto%20circle)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22phases%22)%7B%0A%20%20%20%20const%20ths%3Dp.theta%5Bf%5D%2Crv%3Dp.r_series%3Fp.r_series%5Bf%5D%3A0%3B%0A%20%20%20%20const%20R%3DMath.min(pw%2Cph)*0.38%2Ccxc%3Dpx%2Bpw%2F2%2Ccyc%3Dpy%2Bph%2F2%3B%0A%20%20%20%20ctx.strokeStyle%3D%22%231e2535%22%3Bctx.lineWidth%3D1.5%3Bctx.beginPath()%3Bctx.arc(cxc%2Ccyc%2CR%2C0%2C7)%3Bctx.stroke()%3B%0A%20%20%20%20%2F*%20reference%20circle%20ticks%20at%200%2F90%2F180%2F270%20*%2F%0A%20%20%20%20ctx.fillStyle%3D%22%233a4255%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3B%0A%20%20%20%20for(const%20%5Ba%2Cl%5D%20of%20%5B%5B0%2C%220%22%5D%2C%5BMath.PI%2F2%2C%22%CF%80%2F2%22%5D%2C%5BMath.PI%2C%22%CF%80%22%5D%2C%5B3*Math.PI%2F2%2C%223%CF%80%2F2%22%5D%5D)%7B%0A%20%20%20%20%20%20ctx.fillText(l%2Ccxc%2BMath.cos(a)*(R%2B10*DPR)%2Ccyc-Math.sin(a)*(R%2B10*DPR)%2B3*DPR)%3B%0A%20%20%20%20%7D%0A%20%20%20%20for(const%20th%20of%20ths)%7Bctx.fillStyle%3D%22%237dd3fc%22%3Bctx.beginPath()%3Bctx.arc(cxc%2BR*Math.cos(th)%2Ccyc-R*Math.sin(th)%2C3*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%0A%20%20%20%20%2F*%20mean-field%20vector%20*%2F%0A%20%20%20%20const%20psi%3DMath.atan2(ths.reduce((a%2Ct)%3D%3Ea%2BMath.sin(t)%2C0)%2Cths.reduce((a%2Ct)%3D%3Ea%2BMath.cos(t)%2C0))%3B%0A%20%20%20%20ctx.strokeStyle%3D%22%23f97316%22%3Bctx.lineWidth%3D3*DPR%3Bctx.beginPath()%3Bctx.moveTo(cxc%2Ccyc)%3Bctx.lineTo(cxc%2BR*rv*Math.cos(psi)%2Ccyc-R*rv*Math.sin(psi))%3Bctx.stroke()%3B%0A%20%20%20%20ctx.fillStyle%3D%22%239aa3b0%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(%22r%3D%22%2Brv.toFixed(3)%2Cpx%2B5*DPR%2Cpy%2B14*DPR)%3B%0A%20%20%20%20if(p.title)%7Bctx.fillStyle%3D%22%23b3bac8%22%3Bctx.fillText(p.title%2Cpx%2B5*DPR%2Cpy%2B26*DPR)%3B%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20HEATMAP%20(with%20colorbar)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22heatmap%22%7C%7Cp.kind%3D%3D%3D%22heatmap_pts%22)%7B%0A%20%20%20%20const%20n%3Dp.n%2Cg%3Dp.grid%5Bf%5D%3B%0A%20%20%20%20if(!p._oc)%7Bp._oc%3Ddocument.createElement(%22canvas%22)%3Bp._oc.width%3Dn%3Bp._oc.height%3Dn%3Bp._octx%3Dp._oc.getContext(%222d%22)%3Bp._img%3Dp._octx.createImageData(n%2Cn)%3B%7D%0A%20%20%20%20const%20d%3Dp._img.data%3B%0A%20%20%20%20if(p.kind%3D%3D%3D%22heatmap%22)%7B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3Cn*n%3Bi%2B%2B)%7Bconst%20v%3Dg%5Bi%5D%2F255%3Bd%5B4*i%5D%3DMath.round(255*(0.06%2B0.35*v))%3Bd%5B4*i%2B1%5D%3DMath.round(255*(0.18%2B0.55*v))%3Bd%5B4*i%2B2%5D%3DMath.round(255*(0.25%2B0.55*v))%3Bd%5B4*i%2B3%5D%3D255%3B%7D%0A%20%20%20%20%7D%20else%20%7B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3Cn*n%3Bi%2B%2B)%7Bconst%20v%3Dg%5Bi%5D%2F255%3Bd%5B4*i%5D%3DMath.round(255*(0.18%2B0.60*v))%3Bd%5B4*i%2B1%5D%3DMath.round(255*(0.42-0.08*v))%3Bd%5B4*i%2B2%5D%3DMath.round(255*(0.75-0.42*v))%3Bd%5B4*i%2B3%5D%3D255%3B%7D%0A%20%20%20%20%7D%0A%20%20%20%20p._octx.putImageData(p._img%2C0%2C0)%3B%0A%20%20%20%20const%20cbW%3D16*DPR%2Cs%3DMath.min(pw-cbW-4*DPR%2Cph)*0.9%2Cox%3Dpx%2B(pw-cbW-4*DPR-s)%2F2%2Coy%3Dpy%2B(ph-s)%2F2%3B%0A%20%20%20%20ctx.imageSmoothingEnabled%3Dtrue%3Bctx.drawImage(p._oc%2Cox%2Coy%2Cs%2Cs)%3B%0A%20%20%20%20%2F*%20colorbar%20*%2F%0A%20%20%20%20const%20cbX%3Dox%2Bs%2B4*DPR%2CcbY%3Doy%2CcbH%3Ds%3B%0A%20%20%20%20const%20grad%3Dctx.createLinearGradient(0%2CcbY%2BcbH%2C0%2CcbY)%3B%0A%20%20%20%20if(p.kind%3D%3D%3D%22heatmap%22)%7Bgrad.addColorStop(0%2C%22rgb(15%2C46%2C64)%22)%3Bgrad.addColorStop(0.5%2C%22rgb(26%2C88%2C100)%22)%3Bgrad.addColorStop(1%2C%22rgb(72%2C171%2C156)%22)%3B%7D%0A%20%20%20%20else%7Bgrad.addColorStop(0%2C%22rgb(46%2C107%2C191)%22)%3Bgrad.addColorStop(0.5%2C%22rgb(170%2C142%2C191)%22)%3Bgrad.addColorStop(1%2C%22rgb(237%2C156%2C192)%22)%3B%7D%0A%20%20%20%20ctx.fillStyle%3Dgrad%3Bctx.fillRect(cbX%2CcbY%2CcbW%2CcbH)%3B%0A%20%20%20%20ctx.strokeStyle%3D%22%233a4255%22%3Bctx.lineWidth%3D1%3Bctx.strokeRect(cbX%2CcbY%2CcbW%2CcbH)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3B%0A%20%20%20%20ctx.fillText(%221%22%2CcbX%2BcbW%2B3*DPR%2CcbY%2B8*DPR)%3Bctx.fillText(%220%22%2CcbX%2BcbW%2B3*DPR%2CcbY%2BcbH%2B4*DPR)%3B%0A%20%20%20%20%2F*%20data%20points%20overlay%20*%2F%0A%20%20%20%20if(p.pts)for(const%20pt%20of%20p.pts)%7Bconst%20cx2%3Dox%2Bpt%5B0%5D*s%2Ccy2%3Doy%2Bpt%5B1%5D*s%3Bctx.fillStyle%3Dpt%5B2%5D%3F%22%230f1117%22%3A%22%23ffffff%22%3Bctx.beginPath()%3Bctx.arc(cx2%2Ccy2%2C3*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%0A%20%20%20%20if(p.title)%7Bctx.fillStyle%3D%22%239aa3b0%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(p.title%2Cpx%2Bpw%2F2%2Cpy%2B12*DPR)%3B%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20RNT%20(walker%20paths%20%2B%20MSD%20log-log)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22rnt%22)%7B%0A%20%20%20%20const%20lW%3Dpw*0.46%2CL2%3Dp.L%2Ccxw%3Dpx%2BlW%2F2%2Ccyw%3Dpy%2Bph%2F2%3B%0A%20%20%20%20const%20sc%3DMath.min(lW%2Cph)*0.85%2FL2%3B%0A%20%20%20%20if(p.conf%3E0)%7Bctx.strokeStyle%3D%22%232f3645%22%3Bctx.lineWidth%3D1%3Bctx.beginPath()%3Bctx.arc(cxw%2Ccyw%2Cp.conf*sc%2C0%2C7)%3Bctx.stroke()%3B%7D%0A%20%20%20%20const%20upTo%3DMath.min(f%2B1%2Cp.paths%5B0%5D.length)%3B%0A%20%20%20%20for(let%20w%3D0%3Bw%3Cp.paths.length%3Bw%2B%2B)%7B%0A%20%20%20%20%20%20const%20clr%3D%60hsla(%24%7Bw*47%7D%2C55%25%2C68%25%2C0.6)%60%3Bctx.strokeStyle%3Dclr%3Bctx.lineWidth%3D1%3Bctx.beginPath()%3B%0A%20%20%20%20%20%20for(let%20i%3DMath.max(0%2CupTo-80)%3Bi%3CupTo%3Bi%2B%2B)%7Bconst%5Bwx%2Cwy%5D%3Dp.paths%5Bw%5D%5Bi%5D%3Bconst%20sx%3Dcxw%2Bwx*sc%2Csy%3Dcyw-wy*sc%3Bi%3D%3D%3DMath.max(0%2CupTo-80)%3Fctx.moveTo(sx%2Csy)%3Actx.lineTo(sx%2Csy)%3B%7D%0A%20%20%20%20%20%20ctx.stroke()%3B%0A%20%20%20%20%20%20if(upTo%3E0)%7Bconst%5Bwx%2Cwy%5D%3Dp.paths%5Bw%5D%5BupTo-1%5D%3Bctx.fillStyle%3Dclr%3Bctx.beginPath()%3Bctx.arc(cxw%2Bwx*sc%2Ccyw-wy*sc%2C3*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%0A%20%20%20%20%7D%0A%20%20%20%20ctx.fillStyle%3D%22%235a6880%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(%22walkers%22%2Cpx%2B3*DPR%2Cpy%2B12*DPR)%3B%0A%20%20%20%20%2F*%20MSD%20right%20panel%20with%20proper%20log%20axes%20*%2F%0A%20%20%20%20const%20rPnl%3D%7Bkind%3A%22timeseries%22%2Ctitle%3A%22MSD%20vs%20time%22%2Cxlabel%3A%22time%20t%22%2Cylabel%3A%22%E2%9F%A8r%C2%B2%E2%9F%A9%22%2Clogy%3Atrue%2Cgrow%3Atrue%2Cdot%3A0%2Cys%3A%5Bp.msd_sim%2Cp.msd_th%5D%2Ccolors%3A%5B%22%237dd3fc%22%2C%22%239c3b2c%22%5D%2Clabels%3A%5B%22simulation%22%2C%22theory%22%5D%2Ct%3Ap.t_msd%2Cymin%3AMath.min(...p.msd_sim.filter(v%3D%3Ev%3E0))*0.5%2Cymax%3AMath.max(...p.msd_sim)*2%7D%3B%0A%20%20%20%20drawPanel(rPnl%2Cpi%2Clx%2Blw*0.48%2Cly%2Clw*0.52%2Clh)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%203D%20WAVEPACKET%20HELIX%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22wavepacket_3d%22)%7B%0A%20%20%20%20const%20f_%3DMath.min(frame%2CS.frames-1)%3B%0A%20%20%20%20const%20re%3Dp.psi_re%5Bf_%5D%2Cim%3Dp.psi_im%5Bf_%5D%2Cprob%3Dp.prob%5Bf_%5D%2CV3%3Dp.V3d%3B%0A%20%20%20%20const%20xArr%3Dp.x%2Ccom%3Dp.com_f%5Bf_%5D%2Chalfw%3Dp.halfw_f%5Bf_%5D%3B%0A%20%20%20%20const%20amp%3Dp.amp%2Cpmax%3Dp.pmax%3B%0A%20%20%20%20const%20cxp%3Dpx%2Bpw*0.5%2Ccyp%3Dpy%2Bph*0.52%2Csc3%3DMath.min(pw%2Cph)*0.34%3B%0A%20%20%20%20function%20wxOf(i)%7Breturn%20Math.max(-1.05%2CMath.min(1.05%2C(xArr%5Bi%5D-com)%2Fhalfw))%3B%7D%0A%20%20%20%20function%20Pq(wx%2Cwy%2Cwz)%7Breturn%20P3(wx%2Cwy%2Cwz%2Ccxp%2Ccyp%2Csc3)%3B%7D%0A%20%20%20%20%2F*%20box%20%2B%20grid%20*%2F%0A%20%20%20%20const%20C%3D%5B%5D%3Bfor(const%20sx%20of%5B-1%2C1%5D)for(const%20sy%20of%5B-1%2C1%5D)for(const%20sz%20of%5B-1%2C1%5D)C.push(Pq(sx%2Csy%2Csz))%3B%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(255%2C255%2C255%2C0.09)%22%3Bctx.lineWidth%3D1%3B%0A%20%20%20%20%5B%5B0%2C1%5D%2C%5B0%2C2%5D%2C%5B0%2C4%5D%2C%5B1%2C3%5D%2C%5B1%2C5%5D%2C%5B2%2C3%5D%2C%5B2%2C6%5D%2C%5B3%2C7%5D%2C%5B4%2C5%5D%2C%5B4%2C6%5D%2C%5B5%2C7%5D%2C%5B6%2C7%5D%5D.forEach((%5Ba%2Cb%5D)%3D%3E%7Bctx.beginPath()%3Bctx.moveTo(C%5Ba%5D%5B0%5D%2CC%5Ba%5D%5B1%5D)%3Bctx.lineTo(C%5Bb%5D%5B0%5D%2CC%5Bb%5D%5B1%5D)%3Bctx.stroke()%3B%7D)%3B%0A%20%20%20%20%2F*%20floor%20grid%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(255%2C255%2C255%2C0.07)%22%3B%0A%20%20%20%20for(let%20g%3D-1%3Bg%3C%3D1%3Bg%2B%3D0.5)%7Bconst%20a%3DPq(g%2C0%2C-1)%2Cb%3DPq(g%2C0%2C1)%3Bctx.beginPath()%3Bctx.moveTo(a%5B0%5D%2Ca%5B1%5D)%3Bctx.lineTo(b%5B0%5D%2Cb%5B1%5D)%3Bctx.stroke()%3Bconst%20c%3DPq(-1%2C0%2Cg)%2Cd%3DPq(1%2C0%2Cg)%3Bctx.beginPath()%3Bctx.moveTo(c%5B0%5D%2Cc%5B1%5D)%3Bctx.lineTo(d%5B0%5D%2Cd%5B1%5D)%3Bctx.stroke()%3B%7D%0A%20%20%20%20%2F*%20axis%20labels%20*%2F%0A%20%20%20%20ctx.fillStyle%3D%22%235a6880%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3B%0A%20%20%20%20const%20lx2%3DPq(1.2%2C0%2C0)%2Clim%3DPq(0%2C0%2C1.25)%2Clre%3DPq(-1.1%2C1.3%2C0)%3B%0A%20%20%20%20ctx.fillText(%22x%22%2Clx2%5B0%5D%2Clx2%5B1%5D)%3Bctx.fillText(%22Im(%CF%88)%22%2Clim%5B0%5D%2Clim%5B1%5D)%3Bctx.fillText(%22Re(%CF%88)%22%2Clre%5B0%5D%2Clre%5B1%5D)%3B%0A%20%20%20%20%2F*%20V%20on%20floor%20*%2F%0A%20%20%20%20let%20vmax%3D0%3Bif(V3)for(let%20i%3D0%3Bi%3CV3.length%3Bi%2B%2B)if(Math.abs(V3%5Bi%5D)%3Evmax)vmax%3DMath.abs(V3%5Bi%5D)%3B%0A%20%20%20%20if(vmax%3E1e-9)%7Bctx.strokeStyle%3D%22rgba(0%2C217%2C217%2C0.5)%22%3Bctx.lineWidth%3D1.5*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CxArr.length%3Bi%2B%2B)%7Bconst%20q%3DPq(wxOf(i)%2C0%2CV3%5Bi%5D%2Fvmax*0.65)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%7D%0A%20%20%20%20%2F*%20%7Cpsi%7C%5E2%20(amber)%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22%23f5992e%22%3Bctx.lineWidth%3D2.4*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CxArr.length%3Bi%2B%2B)%7Bconst%20q%3DPq(wxOf(i)%2C0%2Cprob%5Bi%5D%2Fpmax)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20%2F*%20Im(psi)%20green%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(51%2C191%2C128%2C0.9)%22%3Bctx.lineWidth%3D2.2*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CxArr.length%3Bi%2B%2B)%7Bconst%20q%3DPq(wxOf(i)%2Cim%5Bi%5D%2Famp%2C0)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20%2F*%20Re(psi)%20blue%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(90%2C133%2C230%2C0.9)%22%3Bctx.lineWidth%3D2.2*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CxArr.length%3Bi%2B%2B)%7Bconst%20q%3DPq(wxOf(i)%2C0%2Cre%5Bi%5D%2Famp)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20%2F*%20helix%20violet%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(179%2C77%2C242%2C0.92)%22%3Bctx.lineWidth%3D2.5*DPR%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CxArr.length%3Bi%2B%2B)%7Bconst%20q%3DPq(wxOf(i)%2Cim%5Bi%5D%2Famp%2Cre%5Bi%5D%2Famp)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20legend(px%2Cpy%2Cpw%2Cph%2C%5B%5B%22%CF%88%20helix%22%2C%22%23b34df2%22%5D%2C%5B%22Re(%CF%88)%22%2C%22%235a85e6%22%5D%2C%5B%22Im(%CF%88)%22%2C%22%2333bf80%22%5D%2C%5B%22%7C%CF%88%7C%C2%B2%22%2C%22%23f5992e%22%5D%5D)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%234a5568%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(%22drag%20to%20orbit%20%C2%B7%20scroll%20zoom%20%C2%B7%20dblclick%20reset%22%2Cpx%2B5*DPR%2Cpy%2Bph-8*DPR)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20PROB_DENSITY%20(2D%20wavepacket)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22prob_density%22)%7B%0A%20%20%20%20const%20prob%3Dp.prob%5Bf%5D%2CV2%3Dp.V2d%2CxArr%3Dp.x%3B%0A%20%20%20%20const%20dx_nm%3DxArr%5B1%5D-xArr%5B0%5D%3B%20%2F*%20in%20m%3B%20display%20in%20nm%20*%2F%0A%20%20%20%20const%20xnm%3DxArr.map(v%3D%3Ev*1e9)%3B%20%2F*%20%E2%86%92%20nm%20*%2F%0A%20%20%20%20const%20pnl%3D%7Bkind%3A%22timeseries%22%2Ctitle%3Ap.title%7C%7C%22%7C%CF%88(x%2Ct)%7C%C2%B2%20%20(amber)%20%20%C2%B7%20%20V(x)%20scaled%20(cyan)%22%2Cxlabel%3A%22position%20x%20(nm)%22%2Cylabel%3A%22amplitude%22%2Cymin%3A0%2Cymax%3AMath.max(...prob)*1.15%2Cys%3A%5Bprob%5D%2Ccolors%3A%5B%22%23f5992e%22%5D%2Clabels%3A%5B%22%7C%CF%88%7C%C2%B2%22%5D%2Ct%3Axnm%2Cgrow%3Afalse%7D%3B%0A%20%20%20%20%2F*%20inject%20V%20overlay%20before%20rendering%20*%2F%0A%20%20%20%20const%20savedKind%3Dp.kind%3B%0A%20%20%20%20drawPanel(pnl%2Cpi%2Clx%2Cly%2Clw%2Clh)%3B%0A%20%20%20%20%2F*%20now%20overlay%20V%20manually%20*%2F%0A%20%20%20%20const%20xd2%3D%7Bmin%3Axnm%5B0%5D%2Cmax%3Axnm%5Bxnm.length-1%5D%2Clog%3Afalse%7D%3B%0A%20%20%20%20const%20yd2%3D%7Bmin%3A0%2Cmax%3AMath.max(...prob)*1.15%2Clog%3Afalse%7D%3B%0A%20%20%20%20const%20ML%3D46*DPR%2CMR%3D8*DPR%2CMT%3D22*DPR%2CMB%3D28*DPR%3B%0A%20%20%20%20const%20mx2%3DMath.round(lx*cv.width)%2BML%2Cmy2%3DMath.round(ly*cv.height)%2BMT%2Cmw2%3DMath.round(lw*cv.width)-ML-MR%2Cmh2%3DMath.round(lh*cv.height)-MT-MB%3B%0A%20%20%20%20const%20toX2%3Dv%3D%3Emx2%2Bmw2*(v-xd2.min)%2F(xd2.max-xd2.min)%3B%0A%20%20%20%20const%20toY2%3Dv%3D%3Emy2%2Bmh2*(1-v%2Fyd2.max)%3B%0A%20%20%20%20ctx.save()%3Bctx.beginPath()%3Bctx.rect(mx2%2Cmy2%2Cmw2%2Cmh2)%3Bctx.clip()%3B%0A%20%20%20%20if(V2%26%26V2.some(v%3D%3EMath.abs(v)%3E1e-9))%7B%0A%20%20%20%20%20%20let%20vmax%3D0%3Bfor(const%20v%20of%20V2)if(Math.abs(v)%3Evmax)vmax%3DMath.abs(v)%3B%0A%20%20%20%20%20%20const%20vscl%3Dyd2.max*0.5%2Fvmax%3B%0A%20%20%20%20%20%20ctx.strokeStyle%3D%22rgba(0%2C217%2C217%2C0.5)%22%3Bctx.lineWidth%3D1.5*DPR%3Bctx.beginPath()%3B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3Cxnm.length%3Bi%2B%2B)%7Bconst%20X%3DtoX2(xnm%5Bi%5D)%2CY%3DtoY2(Math.max(0%2CV2%5Bi%5D*vscl))%3Bi%3D%3D%3D0%3Fctx.moveTo(X%2CY)%3Actx.lineTo(X%2CY)%3B%7Dctx.stroke()%3B%0A%20%20%20%20%7D%0A%20%20%20%20if(p.markers)for(const%20m%20of%20p.markers)%7Bconst%20X%3DtoX2(m*1e9)%3Bctx.strokeStyle%3D%22rgba(249%2C115%2C22%2C0.55)%22%3Bctx.lineWidth%3D1.1*DPR%3Bctx.setLineDash(%5B3*DPR%2C3*DPR%5D)%3Bctx.beginPath()%3Bctx.moveTo(X%2Cmy2)%3Bctx.lineTo(X%2Cmy2%2Bmh2)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3B%7D%0A%20%20%20%20ctx.restore()%3B%0A%20%20%20%20legend(Math.round(lx*cv.width)%2CMath.round(ly*cv.height)%2CMath.round(lw*cv.width)%2CMath.round(lh*cv.height)%2C%5B%5B%22V(x)%22%2C%22rgba(0%2C217%2C217%2C0.7)%22%5D%2C%5B%22barrier%22%2C%22rgba(249%2C115%2C22%2C0.7)%22%5D%5D)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20CHAIN%20(protein%20folding)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22chain%22)%7B%0A%20%20%20%20const%20xs%3Dp.x%5Bf%5D%2Cys%3Dp.y%5Bf%5D%2Cn%3Dxs.length%3B%0A%20%20%20%20let%20xmin%3DInfinity%2Cxmax%3D-Infinity%2Cymin%3DInfinity%2Cymax%3D-Infinity%3B%0A%20%20%20%20for(let%20i%3D0%3Bi%3Cn%3Bi%2B%2B)%7Bif(xs%5Bi%5D%3Cxmin)xmin%3Dxs%5Bi%5D%3Bif(xs%5Bi%5D%3Exmax)xmax%3Dxs%5Bi%5D%3Bif(ys%5Bi%5D%3Cymin)ymin%3Dys%5Bi%5D%3Bif(ys%5Bi%5D%3Eymax)ymax%3Dys%5Bi%5D%3B%7D%0A%20%20%20%20const%20span%3DMath.max(xmax-xmin%2Cymax-ymin%2C1)*1.4%3B%0A%20%20%20%20const%20cxc%3D(xmin%2Bxmax)%2F2%2Ccyc%3D(ymin%2Bymax)%2F2%3B%0A%20%20%20%20const%20sc%3DMath.min(pw%2Cph)*0.85%2Fspan%2Coxc%3Dpx%2Bpw%2F2%2Coyc%3Dpy%2Bph%2F2%3B%0A%20%20%20%20function%20ts(x%2Cy)%7Breturn%5Boxc%2B(x-cxc)*sc%2Coyc-(y-cyc)*sc%5D%3B%7D%0A%20%20%20%20ctx.strokeStyle%3D%22%234a5568%22%3Bctx.lineWidth%3D2*DPR%3Bctx.beginPath()%3B%0A%20%20%20%20for(let%20i%3D0%3Bi%3Cn-1%3Bi%2B%2B)%7Bconst%5Bax%2Cay%5D%3Dts(xs%5Bi%5D%2Cys%5Bi%5D)%3Bconst%5Bbx%2Cby%5D%3Dts(xs%5Bi%2B1%5D%2Cys%5Bi%2B1%5D)%3Bi%3D%3D%3D0%3Fctx.moveTo(ax%2Cay)%3Actx.lineTo(bx%2Cby)%3B%7Dctx.stroke()%3B%0A%20%20%20%20if(p.native)%7Bctx.strokeStyle%3D%22%239c3b2c%22%3Bctx.lineWidth%3D0.8*DPR%3Bctx.setLineDash(%5B2%2C2%5D)%3Bctx.beginPath()%3Bfor(const%5Ba%2Cb%5Dof%20p.native)%7Bconst%5Bax%2Cay%5D%3Dts(xs%5Ba%5D%2Cys%5Ba%5D)%3Bconst%5Bbx%2Cby%5D%3Dts(xs%5Bb%5D%2Cys%5Bb%5D)%3Bctx.moveTo(ax%2Cay)%3Bctx.lineTo(bx%2Cby)%3B%7Dctx.stroke()%3Bctx.setLineDash(%5B%5D)%3B%7D%0A%20%20%20%20for(let%20i%3D0%3Bi%3Cn%3Bi%2B%2B)%7Bconst%5Bsx%2Csy%5D%3Dts(xs%5Bi%5D%2Cys%5Bi%5D)%3Bctx.fillStyle%3Di%3D%3D%3D0%3F%22%23f97316%22%3Ai%3D%3D%3Dn-1%3F%22%23a855f7%22%3A%22%237dd3fc%22%3Bctx.beginPath()%3Bctx.arc(sx%2Csy%2C4*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%0A%20%20%20%20if(p.title)%7Bctx.fillStyle%3D%22%239aa3b0%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(p.title%2Cpx%2Bpw%2F2%2Cpy%2B12*DPR)%3B%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20FREE%20ENERGY%20F(Q)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22freeenergy%22)%7B%0A%20%20%20%20const%20qb%3Dp.Q_bins%2Cfv%3Dp.F_vals%2Cqn%3Dp.Q_now%3Fp.Q_now%5Bf%5D%3Anull%3B%0A%20%20%20%20const%20ymin%3DMath.min(...fv)%2Cymax%3DMath.max(...fv)%2B0.5%3B%0A%20%20%20%20const%20pnl%3D%7Bkind%3A%22scatter%22%2Ctitle%3Ap.title%7C%7C%22Free%20energy%20F(Q)%22%2Cxlabel%3A%22Q%20(native%20fraction)%22%2Cylabel%3A%22F%20(kT)%22%2Cx%3Aqb%2Cy%3Afv%2Cymin%2Cymax%2Ccolor%3A%22%239c3b2c%22%7D%3B%0A%20%20%20%20drawPanel(pnl%2Cpi%2Clx%2Cly%2Clw%2Clh)%3B%0A%20%20%20%20if(qn!%3D%3Dnull)%7B%2F*%20overlay%20moving%20dot%20*%2F%0A%20%20%20%20%20%20const%20ML%3D46*DPR%2CMT%3D22*DPR%2Cmw2%3DMath.round(lw*cv.width)-ML-8*DPR%2Cmh2%3DMath.round(lh*cv.height)-MT-28*DPR%3B%0A%20%20%20%20%20%20const%20mx2%3DMath.round(lx*cv.width)%2BML%2Cmy2%3DMath.round(ly*cv.height)%2BMT%3B%0A%20%20%20%20%20%20const%20toX2%3Dv%3D%3Emx2%2Bmw2*(v-qb%5B0%5D)%2F(qb%5Bqb.length-1%5D-qb%5B0%5D)%3B%0A%20%20%20%20%20%20const%20toY2%3Dv%3D%3Emy2%2Bmh2*(1-(v-ymin)%2F(ymax-ymin))%3B%0A%20%20%20%20%20%20%2F*%20find%20F%20at%20Q%3Dqn%20by%20interpolation%20*%2F%0A%20%20%20%20%20%20let%20fi%3D0%3Bfor(let%20i%3D0%3Bi%3Cqb.length%3Bi%2B%2B)if(qb%5Bi%5D%3C%3Dqn)fi%3Di%3B%0A%20%20%20%20%20%20const%20fAtQ%3Dfv%5Bfi%5D%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%23f97316%22%3Bctx.beginPath()%3Bctx.arc(toX2(qn)%2CtoY2(fAtQ)%2C5*DPR%2C0%2C7)%3Bctx.fill()%3B%0A%20%20%20%20%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20SWEEP%20CURVE%20(percolation%20S(p)%20%2F%20chi(p))%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22sweep_curve%22)%7B%0A%20%20%20%20const%20data%3Dp.data%2Cpv%3DS.p_arr%7C%7Cp.p_arr%7C%7C%5B%5D%2Cn%3Ddata.length%3B%0A%20%20%20%20const%20ymax%3DMath.max(...data)*1.1%2B0.01%3B%0A%20%20%20%20const%20pnl%3D%7Bkind%3A%22timeseries%22%2Ctitle%3Ap.title%2Cxlabel%3Ap.xlabel%7C%7C%22bond%20prob.%20p%22%2Cylabel%3Ap.ylabel%2Cymin%3A0%2Cymax%2Cys%3A%5Bdata%5D%2Ccolors%3A%5Bp.color%7C%7C%22%237dd3fc%22%5D%2Cgrow%3Atrue%2Ct%3Apv%7D%3B%0A%20%20%20%20drawPanel(pnl%2Cpi%2Clx%2Cly%2Clw%2Clh)%3B%0A%20%20%20%20%2F*%20p_c%20vline%20%2B%20current%20dot%20*%2F%0A%20%20%20%20if(p.pc!%3D%3Dnull%26%26p.pc!%3D%3Dundefined)%7B%0A%20%20%20%20%20%20const%20ML%3D46*DPR%2CMT%3D22*DPR%2Cmw2%3DMath.round(lw*cv.width)-ML-8*DPR%2Cmh2%3DMath.round(lh*cv.height)-MT-28*DPR%3B%0A%20%20%20%20%20%20const%20mx2%3DMath.round(lx*cv.width)%2BML%2Cmy2%3DMath.round(ly*cv.height)%2BMT%3B%0A%20%20%20%20%20%20const%20toX2%3Dv%3D%3Emx2%2Bmw2*(v-pv%5B0%5D)%2F(pv%5Bpv.length-1%5D-pv%5B0%5D)%3B%0A%20%20%20%20%20%20const%20toY2%3Dv%3D%3Emy2%2Bmh2*(1-v%2Fymax)%3B%0A%20%20%20%20%20%20const%20vx%3DtoX2(p.pc)%3Bctx.strokeStyle%3D%22rgba(139%2C148%2C158%2C0.4)%22%3Bctx.lineWidth%3D1%3Bctx.setLineDash(%5B4%2C4%5D)%3Bctx.beginPath()%3Bctx.moveTo(vx%2Cmy2)%3Bctx.lineTo(vx%2Cmy2%2Bmh2)%3Bctx.stroke()%3Bctx.setLineDash(%5B%5D)%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22center%22%3Bctx.fillText(%22p_c%3D%22%2BfmtN(p.pc)%2Cvx%2Cmy2%2B10*DPR)%3B%0A%20%20%20%20%20%20const%20ci%3DMath.min(f%2Cn-1)%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%23f97316%22%3Bctx.beginPath()%3Bctx.arc(toX2(pv%5Bci%5D)%2CtoY2(data%5Bci%5D)%2C4.5*DPR%2C0%2C7)%3Bctx.fill()%3B%0A%20%20%20%20%7D%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20GENERIC%20CURVE%20(CD%20spectrum%2C%20etc.)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22curve%22)%7B%0A%20%20%20%20const%20yy%3Dp.y%5BMath.min(f%2Cp.y.length-1)%5D%2Cymax%3Dp.ymax%7C%7C1%2Cymin%3Dp.ymin%7C%7C-ymax%3B%0A%20%20%20%20const%20wl0%3Dp.wl0%7C%7C190%2Cwl1%3Dp.wl1%7C%7C250%2CN2%3Dyy.length%3B%0A%20%20%20%20const%20pnl%3D%7Bkind%3A%22timeseries%22%2Ctitle%3Ap.title%7C%7Cp.label%3F.%5Bf%5D%7C%7C%22%22%2Cxlabel%3Ap.xlabel%7C%7C%22wavelength%20(nm)%22%2Cylabel%3Ap.ylabel%7C%7C%22CD%22%2Cymin%2Cymax%2Cys%3A%5Byy%5D%2Ccolors%3A%5B%22%237dd3fc%22%5D%2Cgrow%3Afalse%2Ct%3AArray.from(%7Blength%3AN2%7D%2C(_%2Ci)%3D%3Ewl0%2B(wl1-wl0)*i%2F(N2-1))%7D%3B%0A%20%20%20%20drawPanel(pnl%2Cpi%2Clx%2Cly%2Clw%2Clh)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20ATTRACTOR%203D%20(reservoir%20computing)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(p.kind%3D%3D%3D%22attractor_esn%22%7C%7Cp.kind%3D%3D%3D%22attractor%22)%7B%0A%20%20%20%20%2F*%20handled%20inline%20in%20reservoir_computing%20%E2%80%94%20kept%20here%20as%20fallback%20*%2F%0A%20%20%20%20const%20cxp%3Dpx%2Bpw*0.5%2Ccyp%3Dpy%2Bph*0.52%2Csc3%3DMath.min(pw%2Cph)*0.38%3B%0A%20%20%20%20function%20NP(ax%2Cay%2Caz_%2Ci)%7Breturn%20P3((ax%5Bi%5D-S.cx)%2FS.hs%2C(ay%5Bi%5D-S.cy)%2FS.hs%2C(az_%5Bi%5D-S.cz)%2FS.hs%2Ccxp%2Ccyp%2Csc3)%3B%7D%0A%20%20%20%20%2F*%20box%20*%2F%0A%20%20%20%20const%20C%3D%5B%5D%3Bfor(const%20sx%20of%5B-1%2C1%5D)for(const%20sy%20of%5B-1%2C1%5D)for(const%20sz%20of%5B-1%2C1%5D)C.push(P3(sx%2Csy%2Csz%2Ccxp%2Ccyp%2Csc3))%3B%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(255%2C255%2C255%2C0.07)%22%3Bctx.lineWidth%3D1%3B%0A%20%20%20%20%5B%5B0%2C1%5D%2C%5B0%2C2%5D%2C%5B0%2C4%5D%2C%5B1%2C3%5D%2C%5B1%2C5%5D%2C%5B2%2C3%5D%2C%5B2%2C6%5D%2C%5B3%2C7%5D%2C%5B4%2C5%5D%2C%5B4%2C6%5D%2C%5B5%2C7%5D%2C%5B6%2C7%5D%5D.forEach((%5Ba%2Cb%5D)%3D%3E%7Bctx.beginPath()%3Bctx.moveTo(C%5Ba%5D%5B0%5D%2CC%5Ba%5D%5B1%5D)%3Bctx.lineTo(C%5Bb%5D%5B0%5D%2CC%5Bb%5D%5B1%5D)%3Bctx.stroke()%3B%7D)%3B%0A%20%20%20%20const%20TRAIL%3D120%2CDELAY%3D20%3B%0A%20%20%20%20%2F*%20ghost%20*%2F%0A%20%20%20%20ctx.strokeStyle%3D%22rgba(180%2C190%2C210%2C0.07)%22%3Bctx.lineWidth%3D1%3Bctx.beginPath()%3Bfor(let%20i%3D0%3Bi%3CS.ghost_x.length%3Bi%2B%2B)%7Bconst%20q%3DNP(S.ghost_x%2CS.ghost_y%2CS.ghost_z%2Ci)%3Bi%3D%3D%3D0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20%2F*%20true%20*%2F%0A%20%20%20%20const%20t0%3DMath.max(0%2Cf-TRAIL)%3Bctx.strokeStyle%3D%22%237dd3fc%22%3Bctx.lineWidth%3D2*DPR%3Bctx.beginPath()%3Bfor(let%20i%3Dt0%3Bi%3C%3Df%3Bi%2B%2B)%7Bconst%20q%3DNP(S.true_x%2CS.true_y%2CS.true_z%2Ci)%3Bi%3D%3D%3Dt0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3B%0A%20%20%20%20const%20tq%3DNP(S.true_x%2CS.true_y%2CS.true_z%2Cf)%3Bctx.fillStyle%3D%22%237dd3fc%22%3Bctx.beginPath()%3Bctx.arc(tq%5B0%5D%2Ctq%5B1%5D%2C4.5*DPR%2C0%2C7)%3Bctx.fill()%3B%0A%20%20%20%20%2F*%20esn%20*%2F%0A%20%20%20%20const%20fd%3DMath.max(0%2Cf-DELAY)%2Ce0%3DMath.max(0%2Cfd-TRAIL)%3B%0A%20%20%20%20if(fd%3E0)%7Bctx.strokeStyle%3D%22%23f97316%22%3Bctx.lineWidth%3D2*DPR%3Bctx.beginPath()%3Bfor(let%20i%3De0%3Bi%3C%3Dfd%3Bi%2B%2B)%7Bconst%20q%3DNP(S.esn_x%2CS.esn_y%2CS.esn_z%2Ci)%3Bi%3D%3D%3De0%3Fctx.moveTo(q%5B0%5D%2Cq%5B1%5D)%3Actx.lineTo(q%5B0%5D%2Cq%5B1%5D)%3B%7Dctx.stroke()%3Bconst%20eq%3DNP(S.esn_x%2CS.esn_y%2CS.esn_z%2Cfd)%3Bctx.fillStyle%3D%22%23f97316%22%3Bctx.beginPath()%3Bctx.arc(eq%5B0%5D%2Ceq%5B1%5D%2C4.5*DPR%2C0%2C7)%3Bctx.fill()%3B%7D%0A%20%20%20%20legend(px%2Cpy%2Cpw%2Cph%2C%5B%5B%22true%22%2C%22%237dd3fc%22%5D%2C%5B%22ESN%20(20%20frames%20behind)%22%2C%22%23f97316%22%5D%5D)%3B%0A%20%20%20%20ctx.fillStyle%3D%22%234a5568%22%3Bctx.font%3D(9*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(%22drag%20to%20orbit%20%C2%B7%20scroll%20zoom%22%2Cpx%2B5*DPR%2Cpy%2Bph-7*DPR)%3B%0A%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%2F*%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20PASS-THROUGH%20(let%20panel%20handle%20its%20own%20draw)%20%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20%20%20else%20if(typeof%20p._draw%3D%3D%3D%22function%22)%7Bp._draw(ctx%2Cpx%2Cpy%2Cpw%2Cph%2Cf%2CDPR)%3B%7D%0A%0A%20%20%20%20%20%20ctx.restore()%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Controls%20overlay%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20drawControls()%7B%0A%20%20%20%20%20%20const%20W%3Dcv.width%2CH%3Dcv.height%2Cy%3DH-14*DPR%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22rgba(15%2C17%2C23%2C0.75)%22%3Bctx.fillRect(0%2Cy-8*DPR%2CW%2C22*DPR)%3B%0A%20%20%20%20%20%20const%20btn%3Dplaying%3F%22%E2%8F%B8%22%3A%22%E2%96%B6%22%2Cspd%3D%22%C3%97%22%2Bspeed.toFixed(1)%3B%0A%20%20%20%20%20%20ctx.fillStyle%3D%22%238a93a6%22%3Bctx.font%3D(10*DPR)%2B%22px%20monospace%22%3Bctx.textAlign%3D%22left%22%3Bctx.fillText(btn%2B%22%20%20%22%2Bspd%2B%22%20%20frame%20%22%2Bframe%2B%22%2F%22%2BS.frames%2B%22%20%20scroll%3Dzoom%20%20drag%3Dpan%20%20dblclick%3Dreset%22%2C8*DPR%2Cy%2B5*DPR)%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Main%20draw%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20draw()%7B%0A%20%20%20%20%20%20ctx.fillStyle%3DBG%3Bctx.fillRect(0%2C0%2Ccv.width%2Ccv.height)%3B%0A%20%20%20%20%20%20for(let%20i%3D0%3Bi%3CS.panels.length%3Bi%2B%2B)%20try%7BdrawPanel(S.panels%5Bi%5D%2Ci%2C...S.layout%5Bi%5D)%3B%7Dcatch(e)%7B%7D%0A%20%20%20%20%20%20drawControls()%3B%0A%20%20%20%20%7D%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Controls%20(click%20to%20pause%2C%20double-click%20to%20reset%2C%20key%20S%20for%20speed)%20%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20cv.addEventListener(%22click%22%2Ce%3D%3E%7B%0A%20%20%20%20%20%20const%20r%3Dcv.getBoundingClientRect()%2Cy%3De.clientY-r.top%3B%0A%20%20%20%20%20%20if(y%3Ecv.clientHeight-20)%7Bplaying%3D!playing%3B%7D%0A%20%20%20%20%7D)%3B%0A%20%20%20%20addEventListener(%22keydown%22%2Ce%3D%3E%7B%0A%20%20%20%20%20%20if(e.key%3D%3D%3D%22s%22%7C%7Ce.key%3D%3D%3D%22S%22)%7Bspeed%3Dspeed%3C4%3Fspeed%2B0.5%3A0.5%3B%7D%0A%20%20%20%20%20%20if(e.key%3D%3D%3D%22%20%22)%7Bplaying%3D!playing%3B%7D%0A%20%20%20%20%20%20if(e.key%3D%3D%3D%22r%22%7C%7Ce.key%3D%3D%3D%22R%22)%7Bframe%3D0%3BVP.forEach(v%3D%3E%7Bv.zoom%3D1%3Bv.px%3D0%3Bv.py%3D0%3B%7D)%3B%7D%0A%20%20%20%20%7D)%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%203D%20orbit%20events%20(only%20apply%20when%203D%20panel%20hovered)%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20cv.addEventListener(%22mousedown%22%2Ce%3D%3E%7Bconst%20r%3Dcv.getBoundingClientRect()%3Blx3%3De.clientX-r.left%3Bly3%3De.clientY-r.top%3Bdrag3%3Dtrue%3B%7D)%3B%0A%20%20%20%20cv.addEventListener(%22mousemove%22%2Ce%3D%3E%7Bif(!drag3)return%3Bconst%20r%3Dcv.getBoundingClientRect()%3Baz3%2B%3D(e.clientX-r.left-lx3)*0.01%3Bel3-%3D(e.clientY-r.top-ly3)*0.01%3Bel3%3DMath.max(-1.45%2CMath.min(1.45%2Cel3))%3Blx3%3De.clientX-r.left%3Bly3%3De.clientY-r.top%3B%7D)%3B%0A%20%20%20%20addEventListener(%22mouseup%22%2C()%3D%3E%7Bdrag3%3Dfalse%3Bdrag%3Dfalse%3B%7D)%3B%0A%20%20%20%20cv.addEventListener(%22wheel%22%2Ce%3D%3E%7Bzoom3%3DMath.max(0.3%2CMath.min(3%2Czoom3*Math.exp(-e.deltaY*0.001)))%3B%7D%2C%7Bpassive%3Atrue%7D)%3B%0A%0A%20%20%20%20%2F*%20%E2%94%80%E2%94%80%20Animation%20loop%20%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%E2%94%80%20*%2F%0A%20%20%20%20function%20loop(t)%7Bif(playing%26%26t-last%3ES.dt%2Fspeed)%7Bframe%3D(frame%2B1)%25S.frames%3Blast%3Dt%3B%7Dtry%7Bdraw()%3B%7Dcatch(e)%7B%7DrequestAnimationFrame(loop)%3B%7D%0A%20%20%20%20function%20_kick()%7Bfit()%3Bif(cv.width%3E0%26%26cv.height%3E0)%7BrequestAnimationFrame(loop)%3B%7Delse%7BsetTimeout(_kick%2C80)%3B%7D%7D_kick()%3B%0A%0A%20%20%20%20%3C%2Fscript%3E%3C%2Fbody%3E%3C%2Fhtml%3E%22%22%22.replace(%22__SPEC__%22%2C%20json.dumps(_d))%0A%20%20%20%20_st%20%3D%20cd_stats%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22**Recovered**%20helix%20%7B_st%5B'rec'%5D%5B0%5D%3A.0%25%7D%2C%20sheet%20%7B_st%5B'rec'%5D%5B1%5D%3A.0%25%7D%2C%20coil%20%7B_st%5B'rec'%5D%5B2%5D%3A.0%25%7D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22(true%20%7B_st%5B'true'%5D%5B0%5D%3A.0%25%7D%20%2F%20%7B_st%5B'true'%5D%5B1%5D%3A.0%25%7D%20%2F%20%7B_st%5B'true'%5D%5B2%5D%3A.0%25%7D)%20%26nbsp%3B%7C%26nbsp%3B%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22spectral%20R%C2%B2%20%7B_st%5B'r2'%5D%3A.3f%7D%20%26nbsp%3B%7C%26nbsp%3B%20final%20loss%20%7B_st%5B'final_loss'%5D%3A.3f%7D%22%0A%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20mo.iframe(_JS%2C%20height%3D%22470px%22)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20The%20softmax%20output%20guarantees%20the%20prediction%20is%20a%20valid%20composition%2C%20so%20the%20network%20can%20never%0A%20%20%20%20return%20negative%20or%20unnormalised%20fractions.%20The%20physics%20term%20measures%20how%20well%20the%20predicted%0A%20%20%20%20composition%20re-explains%20the%20measured%20spectrum%20through%20the%20fixed%20forward%20model%2C%20which%20is%20what%0A%20%20%20%20keeps%20the%20inversion%20stable%20when%20the%20basis%20curves%20overlap%20and%20noise%20is%20high.%20As%20training%20runs%2C%0A%20%20%20%20the%20reconstructed%20spectrum%20(amber)%20settles%20onto%20the%20measurement%20and%20the%20predicted%20bars%20climb%0A%20%20%20%20toward%20the%20true%20ones.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A