code wiki / (root) / nx_water_activity_test.nx
`nx_water_activity_test.nx` source

↩ module page · 201 lines · 10919 B
1// nx_water_activity_test.nx -- gate for nx_water_activity.  The load-bearing
2// tests are T2 (a real mix lands on the published aw ~0.97 for ice cream mix)
3// and T8/T9, where the hurdle model must NAME the barrier and must refuse to
4// call an ice cream mix safe on water activity -- the reassuring wrong answer.
5// expect_exit: 0   license_tier: ORIGINAL
6import "nx_syscalls.nx"
7import "nx_icecream.nx"
8import "nx_chem_solution.nx"
9import "nx_thermal_process.nx"
10import "nx_water_activity.nx"
11
12func t_puts(s: *u8) -> i64 { var n: i64 = 0; while s[n] != (0 as u8) { n = n + 1 } sys_write(1, s, n); return 0 }
13func t_putn(v: i64) -> i64 { let bb: *u8 = sys_mmap(28); var m: i64 = v; if m < 0 { m = 0 - m; sys_write(1, "-" 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(1, bb, k); return 0 }
14
15func ref_mix() -> *NxIceMix {
  let m: *NxIceMix = nx_ice_mix_new()
  m.fat_g = 120
  m.msnf_g = 100
  m.sugar_g = 150
  m.sugar_mw_q2 = IC_MW_SUCROSE_Q2
  m.sugar_pod = 100
  m.other_solids_g = 30
  m.total_g = 1000
  m.density_q3 = 1100
  return m
26}
27
28func main() -> i64 {
  var pass: i64 = 0
  var total: i64 = 0
31
  // --- T1  pure water is aw = 1.0000 exactly ---
  total = total + 1
  let pure: i64 = wa_from_solutes_q4(1000, 0)
  t_puts("T1 pure water aw x10000 = " as *u8); t_putn(pure); t_puts(" want 10000: " as *u8)
  if pure == 10000 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
37
  // --- T2  a real ice cream mix lands on the published ~0.97 ---
  total = total + 1
  let m: *NxIceMix = ref_mix()
  let aw: i64 = wa_from_mix_q4(m)
  t_puts("T2 ice cream mix aw = " as *u8); t_putn(aw); t_puts(" x10000, published ~0.97 (band 9650..9850): " as *u8)
  var ok2: i64 = 1
  if aw < 9650 { ok2 = 0 }
  if aw > 9850 { ok2 = 0 }
  if ok2 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
47
  // --- T3  aw and freezing point come from the SAME moles: more solute
  //         must both lower aw and deepen the depression ---
  total = total + 1
  let sweet: *NxIceMix = ref_mix()
  sweet.sugar_g = 250
  sweet.other_solids_g = 0
  let aw_s: i64 = wa_from_mix_q4(sweet)
  let fpd_ref: i64 = ic_fpd_milli_c(m)
  let fpd_s: i64 = ic_fpd_milli_c(sweet)
  t_puts("T3 more sugar: aw " as *u8); t_putn(aw); t_puts("->" as *u8); t_putn(aw_s); t_puts(" FPD " as *u8); t_putn(fpd_ref); t_puts("->" as *u8); t_putn(fpd_s); t_puts(" (aw down, FPD up): " as *u8)
  var ok3: i64 = 1
  if aw_s >= aw { ok3 = 0 }
  if fpd_s <= fpd_ref { ok3 = 0 }
  if ok3 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
62
  // --- T4  growth thresholds are ordered as the literature has them ---
  total = total + 1
  let sal: i64 = wa_min_growth_q4(WA_ORG_SALMONELLA)
  let cbot: i64 = wa_min_growth_q4(WA_ORG_CBOT)
  let staph: i64 = wa_min_growth_q4(WA_ORG_STAPH)
  let mold: i64 = wa_min_growth_q4(WA_ORG_MOLD)
  let xero: i64 = wa_min_growth_q4(WA_ORG_XEROPHILE)
  t_puts("T4 min aw salmonella/cbot/staph/mold/xerophile = " as *u8); t_putn(sal); t_puts("/" as *u8); t_putn(cbot); t_puts("/" as *u8); t_putn(staph); t_puts("/" as *u8); t_putn(mold); t_puts("/" as *u8); t_putn(xero); t_puts(" (descending): " as *u8)
  var ok4: i64 = 1
  if cbot >= sal { ok4 = 0 }
  if staph >= cbot { ok4 = 0 }
  if mold >= staph { ok4 = 0 }
  if xero >= mold { ok4 = 0 }
  if ok4 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
77
  // --- T5  THE UNCOMFORTABLE TRUTH: the mix is NOT protected by water ---
  total = total + 1
  let grows_bot: i64 = wa_supports_growth(aw, WA_ORG_CBOT)
  let grows_sal: i64 = wa_supports_growth(aw, WA_ORG_SALMONELLA)
  let nsupp: i64 = wa_organisms_supported(aw)
  let phf: i64 = wa_is_potentially_hazardous(aw)
  t_puts("T5 mix aw supports C.bot=" as *u8); t_putn(grows_bot); t_puts(" salmonella=" as *u8); t_putn(grows_sal); t_puts(" organisms=" as *u8); t_putn(nsupp); t_puts("/6 PHF=" as *u8); t_putn(phf); t_puts(" want 1/1/6/1: " as *u8)
  var ok5: i64 = 1
  if grows_bot != 1 { ok5 = 0 }
  if grows_sal != 1 { ok5 = 0 }
  if nsupp != WA_ORG_N { ok5 = 0 }
  if phf != 1 { ok5 = 0 }
  if ok5 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
91
  // --- T6  a dry product (honey-like, aw 0.60) stops everything ---
  total = total + 1
  let dry: i64 = 6000
  let d_bot: i64 = wa_supports_growth(dry, WA_ORG_CBOT)
  let d_mold: i64 = wa_supports_growth(dry, WA_ORG_MOLD)
  let d_xero: i64 = wa_supports_growth(dry, WA_ORG_XEROPHILE)
  let d_n: i64 = wa_organisms_supported(dry)
  t_puts("T6 aw 0.60: cbot=" as *u8); t_putn(d_bot); t_puts(" mold=" as *u8); t_putn(d_mold); t_puts(" xerophile=" as *u8); t_putn(d_xero); t_puts(" total=" as *u8); t_putn(d_n); t_puts(" want 0/0/0/0: " as *u8)
  if d_n == 0 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
101
  // --- T7  FAIL-CLOSED: no water, bad solute, unknown organism ---
  // An unknown hazard must read as ABLE to grow, never as absent.
  total = total + 1
  let nowater: i64 = wa_from_solutes_q4(0, 1000)
  let negsol: i64 = wa_from_solutes_q4(100, 0 - 5)
  let unk: i64 = wa_supports_growth(9000, 99)
  let inv_grow: i64 = wa_supports_growth(WA_INVALID, WA_ORG_MOLD)
  let inv_phf: i64 = wa_is_potentially_hazardous(WA_INVALID)
  t_puts("T7 no-water=" as *u8); t_putn(nowater); t_puts(" neg-solute=" as *u8); t_putn(negsol); t_puts(" unknown-organism-grows=" as *u8); t_putn(unk); t_puts(" invalid-aw-grows=" as *u8); t_putn(inv_grow); t_puts(" invalid-PHF=" as *u8); t_putn(inv_phf); t_puts(": " as *u8)
  var ok7: i64 = 1
  if nowater != WA_INVALID { ok7 = 0 }
  if negsol != WA_INVALID { ok7 = 0 }
  if unk != 1 { ok7 = 0 }
  if inv_grow != 1 { ok7 = 0 }
  if inv_phf != 1 { ok7 = 0 }
  if ok7 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
118
  // --- T8  THE HURDLE MODEL NAMES THE BARRIER ---
  total = total + 1
  let v_dry: i64 = wa_hurdle_verdict(6600, 6000, 0, 0, 1)
  let v_acid: i64 = wa_hurdle_verdict(3500, 9800, 0, 1, 1)
  let v_retort: i64 = wa_hurdle_verdict(6600, 9800, TP_F0_STANDARD_MS, 0, 1)
  let v_cold: i64 = wa_hurdle_verdict(6600, 9800, 0, 1, 0)
  t_puts("T8 verdicts dry/acid/retort/chilled = " as *u8); t_putn(v_dry); t_puts("/" as *u8); t_putn(v_acid); t_puts("/" as *u8); t_putn(v_retort); t_puts("/" as *u8); t_putn(v_cold); t_puts(" want 1/2/3/4: " as *u8)
  var ok8: i64 = 1
  if v_dry != WA_SAFE_BY_AW { ok8 = 0 }
  if v_acid != WA_SAFE_BY_ACID { ok8 = 0 }
  if v_retort != WA_SAFE_BY_THERMAL { ok8 = 0 }
  if v_cold != WA_SAFE_BY_COLD_CHAIN { ok8 = 0 }
  if ok8 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
132
  // --- T9  ICE CREAM MIX end to end: carried by the COLD CHAIN, and the
  //         model refuses every way of pretending otherwise ---
  total = total + 1
  let real_v: i64 = wa_hurdle_verdict(6600, aw, 0, 1, 0)
  let no_past: i64 = wa_hurdle_verdict(6600, aw, 0, 0, 0)
  let ambient: i64 = wa_hurdle_verdict(6600, aw, 0, 1, 1)
  t_puts("T9 mix chilled+pasteurised=" as *u8); t_putn(real_v); t_puts(" unpasteurised=" as *u8); t_putn(no_past); t_puts(" left ambient=" as *u8); t_putn(ambient); t_puts(" want 4/0/0: " as *u8)
  var ok9: i64 = 1
  if real_v != WA_SAFE_BY_COLD_CHAIN { ok9 = 0 }
  if no_past != WA_UNSAFE { ok9 = 0 }
  if ambient != WA_UNSAFE { ok9 = 0 }
  if ok9 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
145
  // --- T10  the aw hurdle OVERRIDES a low pH and no heat at all ---
  // Milk powder: low-acid, never retorted, ambient -- safe purely on water.
  total = total + 1
  let powder: i64 = wa_hurdle_verdict(6600, 2000, 0, 0, 1)
  let safe_p: i64 = wa_is_safe(6600, 2000, 0, 0, 1)
  t_puts("T10 milk powder aw 0.20 low-acid ambient unheated: verdict=" as *u8); t_putn(powder); t_puts(" is_safe=" as *u8); t_putn(safe_p); t_puts(" want 1/1: " as *u8)
  if powder == WA_SAFE_BY_AW { if safe_p == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) } } else { t_puts("FAIL\n" as *u8) }
153
  // --- T11  the PHF line is exactly 0.85, checked either side ---
  total = total + 1
  let under: i64 = wa_is_potentially_hazardous(8499)
  let at: i64 = wa_is_potentially_hazardous(8500)
  let vu: i64 = wa_hurdle_verdict(6600, 8499, 0, 0, 1)
  let va: i64 = wa_hurdle_verdict(6600, 8500, 0, 0, 1)
  t_puts("T11 aw 0.8499 PHF=" as *u8); t_putn(under); t_puts(" verdict=" as *u8); t_putn(vu); t_puts(" ; aw 0.8500 PHF=" as *u8); t_putn(at); t_puts(" verdict=" as *u8); t_putn(va); t_puts(" want 0/1 then 1/0: " as *u8)
  var ok11: i64 = 1
  if under != 0 { ok11 = 0 }
  if at != 1 { ok11 = 0 }
  if vu != WA_SAFE_BY_AW { ok11 = 0 }
  if va != WA_UNSAFE { ok11 = 0 }
  if ok11 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
167
  // --- T12  inversion: solute needed for a target aw, and the honest
  //          answer when the target is unreachable by adding solute ---
  total = total + 1
  let need: i64 = wa_solute_umol_for_target(600, 8500)
  let rt: i64 = wa_from_solutes_q4(600, need)
  let bad_hi: i64 = wa_solute_umol_for_target(600, 10000)
  let bad_lo: i64 = wa_solute_umol_for_target(600, 0)
  t_puts("T12 solute for aw 0.85 in 600g water = " as *u8); t_putn(need); t_puts(" umol, round-trips to aw=" as *u8); t_putn(rt); t_puts(" ; aw=1.0 -> " as *u8); t_putn(bad_hi); t_puts(" aw=0 -> " as *u8); t_putn(bad_lo); t_puts(": " as *u8)
  var d12: i64 = rt - 8500
  if d12 < 0 { d12 = 0 - d12 }
  var ok12: i64 = 1
  if d12 > 5 { ok12 = 0 }
  if bad_hi != WA_INVALID { ok12 = 0 }
  if bad_lo != WA_INVALID { ok12 = 0 }
  if ok12 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
183
  // --- T13  aw is monotone decreasing in solute (no inversions) ---
  total = total + 1
  var s13: i64 = 0
  var prev13: i64 = 10001
  var ok13: i64 = 1
  while s13 <= 20000000 {
      let a: i64 = wa_from_solutes_q4(600, s13)
      if a > prev13 { ok13 = 0 }
      prev13 = a
      s13 = s13 + 1000000
  }
  t_puts("T13 aw monotone non-increasing in solute (0..20 mol): " as *u8)
  if ok13 == 1 { pass = pass + 1; t_puts("PASS\n" as *u8) } else { t_puts("FAIL\n" as *u8) }
197
  t_puts("WATER-ACTIVITY-GATE passed " as *u8); t_putn(pass); t_puts("/" as *u8); t_putn(total)
  if pass == total { t_puts(" verdict=GREEN\n" as *u8); sys_exit(0); return 0 }
  t_puts(" verdict=RED\n" as *u8); sys_exit(1); return 1
201}
nx_water_activity_test.nx source

`nx_water_activity_test.nx` source
