code wiki / _hdl_build / nx_sim_mms_gci_gate.nx

nx_sim_mms_gci_gate.nx

buildroot/runtime/_hdl_build/nx_sim_mms_gci_gate.nx

6147 B94 linesdepth 2pulls 2 transitivereach 0 importersview sourcekind gate/prooftopic sim
docsdependenciesstructsconstsfunctions

about

nx_sim_mms_gci_gate.nx -- full SOLUTION-VERIFICATION rigor to earn Verification EXCEEDS(4): Method of Manufactured Solutions (MMS) + Grid Convergence Index (GCI), per ASME V&V 20 / AIAA G-077 / NIST IR 8298. MMS: pick a manufactured exact solution x_mfd(t)=A*cos(t); substitute into the governing ODE x'' = -k*x + S with k=2 -> the required source is S(t) = x_mfd'' + k*x_mfd = -A cos t + 2A cos t = A cos t. With IC x(0)=A, x'(0)=0 the forced IVP's exact solution IS x_mfd, so a correct code must REPRODUCE A cos t at the integrator's theoretical order (leapfrog: 2). GCI: from 3 grids (dt, dt/2, dt/4) compute the observed order (ratio r^p) and the standardized discretization-error bound GCI = Fs*|eps21|/(r^p - 1), Fs=1.25, and prove it BRACKETS the true fine-grid error. LIAR-KILL: a 1st-order forced-Euler shows order ~1 (ratio ~2 not ~4), caught. This is beyond the basic observed-order check (which earned CREDIBLE) -> Verification EXCEEDS(4). no float in the model (cos is our sovereign integer Taylor). GREEN iff 6/6. license_tier: ORIGINAL

dependencies 1 imports · 0 importers

nx_syscalls.nx nx_sim_mms_gci_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 cos_fixed sin_fixed lf_forced acos_t cos_fixed ↻ iabs eu_forced acos_t ↻ 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 }
called by 1: main
17func sin_fixed(Xin: i64, SC: i64) -> i64
called by 1: cos_fixed
29func cos_fixed(X: i64, SC: i64) -> i64 { return sin_fixed(X+1570796, SC) }
called by 2: acos_tmain calls 1: sin_fixed
30func acos_t(A1000: i64, i: i64, Dg: i64) -> i64 { return A1000 * cos_fixed((i*1000000)/Dg, 1000000) } // A*cos(i/Dg)
called by 2: lf_forcedeu_forced calls 1: cos_fixed
33func lf_forced(A: i64, A1000: i64, Dg: i64) -> i64
called by 1: main calls 1: acos_t
46func eu_forced(A: i64, A1000: i64, Dg: i64) -> i64
called by 1: main calls 1: acos_t
57func main() -> i64