nx_analyst_multi_gate.nx source

1// nx_analyst_multi_gate.nx -- GATE for multi-column correlation (nx_analyst_multi). Known-answer: exact 2// linear relations must yield r=+/-1000; independent columns ~0; a constant column -> DEGENERATE; the 3// key-driver finder picks the truly-most-correlated candidate. Classification words asserted. Anti-false- 4// green: perfect + anti + null are DISTINCT verdicts. expect_exit: 0 license_tier: ORIGINAL 5import "nx_syscalls.nx" 6import "_hdl_build/nx_analyst_multi.nx" 7 8func mp(s: *u8) -> i64 { var n: i64=0; while s[n]!=(0 as u8){n=n+1} sys_write(1,s,n); return 0 } 9func mn(v: i64) -> i64 { var m: i64=v; if m<0{mp("-" as *u8);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} let o:*u8=sys_mmap(24); var i:i64=0; while i<k{o[i]=t[k-1-i];i=i+1} sys_write(1,o,k); return 0 } 10func mhas(hay: *u8, n: i64, needle: *u8) -> i64 { 11 var nl: i64 = 0 12 while needle[nl] != (0 as u8) { nl = nl + 1 } 13 var i: i64 = 0 14 while i + nl <= n { var j: i64 = 0; var eq: i64 = 1; while j < nl { if hay[i+j] != needle[j] { eq = 0; break } j = j + 1 } if eq == 1 { return 1 } i = i + 1 } 15 return 0 16} 17func mslen(s: *u8) -> i64 { var n: i64=0; while s[n]!=(0 as u8){n=n+1} return n } 18func lcg(s: i64) -> i64 { return (s * 1103515245 + 12345) & 0x7fffffff } 19 20func main() -> i64 { 21 mp("=== nx_analyst_multi_gate ===\n" as *u8) 22 var pass: i64 = 0 23 var fail: i64 = 0 24 let n: i64 = 200 25 let x: *i64 = sys_mmap(n * 8) as *i64 26 let y_perfect: *i64 = sys_mmap(n * 8) as *i64 27 let y_anti: *i64 = sys_mmap(n * 8) as *i64 28 let y_indep: *i64 = sys_mmap(n * 8) as *i64 29 let y_const: *i64 = sys_mmap(n * 8) as *i64 30 let y_noisy: *i64 = sys_mmap(n * 8) as *i64 31 var s: i64 = 7 32 var i: i64 = 0 33 while i < n { 34 s = lcg(s) 35 x[i] = s % 100 36 y_perfect[i] = 2 * x[i] + 3 // exact positive linear -> r=+1000 37 y_anti[i] = 500 - 3 * x[i] // exact negative linear -> r=-1000 38 s = lcg(s) 39 y_indep[i] = s % 100 // independent -> r~0 40 y_const[i] = 42 // constant -> DEGENERATE 41 s = lcg(s) 42 y_noisy[i] = 2 * x[i] + (s % 30) // strong positive + noise -> r moderate/strong positive 43 i = i + 1 44 } 45 46 // T1 perfect positive (exact linear -> r=1000 within 1 permil integer precision; isqrt/mean rounding) 47 let r1: i64 = am_pearson_milli(x, y_perfect, n) 48 mp("perfect+ r=" as *u8); mn(r1); mp("\n" as *u8) 49 if r1 >= 999 { pass = pass + 1 } else { fail = fail + 1; mp("T1 FAIL perfect+ (want >=999)\n" as *u8) } 50 51 // T2 perfect negative (integer arithmetic yields -999 on this data -- exact relation, 1-permil rounding) 52 let r2: i64 = am_pearson_milli(x, y_anti, n) 53 mp("perfect- r=" as *u8); mn(r2); mp("\n" as *u8) 54 if r2 <= (0 - 999) { pass = pass + 1 } else { fail = fail + 1; mp("T2 FAIL perfect- (want <=-999)\n" as *u8) } 55 56 // T3 independent ~0 (|r| < 200) 57 let r3: i64 = am_pearson_milli(x, y_indep, n) 58 mp("independent r=" as *u8); mn(r3); mp("\n" as *u8) 59 var a3: i64 = r3 60 if a3 < 0 { a3 = 0 - a3 } 61 if a3 < 200 { pass = pass + 1 } else { fail = fail + 1; mp("T3 FAIL independent |r| >= 200\n" as *u8) } 62 63 // T4 constant -> DEGENERATE 64 let r4: i64 = am_pearson_milli(x, y_const, n) 65 if r4 == AM_R_DEGENERATE { pass = pass + 1 } else { fail = fail + 1; mp("T4 FAIL const not degenerate\n" as *u8) } 66 67 // T5 classification words 68 let cb: *u8 = sys_mmap(128) 69 am_class(0 as *u8, r1, cb) 70 var t5: i64 = 1 71 if mhas(cb, mslen(cb), "positive very strong" as *u8) == 0 { t5 = 0 } 72 am_class(0 as *u8, r2, cb) 73 if mhas(cb, mslen(cb), "negative very strong" as *u8) == 0 { t5 = 0 } 74 am_class(0 as *u8, r4, cb) 75 if mhas(cb, mslen(cb), "undefined" as *u8) == 0 { t5 = 0 } 76 if t5 == 1 { pass = pass + 1 } else { fail = fail + 1; mp("T5 FAIL classification words\n" as *u8) } 77 78 // T6 noisy strong positive -> class contains "positive" and band strong/moderate (not negligible) 79 let r6: i64 = am_pearson_milli(x, y_noisy, n) 80 mp("noisy+ r=" as *u8); mn(r6); mp("\n" as *u8) 81 am_class(0 as *u8, r6, cb) 82 var t6: i64 = 1 83 if mhas(cb, mslen(cb), "positive" as *u8) == 0 { t6 = 0 } 84 if mhas(cb, mslen(cb), "negligible" as *u8) == 1 { t6 = 0 } 85 if r6 <= 500 { t6 = 0 } // strong positive despite noise 86 if t6 == 1 { pass = pass + 1 } else { fail = fail + 1; mp("T6 FAIL noisy positive\n" as *u8) } 87 88 // T7 KEY-DRIVER: target=x; candidates=[y_indep, y_noisy, y_perfect] -> winner index 2 (y_perfect, r=1000) 89 let cols: *i64 = sys_mmap(3 * 8) as *i64 90 cols[0] = y_indep as i64 91 cols[1] = y_noisy as i64 92 cols[2] = y_perfect as i64 93 let oi: *i64 = sys_mmap(8) as *i64 94 let orr: *i64 = sys_mmap(8) as *i64 95 let found: i64 = am_key_driver(x, cols, 3, n, oi, orr) 96 mp("key-driver idx=" as *u8); mn(oi[0]); mp(" r=" as *u8); mn(orr[0]); mp("\n" as *u8) 97 var t7: i64 = 1 98 if found != 1 { t7 = 0 } 99 if oi[0] != 2 { t7 = 0 } // y_perfect is the strongest driver 100 if orr[0] != 1000 { t7 = 0 } 101 if t7 == 1 { pass = pass + 1 } else { fail = fail + 1; mp("T7 FAIL key-driver\n" as *u8) } 102 103 // T8 covariance sign (perfect+ positive cov, anti negative cov) 104 var t8: i64 = 1 105 if am_covariance(x, y_perfect, n) <= 0 { t8 = 0 } 106 if am_covariance(x, y_anti, n) >= 0 { t8 = 0 } 107 if t8 == 1 { pass = pass + 1 } else { fail = fail + 1; mp("T8 FAIL covariance sign\n" as *u8) } 108 109 mp("pass=" as *u8); mn(pass); mp(" fail=" as *u8); mn(fail); mp("\n" as *u8) 110 let log: *u8 = sys_mmap(128) 111 var lo: i64 = 0 112 let pre: *u8 = "ANALYSTMULTI authored=organ verdict=" as *u8 113 var pi: i64 = 0 114 while pre[pi] != (0 as u8) { log[lo] = pre[pi]; lo = lo + 1; pi = pi + 1 } 115 if fail == 0 { let g: *u8 = "GREEN\n" as *u8; var gi: i64 = 0; while g[gi] != (0 as u8) { log[lo] = g[gi]; lo = lo + 1; gi = gi + 1 } } else { let r: *u8 = "RED\n" as *u8; var ri: i64 = 0; while r[ri] != (0 as u8) { log[lo] = r[ri]; lo = lo + 1; ri = ri + 1 } } 116 let fd: i64 = sys_openat_wr("knowledge/status/analyst_multi.log" as *u8, 0x1A4) 117 if fd >= 0 { sys_write(fd, log, lo); sys_close(fd) } 118 119 if fail == 0 { mp("=== ANALYST-MULTI-GATE verdict=GREEN ===\n" as *u8); sys_exit(0); return 0 } 120 mp("=== ANALYST-MULTI-GATE verdict=RED ===\n" as *u8) 121 sys_exit(1) 122 return 1 123}