nx_step_brepw_gate.nx source

1// nx_step_brepw_gate.nx -- ANALYTIC B-REP WRITE, both capabilities verified through our own gated pipeline: 2// (A) AUTHOR: sw2_box writes a 100x60x40mm analytic B-rep (PLANEs + full edge topology) from parameters -> 3// re-parse -> v2 tessellation must give 6/6 faces, 12 tris, EXACT bbox, exact invariants, mm units. 4// (B) PERSIST: sw2_reemit re-serializes the ENTIRE as1 reachable graph (renumbered) -> re-parse -> full 5// planar+NURBS pipeline must EQUAL the original model (25 planar + 28 round, 1648 tris, order-invariant 6// geometry checksum). expect_exit: 0 license_tier: ORIGINAL 7import "nx_step_brepw.nx" 8import "nx_step_tess2.nx" 9 10func bg2_tooth(name: *u8, pass: i64, fails: *i64) -> i64 { 11 sp_puts(" " as *u8); sp_puts(name); sp_puts(" -> " as *u8) 12 if pass == 1 { sp_puts("PASS\n" as *u8); return 0 } 13 sp_puts("FAIL\n" as *u8) 14 fails[0] = fails[0] + 1 15 return 0 16} 17// build tess ctxs bound to st2 + a fresh mesh; W: [0]=tt [1]=t2 [2]=nt [3]=slab [4]=c [5]=c2 [6]=mesh 18func bg2_ctx(st2: *i64, W: *i64) -> i64 { 19 let tt: *i64 = sys_mmap(128) as *i64 20 let c: *i64 = sys_mmap(128) as *i64 21 let t2: *i64 = sys_mmap(128) as *i64 22 let nt: *i64 = sys_mmap(128) as *i64 23 let c2: *i64 = sys_mmap(128) as *i64 24 let mesh: i64 = sys_mmap(m3_bytes()) as i64 25 m3_init(mesh) 26 var z: i64 = 0 27 while z < 16 { c[z] = 0; c2[z] = 0; z = z + 1 } 28 tt[0] = st2 as i64 29 tt[1] = mesh 30 tt[2] = c as i64 31 tt[3] = sys_mmap(768) as i64 32 tt[4] = 0 33 tt[5] = sys_mmap(8192) as i64 34 tt[6] = sys_mmap(4096) as i64 35 tt[7] = sys_mmap(128) as i64 36 tt[8] = sys_mmap(16) as i64 37 tt[9] = 0 38 t2[0] = st2 as i64 39 t2[1] = c as i64 40 t2[2] = sys_mmap(32) as i64 41 t2[3] = sys_mmap(960) as i64 42 t2[4] = sys_mmap(1024) as i64 43 t2[5] = sys_mmap(128) as i64 44 t2[6] = sys_mmap(512) as i64 45 t2[7] = sys_mmap(256) as i64 46 t2[8] = sys_mmap(8192) as i64 47 t2[9] = sys_mmap(128) as i64 48 t2[10] = sys_mmap(768) as i64 49 t2[11] = sys_mmap(128) as i64 50 t2[12] = sys_mmap(256) as i64 51 nt[0] = st2 as i64 52 nt[1] = mesh 53 nt[2] = c2 as i64 54 nt[3] = sys_mmap(256) as i64 55 nt[4] = sys_mmap(64) as i64 56 nt[5] = sys_mmap(16) as i64 57 nt[6] = sys_mmap(128) as i64 58 nt[7] = sys_mmap(128) as i64 59 W[0] = tt as i64 60 W[1] = t2 as i64 61 W[2] = nt as i64 62 W[3] = sys_mmap(20480) as i64 63 W[4] = c as i64 64 W[5] = c2 as i64 65 W[6] = mesh 66 return 0 67} 68// order-invariant geometry checksum over the mesh triangles 69func bg2_cksum(mesh: i64) -> i64 { 70 let h: *i64 = m3_hdr(mesh) 71 var s: i64 = 0 72 var t: i64 = 0 73 while t < h[1] { 74 let tr: *i64 = m3_tri(mesh, t) 75 var pv: i64 = 0 76 var k: i64 = 0 77 while k < 3 { 78 let v: *i64 = m3_vert(mesh, tr[k]) 79 var m: i64 = 1 80 if k == 1 { m = 5 } 81 if k == 2 { m = 7 } 82 pv = pv + m * (v[0] + 2 * v[1] + 3 * v[2]) 83 k = k + 1 84 } 85 s = s + pv 86 t = t + 1 87 } 88 return s 89} 90// parse a STEP file into fresh tables; returns st ptr (0 fail) 91func bg2_parse(path: *u8, out1: *i64) -> i64 { 92 let buf: *u8 = sys_mmap(8388608) 93 let a1: *i64 = sys_mmap(524288) as *i64 94 let a2: *i64 = sys_mmap(524288) as *i64 95 let a3: *i64 = sys_mmap(524288) as *i64 96 let a4: *i64 = sys_mmap(524288) as *i64 97 let a5: *i64 = sys_mmap(524288) as *i64 98 let st2: *i64 = sys_mmap(256) as *i64 99 let n: i64 = sp_read_file(path, buf, 8388608) 100 if n <= 0 { return 0 } 101 sp_init(st2, buf, n, a1, a2, a3, a4, a5) 102 out1[0] = sp_scan(st2) 103 return st2 as i64 104} 105 106func main() -> i64 { 107 sp_puts("=== nx_step_brepw_gate -- ANALYTIC B-rep write: author + persist, pipeline-verified ===\n" as *u8) 108 let fails: *i64 = sys_mmap(8) as *i64 109 fails[0] = 0 110 let o1: *i64 = sys_mmap(16) as *i64 111 112 // ---- (A) AUTHOR a 100x60x40mm box as analytic B-rep 113 let exp: i64 = sw2_box("knowledge/box_brep.stp" as *u8, 25600, 15360, 10240, "box-100x60x40" as *u8) 114 let stB: i64 = bg2_parse("knowledge/box_brep.stp" as *u8, o1) 115 sp_puts("box: expected entities=" as *u8); sp_putn(exp); sp_puts(" parsed=" as *u8); sp_putn(o1[0]); sp_puts("\n" as *u8) 116 var t1: i64 = 0 117 if stB != 0 { if o1[0] == exp { t1 = 1 } } 118 let ig1: i64 = bg2_tooth("T1 AUTHORED box: entity count matches the writer formula (121)" as *u8, t1, fails) 119 120 let WB: *i64 = sys_mmap(128) as *i64 121 bg2_ctx(stB as *i64, WB) 122 let nb: i64 = sgt2_tessellate(WB[0] as *i64, WB[1] as *i64, WB[3] as *i64) 123 let cB: *i64 = WB[4] as *i64 124 let hB: *i64 = m3_hdr(WB[6]) 125 // bbox 126 var mn0: i64 = 0 127 var mx0: i64 = 0 128 var mn1: i64 = 0 129 var mx1: i64 = 0 130 var mn2: i64 = 0 131 var mx2: i64 = 0 132 if hB[0] > 0 { 133 let v0: *i64 = m3_vert(WB[6], 0) 134 mn0 = v0[0]; mx0 = v0[0]; mn1 = v0[1]; mx1 = v0[1]; mn2 = v0[2]; mx2 = v0[2] 135 var vi: i64 = 1 136 while vi < hB[0] { 137 let v: *i64 = m3_vert(WB[6], vi) 138 if v[0] < mn0 { mn0 = v[0] } 139 if v[0] > mx0 { mx0 = v[0] } 140 if v[1] < mn1 { mn1 = v[1] } 141 if v[1] > mx1 { mx1 = v[1] } 142 if v[2] < mn2 { mn2 = v[2] } 143 if v[2] > mx2 { mx2 = v[2] } 144 vi = vi + 1 145 } 146 } 147 sp_puts(" box pipeline: faces=" as *u8); sp_putn(cB[11]); sp_puts(" tris=" as *u8); sp_putn(hB[1]) 148 sp_puts(" bbox=" as *u8); sp_putn(mx0 - mn0); sp_puts("x" as *u8); sp_putn(mx1 - mn1); sp_puts("x" as *u8); sp_putn(mx2 - mn2) 149 sp_puts(" fx (expect 25600x15360x10240) mismatch=" as *u8); sp_putn(cB[6]); sp_puts("\n" as *u8) 150 var t2t: i64 = 0 151 if cB[11] == 6 { if hB[1] == 12 { if cB[6] == 0 { 152 if mx0 - mn0 == 25600 { if mx1 - mn1 == 15360 { if mx2 - mn2 == 10240 { t2t = 1 } } } 153 } } } 154 let ig2: i64 = bg2_tooth("T2 AUTHORED box reads back: 6/6 planar faces, 12 tris, bbox 100x60x40mm EXACT" as *u8, t2t, fails) 155 156 let mmB: i64 = sgn_unit_mm(stB as *i64) 157 let ig3: i64 = bg2_tooth("T3 authored file DECLARES mm units (self-confirmed)" as *u8, mmB, fails) 158 159 // ---- (B) PERSIST: re-emit the whole as1 reachable graph 160 let stA: i64 = bg2_parse("knowledge/fetched/step_as1.stp" as *u8, o1) 161 if stA == 0 { sp_puts("MISSING as1\nRED\n" as *u8); return 1 } 162 // original baseline 163 let WO: *i64 = sys_mmap(128) as *i64 164 bg2_ctx(stA as *i64, WO) 165 sgt2_tessellate(WO[0] as *i64, WO[1] as *i64, WO[3] as *i64) 166 sgn_tessellate_all(WO[2] as *i64) 167 let cO: *i64 = WO[4] as *i64 168 let c2O: *i64 = WO[5] as *i64 169 let hO: *i64 = m3_hdr(WO[6]) 170 let ckO: i64 = bg2_cksum(WO[6]) 171 // re-emit + re-parse 172 let qn: i64 = sw2_reemit(stA as *i64, "knowledge/as1_reemit.stp" as *u8) 173 let stR: i64 = bg2_parse("knowledge/as1_reemit.stp" as *u8, o1) 174 sp_puts("re-emit: entities=" as *u8); sp_putn(qn); sp_puts(" re-parsed=" as *u8); sp_putn(o1[0]); sp_puts("\n" as *u8) 175 var t4: i64 = 0 176 if qn > 3000 { if stR != 0 { if o1[0] == qn { t4 = 1 } } } 177 let ig4: i64 = bg2_tooth("T4 re-emit: full reachable graph written + re-parsed (count exact)" as *u8, t4, fails) 178 179 let WR: *i64 = sys_mmap(128) as *i64 180 bg2_ctx(stR as *i64, WR) 181 sgt2_tessellate(WR[0] as *i64, WR[1] as *i64, WR[3] as *i64) 182 sgn_tessellate_all(WR[2] as *i64) 183 let cR: *i64 = WR[4] as *i64 184 let c2R: *i64 = WR[5] as *i64 185 let hR: *i64 = m3_hdr(WR[6]) 186 let ckR: i64 = bg2_cksum(WR[6]) 187 sp_puts(" original: planar=" as *u8); sp_putn(cO[11]); sp_puts(" round=" as *u8); sp_putn(c2O[2]) 188 sp_puts(" tris=" as *u8); sp_putn(hO[1]); sp_puts(" cksum=" as *u8); sp_putn(ckO) 189 sp_puts("\n re-emit : planar=" as *u8); sp_putn(cR[11]); sp_puts(" round=" as *u8); sp_putn(c2R[2]) 190 sp_puts(" tris=" as *u8); sp_putn(hR[1]); sp_puts(" cksum=" as *u8); sp_putn(ckR); sp_puts("\n" as *u8) 191 var t5: i64 = 0 192 if cR[11] == cO[11] { if c2R[2] == c2O[2] { if hR[1] == hO[1] { if ckR == ckO { if cR[6] == 0 { t5 = 1 } } } } } 193 let ig5: i64 = bg2_tooth("T5 re-emitted model tessellates IDENTICALLY (25+28 faces, tris, geometry checksum)" as *u8, t5, fails) 194 195 let mmR: i64 = sgn_unit_mm(stR as *i64) 196 let ig6: i64 = bg2_tooth("T6 re-emitted file carries the units block (mm confirmed)" as *u8, mmR, fails) 197 198 // T7 determinism: second re-emit -> same entity count + same byte sum 199 let qn2: i64 = sw2_reemit(stA as *i64, "knowledge/as1_reemit2.stp" as *u8) 200 let bufX: *u8 = sys_mmap(8388608) 201 let nX: i64 = sp_read_file("knowledge/as1_reemit.stp" as *u8, bufX, 8388608) 202 let bufY: *u8 = sys_mmap(8388608) 203 let nY: i64 = sp_read_file("knowledge/as1_reemit2.stp" as *u8, bufY, 8388608) 204 var s1: i64 = 0 205 var s2: i64 = 0 206 var bi: i64 = 0 207 while bi < nX { s1 = s1 + (bufX[bi] as i64); bi = bi + 1 } 208 bi = 0 209 while bi < nY { s2 = s2 + (bufY[bi] as i64); bi = bi + 1 } 210 var t7: i64 = 0 211 if qn2 == qn { if nX == nY { if s1 == s2 { t7 = 1 } } } 212 let ig7: i64 = bg2_tooth("T7 deterministic (byte-identical second re-emit)" as *u8, t7, fails) 213 214 sp_puts("\nP3a WRITE PARITY: (A) AUTHOR analytic B-rep from parameters (planes + full topology, exact dims\n" as *u8) 215 sp_puts("verified through our own reader) + (B) PERSIST any read model (full graph re-serialization incl.\n" as *u8) 216 sp_puts("degree-5 splines + rational Bezier surfaces, pipeline-identical). Faceted write landed earlier.\n" as *u8) 217 sp_puts("Read 53/53 + write author/persist/faceted = INTEROP BOTH WAYS.\n" as *u8) 218 sp_puts("\nfails=" as *u8); sp_putn(fails[0]); sp_puts("\n" as *u8) 219 if fails[0] == 0 { sp_puts("GREEN -- analytic B-rep write 7/7 (author + persist, pipeline-verified)\n" as *u8); return 0 } 220 sp_puts("RED\n" as *u8) 221 return 1 222}