code wiki / _hdl_build / nx_breeder_lab.nx
`nx_breeder_lab.nx` source

↩ module page · 240 lines · 11769 B
1// nx_breeder_lab.nx -- THE GENETICS LAB SCREEN + the tooth that proves the science is VISIBLE.
2// Emits a real game screen: two parents, their Punnett square, and the litter they actually produced --
3// and then MEASURES the rendered litter's visible phenotype ratios against the Punnett prediction.
4// That measurement is the difference between "we have genetics" and "the player can SEE genetics".
5// Educational payload, on screen: genotype notation (Hh), the 4-cell square, the predicted odds, and
6// the observed count -- the exact loop a biology student runs by hand.
7// license_tier: ORIGINAL
8import "nx_syscalls.nx"
9import "nx_game_taming.nx"
10import "nx_creature_render.nx"
11import "nx_game_art.nx"
12import "nx_game_critic.nx"
13import "nx_png_write.nx"
14
15const LW: i64 = 640
16const LH: i64 = 400
17
18func lp(s: *u8) -> i64 { var n: i64 = 0; while s[n] != (0 as u8) { n = n + 1 } sys_write(1, s, n); return 0 }
19func ln2(v: i64) -> i64 {
  let t: *u8 = sys_mmap(32) as *u8
  var m: i64 = v
  var w: i64 = 0
  if m < 0 { t[w] = 45 as u8; w = w + 1; m = 0 - m }
  if m == 0 { t[w] = 48 as u8; sys_write(1, t, w + 1); return 0 }
  let d: *u8 = sys_mmap(32) as *u8
  var k: i64 = 0
  while m > 0 { d[k] = (48 + (m % 10)) as u8; m = m / 10; k = k + 1 }
  var j: i64 = 0
  while j < k { t[w] = d[k - 1 - j]; w = w + 1; j = j + 1 }
  sys_write(1, t, w)
  return 0
32}
33func lnl() -> i64 { lp("\n" as *u8); return 0 }
34
35func lab_png(fb: *i64, path: *u8) -> i64 {
  let rgb: *u8 = sys_mmap(LW*LH*3) as *u8
  var p2: i64 = 0
  let n: i64 = LW*LH
  while p2 < n {
      let v: i64 = fb[p2]
      rgb[p2*3] = (v & 0xff) as u8
      rgb[p2*3+1] = ((v >> 8) & 0xff) as u8
      rgb[p2*3+2] = ((v >> 16) & 0xff) as u8
      p2 = p2 + 1
  }
  return nx_png_write_rgb(path, rgb, LW, LH)
47}
48
49// draw a parent's genotype in real notation, e.g. "Hh" -- the vocabulary the game teaches by using it
50func lab_geno_str(fb: *i64, x: i64, y: i64, g: i64, t: i64, col: i64) -> i64 {
  var cx2: i64 = gr_text(fb, LW, LH, x, y, gx_allele_char(t, gx_a0(g, t)), col)
  cx2 = gr_text(fb, LW, LH, cx2, y, gx_allele_char(t, gx_a1(g, t)), col)
  return cx2
54}
55
56func main() -> i64 {
  let fb: *i64 = sys_mmap(LW*LH*8) as *i64
  art_nebula(fb, LW, LH, 12, 10, 30, 34, 18, 58)
  art_stardust(fb, LW, LH, 150)
  let ink: i64 = gr_pack(238,234,255)
  let dim: i64 = gr_pack(150,146,180)
  let hot: i64 = gr_pack(255,200,110)
63
  gr_rect(fb, LW, LH, 0, 0, LW-1, 24, gr_pack(16,14,34))
  gr_text(fb, LW, LH, 10, 8, "NISHI BREEDERS -- GENETICS LAB" as *u8, ink)
  gr_text(fb, LW, LH, 300, 8, "coat: incomplete dominance (RR crimson / RW rose / WW white)" as *u8, dim)
67
  // ---- the cross under study: two ROSE (heterozygous) parents. Textbook F1 x F1. ----
  let ros: *i64 = sys_mmap(en_bytes(300, COMP_NC)) as *i64
  en_init(ros, 300, COMP_NC)
  var het: i64 = gx_set(0, GX_COAT, 0, 1)          // R W  -> rose
  het = gx_set(het, GX_HORN, 0, 1)                  // H h  -> horned carrier
  het = gx_set(het, GX_WING, 0, 1)                  // W w  -> wingless carrier
  het = gx_set(het, GX_TAIL, 1, 1)
  het = gx_set(het, GX_ELEM, 1, 2)                  // A B  -> codominant AB
76
  let mum: i64 = en_spawn(ros)
  en_set(ros, mum, C_SPEC, 1); en_set(ros, mum, C_IV0, 16); en_set(ros, mum, C_IV1, 16)
  en_set(ros, mum, C_IV2, 16); en_set(ros, mum, C_IV3, 16); en_set(ros, mum, C_LVL, 6)
  en_set(ros, mum, C_XP, 0); en_set(ros, mum, C_NICK, 0)
  en_set(ros, mum, C_GENO, het); en_set(ros, mum, C_SEXG, gx_sexgene(GX_SEX_F, 1, 0))
  en_set(ros, mum, C_POLY, 85); en_set(ros, mum, C_EPI, 0)
  en_set(ros, mum, C_TRUST, 800); en_set(ros, mum, C_WILL, 200); en_set(ros, mum, C_DISP, 800)
84
  let dad: i64 = en_spawn(ros)
  en_set(ros, dad, C_SPEC, 1); en_set(ros, dad, C_IV0, 16); en_set(ros, dad, C_IV1, 16)
  en_set(ros, dad, C_IV2, 16); en_set(ros, dad, C_IV3, 16); en_set(ros, dad, C_LVL, 6)
  en_set(ros, dad, C_XP, 0); en_set(ros, dad, C_NICK, 1)
  en_set(ros, dad, C_GENO, het); en_set(ros, dad, C_SEXG, gx_sexgene(GX_SEX_M, 1, 0))
  en_set(ros, dad, C_POLY, 85); en_set(ros, dad, C_EPI, 0)
  en_set(ros, dad, C_TRUST, 800); en_set(ros, dad, C_WILL, 200); en_set(ros, dad, C_DISP, 800)
92
  // ---- parents panel ----
  gr_rect(fb, LW, LH, 8, 32, 210, 190, gr_pack(20,18,42))
  gr_rect(fb, LW, LH, 8, 32, 210, 34, gr_pack(70,62,120))
  gr_text(fb, LW, LH, 14, 40, "PARENTS" as *u8, ink)
  cr_creature_gx(fb, LW, LH, 62, 110, 36, 1, 16,16,16,16, het, 85, 0)
  cr_creature_gx(fb, LW, LH, 156, 110, 36, 1, 16,16,16,16, het, 85, 0)
  var gx1: i64 = gr_text(fb, LW, LH, 30, 165, "coat " as *u8, dim)
  lab_geno_str(fb, gx1, 165, het, GX_COAT, hot)
  var gx2: i64 = gr_text(fb, LW, LH, 124, 165, "coat " as *u8, dim)
  lab_geno_str(fb, gx2, 165, het, GX_COAT, hot)
  gr_text(fb, LW, LH, 22, 176, "rose (heterozygous)" as *u8, dim)
104
  // ---- PUNNETT SQUARE: computed by gx_punnett, drawn as the real 2x2 ----
  let sq: *i64 = sys_mmap(16 * 8) as *i64
  gx_punnett(het, het, GX_COAT, sq)
  gr_rect(fb, LW, LH, 224, 32, 420, 190, gr_pack(20,18,42))
  gr_rect(fb, LW, LH, 224, 32, 420, 34, gr_pack(70,62,120))
  gr_text(fb, LW, LH, 230, 40, "PUNNETT SQUARE (coat)" as *u8, ink)
  // header alleles
  gr_text(fb, LW, LH, 296, 56, gx_allele_char(GX_COAT, gx_a0(het, GX_COAT)), hot)
  gr_text(fb, LW, LH, 356, 56, gx_allele_char(GX_COAT, gx_a1(het, GX_COAT)), hot)
  gr_text(fb, LW, LH, 262, 84, gx_allele_char(GX_COAT, gx_a0(het, GX_COAT)), hot)
  gr_text(fb, LW, LH, 262, 132, gx_allele_char(GX_COAT, gx_a1(het, GX_COAT)), hot)
  var cell: i64 = 0
  while cell < 4 {
      let ccol: i64 = cell % 2
      let crow: i64 = cell / 2
      let bx: i64 = 284 + ccol*60
      let by: i64 = 72 + crow*48
      gr_rect(fb, LW, LH, bx, by, bx+52, by+40, gr_pack(30,27,58))
      let pa: i64 = sq[cell] / 4
      let pb2: i64 = sq[cell] % 4
      // the cell's own coat swatch = the phenotype that pairing produces
      let ph: i64 = sq[4 + cell]
      gr_rect(fb, LW, LH, bx+4, by+4, bx+48, by+20,
              gr_pack(cr_coat_rgb(ph,0), cr_coat_rgb(ph,1), cr_coat_rgb(ph,2)))
      var tx2: i64 = gr_text(fb, LW, LH, bx+8, by+24, gx_allele_char(GX_COAT, pa), ink)
      gr_text(fb, LW, LH, tx2, by+24, gx_allele_char(GX_COAT, pb2), ink)
      cell = cell + 1
  }
  let odds_rose: i64 = gx_odds_permil(het, het, GX_COAT, 1)
  let odds_crim: i64 = gx_odds_permil(het, het, GX_COAT, 2)
  let odds_white: i64 = gx_odds_permil(het, het, GX_COAT, 0)
  var ox: i64 = gr_text(fb, LW, LH, 230, 170, "predicted  crimson " as *u8, dim)
  ox = gr_num(fb, LW, LH, ox, 170, odds_crim, hot)
  ox = gr_text(fb, LW, LH, ox, 170, "  rose " as *u8, dim)
  ox = gr_num(fb, LW, LH, ox, 170, odds_rose, hot)
  ox = gr_text(fb, LW, LH, ox, 170, "  white " as *u8, dim)
  gr_num(fb, LW, LH, ox, 170, odds_white, hot)
142
  // ---- THE LITTER: actually breed 24 children and render them ----
  let bt: *i64 = sys_mmap((1 + 256*BRD_BEAT_STRIDE) * 8) as *i64
  bt[0] = 0
  let info: *i64 = sys_mmap(8 * 8) as *i64
  gr_rect(fb, LW, LH, 8, 198, 632, 392, gr_pack(18,16,38))
  gr_rect(fb, LW, LH, 8, 198, 632, 200, gr_pack(70,62,120))
  gr_text(fb, LW, LH, 14, 206, "THE LITTER -- 24 children of these two parents" as *u8, ink)
  var nC: i64 = 0
  var nR: i64 = 0
  var nW: i64 = 0
  var nHorn: i64 = 0
  var nWing: i64 = 0
  var kid: i64 = 0
  while kid < 24 {
      let ch: i64 = brd_breed_info(ros, mum, dad, 1, 90210 + kid*31, bt, 256, kid, info)
      if ch > 0 {
          let g: i64 = en_get(ros, ch, C_GENO)
          let ph: i64 = gx_phenotype(g, GX_COAT)
          if ph == 2 { nC = nC + 1 }
          if ph == 1 { nR = nR + 1 }
          if ph == 0 { nW = nW + 1 }
          if gx_phenotype(g, GX_HORN) == 1 { nHorn = nHorn + 1 }
          if gx_phenotype(g, GX_WING) == 0 { nWing = nWing + 1 }
          let col: i64 = kid % 8
          let row: i64 = kid / 8
          cr_creature_gx(fb, LW, LH, 48 + col*76, 246 + row*52, 21,
              en_get(ros, ch, C_SPEC), en_get(ros, ch, C_IV0), en_get(ros, ch, C_IV1),
              en_get(ros, ch, C_IV2), en_get(ros, ch, C_IV3),
              g, en_get(ros, ch, C_POLY), en_get(ros, ch, C_EPI))
      }
      kid = kid + 1
  }
  var sx: i64 = gr_text(fb, LW, LH, 14, 376, "observed  crimson " as *u8, dim)
  sx = gr_num(fb, LW, LH, sx, 376, nC, hot)
  sx = gr_text(fb, LW, LH, sx, 376, "  rose " as *u8, dim)
  sx = gr_num(fb, LW, LH, sx, 376, nR, hot)
  sx = gr_text(fb, LW, LH, sx, 376, "  white " as *u8, dim)
  sx = gr_num(fb, LW, LH, sx, 376, nW, hot)
  sx = gr_text(fb, LW, LH, sx, 376, "   of 24  (1:2:1 predicts 6:12:6)" as *u8, dim)
182
  // ---- the tooth: the VISIBLE ratios must track the Punnett prediction over a big sample ----
  var bC: i64 = 0
  var bR: i64 = 0
  var bW: i64 = 0
  var big: i64 = 0
  while big < 1200 {
      let ch2: i64 = brd_breed_info(ros, mum, dad, 1, 500000 + big*17, bt, 256, big, info)
      if ch2 > 0 {
          let ph2: i64 = gx_expressed(en_get(ros, ch2, C_GENO), en_get(ros, ch2, C_EPI), GX_COAT)
          if ph2 == 2 { bC = bC + 1 }
          if ph2 == 1 { bR = bR + 1 }
          if ph2 == 0 { bW = bW + 1 }
          en_destroy(ros, ch2)
      }
      big = big + 1
  }
  let tot: i64 = bC + bR + bW
  var pC: i64 = 0
  var pR: i64 = 0
  var pW2: i64 = 0
  if tot > 0 { pC = bC*1000/tot; pR = bR*1000/tot; pW2 = bW*1000/tot }
  // ★THE LAW OF LARGE NUMBERS, on screen: a 24-child litter is NOISY (this one came out 8:8:8, not
  // 6:12:6) while 1200 children land on the prediction. Showing both is the honest teaching -- a real
  // breeder's small sample deviates, and that is not the theory being wrong.
  var lx: i64 = gr_text(fb, LW, LH, 14, 386, "over 1200 children: " as *u8, dim)
  lx = gr_num(fb, LW, LH, lx, 386, pC, hot)
  lx = gr_text(fb, LW, LH, lx, 386, " / " as *u8, dim)
  lx = gr_num(fb, LW, LH, lx, 386, pR, hot)
  lx = gr_text(fb, LW, LH, lx, 386, " / " as *u8, dim)
  lx = gr_num(fb, LW, LH, lx, 386, pW2, hot)
  lx = gr_text(fb, LW, LH, lx, 386, " permil -- small samples are noisy, the law holds at scale" as *u8, dim)
214
  lab_png(fb, "knowledge/nx_breeder_lab.png" as *u8)
216
  let cout: *i64 = sys_mmap(GC_N * 8) as *i64
  gc_score(fb, LW, LH, cout)
  let q: i64 = gc_quality(cout)
220
  lp("=== nx_breeder_lab: GENETICS LAB screen emitted ===\n" as *u8)
  lp("  frame knowledge/nx_breeder_lab.png (640x400)\n" as *u8)
  lp("  VISIBLE coat phenotypes over " as *u8); ln2(tot); lp(" rendered children: crimson " as *u8); ln2(pC)
  lp(" / rose " as *u8); ln2(pR); lp(" / white " as *u8); ln2(pW2); lp(" permil\n" as *u8)
  lp("  Punnett predicts                              crimson " as *u8); ln2(odds_crim)
  lp(" / rose " as *u8); ln2(odds_rose); lp(" / white " as *u8); ln2(odds_white); lp(" permil\n" as *u8)
  lp("  litter shown: horned " as *u8); ln2(nHorn); lp("/24, winged " as *u8); ln2(nWing); lp("/24\n" as *u8)
  lp("  critic quality " as *u8); ln2(q); lp(" = " as *u8); lp(gc_verdict(q)); lnl()
  var dC: i64 = pC - odds_crim
  if dC < 0 { dC = 0 - dC }
  var dR: i64 = pR - odds_rose
  if dR < 0 { dR = 0 - dR }
  var dW: i64 = pW2 - odds_white
  if dW < 0 { dW = 0 - dW }
  var ok: i64 = 0
  if dC <= 60 { if dR <= 60 { if dW <= 60 { if tot == 1200 { ok = 1 } } } }
  if ok == 1 { lp("VERDICT GREEN -- what the player SEES matches the Punnett square (max deviation " as *u8); ln2(dC+dR+dW); lp(" permil summed)\n" as *u8); return 0 }
  lp("VERDICT RED -- rendered phenotypes drifted from the prediction\n" as *u8)
  return 1
240}
nx_breeder_lab.nx source

`nx_breeder_lab.nx` source
