`nx_labsci_svc.nx` source

↩ module page · 541 lines · 24538 B
1// nx_labsci_svc.nx -- LAB-SCIENCE SERVICE FACADE.  The second organ built on
2// the nx_service base, over a COMPLETELY DIFFERENT domain from nx_gtm_svc --
3// peptide chemistry, frozen-dessert physics, thermal death kinetics, and
4// water activity.  Its existence is the proof the base generalises: the
5// envelope, escaping, error shape, dispatch and self-description are all
6// inherited unchanged; only the verbs and their handlers differ.
7//
8// Two facades, one base, two unrelated sciences = the architecture scales to
9// all seventeen organs by the same recipe, not just the go-to-market subset.
10//
11// USAGE:  nx_labsci_svc <verb> [args...]
12//   describe
13//   peptide.analyze   <sequence>                 (1-letter code, e.g. DRVYIHPF)
14//   icecream.freeze   <fat_g> <msnf_g> <sugar_g> <sugar_mw_q2> <total_g>
15//   thermal.pasteurize <temp_milli_c> <hold_ms> <highfat 0|1>
16//   thermal.barrier   <ph_milli> <ambient_stable 0|1>
17//   water.activity    <water_g> <solute_umol>
18//   stability.shelf   <ref_days> <ref_temp_c> <target_temp_c> <months> <loss_pct>
19//   stability.mkt     <temp_c>... (2-8 whole Celsius readings, equal duration)
20//   stability.separate <t1> <rh1> <months1> <t2> <rh2> <months2> <t3> <rh3> <months3>
21//
22// Same discipline as the sibling facade: NO new domain logic, every number
23// comes from a gated library, structured errors, fail-closed on bad input.
24//
25// stability.separate is the verb that answers what stability.shelf can only
26// flag.  ICH gives two storage conditions differing in BOTH temperature and
27// humidity, so the activation energy recovered from them is lumped; three
28// conditions break it apart, and the verb REFUSES a design that cannot --
29// reporting the conditioning separately, because a design can be solvable
30// and still useless.
31//
32// The stability verbs report GND_VALIDATED, a tier above the rest of
33// this facade's business-adjacent output, and the reason is specific: their
34// activation energy is RECOVERED from ICH's own equivalence claim and then
35// lands inside published Q10 and Ea bands that were never inputs to the
36// derivation.  That is the same evidence class as the pasteurisation z-value,
37// which is why thermal.* already reports VALIDATED.
38//
39// genealogy_id: nishi_food_science_suite + service_infrastructure
40
41import "nx_syscalls.nx"
42import "nx_service.nx"
43import "nx_grounding.nx"
44import "nx_peptide.nx"
45import "nx_icecream.nx"
46import "nx_thermal_process.nx"
47import "nx_water_activity.nx"
48import "nx_pow10.nx"
49import "nx_shelf_life.nx"
50import "nx_stability.nx"
51import "nx_stability_rh.nx"
52import "nx_ice_recrystal.nx"
53import "nx_ice_distribution.nx"
54const LAB_MAGIC_13000: i64 = 13000
55const LAB_MAGIC_8192: i64 = 8192
56
57const LAB_SVC: *u8 = "labsci" as *u8
58
59// ===== describe =======================================================
60
61func lab_describe(j: *NxJson) -> i64 {
  nsvc_ok_open(j, LAB_SVC, "describe" as *u8)
  nj_kv_str(j, "service" as *u8, "peptide + frozen dessert + thermal + water activity" as *u8)
  nj_comma(j)
  nj_kv_int(j, "verb_count" as *u8, 9)
  nj_comma(j)
  nj_key(j, "verbs" as *u8)
  nj_puts(j, "[" as *u8)
  nj_puts(j, "{\"verb\":\"peptide.analyze\",\"args\":[\"sequence:string\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"icecream.freeze\",\"args\":[\"fat_g:int\",\"msnf_g:int\",\"sugar_g:int\",\"sugar_mw_q2:int\",\"total_g:int\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"thermal.pasteurize\",\"args\":[\"temp_milli_c:int\",\"hold_ms:int\",\"highfat:int\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"thermal.barrier\",\"args\":[\"ph_milli:int\",\"ambient_stable:int\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"water.activity\",\"args\":[\"water_g:int\",\"solute_umol:int\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"stability.shelf\",\"args\":[\"ref_days:int\",\"ref_temp_c:int\",\"target_temp_c:int\",\"months:int\",\"loss_pct:int\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"stability.mkt\",\"args\":[\"temp_c:int...\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"stability.separate\",\"args\":[\"t1\",\"rh1\",\"months1\",\"t2\",\"rh2\",\"months2\",\"t3\",\"rh3\",\"months3\"]}," as *u8)
  nj_puts(j, "{\"verb\":\"icecream.recrystal\",\"args\":[\"fat_g\",\"msnf_g\",\"sugar_g\",\"total_g\",\"t_cold_milli_c\",\"t_warm_milli_c\",\"cycles\",\"days\"]}" as *u8)
  nj_puts(j, "]" as *u8)
  nsvc_ok_close_g(j, GND_DERIVED)
  return 0
81}
82
83// ===== peptide.analyze ================================================
84
85func lab_peptide_analyze(j: *NxJson, seq: *u8) -> i64 {
  if seq[0] == (0 as u8) {
      nsvc_error(j, LAB_SVC, "peptide.analyze" as *u8, NSVC_ERR_BAD_ARGS, "empty sequence" as *u8)
      return 0
  }
  if pep_seq_valid(seq) != 1 {
      nsvc_error(j, LAB_SVC, "peptide.analyze" as *u8, NSVC_ERR_REFUSED, "sequence has a non-standard residue (refused, not approximated)" as *u8)
      return 0
  }
  let mono: i64 = pep_mass_mono_q4(seq)
  nsvc_ok_open(j, LAB_SVC, "peptide.analyze" as *u8)
  nj_kv_str(j, "sequence" as *u8, seq)
  nj_comma(j)
  nj_kv_int(j, "length" as *u8, pep_len(seq))
  nj_comma(j)
  nj_kv_int(j, "monoisotopic_mass_q4" as *u8, mono)
  nj_comma(j)
  nj_kv_int(j, "average_mass_q4" as *u8, pep_mass_avg_q4(seq))
  nj_comma(j)
  nj_kv_int(j, "mh_plus_q4" as *u8, pep_mz_q4(mono, 1))
  nj_comma(j)
  nj_kv_int(j, "m2h_plus_q4" as *u8, pep_mz_q4(mono, 2))
  nj_comma(j)
  nj_kv_int(j, "pi_milli" as *u8, pep_pi_milli(seq))
  nj_comma(j)
  nj_kv_int(j, "gravy_q3" as *u8, pep_gravy_q3(seq))
  nj_comma(j)
  nj_kv_int(j, "ext_coeff_280" as *u8, pep_ext_coeff_280(seq, 0))
  nj_comma(j)
  nj_kv_bool(j, "a280_quantifiable" as *u8, pep_a280_quantifiable(seq, 0))
  // mass/formula are ANCHORED but pI/GRAVY are DERIVED -> weakest link DERIVED
  nsvc_ok_close_g(j, GND_DERIVED)
  return 0
118}
119
120// ===== icecream.freeze ================================================
121
122func lab_icecream_freeze(j: *NxJson, fat: i64, msnf: i64, sugar: i64, mw: i64, total: i64) -> i64 {
  if total <= 0 {
      nsvc_error(j, LAB_SVC, "icecream.freeze" as *u8, NSVC_ERR_BAD_ARGS, "total_g must be positive" as *u8)
      return 0
  }
  let m: *NxIceMix = nx_ice_mix_new()
  m.fat_g = fat
  m.msnf_g = msnf
  m.sugar_g = sugar
  if mw > 0 { m.sugar_mw_q2 = mw }
  m.total_g = total
  let water: i64 = ic_water_g(m)
  if water <= 0 {
      nsvc_error(j, LAB_SVC, "icecream.freeze" as *u8, NSVC_ERR_REFUSED, "solids exceed batch mass (incoherent mix)" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "icecream.freeze" as *u8)
  nj_kv_int(j, "water_g" as *u8, water)
  nj_comma(j)
  nj_kv_int(j, "freezing_point_milli_c" as *u8, ic_freezing_point_milli_c(m))
  nj_comma(j)
  nj_kv_int(j, "frozen_permil_at_minus13c" as *u8, ic_frozen_water_permil(m, 0 - LAB_MAGIC_13000))
  nj_comma(j)
  nj_kv_int(j, "scoop_temp_milli_c" as *u8, ic_scoop_temp_milli_c(m))
  nj_comma(j)
  nj_kv_int(j, "pac_index" as *u8, ic_pac_index(m))
  nj_comma(j)
  nj_kv_bool(j, "sandy_risk" as *u8, ic_sandy_risk(m))
  // freezing point is VALIDATED vs published curves, PAC is DERIVED -> DERIVED
  nsvc_ok_close_g(j, GND_DERIVED)
  return 0
153}
154
155// ===== thermal.pasteurize =============================================
156
157func lab_thermal_pasteurize(j: *NxJson, temp: i64, hold: i64, highfat: i64) -> i64 {
  if hold <= 0 {
      nsvc_error(j, LAB_SVC, "thermal.pasteurize" as *u8, NSVC_ERR_BAD_ARGS, "hold_ms must be positive" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "thermal.pasteurize" as *u8)
  nj_kv_int(j, "temp_milli_c" as *u8, temp)
  nj_comma(j)
  nj_kv_int(j, "hold_ms" as *u8, hold)
  nj_comma(j)
  nj_kv_int(j, "equiv_lethality_ms" as *u8, tp_pasteurization_equiv_ms(temp, hold, highfat))
  nj_comma(j)
  nj_kv_bool(j, "compliant" as *u8, tp_pasteurization_compliant(temp, hold, highfat))
  nj_comma(j)
  nj_kv_int(j, "f0_ms" as *u8, tp_f0_ms(hold, temp))
  // z recovered from two legal PMO schedules + lands in published band
  nsvc_ok_close_g(j, GND_VALIDATED)
  return 0
175}
176
177// ===== thermal.barrier ================================================
178
179func lab_thermal_barrier(j: *NxJson, ph: i64, ambient: i64) -> i64 {
  nsvc_ok_open(j, LAB_SVC, "thermal.barrier" as *u8)
  nj_kv_int(j, "ph_milli" as *u8, ph)
  nj_comma(j)
  nj_kv_bool(j, "high_acid" as *u8, cs_is_high_acid(ph))
  nj_comma(j)
  nj_kv_int(j, "barrier_required" as *u8, tp_barrier_required(ph, ambient))
  // the 4.6 botulinum line is an established regulatory threshold
  nsvc_ok_close_g(j, GND_VALIDATED)
  return 0
189}
190
191// ===== water.activity =================================================
192
193func lab_water_activity(j: *NxJson, water: i64, solute: i64) -> i64 {
  let aw: i64 = wa_from_solutes_q4(water, solute)
  if aw == WA_INVALID {
      nsvc_error(j, LAB_SVC, "water.activity" as *u8, NSVC_ERR_REFUSED, "no water / bad solute -- aw undefined (refused, not reported as dry-and-safe)" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "water.activity" as *u8)
  nj_kv_int(j, "aw_q4" as *u8, aw)
  nj_comma(j)
  nj_kv_bool(j, "potentially_hazardous" as *u8, wa_is_potentially_hazardous(aw))
  nj_comma(j)
  nj_kv_int(j, "organisms_supported" as *u8, wa_organisms_supported(aw))
  nj_comma(j)
  nj_kv_bool(j, "botulinum_can_grow" as *u8, wa_supports_growth(aw, WA_ORG_CBOT))
  // aw is DERIVED (Raoult) but the organism growth thresholds are ASSERTED
  nsvc_ok_close_g(j, GND_ASSERTED)
  return 0
210}
211
212// ===== stability.shelf ================================================
213
214// Projects a shelf life from one storage temperature to another and reports
215// what survives to the end of it.  The temperature move uses the fractional
216// Q10 path, so a non-decade target (the 25 C -> 40 C case the whole category
217// rests on) is answerable instead of refused.
218func lab_stability_shelf(j: *NxJson, ref_days: i64, ref_c: i64, target_c: i64, months: i64, loss_pct: i64) -> i64 {
  var q10: i64 = 0
  var ea: i64 = 0
  var life: i64 = 0
  var retained: i64 = 0
  if ref_days <= 0 {
      nsvc_error(j, LAB_SVC, "stability.shelf" as *u8, NSVC_ERR_BAD_ARGS, "ref_days must be positive" as *u8)
      return 0
  }
  if months <= 0 {
      nsvc_error(j, LAB_SVC, "stability.shelf" as *u8, NSVC_ERR_BAD_ARGS, "months must be positive" as *u8)
      return 0
  }
  q10 = stab_ich_q10_q3()
  ea = stab_ich_ea_j()
  life = sl_at_temp_frac(ref_days, ref_c, q10, target_c)
  if life <= 0 {
      nsvc_error(j, LAB_SVC, "stability.shelf" as *u8, NSVC_ERR_REFUSED, "shelf life not computable for this temperature pair (refused, not extrapolated)" as *u8)
      return 0
  }
  retained = stab_potency_retained_permil(months, months, loss_pct)
  if retained == STAB_INVALID {
      nsvc_error(j, LAB_SVC, "stability.shelf" as *u8, NSVC_ERR_BAD_ARGS, "loss_pct must be in 1..99" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "stability.shelf" as *u8)
  nj_kv_int(j, "shelf_days_at_target" as *u8, life)
  nj_comma(j)
  nj_kv_int(j, "q10_q3_recovered" as *u8, q10)
  nj_comma(j)
  nj_kv_int(j, "activation_energy_j_per_mol" as *u8, ea)
  nj_comma(j)
  nj_kv_int(j, "potency_retained_permil" as *u8, retained)
  nj_comma(j)
  nj_kv_int(j, "overage_needed_permil" as *u8, stab_overage_permil(retained))
  nj_comma(j)
  nj_kv_bool(j, "significant_change" as *u8, stab_is_significant_change(retained))
  nj_comma(j)
  // A hard limit the consumer must not be able to route around.
  nj_kv_bool(j, "accelerated_alone_sets_shelf_life" as *u8, stab_accelerated_alone_sets_shelf_life())
  nj_comma(j)
  // ICH's two conditions differ in humidity as well as temperature, so the
  // recovered figure is an EFFECTIVE Ea. Surfaced at the wire rather than
  // buried in a header, because the number is only safe along that path.
  nj_kv_bool(j, "activation_energy_is_lumped_with_humidity" as *u8, stab_ea_is_lumped_with_humidity())
  nj_comma(j)
  nj_kv_str(j, "valid_along" as *u8, "the ICH storage path (60-75% RH); off-path use over-states degradation for a drier product" as *u8)
  nsvc_ok_close_g(j, GND_VALIDATED)
  return 0
267}
268
269// ===== stability.mkt ==================================================
270
271// Mean kinetic temperature of an equal-duration temperature history, beside
272// the arithmetic mean it is so often confused with.  Reporting both in one
273// response is the point: the gap between them IS the finding.
274func lab_stability_mkt(j: *NxJson, argc: i64, argv: *i64) -> i64 {
  var n: i64 = 0
  var i: i64 = 0
  var ea: i64 = 0
  var mkt: i64 = 0
  var mean: i64 = 0
  let temps: *i64 = sys_mmap(16 * 8)
  let hours: *i64 = sys_mmap(16 * 8)
  while i < 8 {
      if argc > i + 2 {
          temps[n] = nsvc_arg_int(argc, argv, i + 2)
          hours[n] = 1
          n = n + 1
      }
      i = i + 1
  }
  if n < 2 {
      nsvc_error(j, LAB_SVC, "stability.mkt" as *u8, NSVC_ERR_BAD_ARGS, "give 2..8 whole-Celsius readings" as *u8)
      return 0
  }
  ea = stab_ich_ea_j()
  mkt = stab_mkt_centi_kelvin(temps, hours, n, ea)
  mean = stab_mean_centi_kelvin(temps, hours, n)
  if mkt == STAB_INVALID {
      nsvc_error(j, LAB_SVC, "stability.mkt" as *u8, NSVC_ERR_REFUSED, "history not computable (refused, not averaged)" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "stability.mkt" as *u8)
  nj_kv_int(j, "readings" as *u8, n)
  nj_comma(j)
  nj_kv_int(j, "mkt_centi_kelvin" as *u8, mkt)
  nj_comma(j)
  nj_kv_int(j, "mkt_centi_c" as *u8, mkt - STAB_KELVIN_OFF_CENTI)
  nj_comma(j)
  nj_kv_int(j, "arithmetic_mean_centi_c" as *u8, mean - STAB_KELVIN_OFF_CENTI)
  nj_comma(j)
  nj_kv_int(j, "understatement_centi_c" as *u8, mkt - mean)
  nj_comma(j)
  nj_kv_int(j, "activation_energy_j_per_mol" as *u8, ea)
  nj_comma(j)
  nj_kv_str(j, "note" as *u8, "degradation is convex in temperature, so the arithmetic mean UNDER-states it" as *u8)
  nsvc_ok_close_g(j, GND_VALIDATED)
  return 0
317}
318
319// ===== stability.separate =============================================
320
321// Separates the temperature axis from the moisture axis given three real
322// storage conditions.  This is the verb that answers what stability.shelf
323// can only flag: the lumped parameter is only lumped because ICH gives two
324// conditions, and three break it apart.
325//
326// REFUSES rather than guesses on a design that cannot support the split,
327// and reports the CONDITIONING separately -- a design can be solvable and
328// still useless, which is the failure mode that produces a confident,
329// specific, meaningless activation energy.
330func lab_stability_separate(j: *NxJson, argc: i64, argv: *i64) -> i64 {
  var t1: i64 = 0
  var r1: i64 = 0
  var m1: i64 = 0
  var t2: i64 = 0
  var r2: i64 = 0
  var m2: i64 = 0
  var t3: i64 = 0
  var r3: i64 = 0
  var m3: i64 = 0
  var rc: i64 = 0
  var cond: i64 = 0
  let out: *i64 = sys_mmap(8 * 4)
  if argc < 11 {
      nsvc_error(j, LAB_SVC, "stability.separate" as *u8, NSVC_ERR_BAD_ARGS, "need 9 args: t1 rh1 months1 t2 rh2 months2 t3 rh3 months3" as *u8)
      return 0
  }
  t1 = nsvc_arg_int(argc, argv, 2); r1 = nsvc_arg_int(argc, argv, 3); m1 = nsvc_arg_int(argc, argv, 4)
  t2 = nsvc_arg_int(argc, argv, 5); r2 = nsvc_arg_int(argc, argv, 6); m2 = nsvc_arg_int(argc, argv, 7)
  t3 = nsvc_arg_int(argc, argv, 8); r3 = nsvc_arg_int(argc, argv, 9); m3 = nsvc_arg_int(argc, argv, 10)
  cond = srh_design_conditioning_permil(t1, r1, t2, r2, t3, r3)
  rc = srh_separate(t1, r1, m1, t2, r2, m2, t3, r3, m3, out)
  if rc == SRH_UNDERDETERMINED {
      nsvc_error(j, LAB_SVC, "stability.separate" as *u8, NSVC_ERR_REFUSED, "these three conditions cannot separate temperature from humidity (singular design -- one axis never varies, or the points are collinear in 1/T vs RH)" as *u8)
      return 0
  }
  if rc != 0 {
      nsvc_error(j, LAB_SVC, "stability.separate" as *u8, NSVC_ERR_BAD_ARGS, "lifetimes must all be positive" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "stability.separate" as *u8)
  nj_kv_int(j, "thermal_activation_energy_j_per_mol" as *u8, srh_ea_from_e10(out[0]))
  nj_comma(j)
  nj_kv_int(j, "moisture_sensitivity_b10_milli_per_pct_rh" as *u8, out[1])
  nj_comma(j)
  nj_kv_bool(j, "moisture_sensitivity_in_published_band" as *u8, srh_b10_in_published_band(out[1]))
  nj_comma(j)
  nj_kv_int(j, "lumped_activation_energy_j_per_mol" as *u8, stab_ich_ea_j())
  nj_comma(j)
  nj_kv_int(j, "design_conditioning_permil" as *u8, cond)
  nj_comma(j)
  // Solvable and usable are different questions, and only one of them is
  // answered by the determinant being non-zero.
  nj_kv_bool(j, "design_is_well_conditioned" as *u8, srh_design_is_well_conditioned(t1, r1, t2, r2, t3, r3))
  nj_comma(j)
  nj_kv_bool(j, "ich_pair_alone_could_separate" as *u8, srh_ich_pair_alone_is_separable())
  nsvc_ok_close_g(j, GND_VALIDATED)
  return 0
378}
379
380// ===== icecream.recrystal =============================================
381
382// When does this ice cream go coarse?  Not a chemical question: frozen
383// dessert fails by ice crystal growth, so this verb answers in crystal
384// diameter and texture grade rather than in potency.
385//
386// It reports the steady-storage and cycled outcomes SIDE BY SIDE from the
387// same starting point, because the gap between them is the whole finding --
388// two histories at the same mean temperature do not give the same product,
389// and a consumer that only saw one number could not tell.
390func lab_icecream_recrystal(j: *NxJson, fat: i64, msnf: i64, sugar: i64, total: i64, tcold: i64, twarm: i64, cycles: i64, days: i64) -> i64 {
  if total <= 0 {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_BAD_ARGS, "total_g must be positive" as *u8)
      return 0
  }
  if days < 0 {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_BAD_ARGS, "days must not be negative" as *u8)
      return 0
  }
  let m: *NxIceMix = nx_ice_mix_new()
  m.fat_g = fat
  m.msnf_g = msnf
  m.sugar_g = sugar
  m.total_g = total
  if ic_water_g(m) <= 0 {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_REFUSED, "solids exceed batch mass (incoherent mix)" as *u8)
      return 0
  }
  let cyc: i64 = ir_cycled_fraction_permil(m, tcold, twarm)
  if cyc == IR_INVALID {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_REFUSED, "frozen fraction not computable at these temperatures (refused, not extrapolated)" as *u8)
      return 0
  }
  let k: i64 = IR_K_DEFAULT
  // Use the coupling DERIVED from a crystal size distribution rather than
  // the asserted constant. Measurement showed the assertion optimistic by
  // roughly 85%: it implicitly equates the number fraction of crystals lost
  // with the mass fraction melted, and small crystals are numerous and
  // light. Reporting the optimistic figure would under-state damage, which
  // is the direction that calls a defective product acceptable.
  let dist: *IceDist = dist_reference()
  let growth: i64 = dist_volume_growth_permil(dist, cyc)
  if growth == DIST_INVALID {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_REFUSED, "cycled fraction outside the range the distribution can resolve (refused, not extrapolated)" as *u8)
      return 0
  }
  let steady: i64 = ir_size_isothermal(IR_FRESH_D, k, days)
  let cycled: i64 = ir_size_combined_growth(IR_FRESH_D, k, days, growth, cycles)
  let optimistic: i64 = ir_size_combined(IR_FRESH_D, k, days, cyc, cycles)
  if cycled == IR_INVALID {
      nsvc_error(j, LAB_SVC, "icecream.recrystal" as *u8, NSVC_ERR_REFUSED, "crystal growth exceeded the model's representable range" as *u8)
      return 0
  }
  nsvc_ok_open(j, LAB_SVC, "icecream.recrystal" as *u8)
  nj_kv_int(j, "cycled_ice_fraction_permil" as *u8, cyc)
  nj_comma(j)
  nj_kv_int(j, "steady_diameter_deci_um" as *u8, steady)
  nj_comma(j)
  nj_kv_int(j, "cycled_diameter_deci_um" as *u8, cycled)
  nj_comma(j)
  nj_kv_int(j, "steady_texture_grade" as *u8, ir_texture_grade(steady))
  nj_comma(j)
  nj_kv_int(j, "cycled_texture_grade" as *u8, ir_texture_grade(cycled))
  nj_comma(j)
  nj_kv_int(j, "per_cycle_volume_growth_permil" as *u8, growth)
  nj_comma(j)
  nj_kv_int(j, "diameter_if_coupling_were_asserted" as *u8, optimistic)
  nj_comma(j)
  nj_kv_bool(j, "coupling_derived_from_distribution" as *u8, 1)
  nj_comma(j)
  nj_kv_bool(j, "distribution_is_reference_not_measured" as *u8, dist_profile_is_reference())
  nj_comma(j)
  // Say which numbers are load-bearing assumptions rather than physics.
  nj_kv_bool(j, "rate_constant_is_recovered_from_data" as *u8, ir_k_is_recovered(k))
  nj_comma(j)
  nj_kv_bool(j, "geometric_coupling_is_asserted" as *u8, 0)
  nj_comma(j)
  nj_kv_bool(j, "mkt_describes_this_mechanism" as *u8, ir_mkt_describes_recrystallisation())
  nj_comma(j)
  nj_kv_str(j, "note" as *u8, "recrystallisation is driven by the temperature OSCILLATION, so no single effective temperature (including MKT) can represent it" as *u8)
  // The cycled ice fraction is DERIVED from the validated freezing curve,
  // but the geometric coupling and the default rate constant are asserted,
  // and weakest-link governs.
  nsvc_ok_close_g(j, GND_ASSERTED)
  return 0
465}
466
467// ===== argv join (shared shape with the sibling facade) ===============
468
469func lab_join_args(argc: i64, argv: *i64, start: i64, out: *u8, cap: i64) -> i64 {
  var w: i64 = 0
  var k: i64 = start
  while k < argc {
      if w > 0 {
          if w < cap { out[w] = 32 as u8; w = w + 1 }
      }
      let s: *u8 = argv[k] as *u8
      var i: i64 = 0
      while s[i] != (0 as u8) {
          if w < cap { out[w] = s[i]; w = w + 1 }
          i = i + 1
      }
      k = k + 1
  }
  if w < cap { out[w] = 0 as u8 } else { out[cap - 1] = 0 as u8 }
  return w
486}
487
488// ===== dispatch =======================================================
489
490func lab_dispatch(argc: i64, argv: *i64, j: *NxJson) -> i64 {
  let verb: *u8 = nsvc_arg(argc, argv, 1)
  if verb[0] == (0 as u8) { lab_describe(j); return 0 }
  if nsvc_streq(verb, "describe" as *u8) == 1 { lab_describe(j); return 0 }
  if nsvc_streq(verb, "peptide.analyze" as *u8) == 1 {
      let sbuf: *u8 = sys_mmap(512)
      lab_join_args(argc, argv, 2, sbuf, 512)
      lab_peptide_analyze(j, sbuf)
      return 0
  }
  if nsvc_streq(verb, "icecream.freeze" as *u8) == 1 {
      lab_icecream_freeze(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3), nsvc_arg_int(argc, argv, 4), nsvc_arg_int(argc, argv, 5), nsvc_arg_int(argc, argv, 6))
      return 0
  }
  if nsvc_streq(verb, "thermal.pasteurize" as *u8) == 1 {
      lab_thermal_pasteurize(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3), nsvc_arg_int(argc, argv, 4))
      return 0
  }
  if nsvc_streq(verb, "thermal.barrier" as *u8) == 1 {
      lab_thermal_barrier(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3))
      return 0
  }
  if nsvc_streq(verb, "stability.shelf" as *u8) == 1 {
      lab_stability_shelf(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3), nsvc_arg_int(argc, argv, 4), nsvc_arg_int(argc, argv, 5), nsvc_arg_int(argc, argv, 6))
      return 0
  }
  if nsvc_streq(verb, "stability.mkt" as *u8) == 1 {
      lab_stability_mkt(j, argc, argv)
      return 0
  }
  if nsvc_streq(verb, "stability.separate" as *u8) == 1 {
      lab_stability_separate(j, argc, argv)
      return 0
  }
  if nsvc_streq(verb, "icecream.recrystal" as *u8) == 1 {
      lab_icecream_recrystal(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3), nsvc_arg_int(argc, argv, 4), nsvc_arg_int(argc, argv, 5), nsvc_arg_int(argc, argv, 6), nsvc_arg_int(argc, argv, 7), nsvc_arg_int(argc, argv, 8), nsvc_arg_int(argc, argv, 9))
      return 0
  }
  if nsvc_streq(verb, "water.activity" as *u8) == 1 {
      lab_water_activity(j, nsvc_arg_int(argc, argv, 2), nsvc_arg_int(argc, argv, 3))
      return 0
  }
  nsvc_error(j, LAB_SVC, verb, NSVC_ERR_UNKNOWN_VERB, "unknown verb; call describe for the catalog" as *u8)
  return 1
534}
535
536func main(argc: i64, argv: *i64) -> i64 {
  let j: *NxJson = nx_json_new(LAB_MAGIC_8192)
  let rc: i64 = lab_dispatch(argc, argv, j)
  nj_flush(j)
  return rc
541}
nx_labsci_svc.nx source

`nx_labsci_svc.nx` source
