code wiki / _hdl_build / nx_research_dynamics_gate.nx

nx_research_dynamics_gate.nx

buildroot/runtime/_hdl_build/nx_research_dynamics_gate.nx

5604 B94 linesdepth 2pulls 2 transitivereach 0 importersview sourcekind gate/prooftopic research
docsdependenciesstructsconstsfunctions

about

nx_research_dynamics_gate.nx -- TIME-STEPPING dynamics simulation, sovereign + fixed-point (no float), with a conservation-law verifier (operator NEXT: "time-stepping dynamics sims"). Simulates a harmonic oscillator (a = -x, unit mass/stiffness) two ways and lets the verifier tell them apart by ENERGY conservation: LEAPFROG (symplectic): v_half=v - x/(2D); x'=x + v_half/D; v'=v_half - x'/(2D) -- energy stays BOUNDED. FORWARD-EULER (naive): x'=x + v/D; v'=v - x/D -- energy DRIFTS upward. The well-known physics: a symplectic integrator conserves energy to bounded error over arbitrarily many steps; Euler injects energy and the orbit spirals out. The verifier measures max energy deviation; LIAR-KILL: Euler produces a perfectly plausible trajectory but FAILS the conservation bar leapfrog passes -- a drifting integrator cannot be rubber-stamped as energy-conserving. Fixed-point integers, deterministic. GREEN iff 6/6. license_tier: ORIGINAL

dependencies 1 imports · 0 importers

nx_syscalls.nx nx_research_dynamics_gate.nx

imports: nx_syscalls.nx

imported by: nobody (leaf or entry point)

call flow from main pre-order; caps 40 nodes / depth 6 declared; ↻ = already shown

main sys_mmap g_w sys_write sim_leapfrog iabs sim_euler iabs ↻ permille g_n g_w ↻ sys_mmap ↻ sys_write ↻ g_row g_w ↻ sys_exit

structs

none

consts

none

functions

13func g_w(s: *u8) -> i64 { var n: i64=0; while s[n]!=(0 as u8){n=n+1} sys_write(1,s,n); return 0 }
called by 3: g_ng_rowmain calls 1: sys_write
14func g_n(v: i64) -> i64 { var m: i64=v; if m<0{g_w("-");m=0-m} let t:*u8=sys_mmap(24); var k:i64=0; if m==0{t[0]=48 as u8;k=1}; while m>0{t[k]=(48+(m%10)) as u8;m=m/10;k=k+1}; var i:i64=0; let o:*u8=sys_mmap(24); while i<k{o[i]=t[k-1-i];i=i+1}; sys_write(1,o,k); return 0 }
called by 1: main calls 3: g_wsys_mmapsys_write
15func g_row(id: *u8, ok: i64, pass: *i64) -> i64 { g_w(" "); g_w(id); g_w(": "); if ok==1 { g_w("OK\n"); pass[0]=pass[0]+1 } else { g_w("FAIL\n") } return 0 }
called by 1: main calls 1: g_w
16func iabs(x: i64) -> i64 { if x<0 { return 0-x } return x }
19func sim_leapfrog(x0: i64, v0: i64, D: i64, N: i64, out: *i64) -> i64
called by 1: main calls 1: iabs
38func sim_euler(x0: i64, v0: i64, D: i64, N: i64, out: *i64) -> i64
called by 1: main calls 1: iabs
56func permille(dev: i64, E0: i64) -> i64 { return dev*1000/E0 } // dev*1000 < 2^63 for our magnitudes
called by 1: main
58func main() -> i64