code wiki / _hdl_build / _f64_gamma_refl_gate.nx

_f64_gamma_refl_gate.nx

buildroot/runtime/_hdl_build/_f64_gamma_refl_gate.nx

2823 B55 linesdepth 8pulls 17 transitivereach 0 importersview sourcekind gate/prooftopic f64
docsdependenciesstructsconstsfunctions

about

_f64_gamma_refl_gate.nx -- ULP gate for the gamma reflection kernel (0<x<0.5) vs the full-domain sovereign oracle bf_gamma_f64 (which does Stirling+shift+REFLECTION at 120 bits, reflection-product proven 40/40). x = k/64, k=1..31 (0.0156..0.484). The reflection composes Gamma(1-x) (<=7 ulp rung) so the honest bar is <=16 ulp (matches the gamma rung-1 relaxed ship); the MEASURED max is logged either way. Durable: GAMMAREFL-GATE -> knowledge/status/math_engine.log. license_tier: ORIGINAL

dependencies 5 imports · 0 importers

nx_syscalls.nx nx_bigfloat120.nx nx_bigfloat120_div.nx nx_bigfloat120_gamma.nx _pe_f64gamma_refl.nx _f64_gamma_refl_gate.nx

imports: nx_syscalls.nxnx_bigfloat120.nxnx_bigfloat120_div.nxnx_bigfloat120_gamma.nx_pe_f64gamma_refl.nx

imported by: nobody (leaf or entry point)

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

main gr_p sys_write bf_new sys_mmap bf_set_int bf_div_small bf_is_zero bf_copy bf_norm bf_to_f64 bf_is_zero ↻ bf_gamma_f64 bf_new ↻ bf_set_f64 bfg_gamma_pos bf_new ↻ bf_set_int ↻ bf_copy ↻ bf_cmp bf_is_zero ↻ bf_mul bf_is_zero ↻ bf_copy ↻ sys_mmap ↻ bf_add bf_is_zero ↻ bf_copy ↻ sys_mmap ↻ bf_shr2 bf_norm ↻ bfg_lngamma_big bf_new ↻ bf_set_int ↻ bf_copy ↻ bf_sub bf_is_zero ↻ bf_copy ↻ sys_mmap ↻ bf_shr2 ↻

structs

none

consts

none

functions

14func gr_p(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 1: main calls 1: sys_write
15func gr_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 }
called by 1: main calls 1: sys_write
16func gr_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 }
called by 1: main calls 2: sys_mmapsys_write
18func gr_ulp(a: i64, b: i64) -> i64 { if a >= b { return a - b } return b - a }
called by 1: main
20func main() -> i64
calls 14: gr_pbf_newbf_set_intbf_div_smallbf_to_f64bf_gamma_f64+8