nx_geo_buffer_gate.nx source

1// nx_geo_buffer_gate.nx -- GATE for GEO-017 buffer/offset measure. Polygon = square side 4e6 2// (geo_area2 = 3.2e13, perimeter = 16e6). Proves the convex Minkowski-disk buffer formula: 3// r=0 -> buffer_area2 == area2(P) == 3.2e13 (no buffer) 4// perim == 16e6 (sum of 4 edges of 4e6) 5// r=1e6 -> buffer_area2 == 3.2e13 + 2*16e6*1e6 + 710*1e12/113 6// == 32000000000000 + 32000000000000 + 6283185840707 == 70283185840707 (integer + 355/113 pi) 7// growth -> buffer(r=1e6) > buffer(r=0) 8// 9// Evidence -> knowledge/status/geo_buffer.log (GEOBUFGATE authored=organ ... verdict=GREEN). 10// license_tier: ORIGINAL 11import "nx_geo_buffer.nx" 12import "nx_syscalls.nx" 13 14const GB_LOG: *u8 = "knowledge/status/geo_buffer.log" 15 16func gb_w(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 } 17func gb_wn(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;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 } 18 19func gb_set(p: *i64, i: i64, lat: i64, lon: i64) -> i64 { p[i * 2] = lat; p[i * 2 + 1] = lon; return 0 } 20func gb_square(p: *i64, lo: i64, hi: i64) -> i64 { gb_set(p, 0, lo, lo); gb_set(p, 1, lo, hi); gb_set(p, 2, hi, hi); gb_set(p, 3, hi, lo); return 0 } 21 22func gb_emit(fd: i64, per: i64, b0: i64, br: i64, ok: i64) -> i64 { 23 gb_w(fd, "GEOBUFGATE authored=organ method=convex-minkowski-disk perimeter=" as *u8); gb_wn(fd, per) 24 gb_w(fd, " buffer0_area2=" as *u8); gb_wn(fd, b0) 25 gb_w(fd, " buffer_r1e6_area2=" as *u8); gb_wn(fd, br) 26 if ok == 1 { gb_w(fd, " verdict=GREEN\n" as *u8) } else { gb_w(fd, " verdict=RED\n" as *u8) } 27 return 0 28} 29 30func main() -> i64 { 31 let p: *i64 = sys_mmap(8 * 8) as *i64 32 gb_square(p, 0, 4000000) 33 34 let per: i64 = geo_perimeter(p, 4) 35 let b0: i64 = geo_buffer_area2(p, 4, 0) 36 let br: i64 = geo_buffer_area2(p, 4, 1000000) 37 38 var ok: i64 = 1 39 if per != 16000000 { ok = 0 } 40 if b0 != 32000000000000 { ok = 0 } 41 if br != 70283185840707 { ok = 0 } 42 if br <= b0 { ok = 0 } 43 44 gb_emit(1, per, b0, br, ok) 45 let lf: i64 = sys_openat_append(GB_LOG, 420) 46 if lf >= 0 { gb_emit(lf, per, b0, br, ok); sys_close(lf) } 47 48 if ok == 1 { return 0 } 49 return 1 50}