_f64_erfc_gate.nx source

1// _f64_erfc_gate.nx -- ULP gate for the f64 erfc kernel vs the sovereign 120-bit erfc 2// oracle (ME2-ERF-002). For x = k/16, k=28..96 (1.75 .. 6.0) compares nx_f64_erfc bits 3// against the bf_erfc oracle rounded to f64. RUNG-1 BAR: max ULP <= 4 (matches erf). 4// Durable: F64ERFC-GATE + ERFCGATE rows -> knowledge/status/math_engine.log. 5// license_tier: ORIGINAL 6import "nx_syscalls.nx" 7import "nx_bigfloat120.nx" 8import "nx_bigfloat120_div.nx" 9import "nx_bigfloat120_trig.nx" 10import "nx_bigfloat120_erfc.nx" 11import "_pe_f64erfc.nx" 12 13func fg_p(s: *u8) -> i64 { var n: i64=0; while s[n]!=(0 as u8){n=n+1} sys_write(1,s,n); return 0 } 14func fg_f(fd: i64, s: *u8) -> i64 { var n: i64=0; while s[n]!=(0 as u8){n=n+1} sys_write(fd,s,n); return 0 } 15func fg_n(fd: i64, v: i64) -> i64 { let bb: *u8=sys_mmap(28); var m: i64=v; if m<0{m=0-m; sys_write(fd,"-" as *u8,1)}; let t: *u8=sys_mmap(28); 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; while i<k{bb[i]=t[k-1-i];i=i+1}; sys_write(fd,bb,k); return 0 } 16 17func fg_ulp(a: i64, b: i64) -> i64 { if a >= b { return a - b } return b - a } 18 19func main() -> i64 { 20 fg_p("=== F64-ERFC GATE: f64 kernel vs sovereign 120-bit erfc oracle (1.75..6.0) ===\n" as *u8) 21 let kb: *i64 = bf_new() 22 let xb: *i64 = bf_new() 23 let oc: *i64 = bf_new() 24 var maxulp: i64 = 0 25 var worst: i64 = 0 26 var npts: i64 = 0 27 var k: i64 = 28 28 while k <= 96 { 29 bf_set_int(kb, k) 30 bf_div_small(xb, kb, 16) 31 let xbits: i64 = bf_to_f64(xb, 0, 0) 32 bf_erfc(oc, xb) 33 let want: i64 = bf_to_f64(oc, 0, 0) 34 let got: i64 = nx_f64_erfc(xbits) 35 let u: i64 = fg_ulp(got, want) 36 if u > maxulp { maxulp = u; worst = xbits } 37 npts = npts + 1 38 k = k + 1 39 } 40 let lfd: i64 = sys_openat_append("knowledge/status/math_engine.log" as *u8, 0x1a4) 41 if lfd >= 0 { 42 fg_f(lfd, "F64ERFC-GATE epoch=" as *u8); fg_n(lfd, sys_now_realtime_sec()) 43 fg_f(lfd, " points=" as *u8); fg_n(lfd, npts) 44 fg_f(lfd, " max_ulp=" as *u8); fg_n(lfd, maxulp) 45 fg_f(lfd, " worst_x_bits=" as *u8); fg_n(lfd, worst) 46 fg_f(lfd, " domain=1.75..6.0 oracle=nx_bigfloat120_erfc(SOVEREIGN,erf+erfc=1-proven)" as *u8) 47 if maxulp <= 4 { fg_f(lfd, " verdict=GREEN\n" as *u8) } else { fg_f(lfd, " verdict=RED\n" as *u8) } 48 // canonical ch7-completion anchor row 49 fg_f(lfd, "ERFCGATE epoch=" as *u8); fg_n(lfd, sys_now_realtime_sec()) 50 fg_f(lfd, " max_ulp=" as *u8); fg_n(lfd, maxulp) 51 if maxulp <= 4 { fg_f(lfd, " verdict=GREEN\n" as *u8) } else { fg_f(lfd, " verdict=RED\n" as *u8) } 52 sys_close(lfd) 53 } 54 fg_p(" points=" as *u8); fg_n(1, npts) 55 fg_p(" max_ulp=" as *u8); fg_n(1, maxulp) 56 if maxulp <= 4 { fg_p(" F64ERFC-GATE: GREEN (rung-1 bar ulp<=4 held vs sovereign oracle)\n" as *u8); sys_exit(0); return 0 } 57 fg_p(" F64ERFC-GATE: RED (bar ulp<=4 broken -- measured gap named, kernel needs refinement)\n" as *u8) 58 sys_exit(1) 59 return 1 60}