nx_geo_query_gate.nx source

1// nx_geo_query_gate.nx -- GATE for GEO-009 spatial query. Proves: 2// nearest : among {LA, Boston, Philadelphia}, nearest to NYC is Philadelphia (index 2) 3// within : within 200 km of NYC = {Philadelphia} only (Boston ~306 km out, LA out) 4// bearing : due N=0, due E=90, due S=180, due W=270, NE@equator=45 (the compass quadrants) 5// 6// Evidence -> knowledge/status/geo_query.log (GEOQUERYGATE authored=organ ... verdict=GREEN). 7// license_tier: ORIGINAL 8import "nx_geo_query.nx" 9import "nx_syscalls.nx" 10 11const GQ_LOG: *u8 = "knowledge/status/geo_query.log" 12 13func gq_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 } 14func gq_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 } 15 16func gq_near(a: i64, b: i64, tol: i64) -> i64 { var d: i64 = a - b; if d < 0 { d = 0 - d } if d <= tol { return 1 } return 0 } 17 18func gq_set(p: *i64, i: i64, lat: i64, lon: i64) -> i64 { p[i * 2] = lat; p[i * 2 + 1] = lon; return 0 } 19 20func gq_emit(fd: i64, near: i64, wcount: i64, widx: i64, bN: i64, bE: i64, bS: i64, bW: i64, bNE: i64, ok: i64) -> i64 { 21 gq_w(fd, "GEOQUERYGATE authored=organ nearest_idx=" as *u8); gq_wn(fd, near) 22 gq_w(fd, " within200km_count=" as *u8); gq_wn(fd, wcount) 23 gq_w(fd, " within_idx0=" as *u8); gq_wn(fd, widx) 24 gq_w(fd, " bearing_N=" as *u8); gq_wn(fd, bN) 25 gq_w(fd, " bearing_E=" as *u8); gq_wn(fd, bE) 26 gq_w(fd, " bearing_S=" as *u8); gq_wn(fd, bS) 27 gq_w(fd, " bearing_W=" as *u8); gq_wn(fd, bW) 28 gq_w(fd, " bearing_NE=" as *u8); gq_wn(fd, bNE) 29 if ok == 1 { gq_w(fd, " verdict=GREEN\n" as *u8) } else { gq_w(fd, " verdict=RED\n" as *u8) } 30 return 0 31} 32 33func main() -> i64 { 34 // candidates: 0=LA, 1=Boston, 2=Philadelphia (microdegrees). 35 let pts: *i64 = sys_mmap(8 * 6) as *i64 36 gq_set(pts, 0, 34052200, 0 - 118243700) // LA 37 gq_set(pts, 1, 42360100, 0 - 71058900) // Boston 38 gq_set(pts, 2, 39952600, 0 - 75165200) // Philadelphia 39 let qlat: i64 = 40712800 // NYC 40 let qlon: i64 = 0 - 74006000 41 42 let near: i64 = geo_nearest(qlat, qlon, pts, 3) // expect 2 (Philadelphia) 43 let widx: *i64 = sys_mmap(8 * 3) as *i64 44 let wcount: i64 = geo_within(qlat, qlon, 200000, pts, 3, widx) // expect 1 (Philadelphia) 45 var w0: i64 = 0 - 1 46 if wcount >= 1 { w0 = widx[0] } 47 48 // bearings from origin (equator) for clean quadrants. 49 let bN: i64 = geo_bearing(0, 0, 1000000, 0) // due north -> 0 50 let bE: i64 = geo_bearing(0, 0, 0, 1000000) // due east -> 90 51 let bS: i64 = geo_bearing(0, 0, 0 - 1000000, 0) // due south -> 180 52 let bW: i64 = geo_bearing(0, 0, 0, 0 - 1000000) // due west -> 270 53 let bNE: i64 = geo_bearing(0, 0, 1000000, 1000000) // NE @ equator -> 45 54 55 var ok: i64 = 1 56 if near != 2 { ok = 0 } 57 if wcount != 1 { ok = 0 } 58 if w0 != 2 { ok = 0 } 59 if gq_near(bN, 0, 2) != 1 { ok = 0 } 60 if gq_near(bE, 90, 2) != 1 { ok = 0 } 61 if gq_near(bS, 180, 2) != 1 { ok = 0 } 62 if gq_near(bW, 270, 2) != 1 { ok = 0 } 63 if gq_near(bNE, 45, 2) != 1 { ok = 0 } 64 65 gq_emit(1, near, wcount, w0, bN, bE, bS, bW, bNE, ok) 66 let lf: i64 = sys_openat_append(GQ_LOG, 420) 67 if lf >= 0 { gq_emit(lf, near, wcount, w0, bN, bE, bS, bW, bNE, ok); sys_close(lf) } 68 69 if ok == 1 { return 0 } 70 return 1 71}