code wiki / _hdl_build / _f64_gamma_refl_neg_gate.nx

_f64_gamma_refl_neg_gate.nx

buildroot/runtime/_hdl_build/_f64_gamma_refl_neg_gate.nx

3064 B56 linesdepth 8pulls 17 transitivereach 0 importersview sourcekind gate/prooftopic f64
docsdependenciesstructsconstsfunctions

about

_f64_gamma_refl_neg_gate.nx -- proves the SAME reflection kernel (_pe_f64gamma_refl) is correct for NEGATIVE x, so gamma is full-domain across the reals minus the non-positive-integer poles. Grid x = -(2j+1)/4 (quarter-points -0.25..-9.25, far from the integer poles where sin(pi*x)=0). Oracle = the full-domain bf_gamma_f64 (Stirling +shift+reflection, 120-bit). Bar <= 16 ulp (composes Gamma(1-x) <=7ulp + sin + div). Durable: GAMMAREFLNEG-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_neg_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 gn_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 gn_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 gn_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 gn_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
17func gn_ulp(a: i64, b: i64) -> i64 { if a >= b { return a - b } return b - a }
called by 1: main
19func main() -> i64
calls 14: gn_pbf_newbf_set_intbf_div_smallbf_to_f64bf_gamma_f64+8