From c948f4df7636aca7279e04d30d763ec87cd618f1 Mon Sep 17 00:00:00 2001
From: Claude <claude@anthropic.com>
Date: Mon, 19 Oct 2026 12:49:20 +0700
Subject: [PATCH] analog: PSP103 drain and gate charge tables (tt, ss), the
 charge-aware table solver, the 500 fF inverter within 0.02 % of ngspice
 (Closes #10880)

M3c of #10486, the charge half of #10777 (Part of #10777).

- specs/analog/decks/psp103_q_tt.cir, psp103_q_ss.cir: the drain and gate charge of sg13_lv_nmos / sg13_lv_pmos on the M3b grid, taken from the terminal currents of one backward Euler step (PSP103 reports no terminal charge), with the current tables and ngspice's own op-point capacitances printed in the same run.
- mos_psp103_data.t27: TT_QD_N TT_QG_N TT_QD_P TT_QG_P and the SS_ four, the 401-sample ngspice backward Euler waveform at 500 fF (INV500_BE, printed by psp103_inv_be_500f.cir), pins.
- mos_table.t27: the solver generalized with the two drain-charge tables and the previous step's gate (the *_q functions); the current-only functions are those with zero charge tables, nothing is copied, M3a's pinned run is untouched. Three tests: zero charge is the old path, a plane of charge gives the closed form of the gate-edge kick, the mirror holds with four tables swapped.
- mos_psp103.t27: tt inverter at 500 fF and dt = 100 ps against ngspice backward Euler on the same grid, tpLH 0.020 % and tpHL 0.0041 % (3.09 % and 3.63 % without charge); the three errors apart; the ss charge tables in the tt run fail the check; the table slopes agree with ngspice's op-point capacitances within 0.3 %; the weak-inversion 1 % bound stated with its reasons (follow-up #10881). Corners ff, fs, sf: #10882.

Part of #10486

Co-Authored-By: Claude Sonnet 5.5 <noreply@anthropic.com>
---
 .trinity/seals/analog_AnalogMosPsp103.json    |  20 +-
 .../seals/analog_AnalogMosPsp103Data.json     |  16 +-
 .trinity/seals/analog_AnalogMosTable.json     |  16 +-
 .../seals/analog_AnalogSarComparator.json     |   4 +-
 .trinity/seals/analog_AnalogSarSwitch.json    |   4 +-
 specs/analog/decks/psp103_inv_be_500f.cir     |  11 +
 specs/analog/decks/psp103_q_ss.cir            | 140 +++++
 specs/analog/decks/psp103_q_tt.cir            | 140 +++++
 specs/analog/mos_psp103.t27                   | 414 ++++++++++++-
 specs/analog/mos_psp103_data.t27              | 555 +++++++++++++++++-
 specs/analog/mos_table.t27                    | 268 ++++++++-
 11 files changed, 1517 insertions(+), 71 deletions(-)
 create mode 100644 specs/analog/decks/psp103_q_ss.cir
 create mode 100644 specs/analog/decks/psp103_q_tt.cir

diff --git a/.trinity/seals/analog_AnalogMosPsp103.json b/.trinity/seals/analog_AnalogMosPsp103.json
index 147dcf14ba..d265286ba3 100644
--- a/.trinity/seals/analog_AnalogMosPsp103.json
+++ b/.trinity/seals/analog_AnalogMosPsp103.json
@@ -7,31 +7,31 @@
     },
     {
       "spec": "specs/analog/mos_psp103_data.t27",
-      "spec_hash": "sha256:d5f3244ff1936f2b96e95b06e022d6b42b489ac49e025ec99486d6d36e6c6580"
+      "spec_hash": "sha256:be100dec66ee80de6cf520f79038e7668bc1171102db1fefe6a4c74559d029de"
     },
     {
       "spec": "specs/analog/mos_table.t27",
-      "spec_hash": "sha256:9c0d1537ed5d8d60c4bdcb67d1e4c73ef80270e8cab224cb089c5b568fe11182"
+      "spec_hash": "sha256:16d28f10fa8663926dddb4bcced981cdc4cd8f38a6b748e4b7688f2e56bbab6c"
     }
   ],
   "config": "gen=c,rust,verilog,zig;tests=zig",
-  "gen_hash_c": "sha256:d7a49c643703a544723718186df74850018c4bf067fdb06ca9c860d775e069e8",
-  "gen_hash_rust": "sha256:32ee0d51345b936e45ff10956b90fdcafb6106bca46c2700de2bf43bf4407119",
-  "gen_hash_verilog": "sha256:0e3d19e85e67353129e8213c6135a5d1b518de4d7502c2d92512695b92de832e",
-  "gen_hash_zig": "sha256:445a5f05e794204ad3e4c198735c6fb1f4ab816e26fd36d3c408bba83168164f",
+  "gen_hash_c": "sha256:a5aeefd7abe3aee47919f32153161040121a31fba6fbfc6e3517fe1796191a8a",
+  "gen_hash_rust": "sha256:0660934f0f63c77f8b204796f1c8e3f5db1d35603797d7ddf912f66a2039d4a2",
+  "gen_hash_verilog": "sha256:06e55f455e71a7b274d2b8e766028a05ec666c1814b4ed64f8e7ad381706e61e",
+  "gen_hash_zig": "sha256:460548424dd136ec773917ebc3dd3c6ad31561cbde36cabf8a69266c2e8c7b6c",
   "module": "AnalogMosPsp103",
   "ring": 12,
   "seal_schema": 2,
-  "sealed_at": "2026-10-19T02:10:22Z",
+  "sealed_at": "2026-10-19T05:34:01Z",
   "sealed_by": "t27c-bootstrap@0.5.2",
-  "spec_hash": "sha256:c5bd55a847953654fbef792dcfc61da7f6434b8023341df9ab4a53e5d082bbe0",
+  "spec_hash": "sha256:811617266484683542dd1723db77687ad7fdf7f5a665953f5e6a920a7bdf0118",
   "spec_path": "specs/analog/mos_psp103.t27",
   "tests": {
     "failed": 0,
     "failing": [],
     "forced": false,
-    "passed": 11,
-    "total": 11,
+    "passed": 17,
+    "total": 17,
     "vacuous": 0
   },
   "toolchain": {
diff --git a/.trinity/seals/analog_AnalogMosPsp103Data.json b/.trinity/seals/analog_AnalogMosPsp103Data.json
index 06125dea5e..d1dac009ec 100644
--- a/.trinity/seals/analog_AnalogMosPsp103Data.json
+++ b/.trinity/seals/analog_AnalogMosPsp103Data.json
@@ -7,23 +7,23 @@
     }
   ],
   "config": "gen=c,rust,verilog,zig;tests=zig",
-  "gen_hash_c": "sha256:d0135ebcf91901c442dc8a0c3e3a431f50a57cbc563d430a34e45e68392c9dfe",
-  "gen_hash_rust": "sha256:9dad368c2549f028ae0828213c7231a0e8ff613988ff42ce7d3aac6fcf9c4e52",
-  "gen_hash_verilog": "sha256:52f448d5fe93db1de3faa5fa13cb881486ec67a2ed17d1db475cd4d8d222ac82",
-  "gen_hash_zig": "sha256:1bdd4c520cb0864ffa3c4d2f7ce80418e7a0752e37610b81b4a000c528c413f9",
+  "gen_hash_c": "sha256:be16cbf1fbac3021d77530b1e86f19bfacf626eb7cfe834900c984331491c060",
+  "gen_hash_rust": "sha256:cbd4b0ad5ced378295e73bec1bcc69ff8270b99f4a4e27064a6aa2b0aa02c222",
+  "gen_hash_verilog": "sha256:f062ed3fc219deec9fa254d48b0a8e3ac8de0d2c525810154d1f25ab20a95162",
+  "gen_hash_zig": "sha256:03d50a36f22f40daf936ade0b764ae11b4142c3f667bf1ab640007eb4fdbcb45",
   "module": "AnalogMosPsp103Data",
   "ring": 12,
   "seal_schema": 2,
-  "sealed_at": "2026-10-19T02:10:08Z",
+  "sealed_at": "2026-10-19T05:33:43Z",
   "sealed_by": "t27c-bootstrap@0.5.2",
-  "spec_hash": "sha256:d5f3244ff1936f2b96e95b06e022d6b42b489ac49e025ec99486d6d36e6c6580",
+  "spec_hash": "sha256:be100dec66ee80de6cf520f79038e7668bc1171102db1fefe6a4c74559d029de",
   "spec_path": "specs/analog/mos_psp103_data.t27",
   "tests": {
     "failed": 0,
     "failing": [],
     "forced": false,
-    "passed": 4,
-    "total": 4,
+    "passed": 7,
+    "total": 7,
     "vacuous": 0
   },
   "toolchain": {
diff --git a/.trinity/seals/analog_AnalogMosTable.json b/.trinity/seals/analog_AnalogMosTable.json
index 5b1d2673a8..de7669cec8 100644
--- a/.trinity/seals/analog_AnalogMosTable.json
+++ b/.trinity/seals/analog_AnalogMosTable.json
@@ -7,23 +7,23 @@
     }
   ],
   "config": "gen=c,rust,verilog,zig;tests=zig",
-  "gen_hash_c": "sha256:3de4432818944eab74fae03952b4b4cffd104b00b9cda80dd73bb6ed6e819f0c",
-  "gen_hash_rust": "sha256:775c36a7f0370d1bc3dd0752f72f435e3fbc6e4972ad9b5c77efc35fb2784750",
-  "gen_hash_verilog": "sha256:0f3b42817e930cb80890ec45f5ecb5ecdf2a73bb764d029149ab8519dfd43dfc",
-  "gen_hash_zig": "sha256:05393568aa75f13169560fb6fd1118d46141ae46931458fc181ac12ba42ca196",
+  "gen_hash_c": "sha256:e0ca09255f83209260afff210947c8af851f5412f9aa7a40f7f5b148a29ac739",
+  "gen_hash_rust": "sha256:dd4d433feaf2d937fc3acb769f6a9a8febc6d7cc58363dbc061ad9922eba3eb6",
+  "gen_hash_verilog": "sha256:152706909b6dd98a9fda6242ef666dfa85a03c0a6e4d37564088d2e43da3aafa",
+  "gen_hash_zig": "sha256:652d1bfbd355d08e8eca4a3f13a4957322006510c310d796111db8aac65ef1d0",
   "module": "AnalogMosTable",
   "ring": 12,
   "seal_schema": 2,
-  "sealed_at": "2026-10-19T02:09:56Z",
+  "sealed_at": "2026-10-19T05:33:16Z",
   "sealed_by": "t27c-bootstrap@0.5.2",
-  "spec_hash": "sha256:9c0d1537ed5d8d60c4bdcb67d1e4c73ef80270e8cab224cb089c5b568fe11182",
+  "spec_hash": "sha256:16d28f10fa8663926dddb4bcced981cdc4cd8f38a6b748e4b7688f2e56bbab6c",
   "spec_path": "specs/analog/mos_table.t27",
   "tests": {
     "failed": 0,
     "failing": [],
     "forced": false,
-    "passed": 14,
-    "total": 14,
+    "passed": 17,
+    "total": 17,
     "vacuous": 0
   },
   "toolchain": {
diff --git a/.trinity/seals/analog_AnalogSarComparator.json b/.trinity/seals/analog_AnalogSarComparator.json
index df7ea1875f..5a33e5f807 100644
--- a/.trinity/seals/analog_AnalogSarComparator.json
+++ b/.trinity/seals/analog_AnalogSarComparator.json
@@ -7,7 +7,7 @@
     },
     {
       "spec": "specs/analog/mos_table.t27",
-      "spec_hash": "sha256:9c0d1537ed5d8d60c4bdcb67d1e4c73ef80270e8cab224cb089c5b568fe11182"
+      "spec_hash": "sha256:16d28f10fa8663926dddb4bcced981cdc4cd8f38a6b748e4b7688f2e56bbab6c"
     }
   ],
   "config": "gen=c,rust,verilog,zig;tests=zig",
@@ -18,7 +18,7 @@
   "module": "AnalogSarComparator",
   "ring": 12,
   "seal_schema": 2,
-  "sealed_at": "2026-10-19T02:37:52Z",
+  "sealed_at": "2026-10-19T05:38:44Z",
   "sealed_by": "t27c-bootstrap@0.5.2",
   "spec_hash": "sha256:b88d05d5d17937024580ab40e80b698b63962cab531164d4a94a5dddd05a72a8",
   "spec_path": "specs/analog/sar_comparator.t27",
diff --git a/.trinity/seals/analog_AnalogSarSwitch.json b/.trinity/seals/analog_AnalogSarSwitch.json
index 3f8d0787ed..623e9ac7f9 100644
--- a/.trinity/seals/analog_AnalogSarSwitch.json
+++ b/.trinity/seals/analog_AnalogSarSwitch.json
@@ -7,7 +7,7 @@
     },
     {
       "spec": "specs/analog/mos_table.t27",
-      "spec_hash": "sha256:9c0d1537ed5d8d60c4bdcb67d1e4c73ef80270e8cab224cb089c5b568fe11182"
+      "spec_hash": "sha256:16d28f10fa8663926dddb4bcced981cdc4cd8f38a6b748e4b7688f2e56bbab6c"
     },
     {
       "spec": "specs/analog/sar_cdac.t27",
@@ -26,7 +26,7 @@
   "module": "AnalogSarSwitch",
   "ring": 12,
   "seal_schema": 2,
-  "sealed_at": "2026-10-19T03:59:52Z",
+  "sealed_at": "2026-10-19T05:48:24Z",
   "sealed_by": "t27c-bootstrap@0.5.2",
   "spec_hash": "sha256:122726502a51cb4c8980674b819ee3251bf8ee69a5586d6cd44a5eafcf284cb4",
   "spec_path": "specs/analog/sar_switch.t27",
diff --git a/specs/analog/decks/psp103_inv_be_500f.cir b/specs/analog/decks/psp103_inv_be_500f.cir
index ddb03914c7..104c40c34a 100644
--- a/specs/analog/decks/psp103_inv_be_500f.cir
+++ b/specs/analog/decks/psp103_inv_be_500f.cir
@@ -11,6 +11,8 @@ IHP SG13G2 inverter, 500 fF, backward Euler on a fixed 100 ps grid, corner mos_t
 * at (k + 2) dt; step 0 starts at t = dt = 100 ps and step 200 at 201 dt = 20.1 ns.
 * What it prints (ngspice .meas, 6 digits, in s): tplh and tphl, the times at which the output crosses 0.75 V rising (after step 0's edge) and
 * falling (after step 200's edge); the delays are these minus 100 ps and minus 20.1 ns.
+* M3c adds the waveform: between "BEGIN tran_be" and "END tran_be", 401 rows "<index> TAB <time> TAB <y>", y the output in units of 2^-24 V at
+* t = 100 ps, 200 ps, ... 40.1 ns (row 0 is the circuit at rest before step 0, row k + 1 the output after step k), as psp103_inv_be.cir does at 2 ns.
 .lib $IHP_PDK/ihp-sg13g2/libs.tech/ngspice/models/cornerMOSlv.lib mos_tt
 VDD vdd 0 1.5
 Bin in 0 V = 1.5*(1-min(max(nint(time/1e-10)-1,0),1)+min(max(nint(time/1e-10)-201,0),1))
@@ -26,5 +28,14 @@ CL out 0 500f
 pre_osdi $IHP_OSDI/psp103.osdi
 pre_osdi $IHP_OSDI/psp103_nqs.osdi
 run
+* M3c (#10777): after the run, which prints the two .meas numbers as before, the waveform is printed the way psp103_inv_be.cir prints it.
+set numdgt=10
+set noprintscale
+let y = nint(v(out) * 16777216)
+let ts = time[8,408]
+let ys = y[8,408]
+echo BEGIN tran_be
+print ts ys
+echo END tran_be
 .endc
 .end
diff --git a/specs/analog/decks/psp103_q_ss.cir b/specs/analog/decks/psp103_q_ss.cir
new file mode 100644
index 0000000000..cd4971863c
--- /dev/null
+++ b/specs/analog/decks/psp103_q_ss.cir
@@ -0,0 +1,140 @@
+IHP SG13G2 PSP103 drain and gate charge tables, corner mos_ss
+* specs/analog/decks/psp103_q_ss.cir -- charge table data for specs/analog/mos_psp103_data.t27 (M3c, #10777, epic #10486). A netlist plus a
+* .control block that prints the numbers. The same deck exists for corner mos_tt (psp103_q_tt.cir); the two differ in the title and the .lib section.
+* The IHP corner file is read at run time and never copied here.
+* Model: IHP-Open-PDK commit 21133f219b1a6da10d803e345761e77415ceff40, ihp-sg13g2/libs.tech/ngspice/models/cornerMOSlv.lib (PSP 103.6,
+* revision 200310), section mos_ss; the model code is libs.tech/verilog-a/psp103/psp103.va of the same commit, compiled with
+* OpenVAF-Reloaded v24.0.2mob to psp103.osdi (and psp103_nqs.osdi). Temperature 27 C (TNOM 27 C). Devices: sg13_lv_nmos w=1.0u l=0.45u
+* ng=1 m=1 and sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1, exactly the devices of psp103_tab_ss.cir (the same netlist lines), on the same
+* 25 x 25 grid of 62.5 mV over 0 .. 1.5 V.
+* Solver options: reltol 1e-6, vntol 1e-9, abstol 1e-14, as psp103_tab_ss.cir.
+* Paths are NOT in this file. ngspice reads two environment variables: IHP_PDK, the IHP-Open-PDK checkout, and IHP_OSDI, the directory that
+* holds the compiled psp103.osdi and psp103_nqs.osdi.
+* Run: IHP_PDK=<checkout> IHP_OSDI=<osdi dir> ngspice -b psp103_q_ss.cir > psp103_q_ss.out    (ngspice-47)
+* WHICH VARIABLES. PSP103 reports no terminal charge as an output variable: `show` on the OSDI instance lists currents (ids, ige, ...), conductances
+* (gm, gmb, gds) and capacitances (cdd, cdg, cgg, cgd, ..., plus the overlap cgsol, cgdol and the junction cjs, cjd, each reported apart), and
+* no q*. The charge is therefore taken from the terminal CURRENTS of a backward Euler step, which is what ngspice itself integrates: for a step
+* of dt from state A to state B the current into a terminal is i(B) = Idc(B) + (Q(B) - Q(A)) / dt, with Q the terminal's whole charge
+* (channel, overlap, fringe and junction parts together), so Q(B) - Q(A) = (i(B) - Idc(B)) dt. The drive is a staircase: node k of the grid
+* (ix = k / 25 along Vgs, iy = k mod 25 along Vds) is visited from the origin A = (0, 0) and the bias goes back to A on the next step. The
+* return step is the one measured: its end state is the origin, where the device is off and Idc(A) is the current of the first step, so the
+* subtraction is of a current of picoamps, not of the 400 uA of a node in strong inversion. Then
+*     Q(B) - Q(A) = (i_A - i_return) dt        with i_A the terminal current at the origin before the first jump and i_return at the end of the return step.
+* The step is dt = 10 ns, long against the gate resistance of the model (about 29 ohm times a few fF is 0.1 ps), so the gate network has settled
+* and the charge is the quasi-static one. (At dt = 100 ps the settling is 1 part in 1000 and the outgoing step, measured against the DC current of the
+* node, shows it as a current deficit of tau gm dV, which is why the outgoing step is not used.)
+* What it prints. Between "BEGIN name" and "END name" every row "<index> TAB <value>" is one number, in sweep order (index = ix * 25 + iy). The
+* values are pasted into the spec in that order, none changed.
+*   grid_n, grid_p  Id of the nmos at Vgs = ix * 62.5 mV, Vds = iy * 62.5 mV and of the pmos at Vsg, Vsd, exactly as psp103_tab_ss.cir prints them
+*           (units of 2^-24 V, dt/C = 200 V/A): the current tables and the charge tables come from one run of one netlist. Negative values held at 0.
+*   op_n    cdd, cdg, cgg, cgsol, cgdol, cjd of the nmos at Vgs = 0.9 V, Vds = 0.75 V in aF, as ngspice reports them: the channel capacitances
+*           (cdd = dQd/dVd, cdg = dQd/dVg, cgg = dQg/dVg of the intrinsic charges), the gate-source and gate-drain overlaps and the drain junction, each apart.
+*   qd_n    the DRAIN charge of the nmos. Unit: Q / 500 fF in units of 2^-24 V, that is the voltage a 500 fF node moves when it takes the charge, so
+*           one unit is 500 fF * 59.6 nV = 2.98e-20 C (0.19 electron). Entry = nint(2e12 * 2^24 * (Q(B) - Q(A))) + 2^21: the charge relative to the
+*           origin Vgs = Vds = 0, plus 2^21 so that no entry is negative. The offset is arbitrary; only differences are physical.
+*   qg_n    the GATE charge of the nmos, the same unit and offset.
+*   qd_p    the DRAIN charge of the pmos at Vsg = ix * 62.5 mV, Vsd = iy * 62.5 mV (source and bulk on the 1.5 V rail), NEGATED: entry =
+*           2^21 - nint(2e12 * 2^24 * (Q(B) - Q(A))), so that it grows with Vsd as the nmos' does (capacitances are positive: the drain charge grows with
+*           the drain voltage, which falls as Vsd grows).
+*   qg_p    the GATE charge of the pmos, negated the same way (grows with Vsg).
+*   tchk    the largest distance of an ngspice time point used from the grid k * 10 ns, in s: 0 when every row used is a step end.
+* nmos tran device: drain and gate on B sources through zero-volt current meters; the source and bulk at 0 V
+.lib $IHP_PDK/ihp-sg13g2/libs.tech/ngspice/models/cornerMOSlv.lib mos_ss
+.func kst() {nint(time/1e-8)}
+.func nd() {max(floor(kst()/2)-1,0)}
+.func act() {(time<1e-3)*(kst()>=2)*((kst()-2*floor(kst()/2))==0)}
+.func ixf() {floor(nd()/25)}
+.func iyf() {nd()-25*floor(nd()/25)}
+* (act() is 0 in a DC sweep: ngspice puts the swept value into `time`, which is at least 62.5 mV there and below 13 us in the transient.)
+BDN tdn 0 V = act()*0.0625*iyf()
+BGN tgn 0 V = act()*0.0625*ixf()
+VMDN tdn tdnn 0
+VMGN tgn tgnn 0
+XTN tdnn tgnn 0 0 sg13_lv_nmos w=1.0u l=0.45u ng=1 m=1
+* pmos tran device: source and bulk at the 1.5 V rail; the origin is Vsg = Vsd = 0
+VTSP tsp 0 1.5
+BDP tdp 0 V = 1.5-act()*0.0625*iyf()
+BGP tgp 0 V = 1.5-act()*0.0625*ixf()
+VMDP tdp tdpp 0
+VMGP tgp tgpp 0
+XTP tdpp tgpp tsp tsp sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1
+* nmos table device: drain VDN, gate VGN, source and bulk at 0 V (as psp103_tab_ss.cir)
+VDN dn 0 0
+VGN gn 0 0
+XMN dn gn 0 0 sg13_lv_nmos w=1.0u l=0.45u ng=1 m=1
+* pmos table device: source and bulk at the 1.5 V rail, drain VDP and gate VGP, both swept downward from 1.5 V
+VSP sp 0 1.5
+VDP dp 0 1.5
+VGP gp 0 1.5
+XMP dp gp sp sp sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1
+.temp 27
+.options reltol=1e-6 vntol=1e-9 abstol=1e-14
+.options method=gear maxord=1 trtol=1e7 itl4=200
+.control
+pre_osdi $IHP_OSDI/psp103.osdi
+pre_osdi $IHP_OSDI/psp103_nqs.osdi
+set numdgt=10
+set noprintscale
+echo BEGIN grid_n
+dc VDN 0 1.5 0.0625 VGN 0 1.5 0.0625
+let u = nint(-i(VDN) * 200 * 16777216)
+let u = (u + abs(u)) / 2
+print u
+echo END grid_n
+echo BEGIN grid_p
+dc VDP 1.5 0 -0.0625 VGP 1.5 0 -0.0625
+let u = nint(i(VDP) * 200 * 16777216)
+let u = (u + abs(u)) / 2
+print u
+echo END grid_p
+* the capacitances ngspice itself reports at one operating point of the nmos, Vgs = 0.9 V, Vds = 0.75 V, in aF (nint(C * 1e18)): the cross-check of
+* the charge tables, which are made from currents and never see these
+alter VDN = 0.75
+alter VGN = 0.9
+op
+let c_dd = nint(@n.xmn.nsg13_lv_nmos[cdd] * 1e18)
+let c_dg = nint(@n.xmn.nsg13_lv_nmos[cdg] * 1e18)
+let c_gg = nint(@n.xmn.nsg13_lv_nmos[cgg] * 1e18)
+let c_gsol = nint(@n.xmn.nsg13_lv_nmos[cgsol] * 1e18)
+let c_gdol = nint(@n.xmn.nsg13_lv_nmos[cgdol] * 1e18)
+let c_jd = nint(@n.xmn.nsg13_lv_nmos[cjd] * 1e18)
+echo BEGIN op_n
+print c_dd c_dg c_gg c_gsol c_gdol c_jd
+echo END op_n
+tran 7.8125n 12.54u 0 10n
+let dn0 = i(vmdn)[8]
+let gn0 = i(vmgn)[8]
+let dp0 = i(vmdp)[8]
+let gp0 = i(vmgp)[8]
+let qdn = unitvec(625)
+let qgn = unitvec(625)
+let qdp = unitvec(625)
+let qgp = unitvec(625)
+let tchk = abs(time[8] - 1e-8)
+let k = 0
+while k < 625
+  let r = 2*k + 10
+  let tchk = max(tchk, abs(time[r] - (2*k+3)*1e-8))
+  let qdn[k] = nint((dn0 - i(vmdn)[r]) * 1e-8 * 2e12 * 16777216) + 2097152
+  let qgn[k] = nint((gn0 - i(vmgn)[r]) * 1e-8 * 2e12 * 16777216) + 2097152
+  let qdp[k] = 2097152 - nint((dp0 - i(vmdp)[r]) * 1e-8 * 2e12 * 16777216)
+  let qgp[k] = 2097152 - nint((gp0 - i(vmgp)[r]) * 1e-8 * 2e12 * 16777216)
+  let k = k + 1
+end
+echo BEGIN qd_n
+print qdn
+echo END qd_n
+echo BEGIN qg_n
+print qgn
+echo END qg_n
+echo BEGIN qd_p
+print qdp
+echo END qd_p
+echo BEGIN qg_p
+print qgp
+echo END qg_p
+echo BEGIN tchk
+print tchk
+echo END tchk
+.endc
+.end
diff --git a/specs/analog/decks/psp103_q_tt.cir b/specs/analog/decks/psp103_q_tt.cir
new file mode 100644
index 0000000000..552bbafe6c
--- /dev/null
+++ b/specs/analog/decks/psp103_q_tt.cir
@@ -0,0 +1,140 @@
+IHP SG13G2 PSP103 drain and gate charge tables, corner mos_tt
+* specs/analog/decks/psp103_q_tt.cir -- charge table data for specs/analog/mos_psp103_data.t27 (M3c, #10777, epic #10486). A netlist plus a
+* .control block that prints the numbers. The same deck exists for corner mos_ss (psp103_q_ss.cir); the two differ in the title and the .lib section.
+* The IHP corner file is read at run time and never copied here.
+* Model: IHP-Open-PDK commit 21133f219b1a6da10d803e345761e77415ceff40, ihp-sg13g2/libs.tech/ngspice/models/cornerMOSlv.lib (PSP 103.6,
+* revision 200310), section mos_tt; the model code is libs.tech/verilog-a/psp103/psp103.va of the same commit, compiled with
+* OpenVAF-Reloaded v24.0.2mob to psp103.osdi (and psp103_nqs.osdi). Temperature 27 C (TNOM 27 C). Devices: sg13_lv_nmos w=1.0u l=0.45u
+* ng=1 m=1 and sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1, exactly the devices of psp103_tab_tt.cir (the same netlist lines), on the same
+* 25 x 25 grid of 62.5 mV over 0 .. 1.5 V.
+* Solver options: reltol 1e-6, vntol 1e-9, abstol 1e-14, as psp103_tab_tt.cir.
+* Paths are NOT in this file. ngspice reads two environment variables: IHP_PDK, the IHP-Open-PDK checkout, and IHP_OSDI, the directory that
+* holds the compiled psp103.osdi and psp103_nqs.osdi.
+* Run: IHP_PDK=<checkout> IHP_OSDI=<osdi dir> ngspice -b psp103_q_tt.cir > psp103_q_tt.out    (ngspice-47)
+* WHICH VARIABLES. PSP103 reports no terminal charge as an output variable: `show` on the OSDI instance lists currents (ids, ige, ...), conductances
+* (gm, gmb, gds) and capacitances (cdd, cdg, cgg, cgd, ..., plus the overlap cgsol, cgdol and the junction cjs, cjd, each reported apart), and
+* no q*. The charge is therefore taken from the terminal CURRENTS of a backward Euler step, which is what ngspice itself integrates: for a step
+* of dt from state A to state B the current into a terminal is i(B) = Idc(B) + (Q(B) - Q(A)) / dt, with Q the terminal's whole charge
+* (channel, overlap, fringe and junction parts together), so Q(B) - Q(A) = (i(B) - Idc(B)) dt. The drive is a staircase: node k of the grid
+* (ix = k / 25 along Vgs, iy = k mod 25 along Vds) is visited from the origin A = (0, 0) and the bias goes back to A on the next step. The
+* return step is the one measured: its end state is the origin, where the device is off and Idc(A) is the current of the first step, so the
+* subtraction is of a current of picoamps, not of the 400 uA of a node in strong inversion. Then
+*     Q(B) - Q(A) = (i_A - i_return) dt        with i_A the terminal current at the origin before the first jump and i_return at the end of the return step.
+* The step is dt = 10 ns, long against the gate resistance of the model (about 29 ohm times a few fF is 0.1 ps), so the gate network has settled
+* and the charge is the quasi-static one. (At dt = 100 ps the settling is 1 part in 1000 and the outgoing step, measured against the DC current of the
+* node, shows it as a current deficit of tau gm dV, which is why the outgoing step is not used.)
+* What it prints. Between "BEGIN name" and "END name" every row "<index> TAB <value>" is one number, in sweep order (index = ix * 25 + iy). The
+* values are pasted into the spec in that order, none changed.
+*   grid_n, grid_p  Id of the nmos at Vgs = ix * 62.5 mV, Vds = iy * 62.5 mV and of the pmos at Vsg, Vsd, exactly as psp103_tab_tt.cir prints them
+*           (units of 2^-24 V, dt/C = 200 V/A): the current tables and the charge tables come from one run of one netlist. Negative values held at 0.
+*   op_n    cdd, cdg, cgg, cgsol, cgdol, cjd of the nmos at Vgs = 0.9 V, Vds = 0.75 V in aF, as ngspice reports them: the channel capacitances
+*           (cdd = dQd/dVd, cdg = dQd/dVg, cgg = dQg/dVg of the intrinsic charges), the gate-source and gate-drain overlaps and the drain junction, each apart.
+*   qd_n    the DRAIN charge of the nmos. Unit: Q / 500 fF in units of 2^-24 V, that is the voltage a 500 fF node moves when it takes the charge, so
+*           one unit is 500 fF * 59.6 nV = 2.98e-20 C (0.19 electron). Entry = nint(2e12 * 2^24 * (Q(B) - Q(A))) + 2^21: the charge relative to the
+*           origin Vgs = Vds = 0, plus 2^21 so that no entry is negative. The offset is arbitrary; only differences are physical.
+*   qg_n    the GATE charge of the nmos, the same unit and offset.
+*   qd_p    the DRAIN charge of the pmos at Vsg = ix * 62.5 mV, Vsd = iy * 62.5 mV (source and bulk on the 1.5 V rail), NEGATED: entry =
+*           2^21 - nint(2e12 * 2^24 * (Q(B) - Q(A))), so that it grows with Vsd as the nmos' does (capacitances are positive: the drain charge grows with
+*           the drain voltage, which falls as Vsd grows).
+*   qg_p    the GATE charge of the pmos, negated the same way (grows with Vsg).
+*   tchk    the largest distance of an ngspice time point used from the grid k * 10 ns, in s: 0 when every row used is a step end.
+* nmos tran device: drain and gate on B sources through zero-volt current meters; the source and bulk at 0 V
+.lib $IHP_PDK/ihp-sg13g2/libs.tech/ngspice/models/cornerMOSlv.lib mos_tt
+.func kst() {nint(time/1e-8)}
+.func nd() {max(floor(kst()/2)-1,0)}
+.func act() {(time<1e-3)*(kst()>=2)*((kst()-2*floor(kst()/2))==0)}
+.func ixf() {floor(nd()/25)}
+.func iyf() {nd()-25*floor(nd()/25)}
+* (act() is 0 in a DC sweep: ngspice puts the swept value into `time`, which is at least 62.5 mV there and below 13 us in the transient.)
+BDN tdn 0 V = act()*0.0625*iyf()
+BGN tgn 0 V = act()*0.0625*ixf()
+VMDN tdn tdnn 0
+VMGN tgn tgnn 0
+XTN tdnn tgnn 0 0 sg13_lv_nmos w=1.0u l=0.45u ng=1 m=1
+* pmos tran device: source and bulk at the 1.5 V rail; the origin is Vsg = Vsd = 0
+VTSP tsp 0 1.5
+BDP tdp 0 V = 1.5-act()*0.0625*iyf()
+BGP tgp 0 V = 1.5-act()*0.0625*ixf()
+VMDP tdp tdpp 0
+VMGP tgp tgpp 0
+XTP tdpp tgpp tsp tsp sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1
+* nmos table device: drain VDN, gate VGN, source and bulk at 0 V (as psp103_tab_tt.cir)
+VDN dn 0 0
+VGN gn 0 0
+XMN dn gn 0 0 sg13_lv_nmos w=1.0u l=0.45u ng=1 m=1
+* pmos table device: source and bulk at the 1.5 V rail, drain VDP and gate VGP, both swept downward from 1.5 V
+VSP sp 0 1.5
+VDP dp 0 1.5
+VGP gp 0 1.5
+XMP dp gp sp sp sg13_lv_pmos w=3.0u l=0.45u ng=1 m=1
+.temp 27
+.options reltol=1e-6 vntol=1e-9 abstol=1e-14
+.options method=gear maxord=1 trtol=1e7 itl4=200
+.control
+pre_osdi $IHP_OSDI/psp103.osdi
+pre_osdi $IHP_OSDI/psp103_nqs.osdi
+set numdgt=10
+set noprintscale
+echo BEGIN grid_n
+dc VDN 0 1.5 0.0625 VGN 0 1.5 0.0625
+let u = nint(-i(VDN) * 200 * 16777216)
+let u = (u + abs(u)) / 2
+print u
+echo END grid_n
+echo BEGIN grid_p
+dc VDP 1.5 0 -0.0625 VGP 1.5 0 -0.0625
+let u = nint(i(VDP) * 200 * 16777216)
+let u = (u + abs(u)) / 2
+print u
+echo END grid_p
+* the capacitances ngspice itself reports at one operating point of the nmos, Vgs = 0.9 V, Vds = 0.75 V, in aF (nint(C * 1e18)): the cross-check of
+* the charge tables, which are made from currents and never see these
+alter VDN = 0.75
+alter VGN = 0.9
+op
+let c_dd = nint(@n.xmn.nsg13_lv_nmos[cdd] * 1e18)
+let c_dg = nint(@n.xmn.nsg13_lv_nmos[cdg] * 1e18)
+let c_gg = nint(@n.xmn.nsg13_lv_nmos[cgg] * 1e18)
+let c_gsol = nint(@n.xmn.nsg13_lv_nmos[cgsol] * 1e18)
+let c_gdol = nint(@n.xmn.nsg13_lv_nmos[cgdol] * 1e18)
+let c_jd = nint(@n.xmn.nsg13_lv_nmos[cjd] * 1e18)
+echo BEGIN op_n
+print c_dd c_dg c_gg c_gsol c_gdol c_jd
+echo END op_n
+tran 7.8125n 12.54u 0 10n
+let dn0 = i(vmdn)[8]
+let gn0 = i(vmgn)[8]
+let dp0 = i(vmdp)[8]
+let gp0 = i(vmgp)[8]
+let qdn = unitvec(625)
+let qgn = unitvec(625)
+let qdp = unitvec(625)
+let qgp = unitvec(625)
+let tchk = abs(time[8] - 1e-8)
+let k = 0
+while k < 625
+  let r = 2*k + 10
+  let tchk = max(tchk, abs(time[r] - (2*k+3)*1e-8))
+  let qdn[k] = nint((dn0 - i(vmdn)[r]) * 1e-8 * 2e12 * 16777216) + 2097152
+  let qgn[k] = nint((gn0 - i(vmgn)[r]) * 1e-8 * 2e12 * 16777216) + 2097152
+  let qdp[k] = 2097152 - nint((dp0 - i(vmdp)[r]) * 1e-8 * 2e12 * 16777216)
+  let qgp[k] = 2097152 - nint((gp0 - i(vmgp)[r]) * 1e-8 * 2e12 * 16777216)
+  let k = k + 1
+end
+echo BEGIN qd_n
+print qdn
+echo END qd_n
+echo BEGIN qg_n
+print qgn
+echo END qg_n
+echo BEGIN qd_p
+print qdp
+echo END qd_p
+echo BEGIN qg_p
+print qgp
+echo END qg_p
+echo BEGIN tchk
+print tchk
+echo END tchk
+.endc
+.end
diff --git a/specs/analog/mos_psp103.t27 b/specs/analog/mos_psp103.t27
index 2f76eeb226..79da6016d4 100644
--- a/specs/analog/mos_psp103.t27
+++ b/specs/analog/mos_psp103.t27
@@ -1,6 +1,6 @@
 // SPDX-License-Identifier: Apache-2.0
 ; specs/analog/mos_psp103.t27 -- the CMOS inverter on the IHP SG13G2 PSP103 tables, per corner, judged by ngspice
-; M3b of the analog emulation epic (#10486, issue #10516). M3a (specs/analog/mos_table.t27) solved a CMOS inverter whose MOS currents
+; M3b of the analog emulation epic (#10486, issue #10516); M3c (#10777, sub-issue #10880) adds the terminal charge, THE CHARGE below. M3a (specs/analog/mos_table.t27) solved a CMOS inverter whose MOS currents
 ; were read from a table made of the square law. This file runs the same solver, unchanged but for taking the NMOS table and the
 ; PMOS table as two arguments, on tables made from the foundry's model: Id(Vgs, Vds) of sg13_lv_nmos (W 1.0 um) and sg13_lv_pmos
 ; (W 3.0 um), L 0.45 um, PSP 103.6 through ngspice-47 and the OSDI interface, IHP-Open-PDK commit 21133f219b1a6da10d803e345761e77415ceff40,
@@ -58,13 +58,43 @@
 ; THE MISSING CHARGE. At 500 fF and dt = 100 ps the model is 3.09 % and 3.63 % faster than ngspice backward Euler: PSP103's device
 ; capacitances (gate-drain overlap and channel, drain junctions) are now a few percent of the load and the table has none. The gap
 ; grows 21.8 and 21.1 times for a load 20 times smaller: it scales with 1 / C. So it is the capacitance, not the table, and nearly all of
-; the 0.14 % and 0.17 % at 10 pF is the same missing charge. M3b tabulates current only; terminal charges are the next step.
+; the 0.14 % and 0.17 % at 10 pF is the same missing charge. M3b tabulated current only; M3c tabulates the charge, next paragraph.
+; -
+; THE CHARGE (M3c, issue #10777). mos_psp103_data.t27 holds the DRAIN charge of each device on the same grid and corner, TT_QD_N and TT_QD_P (and
+; the gate charge TT_QG_N, TT_QG_P, for the first circuit that has a gate node), and the same four for ss. specs/analog/decks/psp103_q_tt.cir and
+; psp103_q_ss.cir print them: PSP103 reports no terminal charge, so the deck takes it from the terminal currents of one backward Euler step, which is
+; what ngspice integrates, and prints the current tables in the same run (they are TT_N, TT_P, SS_N, SS_P byte for byte). One unit is the voltage a
+; 500 fF node moves; the PMOS tables are negated so that every drain table grows with Vds. The solver is M3a's, generalized in mos_table.t27 (the
+; `*_q` functions; the current-only ones are the same functions with zero charge tables, so M3a's pinned run is untouched). The output node holds
+; C y + Qd_n - Qd_p, so the backward Euler step is
+;   (y - y0) + (qn(A, y) - qn(A0, y0)) - (qp(Vdd - A, Vdd - y) - qp(Vdd - A0, Vdd - y0)) + In - Ip = 0
+; with A0 the gate of the previous step: a gate edge is a jump of charge (the Miller kick), not a current. The Newton derivative adds the Vds slope of
+; each charge table to the two current slopes, each read off its own table (FR-002; the slope is the interpolant's movement across the cell, a test).
+; Against ngspice backward Euler on the same grid, tt, 500 fF, dt = 100 ps, 400 steps (the three errors apart):
+;   table + model  tpLH 13.4379 steps against ngspice's 13.4352: 0.020 % (3.09 % without charge); tpHL 10.1826 against 10.1830: 0.0041 % (3.63 %).
+;                  The waveform's worst of 401 samples is 6693 units (0.40 mV, 0.027 % of Vdd), 438317 (26 mV) without charge; after step 0 the
+;                  fixed-point run is 27 units (1.6 uV) from ngspice's sample, after the falling edge 2566 (0.15 mV).
+;   arithmetic     the fixed-point run against the f64 run of the same iteration: 4.03 units (0.24 uV), no delay (< 1e-5 step). The last Newton
+;                  residual is at most 2 units at every step (the bound, reached); K = 1 leaves 11636.
+;   discretization ngspice backward Euler at 100 ps against ngspice converged: the delays are 0.52 % and 0.67 % longer, which the charge does not
+;                  touch. The model with charge is 0.54 % and 0.67 % from the converged run: now all of it is the scheme's.
+; The tables never saw ngspice's own capacitance numbers, and agree with them: at Vgs = 0.9 V, Vds = 0.75 V the slopes of TT_QD_N and TT_QG_N are within
+; 0.16 %, 0.29 % and 0.19 % of cdd + cjd + cgdol, -(cdg + cgdol) and cgg + cgsol + cgdol as ngspice reports them (the deck's op_n block).
+; THE CHARGE'S NEGATIVE CONTROL. The ss charge tables in the tt run differ from tt's by 4.3 .. 5.3 % of their range, which moves the delays by only
+; 0.13 % and 0.17 % (and the waveform by 24012 units, 1.4 mV): the issue's 1 % bar cannot see it, so the check holds 0.05 % and 6700 units, and refuses
+; the ss charge by 2.5, 3.3 and 3.6 times. The charge missing (62 and 73 times over), the N and P charge tables in each other's slots (11 to 14
+; times), and either device's charge alone fail too.
+; THE WEAK INVERSION IS NOT MET (part 3 of #10777, follow-up #10881). 1 % with a floor of 1e-9 A at the midpoints and R2 points stays at the measured
+; 0.55 .. 0.73, and a test says why no read of these tables can do it: a table word is 298 pA, so rounding alone is up to 149 pA, 15 % of the floor
+; where 1 % allows 10 pA; and the NMOS current grows 6.5 .. 7 times per 62.5 mV cell, where a chord is more than 41 % above the exponential at the
+; middle. The reader that meets it is a log-domain table in finer units with an integer exp, a different reader: #10881.
 ; -
 ; THE NEGATIVE CONTROLS, tests that pass by seeing the failure. `check_tt` holds a pair of tables to the tt bounds (Id at the 1000
 ; R2 points of each device, the DC transfer and the 10 pF transient against ngspice tt). The tt pair passes. The ss pair used where
 ; tt is expected fails every measure, by a factor of 36 to 830 over its bound: 59 uA at a held-out point, 283 mV of DC output, 10.3 mV of
 ; trip point, 112 mV of waveform, 6.5 % and 16 % of delay. Each of ff, fs and sf used as tt fails too. The pair in the wrong slots
-; (PMOS table for the NMOS and the reverse) fails by factors of 100 to 10000. K = 1 breaks the residual bound by 4000.
+; (PMOS table for the NMOS and the reverse) fails by factors of 100 to 10000. K = 1 breaks the residual bound by 4000. The charge has its own
+; control, above.
 ; -
 ; THE HARDWARE. Everything the board would run is M3a's: `board_step` takes bool and unsigned words and no float, with the NMOS
 ; table and the PMOS table as two slices that a board would load at run time (a table is 625 words, each under 2^21): one bitstream per
@@ -77,13 +107,17 @@
 ; every tt held-out point and within 0.57 % where the current is at least 100 uA; the inverter built from the tt pair, solved by M3a's
 ; fixed three Newton iterations (last residual at most 1 unit), follows ngspice's backward Euler within 1.29 mV and 0.18 % of the delay
 ; at 10 pF, and ngspice's converged run within 0.37 % and 0.51 % of the delay; its DC transfer is within 5.5 mV and its trip point
-; within 12 uV of ngspice's, and in every corner within 6.6 mV and 0.27 mV; and the ss pair, or any other, used as tt is refused.
+; within 12 uV of ngspice's, and in every corner within 6.6 mV and 0.27 mV; and the ss pair, or any other, used as tt is refused. With the tt
+; drain-charge tables (read by the same reader, their slopes from the same nodes) the 500 fF inverter follows ngspice's backward Euler on the same
+; grid within 0.020 % and 0.0041 % of the two 50 % delays (the issue asked for 1 %) and 0.40 mV, three Newton iterations with a last residual of
+; at most 2 units, and the ss charge tables used as tt are refused.
 ; WHAT IT DOES NOT SAY. That this is PSP103: it is a TABLE MODEL of it, bilinear, at 27 C, for two device sizes, no body effect (bulk at
 ; the source), no Vds < 0, no Vgs < 0 (the table's range is 0 .. 1.5 V), nothing about another W, L, m, temperature, the high-voltage
 ; devices, mismatch or Monte Carlo. That the relative-error proposal of the issue (1 %, floor 1e-9 A) is met: it is not (see above); the
-; bound is the measured one. Nothing about charge: with 500 fF the model is 3 .. 4 % off, and a circuit whose loads are fF (the
-; golden SAR's comparator) needs charge tables. The transient is checked for tt only, on two circuits of one dt / C and one input (two
-; levels, an ideal step); the other corners are checked on Id and on the DC transfer. Nothing about a board (no bitstream, BRAM, LUT, DSP,
+; bound is the measured one (#10881). Charge: the drain charge of the tt and ss corners only (ff, fs, sf: #10882), for a 500 fF node and a
+; 100 ps step; the gate charge is tabulated and no circuit here reads it; the source and bulk charge is not tabulated (the source and bulk are
+; on the rails); the 10 pF run has no charge in it (its device capacitance is 0.2 % of the load). The transient is checked for tt only, on two
+; circuits of one dt / C and one input (two levels, an ideal step); the other corners are checked on Id and on the DC transfer. Nothing about a board (no bitstream, BRAM, LUT, DSP,
 ; Fmax, UART line), nothing about power, and no speed claim against SPICE (`may_claim_speedup` in bench.t27 says no and stays no).
 ; The uberWAVE deck itself is not run here, only its inverter's sizes are used. ASCII only, English only.
 ; phi^2 + 1/phi^2 = 3 | TRINITY
@@ -91,8 +125,8 @@
 module AnalogMosPsp103;
 
 use analog::linear::{fabs64};
-use analog::mos_table::{CELL, NCELL, NODES, VDD, K_NEWTON, RES_BOUND, RUN_STEPS, lim, id_of, gm_cell, gds_cell, mag_diff, monotone, fmax64, run_fix, run_tabf, max_of, board_step, inv_step};
-use analog::mos_psp103_data::{R2_A1, R2_A2, R2_TT_COUNT, CORNER_FF, CORNER_SS, CORNER_FS, CORNER_SF, TT_N, TT_P, SS_N, SS_P, TT_INV, TT_MID_N, TT_MID_P, TT_R2_N, TT_R2_P, INV10_BE, INV10_FINE, INV10_FINE_TPLH_FS, INV10_FINE_TPHL_FS, INV500_BE_TPLH_FS, INV500_BE_TPHL_FS, INV500_FINE_TPLH_FS, INV500_FINE_TPHL_FS, table_n, table_p, inv_dc, r2_n, r2_p, r2_count};
+use analog::mos_table::{CELL, NCELL, NODES, VDD, K_NEWTON, RES_BOUND, RUN_STEPS, lim, id_of, gm_cell, gds_cell, mag_diff, monotone, monotone_along, fmax64, run_fix, run_fix_q, run_tabf, run_tabf_q, max_of, board_step, inv_step, inv_step_q, g_resid_q, NOQ};
+use analog::mos_psp103_data::{R2_A1, R2_A2, R2_TT_COUNT, CORNER_FF, CORNER_SS, CORNER_FS, CORNER_SF, TT_N, TT_P, SS_N, SS_P, TT_INV, TT_MID_N, TT_MID_P, TT_R2_N, TT_R2_P, INV10_BE, INV10_FINE, INV10_FINE_TPLH_FS, INV10_FINE_TPHL_FS, INV500_BE_TPLH_FS, INV500_BE_TPHL_FS, INV500_FINE_TPLH_FS, INV500_FINE_TPHL_FS, INV500_BE, TT_QD_N, TT_QG_N, TT_QD_P, TT_QG_P, SS_QD_N, SS_QG_N, SS_QD_P, SS_QG_P, OP_CDD_AF, OP_CDG_AF, OP_CGG_AF, OP_CGSOL_AF, OP_CGDOL_AF, OP_CJD_AF, table_n, table_p, inv_dc, r2_n, r2_p, r2_count};
 
 pub const KIND : str = "analog-cell-check";
 pub const ID : str = "analog/mos_psp103";
@@ -134,9 +168,21 @@ pub const B_DELAY_REL : f64 = 0.0018;
 pub const B_DC_ALL_UNITS : f64 = 112000.0;
 pub const B_DC_ALL_VM : f64 = 0.0085;
 
+/// The 500 fF check, tt, against ngspice backward Euler on the same grid: each 50 % delay within 0.05 % (the model with its own charge is 0.020 %
+/// and 0.0041 % off), the waveform within 6700 units (0.40 mV) at all 401 samples. The issue asks for 1 %; this is the bound the tests hold,
+/// 20 times tighter, so that a wrong corner's charge (0.13 % and 0.17 %) is seen.
+pub const B_Q_DELAY_REL : f64 = 0.0005;
+pub const B_Q_WAVE_UNITS : f64 = 6700.0;
+/// The issue's bar for the delays: 1 %.
+pub const B_ISSUE_DELAY_REL : f64 = 0.01;
+/// aF per (table unit per cell): a charge table's slope read as a capacitance. 500 fF * 2^-24 V per unit, over 62.5 mV (2^20 units a cell).
+pub const AF_PER_UNIT_CELL : f64 = 0.476837158203125;
+
 ; --- what the tt run must give, derived by the tests below -----------------------------------------------
 /// The fold of the 400 samples of the fixed-point tt run (K_NEWTON iterations, from the DC rest value).
 pub const EXP_S_TT : u64 = 0xD9E25A3AAF1D8068;
+/// The same with the tt drain-charge tables, the circuit at 500 fF / 100 ps, the gate at Vdd before step 0.
+pub const EXP_S_TT_Q : u64 = 0x215DC5EAC4C48576;
 
 ; --- the held-out points ----------------------------------------------------------------------------
 
@@ -320,6 +366,81 @@ pub fn check_ok(id_abs: f64, id_b4: f64, dc_worst: f64, dc_vm: f64, tr_max: f64,
     return true;
 }
 
+; --- the 500 fF check with charge (M3c) -------------------------------------------------------------
+
+pub struct CheckQ { up_rel: f64, dn_rel: f64, wave: f64 }
+
+/// The 500 fF check, M3c: the f64 model of the inverter built from the current tables tn, tp and the drain-charge tables qn, qp
+/// (the gate before step 0 at Vdd) against ngspice's tt backward Euler at 500 fF and 100 ps on the same grid: the relative gap of
+/// each 50 % delay to ngspice's own .meas crossing, and the largest |dY| of the 401 samples in units.
+pub fn check_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32) -> CheckQ {
+    var outf : [400]f64 = [0.0] ** 400;
+    var fl : [401]f64 = [0.0] ** 401;
+    var be : [401]f64 = [0.0] ** 401;
+    const rest : u32 = dc_out(tn, tp, VDD);
+    run_tabf_q(tn, tp, qn, qp, K_NEWTON, rest as f64, VDD, &outf);
+    shift_run(rest as f64, &outf, &fl);
+    lift_u32(&INV500_BE, 401, &be);
+    return CheckQ { up_rel: rel_gap(delay_up(&fl), steps_of_fs(INV500_BE_TPLH_FS, EDGE1_500F_FS, DT_500F_FS)), dn_rel: rel_gap(delay_dn(&fl), steps_of_fs(INV500_BE_TPHL_FS, EDGE2_500F_FS, DT_500F_FS)), wave: max_diff(&fl, &be, 401) };
+}
+
+/// Is every measure of a 500 fF check inside the bounds?
+pub fn check_q_ok(up_rel: f64, dn_rel: f64, wave: f64) -> bool {
+    if (up_rel > B_Q_DELAY_REL) { return false; }
+    if (dn_rel > B_Q_DELAY_REL) { return false; }
+    if (wave > B_Q_WAVE_UNITS) { return false; }
+    return true;
+}
+
+/// How many of the 576 cells of a table have the interpolant's movement across the cell equal to the cell's slope, within the 2 units
+/// of the two reads and the slope's rounding: ax 1 across a cell of Vds against `gds_cell`, ax 0 across a cell of Vgs against
+/// `gm_cell`. For a table nondecreasing along that axis. The slope is read off the same nodes as the value (FR-002).
+pub fn slope_cells_ok(t: []const u32, ax: u32) -> u32 {
+    var ok : u32 = 0;
+    var cx : u32 = 0;
+    while (cx < NCELL) {
+        const fx : u32 = (cx * 977 + 31) % CELL;
+        var cy : u32 = 0;
+        while (cy < NCELL) {
+            const fy : u32 = (cy * 4093 + 17) % CELL;
+            if (ax == 1) {
+                const dy : u32 = id_of(t, cx * CELL + fx, (cy + 1) * CELL) - id_of(t, cx * CELL + fx, cy * CELL);
+                if (mag_diff(dy, gds_cell(t, cx * CELL + fx, cy * CELL + 12345)) <= 2) { ok = ok + 1; }
+            } else {
+                const dx : u32 = id_of(t, (cx + 1) * CELL, cy * CELL + fy) - id_of(t, cx * CELL, cy * CELL + fy);
+                if (mag_diff(dx, gm_cell(t, cx * CELL + 12345, cy * CELL + fy)) <= 2) { ok = ok + 1; }
+            }
+            cy = cy + 1;
+        }
+        cx = cx + 1;
+    }
+    return ok;
+}
+
+/// A charge table's slope across one cell as a capacitance in aF: the node values a (at the lower node) and b (at the upper), as the
+/// signed difference b - a in table units (a charge table that falls with the axis gives a negative capacitance).
+pub fn cap_af(a: u32, b: u32) -> f64 {
+    return (((b as f64) - (a as f64)) * AF_PER_UNIT_CELL);
+}
+
+/// The largest |a[i] - b[i]| over the 625 entries of two tables, as a fraction of the range (largest minus smallest entry) of a.
+pub fn table_gap_of_range(a: []const u32, b: []const u32) -> f64 {
+    var lo : u32 = a[0];
+    var hi : u32 = a[0];
+    var gap : u32 = 0;
+    var i : u32 = 0;
+    while (i < 625) {
+        const x : u32 = a[i as usize];
+        const y : u32 = b[i as usize];
+        if (x < lo) { lo = x; }
+        if (x > hi) { hi = x; }
+        const d : u32 = mag_diff(x, y);
+        if (d > gap) { gap = d; }
+        i = i + 1;
+    }
+    return (gap as f64) / ((hi - lo) as f64);
+}
+
 ; --- the tests ----------------------------------------------------------------------------------
 
 test the_units_the_grid_and_the_held_out_points_agree_with_m3a_and_the_decks {
@@ -800,3 +921,278 @@ test one_newton_iteration_breaks_the_residual_bound_on_these_tables {
     assert board_step(&TT_N, &TT_P, 9000000, 0) > 9000000;
     assert s2 != s1;
 }
+
+test the_charge_tables_are_capacitances_read_off_one_table_and_agree_with_ngspices_own_op_point {
+    // THE SHAPE. A drain table grows with Vds and a gate table with Vgs (capacitances are positive), in both corners, for both
+    // devices (the PMOS tables are stored negated, so they grow with Vsd and Vsg). And the other axis does NOT grow: the drain
+    // charge falls as the gate rises (the gate-drain coupling, the Miller kick), the gate charge falls as the drain rises.
+    assert(monotone_along(&TT_QD_N, 1));
+    assert(monotone_along(&TT_QD_P, 1));
+    assert(monotone_along(&TT_QG_N, 0));
+    assert(monotone_along(&TT_QG_P, 0));
+    assert(monotone_along(&SS_QD_N, 1));
+    assert(monotone_along(&SS_QD_P, 1));
+    assert(monotone_along(&SS_QG_N, 0));
+    assert(monotone_along(&SS_QG_P, 0));
+    assert(monotone_along(&TT_QD_N, 0) == false);
+    assert(monotone_along(&TT_QD_P, 0) == false);
+    assert(monotone_along(&TT_QG_N, 1) == false);
+    assert(monotone_along(&TT_QG_P, 1) == false);
+    // FR-002. The slope the Newton step uses is read off the same table as the charge: across every one of the 576 cells of each of
+    // the eight tables the interpolant moves by exactly the cell's slope (to the 2 units of the reads), 4608 cell checks.
+    assert_eq(slope_cells_ok(&TT_QD_N, 1), 576);
+    assert_eq(slope_cells_ok(&TT_QD_P, 1), 576);
+    assert_eq(slope_cells_ok(&TT_QG_N, 0), 576);
+    assert_eq(slope_cells_ok(&TT_QG_P, 0), 576);
+    assert_eq(slope_cells_ok(&SS_QD_N, 1), 576);
+    assert_eq(slope_cells_ok(&SS_QD_P, 1), 576);
+    assert_eq(slope_cells_ok(&SS_QG_N, 0), 576);
+    assert_eq(slope_cells_ok(&SS_QG_P, 0), 576);
+    // SIZES. The drain capacitance of the NMOS at zero bias is 1.01 fF (overlap and junction), its gate capacitance in strong inversion
+    // 6.1 fF (the oxide, 0.45 um x 1 um), and the PMOS, three times as wide, 17 fF: a ratio of 2.8.
+    const cdd0 : f64 = cap_af(TT_QD_N[0], TT_QD_N[1]);
+    assert cdd0 > 1000.0;
+    assert cdd0 < 1025.0;
+    const cgn : f64 = cap_af(TT_QG_N[19 * 25 + 1], TT_QG_N[20 * 25 + 1]);
+    const cgp : f64 = cap_af(TT_QG_P[19 * 25 + 1], TT_QG_P[20 * 25 + 1]);
+    assert cgn > 6000.0;
+    assert cgn < 6200.0;
+    assert cgp / cgn > 2.7;
+    assert cgp / cgn < 2.95;
+    // AN INDEPENDENT CHECK against what ngspice reports itself. At the NMOS operating point Vgs = 0.9 V, Vds = 0.75 V the model gives
+    // the channel capacitances cdd, cdg, cgg and, apart, the overlaps cgsol, cgdol and the junction cjd (aF, pasted from the deck's
+    // op_n block). The drain terminal's capacitances are cdd + cjd + cgdol to the drain, -(cdg + cgdol) to the gate; the gate's own is
+    // cgg + cgsol + cgdol. The tables were made from currents and never saw these numbers; their slopes at that point are within 1 %.
+    const vg : u32 = 15099494;
+    const vd : u32 = 12 * CELL;
+    const sl : u32 = gds_cell(&TT_QD_N, vg, vd - 1000);
+    const sh : u32 = gds_cell(&TT_QD_N, vg, vd + 1000);
+    const cdd_t : f64 = ((sl + sh) as f64) * 0.5 * AF_PER_UNIT_CELL;
+    const cdg_t : f64 = 0.0 - cap_af(TT_QD_N[14 * 25 + 12], TT_QD_N[15 * 25 + 12]);
+    const cgg_t : f64 = cap_af(TT_QG_N[14 * 25 + 12], TT_QG_N[15 * 25 + 12]);
+    const cdd_o : f64 = ((OP_CDD_AF + OP_CJD_AF + OP_CGDOL_AF) as f64);
+    const cdg_o : f64 = ((OP_CDG_AF + OP_CGDOL_AF) as f64);
+    const cgg_o : f64 = ((OP_CGG_AF + OP_CGSOL_AF + OP_CGDOL_AF) as f64);
+    assert rel_gap(cdd_t, cdd_o) < 0.01;
+    assert rel_gap(cdg_t, cdg_o) < 0.01;
+    assert rel_gap(cgg_t, cgg_o) < 0.01;
+    // 0.16 %, 0.29 %, 0.19 %: both sides.
+    assert rel_gap(cdd_t, cdd_o) > 0.001;
+    assert rel_gap(cdg_t, cdg_o) > 0.002;
+    assert rel_gap(cgg_t, cgg_o) > 0.001;
+    // NEGATIVE CONTROL: the ss table's gate-drain coupling at the same point is 2.6 % off tt's, which the 1 % sees.
+    const cdg_s : f64 = 0.0 - cap_af(SS_QD_N[14 * 25 + 12], SS_QD_N[15 * 25 + 12]);
+    assert rel_gap(cdg_s, cdg_o) > 0.02;
+    assert rel_gap(cdg_s, cdg_o) < 0.03;
+}
+
+test the_tt_inverter_at_500_fF_with_charge_is_within_1_percent_and_within_the_measured_bounds_of_ngspice_backward_euler {
+    var out : [400]u32 = [0] ** 400;
+    var res : [400]u32 = [0] ** 400;
+    var outf : [400]f64 = [0.0] ** 400;
+    var tmp : [400]f64 = [0.0] ** 400;
+    var fix : [401]f64 = [0.0] ** 401;
+    var fl : [401]f64 = [0.0] ** 401;
+    var be : [401]f64 = [0.0] ** 401;
+    const rest : u32 = dc_out(&TT_N, &TT_P, VDD);
+    const s : u64 = run_fix_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P, K_NEWTON, rest, VDD, false, &out, &res);
+    run_tabf_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P, K_NEWTON, rest as f64, VDD, &outf);
+    lift_u32(&out, 400, &tmp);
+    shift_run(rest as f64, &tmp, &fix);
+    shift_run(rest as f64, &outf, &fl);
+    lift_u32(&INV500_BE, 401, &be);
+    const up_be5 : f64 = steps_of_fs(INV500_BE_TPLH_FS, EDGE1_500F_FS, DT_500F_FS);
+    const dn_be5 : f64 = steps_of_fs(INV500_BE_TPHL_FS, EDGE2_500F_FS, DT_500F_FS);
+    const up_fix : f64 = delay_up(&fix);
+    const dn_fix : f64 = delay_dn(&fix);
+    // The run: from the DC rest value (3 units), 400 steps, the fold pinned; the last Newton residual is at most 2 units (119 nV) at
+    // every step, at the bound.
+    assert_eq(rest, 3);
+    assert_eq(s, EXP_S_TT_Q);
+    assert max_of(&res) <= RES_BOUND;
+    assert_eq(max_of(&res), 2);
+    // THE ISSUE'S BAR: both 50 % delays within 1 % of ngspice backward Euler on the same grid (3.09 % and 3.63 % without charge).
+    assert rel_gap(up_fix, up_be5) < B_ISSUE_DELAY_REL;
+    assert rel_gap(dn_fix, dn_be5) < B_ISSUE_DELAY_REL;
+    // What is measured, both sides: tpLH 13.4379 steps against ngspice's 13.4352 (0.020 %, 0.27 ps longer), tpHL 10.1826 against
+    // 10.1830 (0.0041 %, 0.04 ps shorter).
+    assert rel_gap(up_fix, up_be5) < 0.00021;
+    assert rel_gap(up_fix, up_be5) > 0.00019;
+    assert rel_gap(dn_fix, dn_be5) < 0.000044;
+    assert rel_gap(dn_fix, dn_be5) > 0.000039;
+    assert up_fix > up_be5;
+    assert dn_fix < dn_be5;
+    // ngspice's two reads of its own delay agree: the crossing read off the 100 ps samples and its .meas, to 0.0001 step (10 fs).
+    assert fabs64(delay_up(&be) - up_be5) < 0.0001;
+    assert fabs64(delay_dn(&be) - dn_be5) < 0.0001;
+    // The waveform: the f64 model against ngspice's 401 samples, the largest |dY| is 6693 units (0.40 mV, 0.027 % of Vdd), bound 6700.
+    const wave : f64 = max_diff(&fl, &be, 401);
+    assert wave < B_Q_WAVE_UNITS;
+    assert wave > 6680.0;
+    // The two gate edges, where the charge does its work. After step 0 the output is 41 mV off the rail: the fixed-point run lands 27
+    // units (1.6 uV) from ngspice's sample; after step 200 (the falling edge) 2566 units (0.15 mV).
+    assert mag_diff(out[0], INV500_BE[1]) < 30;
+    assert mag_diff(out[0], INV500_BE[1]) > 20;
+    assert mag_diff(out[200], INV500_BE[201]) < 2600;
+    assert mag_diff(out[200], INV500_BE[201]) > 2500;
+    // ARITHMETIC: the fixed-point run against the f64 run of the same iteration, 4.03 units (0.24 uV), no delay (< 1e-5 step).
+    const arith : f64 = max_diff(&fix, &fl, 401);
+    assert arith < 4.1;
+    assert arith > 3.95;
+    assert fabs64(up_fix - delay_up(&fl)) < 0.00001;
+    assert fabs64(dn_fix - delay_dn(&fl)) < 0.00001;
+}
+
+test the_three_errors_at_500_fF_kept_apart_and_what_the_charge_removed {
+    const wq : CheckQ = check_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P);
+    const nq : CheckQ = check_q(&TT_N, &TT_P, &NOQ, &NOQ);
+    // TABLE + MODEL (the f64 model against ngspice backward Euler on the same grid): with the charge 0.020 % and 0.0041 % of the
+    // delay and 6693 units on the waveform; without it 3.09 % and 3.63 % and 438317 units (26 mV, 1.7 % of Vdd).
+    assert wq.up_rel < 0.00021;
+    assert wq.dn_rel < 0.000044;
+    assert nq.up_rel > 0.0307;
+    assert nq.up_rel < 0.0311;
+    assert nq.dn_rel > 0.0361;
+    assert nq.dn_rel < 0.0364;
+    assert nq.wave > 438000.0;
+    assert nq.wave < 438700.0;
+    // The charge removes 99.4 % of the delay gap in tpLH and 99.9 % in tpHL, and 98.5 % of the waveform's.
+    assert wq.up_rel * 100.0 < nq.up_rel;
+    assert wq.dn_rel * 500.0 < nq.dn_rel;
+    assert wq.wave * 60.0 < nq.wave;
+    // DISCRETIZATION: ngspice backward Euler at 100 ps against ngspice converged at 500 fF, which the charge does not touch: the
+    // delays are 0.52 % and 0.67 % longer, as at 10 pF.
+    const up_be5 : f64 = steps_of_fs(INV500_BE_TPLH_FS, EDGE1_500F_FS, DT_500F_FS);
+    const dn_be5 : f64 = steps_of_fs(INV500_BE_TPHL_FS, EDGE2_500F_FS, DT_500F_FS);
+    const up_fine5 : f64 = steps_of_fs(INV500_FINE_TPLH_FS, EDGE1_500F_FS, DT_500F_FS);
+    const dn_fine5 : f64 = steps_of_fs(INV500_FINE_TPHL_FS, EDGE2_500F_FS, DT_500F_FS);
+    assert rel_gap(up_be5, up_fine5) > 0.0050;
+    assert rel_gap(up_be5, up_fine5) < 0.0054;
+    assert rel_gap(dn_be5, dn_fine5) > 0.0065;
+    assert rel_gap(dn_be5, dn_fine5) < 0.0069;
+    // All three together, the fixed-point model against ngspice's converged run: 0.54 % and 0.67 % of the delay; the discretization is
+    // now the whole of it (ngspice's own scheme at 100 ps), and the model with charge follows that scheme to 0.02 %.
+    var out : [400]u32 = [0] ** 400;
+    var res : [400]u32 = [0] ** 400;
+    var tmp : [400]f64 = [0.0] ** 400;
+    var fix : [401]f64 = [0.0] ** 401;
+    const rest : u32 = dc_out(&TT_N, &TT_P, VDD);
+    const s : u64 = run_fix_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P, K_NEWTON, rest, VDD, false, &out, &res);
+    lift_u32(&out, 400, &tmp);
+    shift_run(rest as f64, &tmp, &fix);
+    assert rel_gap(delay_up(&fix), up_fine5) > 0.0053;
+    assert rel_gap(delay_up(&fix), up_fine5) < 0.0056;
+    assert rel_gap(delay_dn(&fix), dn_fine5) > 0.0065;
+    assert rel_gap(delay_dn(&fix), dn_fine5) < 0.0069;
+    // K = 1 breaks the residual bound with charge too (11636 units); K = 2 already meets it.
+    var out1 : [400]u32 = [0] ** 400;
+    var res1 : [400]u32 = [0] ** 400;
+    const s1 : u64 = run_fix_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P, 1, rest, VDD, false, &out1, &res1);
+    assert max_of(&res1) > 10000;
+    assert s1 != s;
+}
+
+test the_ss_charge_tables_in_the_tt_run_fail_the_500_fF_check_and_so_do_the_other_wrong_charges {
+    // THE NEGATIVE CONTROL. The check is the tt 500 fF check: both delays and the waveform against ngspice tt backward Euler. The tt
+    // current tables with the tt charge pass it; with the ss charge tables (the same tables of the slower corner, 4.3 .. 5.3 % smaller
+    // in range) it fails every measure. The issue's bar of 1 % would not see the ss charge (0.13 % and 0.17 %): the check holds 0.05 %.
+    const g : CheckQ = check_q(&TT_N, &TT_P, &TT_QD_N, &TT_QD_P);
+    assert check_q_ok(g.up_rel, g.dn_rel, g.wave);
+    const b : CheckQ = check_q(&TT_N, &TT_P, &SS_QD_N, &SS_QD_P);
+    assert check_q_ok(b.up_rel, b.dn_rel, b.wave) == false;
+    assert b.up_rel > 2.4 * B_Q_DELAY_REL;
+    assert b.dn_rel > 3.2 * B_Q_DELAY_REL;
+    assert b.wave > 3.5 * B_Q_WAVE_UNITS;
+    // 0.127 % and 0.167 % of the delay (0.17 ps and 0.17 ps) and 24012 units (1.43 mV) of waveform, 6 .. 40 times the tt's own error.
+    assert b.up_rel > 0.0012;
+    assert b.up_rel < 0.0013;
+    assert b.dn_rel > 0.0016;
+    assert b.dn_rel < 0.0017;
+    assert b.wave > 24000.0;
+    assert b.wave < 24100.0;
+    assert b.up_rel > 6.0 * g.up_rel;
+    assert b.dn_rel > 40.0 * g.dn_rel;
+    assert b.up_rel < B_ISSUE_DELAY_REL;
+    assert b.dn_rel < B_ISSUE_DELAY_REL;
+    // The charge missing altogether fails (3.09 % and 3.63 %), the NMOS and PMOS charge tables in each other's slots fail (0.59 %
+    // and 0.57 %, 91415 units), and either device's charge alone is not enough (the NMOS's: 2.0 % and 2.9 %; the PMOS's: 1.0 % and 0.8 %).
+    const nq : CheckQ = check_q(&TT_N, &TT_P, &NOQ, &NOQ);
+    assert check_q_ok(nq.up_rel, nq.dn_rel, nq.wave) == false;
+    assert nq.up_rel > 60.0 * B_Q_DELAY_REL;
+    const sw : CheckQ = check_q(&TT_N, &TT_P, &TT_QD_P, &TT_QD_N);
+    assert check_q_ok(sw.up_rel, sw.dn_rel, sw.wave) == false;
+    assert sw.up_rel > 10.0 * B_Q_DELAY_REL;
+    assert sw.dn_rel > 10.0 * B_Q_DELAY_REL;
+    assert sw.wave > 13.0 * B_Q_WAVE_UNITS;
+    const only_n : CheckQ = check_q(&TT_N, &TT_P, &TT_QD_N, &NOQ);
+    assert check_q_ok(only_n.up_rel, only_n.dn_rel, only_n.wave) == false;
+    assert only_n.up_rel > 30.0 * B_Q_DELAY_REL;
+    const only_p : CheckQ = check_q(&TT_N, &TT_P, &NOQ, &TT_QD_P);
+    assert check_q_ok(only_p.up_rel, only_p.dn_rel, only_p.wave) == false;
+    assert only_p.up_rel > 15.0 * B_Q_DELAY_REL;
+    // The other direction: the ss current tables with the tt charge fail the check too (the current error dominates).
+    const wrong_i : CheckQ = check_q(&SS_N, &SS_P, &TT_QD_N, &TT_QD_P);
+    assert check_q_ok(wrong_i.up_rel, wrong_i.dn_rel, wrong_i.wave) == false;
+}
+
+test the_issue_s_weak_inversion_one_percent_is_out_of_reach_of_a_chord_over_these_tables_and_why {
+    // PART 3 OF #10777, the part that is a bound and not a fix. The issue proposes at most 1 % relative error with a floor of 1e-9 A at the
+    // 576 midpoints and 1000 R2 points. M3b measured 0.55 .. 0.73 (test the_tt_tables_at_the_held_out_points_against_ngspice, above:
+    // rel9 between 0.5 and 0.75). Two facts, both of the tables and not of the solver, say why no read of THESE tables can do better.
+    // (a) ROUNDING. A table word is 298 pA; half a word, the rounding error of a node, is 149 pA, 15 % of the 1 nA floor, where the
+    // criterion allows 10 pA. Below a few nA the table itself is 10 .. 100 % off before any interpolation.
+    const half : f64 = amps_of(0.5);
+    assert half > 0.00000000014;
+    assert half < 0.00000000015;
+    assert half / 0.000000001 > 0.14;
+    assert half > 10.0 * (0.01 * 0.000000001);
+    // (b) GROWTH. In weak inversion the NMOS current grows between 6.5 and 7 times per 62.5 mV cell (ngspice's own currents at three
+    // neighbouring cell midpoints, Vds = 1.4 V, from 0.6 nA to 200 nA); a straight line between two nodes whose ratio is r is above the
+    // exponential at the middle by (1 + r) / (2 sqrt(r)), and ((1 + r) / 2)^2 > 2 r means by more than 41 %.
+    var i : u32 = 0;
+    var checked : u32 = 0;
+    while (i < 3) {
+        const a : f64 = (TT_MID_N[(i * 24 + 20) as usize]) as f64;
+        const b : f64 = (TT_MID_N[((i + 1) * 24 + 20) as usize]) as f64;
+        const r : f64 = b / a;
+        assert r > 6.5;
+        assert r < 7.0;
+        assert ((1.0 + r) * 0.5) * ((1.0 + r) * 0.5) > 2.0 * r;
+        checked = checked + 1;
+        i = i + 1;
+    }
+    assert_eq(checked, 3);
+    // (c) WHAT WOULD DO IT. Tabulating the current in units of 0.3 pA would put 1 % of the floor at 33 words, but the largest word, 402 uA,
+    // is then 1.3e9, and the reader's u64 interpolation needs every word under 2^23: the read would need 128-bit sums or a table in the log
+    // domain (an integer exp). That is a different reader, not a finer grid: the next issue.
+    var wide : f64 = 0.000402305 / 0.0000000000003;
+    var bound : f64 = 8388608.0;
+    assert wide > 100.0 * bound;
+    // And on the grid there is no free lunch: the midpoint chord error of the issue's measure is above 0.5 for the NMOS and above 0.6
+    // for the PMOS (0.51 .. 0.62 and 0.64 .. 0.73 over midpoints and R2 points).
+    const mn : Hold = hold_stats(&TT_N, &TT_MID_N, 0, 576);
+    const mp : Hold = hold_stats(&TT_P, &TT_MID_P, 0, 576);
+    assert mn.rel9 > 0.5;
+    assert mp.rel9 > 0.6;
+    assert mn.rel9 > 50.0 * 0.01;
+}
+
+test the_ss_charge_tables_differ_from_tt_by_four_to_five_percent_of_range {
+    // Why the negative control of the charge is the delay at 0.05 % and not a coarse count: the corners' charge tables are alike. The largest
+    // |tt - ss| entry is 4.7 % (NMOS drain), 5.3 % (NMOS gate), 4.3 % (PMOS drain) and 4.6 % (PMOS gate) of the tt table's range, so
+    // the delays move by 0.13 % and 0.17 %.
+    const dn : f64 = table_gap_of_range(&TT_QD_N, &SS_QD_N);
+    const gn : f64 = table_gap_of_range(&TT_QG_N, &SS_QG_N);
+    const dp : f64 = table_gap_of_range(&TT_QD_P, &SS_QD_P);
+    const gp : f64 = table_gap_of_range(&TT_QG_P, &SS_QG_P);
+    assert dn > 0.046;
+    assert dn < 0.048;
+    assert gn > 0.052;
+    assert gn < 0.054;
+    assert dp > 0.043;
+    assert dp < 0.044;
+    assert gp > 0.045;
+    assert gp < 0.047;
+    // A table against itself is 0.
+    assert table_gap_of_range(&TT_QD_N, &TT_QD_N) < 0.0000001;
+}
diff --git a/specs/analog/mos_psp103_data.t27 b/specs/analog/mos_psp103_data.t27
index 5ad6a93cac..1061dfbf90 100644
--- a/specs/analog/mos_psp103_data.t27
+++ b/specs/analog/mos_psp103_data.t27
@@ -1,6 +1,6 @@
 // SPDX-License-Identifier: Apache-2.0
 ; specs/analog/mos_psp103_data.t27 -- the IHP SG13G2 PSP103 current tables, per corner, from ngspice: data and nothing else
-; M3b of the analog emulation epic (#10486, issue #10516). M3a (specs/analog/mos_table.t27) reads a 25 x 25 table by exact-integer
+; M3b of the analog emulation epic (#10486, issue #10516); M3c (#10777, sub-issue #10880) adds the charge tables, THE CHARGE TABLES below. M3a (specs/analog/mos_table.t27) reads a 25 x 25 table by exact-integer
 ; bilinear interpolation; this file holds the tables it reads when the device is the foundry's, Id(Vgs, Vds) of sg13_lv_nmos
 ; (W 1.0 um, L 0.45 um) and sg13_lv_pmos (W 3.0 um, L 0.45 um) for the five corners mos_tt, mos_ff, mos_ss, mos_fs and mos_sf, and the
 ; numbers the tests judge them by (held-out currents, the inverter's DC transfer, the inverter's transient). The functions are in
@@ -13,7 +13,9 @@
 ;   psp103_tab_tt.cir, _ff, _ss, _fs, _sf   the two tables, the inverter's DC transfer, the held-out currents (the tt deck also prints the
 ;                                         576 cell midpoints and 1000 R2 points; the others 200 R2 points)
 ;   psp103_inv_be.cir, psp103_inv_fine.cir   the inverter transient at 10 pF, backward Euler on a fixed 2 ns grid and a converged run
-;   psp103_inv_be_500f.cir, psp103_inv_fine_500f.cir   the same at 500 fF and 100 ps: two crossing times each
+;   psp103_inv_be_500f.cir, psp103_inv_fine_500f.cir   the same at 500 fF and 100 ps: two crossing times each (the backward Euler deck
+;                                         also prints the 401-sample waveform, INV500_BE)
+;   psp103_q_tt.cir, psp103_q_ss.cir   (M3c) the drain and gate charge tables, the current tables again, and the capacitances ngspice reports at one point
 ; Each is run as `IHP_PDK=<checkout> IHP_OSDI=<dir with psp103.osdi> ngspice -b <deck>`; no deck holds a path. ngspice rounds each current
 ; to a table unit itself (nint in the deck) and prints it between "BEGIN name" and "END name"; the rows were wrapped onto lines
 ; below in print order, none changed or reordered. The pins (order-sensitive folds, as M3a pins its run) fail on any changed number.
@@ -31,13 +33,29 @@
 ; Vgs (Vsg) = 1.5 frac(k * 0.7548776662466927) and Vds (Vsd) = 1.5 frac(k * 0.5698402909980532), 1/plastic number and its square; the
 ; other corners print the first 200 of those points. Each pair of devices has its own set of currents, the same coordinates.
 ; -
+; THE CHARGE TABLES (M3c, issue #10777). Two more decks, specs/analog/decks/psp103_q_tt.cir and psp103_q_ss.cir (corners mos_tt and
+; mos_ss), print the DRAIN and GATE charge of each device on the same 25 x 25 grid: TT_QD_N TT_QG_N TT_QD_P TT_QG_P and the SS_ four.
+; PSP103 reports no terminal charge as an output variable (`show` lists currents, conductances and the capacitances cdd, cdg, cgg, ...
+; with the overlap cgdol and junction cjd apart), so the deck takes the whole charge of a terminal, channel plus overlap plus fringe plus
+; junction, from the terminal CURRENT of one backward Euler step, which is the quantity ngspice integrates: Q(B) - Q(A) = (i_A - i_return) dt
+; for a step to the grid node B and back to the origin A = (Vgs, Vds) = (0, 0), dt = 10 ns, where the device is off. The same deck run
+; prints the current tables (grid_n, grid_p), which are byte for byte TT_N, TT_P, SS_N, SS_P: the current tables and the charge tables come
+; from one run of one netlist. A charge entry is Q / 500 fF in units of 2^-24 V (the voltage a 500 fF node moves when it takes the
+; charge; one unit is 2.98e-20 C), relative to the origin plus 2^21 so that no entry is negative (only differences are physical); the PMOS
+; tables are NEGATED, so that every drain table grows with Vds and every gate table with Vgs, as capacitances are positive. The scale
+; ties a charge table to C = 500 fF, as a current table is tied to dt / C: for another load the entries scale by 500 fF / C.
+; The tt and ss charge tables differ by 4.3 .. 5.3 % of their range, which is why the negative control in mos_psp103.t27 is the delay and
+; not a count of volts.
+; -
 ; WHAT A PASS SAYS. The numbers in this file are the numbers the decks printed, in order (the pins), the five corners are in the order
 ; the foundry names them (ff above tt above ss in drive current for both devices; fs a fast NMOS with a slow PMOS, sf the reverse), and
 ; the units are the ones the header states.
 ; WHAT IT DOES NOT SAY. That these are the tables a board loads: the board gets the same words, but no bitstream, BRAM or UART line is
 ; in this milestone. That the numbers are right: they are ngspice's, at its solver tolerances (reltol 1e-6, vntol 1e-9, abstol 1e-14),
 ; for the two device sizes and 27 C, nothing about another W, L, m, temperature, the high-voltage devices, mismatch or Monte Carlo.
-; That PSP103 is a physical transistor: it is the foundry's model of one. The charge (capacitance) of the devices is not tabulated.
+; That PSP103 is a physical transistor: it is the foundry's model of one. The drain and gate charge is tabulated for the corners tt and ss only
+; (ff, fs, sf have none), the source and bulk charge not at all (a source or bulk off the rail is not a circuit this file serves), and
+; a charge table is for C = 500 fF.
 ; ASCII only, English only.
 ; phi^2 + 1/phi^2 = 3 | TRINITY
 
@@ -74,6 +92,17 @@ pub const R2_A2 : f64 = 0.5698402909980532;
 /// The R2 points given per corner: 1000 for tt, 200 for the others.
 pub const R2_TT_COUNT : u32 = 1000;
 pub const R2_OTHER_COUNT : u32 = 200;
+/// The charge tables: an entry is Q / C_REF in units of 2^-24 V, relative to the grid origin, plus Q_OFFSET (so no entry is negative).
+pub const C_REF_FF : u32 = 500;
+pub const Q_OFFSET : u32 = 2097152;
+/// The capacitances ngspice itself reports at the NMOS operating point Vgs = 0.9 V, Vds = 0.75 V, in aF (deck psp103_q_tt.cir, block op_n):
+/// the channel capacitances cdd, cdg, cgg and, each apart, the gate-source and gate-drain overlaps and the drain junction.
+pub const OP_CDD_AF : u32 = 81;
+pub const OP_CDG_AF : u32 = 1398;
+pub const OP_CGG_AF : u32 = 3483;
+pub const OP_CGSOL_AF : u32 = 641;
+pub const OP_CGDOL_AF : u32 = 615;
+pub const OP_CJD_AF : u32 = 328;
 
 ; --- order-sensitive folds of the pasted numbers (fold of linear.t27): a changed number changes a pin ---------
 pub const PIN_TT : u64 = 0x814310C06BAD1CA9;
@@ -84,6 +113,11 @@ pub const PIN_SF : u64 = 0x13F09713FAAE89F6;
 pub const PIN_TT_HELD : u64 = 0x4D63964ABAF11967;
 pub const PIN_OTHER_HELD : u64 = 0x63623E6E75519E74;
 pub const PIN_TRAN : u64 = 0x859D61FDD24C4A10;
+/// The charge tables of corners tt and ss, in the order qd_n, qg_n, qd_p, qg_p (the four tables of a corner folded one after another).
+/// The 500 fF waveform INV500_BE.
+pub const PIN_TRAN500 : u64 = 0xB35700152266D9B9;
+pub const PIN_TT_Q : u64 = 0xFE5696611BC5AB3B;
+pub const PIN_SS_Q : u64 = 0x06627578C483E909;
 
 ; --- the numbers --------------------------------------------------------------------------------
 ; --- corner mos_tt ---------------------------------------------------------------------------
@@ -1069,6 +1103,425 @@ pub const INV500_BE_TPHL_FS : u32 = 21118300;
 pub const INV500_FINE_TPLH_FS : u32 = 1436550;
 pub const INV500_FINE_TPHL_FS : u32 = 21111500;
 
+/// ngspice backward Euler on the 100 ps grid, 500 fF: output in units of 2^-24 V at t = 100 ps (row 0, at rest), then after steps 0..399 (deck psp103_inv_be_500f.cir, block tran_be).
+pub const INV500_BE : [401]u32 = [
+    11, 689206, 1688379, 2683758, 3674715, 4660554, 5640416, 6613237, 7577703, 8532198, 9474757, 10403029, 11314268, 12205373,
+    13072974, 13913588, 14723817, 15500554, 16241174, 16943663, 17606687, 18229588, 18812333, 19355431, 19859838, 20326861, 20758070, 21155219,
+    21520180, 21854886, 22161291, 22441331, 22696901, 22929832, 23141878, 23334705, 23509886, 23668900, 23813125, 23943845, 24062252, 24169443,
+    24266432, 24354151, 24433452, 24505118, 24569861, 24628334, 24681130, 24728787, 24771799, 24810608, 24845621, 24877204, 24905689, 24931375,
+    24954537, 24975419, 24994245, 25011215, 25026512, 25040298, 25052724, 25063922, 25074013, 25083107, 25091301, 25098684, 25105337, 25111331,
+    25116732, 25121597, 25125981, 25129930, 25133488, 25136694, 25139582, 25142183, 25144527, 25146638, 25148540, 25150253, 25151797, 25153187,
+    25154439, 25155568, 25156584, 25157499, 25158324, 25159067, 25159736, 25160339, 25160882, 25161371, 25161812, 25162209, 25162566, 25162888,
+    25163178, 25163440, 25163675, 25163887, 25164078, 25164250, 25164405, 25164545, 25164670, 25164784, 25164886, 25164978, 25165060, 25165135,
+    25165202, 25165263, 25165317, 25165366, 25165410, 25165450, 25165486, 25165518, 25165548, 25165574, 25165597, 25165619, 25165638, 25165655,
+    25165671, 25165685, 25165697, 25165709, 25165719, 25165728, 25165736, 25165744, 25165751, 25165757, 25165762, 25165767, 25165772, 25165776,
+    25165779, 25165782, 25165785, 25165788, 25165790, 25165792, 25165794, 25165796, 25165798, 25165799, 25165800, 25165802, 25165803, 25165803,
+    25165804, 25165805, 25165806, 25165806, 25165807, 25165807, 25165808, 25165808, 25165809, 25165809, 25165809, 25165809, 25165810, 25165810,
+    25165810, 25165810, 25165810, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811, 25165811,
+    25165811, 25165811, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812, 25165812,
+    25165812, 25165812, 25165812, 25165812, 25165812, 24235796, 22909274, 21591916, 20285047, 18990196, 17709362, 16445188, 15201183, 13981967,
+    12793447, 11642816, 10538225, 9488083, 8500108, 7580413, 6732891, 5959040, 5258154, 4627741, 4064025, 3562407, 3117856, 2725209,
+    2379377, 2075490, 1808981, 1575635, 1371605, 1193415, 1037945, 902409, 784334, 681532, 592072, 514258, 446597, 387784,
+    336676, 292273, 253703, 220206, 191119, 165865, 143941, 124909, 108390, 94053, 81610, 70812, 61441, 53309,
+    46254, 40131, 34819, 30210, 26211, 22741, 19731, 17119, 14853, 12887, 11182, 9702, 8419, 7305,
+    6339, 5500, 4773, 4142, 3595, 3120, 2708, 2351, 2041, 1772, 1539, 1336, 1161, 1008,
+    876, 762, 662, 576, 501, 436, 380, 331, 288, 252, 220, 192, 168, 147,
+    129, 113, 100, 88, 78, 69, 61, 54, 49, 44, 39, 35, 32, 29,
+    27, 25, 23, 21, 20, 19, 18, 17, 16, 15, 15, 14, 14, 13,
+    13, 13, 12, 12, 12, 12, 12, 11, 11, 11, 11, 11, 11, 11,
+    11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11,
+    11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11,
+    11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11,
+    11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11,
+    11, 11, 11, 11, 11, 11, 11, 11, 11
+];
+
+; --- the charge tables, corners mos_tt and mos_ss -----------------------------------------------------
+/// QD_N, corner mos_tt: NMOS drain charge, Vgs = ix * 62.5 mV, Vds = iy * 62.5 mV, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_tt.cir, block qd_n).
+pub const TT_QD_N : [625]u32 = [
+    2097152, 2099276, 2101369, 2103434, 2105472, 2107485, 2109473, 2111439, 2113383, 2115305, 2117207, 2119090, 2120953, 2122798,
+    2124625, 2126435, 2128229, 2130006, 2131768, 2133515, 2135248, 2136966, 2138671, 2140362, 2142041, 2095875, 2098007, 2100109,
+    2102182, 2104228, 2106250, 2108247, 2110222, 2112175, 2114106, 2116018, 2117909, 2119782, 2121636, 2123472, 2125291, 2127093,
+    2128879, 2130649, 2132405, 2134145, 2135871, 2137583, 2139282, 2140968, 2094581, 2096727, 2098837, 2100918, 2102973, 2105003,
+    2107009, 2108993, 2110955, 2112896, 2114816, 2116717, 2118599, 2120462, 2122307, 2124135, 2125946, 2127741, 2129519, 2131283,
+    2133031, 2134765, 2136485, 2138192, 2139885, 2093224, 2095416, 2097537, 2099626, 2101689, 2103728, 2105743, 2107736, 2109707,
+    2111657, 2113587, 2115497, 2117388, 2119260, 2121114, 2122951, 2124771, 2126575, 2128362, 2130134, 2131891, 2133634, 2135362,
+    2137076, 2138777, 2091539, 2093955, 2096109, 2098210, 2100282, 2102330, 2104354, 2106356, 2108337, 2110296, 2112236, 2114156,
+    2116056, 2117938, 2119803, 2121649, 2123479, 2125291, 2127088, 2128869, 2130635, 2132386, 2134122, 2135845, 2137554, 2088850,
+    2091937, 2094234, 2096363, 2098450, 2100510, 2102547, 2104560, 2106552, 2108523, 2110474, 2112405, 2114317, 2116211, 2118086,
+    2119943, 2121784, 2123607, 2125414, 2127206, 2128982, 2130742, 2132488, 2134220, 2135939, 2084917, 2088868, 2091611, 2093839,
+    2095959, 2098040, 2100094, 2102123, 2104131, 2106116, 2108081, 2110027, 2111952, 2113859, 2115748, 2117619, 2119472, 2121308,
+    2123127, 2124931, 2126719, 2128491, 2130248, 2131991, 2133720, 2080177, 2084746, 2088197, 2090716, 2092925, 2095046, 2097126,
+    2099177, 2101203, 2103207, 2105189, 2107151, 2109093, 2111015, 2112919, 2114805, 2116672, 2118523, 2120356, 2122173, 2123974,
+    2125759, 2127530, 2129285, 2131026, 2074982, 2079923, 2084017, 2087059, 2089482, 2091689, 2093815, 2095897, 2097949, 2099974,
+    2101976, 2103957, 2105917, 2107856, 2109777, 2111679, 2113562, 2115428, 2117277, 2119108, 2120923, 2122723, 2124506, 2126275,
+    2128029, 2069507, 2074681, 2079231, 2082853, 2085653, 2088033, 2090244, 2092376, 2094463, 2096516, 2098542, 2100545, 2102524,
+    2104483, 2106422, 2108341, 2110241, 2112124, 2113988, 2115835, 2117666, 2119479, 2121277, 2123060, 2124827, 2063843, 2069176,
+    2074032, 2078143, 2081409, 2084073, 2086431, 2088645, 2090784, 2092876, 2094932, 2096960, 2098963, 2100943, 2102902, 2104840,
+    2106759, 2108658, 2110539, 2112403, 2114249, 2116078, 2117891, 2119688, 2121470, 2058043, 2063489, 2068557, 2073037, 2076756,
+    2079783, 2082364, 2084706, 2086924, 2089069, 2091165, 2093225, 2095255, 2097260, 2099240, 2101199, 2103137, 2105055, 2106954,
+    2108835, 2110698, 2112543, 2114372, 2116184, 2117980, 2052138, 2057670, 2062891, 2067642, 2071744, 2075157, 2078024, 2080548,
+    2082878, 2085097, 2087246, 2089347, 2091410, 2093443, 2095449, 2097430, 2099390, 2101328, 2103246, 2105145, 2107025, 2108887,
+    2110732, 2112560, 2114372, 2046150, 2051750, 2057086, 2062037, 2066447, 2070218, 2073403, 2076157, 2078639, 2080958, 2083177,
+    2085329, 2087432, 2089498, 2091533, 2093541, 2095523, 2097483, 2099422, 2101340, 2103239, 2105119, 2106981, 2108825, 2110653,
+    2040092, 2045747, 2051173, 2056277, 2060927, 2065011, 2068513, 2071529, 2074199, 2076647, 2078955, 2081172, 2083325, 2085431,
+    2087499, 2089536, 2091544, 2093528, 2095488, 2097427, 2099346, 2101244, 2103124, 2104986, 2106830, 2033978, 2039678, 2045176,
+    2050399, 2055239, 2059583, 2063380, 2066668, 2069555, 2072161, 2074581, 2076878, 2079091, 2081244, 2083351, 2085420, 2087457,
+    2089467, 2091451, 2093412, 2095351, 2097270, 2099168, 2101048, 2102909, 2027815, 2033553, 2039110, 2044428, 2049418, 2053977,
+    2058035, 2061590, 2064710, 2067498, 2070051, 2072446, 2074732, 2076940, 2079092, 2081198, 2083268, 2085306, 2087315, 2089300,
+    2091261, 2093200, 2095119, 2097017, 2098897, 2021610, 2027381, 2032988, 2038383, 2043493, 2048229, 2052513, 2056317, 2059673,
+    2062659, 2065367, 2067877, 2070249, 2072523, 2074726, 2076875, 2078980, 2081049, 2083087, 2085097, 2087082, 2089043, 2090982,
+    2092900, 2094798, 2015369, 2021169, 2026817, 2032276, 2037485, 2042366, 2046843, 2050873, 2054458, 2057650, 2060528, 2063171,
+    2065643, 2067995, 2070258, 2072455, 2074600, 2076704, 2078772, 2080809, 2082818, 2084802, 2086763, 2088702, 2090619, 2009096,
+    2014921, 2020604, 2026117, 2031408, 2036410, 2041052, 2045280, 2049081, 2052481, 2055539, 2058329, 2060917, 2063357, 2065690,
+    2067942, 2070132, 2072273, 2074374, 2076440, 2078476, 2080485, 2082468, 2084428, 2086366, 2002796, 2008642, 2014357, 2019915,
+    2025275, 2030378, 2035159, 2039563, 2043561, 2047162, 2050405, 2053353, 2056069, 2058611, 2061023, 2063338, 2065579, 2067763,
+    2069899, 2071997, 2074061, 2076095, 2078102, 2080085, 2082044, 1996472, 2002337, 2008078, 2013676, 2019094, 2024282, 2029182,
+    2033738, 2037916, 2041707, 2045135, 2048248, 2051103, 2053758, 2056259, 2058645, 2060944, 2063175, 2065351, 2067483, 2069577,
+    2071638, 2073671, 2075676, 2077657, 1990126, 1996008, 2001773, 2007405, 2012873, 2018133, 2023134, 2027824, 2032162, 2036130,
+    2039737, 2043017, 2046020, 2048797, 2051399, 2053866, 2056229, 2058512, 2060732, 2062902, 2065028, 2067119, 2069177, 2071207,
+    2073211, 1983762, 1989658, 1995444, 2001106, 2006617, 2011939, 2017027, 2021833, 2026313, 2030444, 2034222, 2037669, 2040823,
+    2043732, 2046443, 2048999, 2051435, 2053777, 2056046, 2058256, 2060419, 2062540, 2064626, 2066682, 2068709, 1977380, 1983289,
+    1989094, 1994782, 2000330, 2005707, 2010870, 2015776, 2020383, 2024661, 2028600, 2032209, 2035517, 2038563, 2041392, 2044046,
+    2046562, 2048971, 2051294, 2053549, 2055750, 2057905, 2060021, 2062103, 2064155
+];
+/// QG_N, corner mos_tt: NMOS gate charge, same axes, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_tt.cir, block qg_n).
+pub const TT_QG_N : [625]u32 = [
+    2097152, 2095883, 2094623, 2093371, 2092128, 2090893, 2089668, 2088451, 2087244, 2086046, 2084857, 2083677, 2082506, 2081345,
+    2080192, 2079048, 2077913, 2076786, 2075668, 2074557, 2073455, 2072360, 2071273, 2070194, 2069122, 2101014, 2099736, 2098468,
+    2097207, 2095955, 2094712, 2093478, 2092252, 2091035, 2089828, 2088630, 2087441, 2086261, 2085091, 2083929, 2082776, 2081632,
+    2080497, 2079370, 2078252, 2077142, 2076039, 2074945, 2073858, 2072778, 2104866, 2103573, 2102296, 2101028, 2099768, 2098516,
+    2097272, 2096038, 2094812, 2093596, 2092388, 2091190, 2090001, 2088821, 2087651, 2086489, 2085336, 2084192, 2083057, 2081930,
+    2080812, 2079702, 2078599, 2077505, 2076418, 2108799, 2107447, 2106157, 2104880, 2103612, 2102351, 2101099, 2099855, 2098621,
+    2097395, 2096178, 2094971, 2093772, 2092583, 2091403, 2090233, 2089071, 2087918, 2086774, 2085639, 2084512, 2083393, 2082283,
+    2081181, 2080086, 2113288, 2111638, 2110308, 2109018, 2107740, 2106470, 2105208, 2103955, 2102710, 2101475, 2100248, 2099031,
+    2097823, 2096624, 2095434, 2094253, 2093082, 2091920, 2090766, 2089622, 2088486, 2087359, 2086240, 2085129, 2084026, 2119557,
+    2117020, 2115507, 2114183, 2112888, 2111604, 2110329, 2109064, 2107807, 2106559, 2105320, 2104091, 2102871, 2101660, 2100458,
+    2099266, 2098084, 2096910, 2095746, 2094591, 2093444, 2092307, 2091178, 2090057, 2088945, 2128052, 2124375, 2122293, 2120848,
+    2119513, 2118205, 2116911, 2115628, 2114354, 2113090, 2111836, 2110591, 2109356, 2108130, 2106914, 2105708, 2104512, 2103325,
+    2102147, 2100979, 2099820, 2098670, 2097529, 2096396, 2095273, 2137999, 2133492, 2130499, 2128690, 2127248, 2125893, 2124569,
+    2123261, 2121966, 2120682, 2119408, 2118145, 2116893, 2115650, 2114418, 2113195, 2111983, 2110781, 2109588, 2108405, 2107232,
+    2106068, 2104913, 2103768, 2102631, 2148768, 2143751, 2139908, 2137437, 2135731, 2134272, 2132893, 2131549, 2130226, 2128917,
+    2127622, 2126338, 2125066, 2123804, 2122553, 2121313, 2120084, 2118865, 2117656, 2116457, 2115268, 2114089, 2112920, 2111760,
+    2110610, 2160041, 2154700, 2150240, 2147017, 2144838, 2143166, 2141686, 2140283, 2138918, 2137578, 2136255, 2134947, 2133652,
+    2132370, 2131099, 2129840, 2128592, 2127355, 2126129, 2124913, 2123708, 2122513, 2121329, 2120154, 2118989, 2171652, 2166092,
+    2161207, 2157335, 2154559, 2152534, 2150874, 2149372, 2147946, 2146561, 2145204, 2143866, 2142545, 2141239, 2139947, 2138667,
+    2137399, 2136143, 2134898, 2133665, 2132443, 2131231, 2130030, 2128839, 2127659, 2183511, 2177791, 2172609, 2168236, 2164870,
+    2162384, 2160446, 2158789, 2157268, 2155821, 2154417, 2153043, 2151691, 2150357, 2149040, 2147737, 2146448, 2145171, 2143907,
+    2142655, 2141414, 2140185, 2138968, 2137761, 2136564, 2195558, 2189717, 2184318, 2179571, 2175694, 2172711, 2170413, 2168531,
+    2166873, 2165337, 2163871, 2162449, 2161060, 2159694, 2158348, 2157020, 2155707, 2154408, 2153123, 2151851, 2150591, 2149344,
+    2148108, 2146884, 2145671, 2207755, 2201820, 2196255, 2191227, 2186933, 2183479, 2180775, 2178605, 2176760, 2175102, 2173553,
+    2172072, 2170635, 2169231, 2167852, 2166495, 2165156, 2163833, 2162525, 2161232, 2159952, 2158686, 2157431, 2156189, 2154959,
+    2220074, 2214062, 2208369, 2203124, 2198500, 2194630, 2191514, 2189011, 2186932, 2185116, 2183459, 2181900, 2180405, 2178955,
+    2177539, 2176148, 2174780, 2173431, 2172099, 2170783, 2169482, 2168194, 2166921, 2165660, 2164412, 2232495, 2226420, 2220623,
+    2215208, 2210321, 2206098, 2202594, 2199741, 2197388, 2195376, 2193583, 2191928, 2190362, 2188857, 2187395, 2185968, 2184567,
+    2183189, 2181830, 2180490, 2179166, 2177857, 2176563, 2175283, 2174016, 2245002, 2238873, 2232992, 2227441, 2222342, 2217825,
+    2213972, 2210772, 2208122, 2205883, 2203924, 2202150, 2200497, 2198926, 2197413, 2195943, 2194505, 2193095, 2191708, 2190341,
+    2188993, 2187662, 2186346, 2185045, 2183759, 2257582, 2251408, 2245457, 2239794, 2234523, 2229762, 2225604, 2222075, 2219121,
+    2216630, 2214479, 2212564, 2210807, 2209157, 2207583, 2206063, 2204585, 2203140, 2201721, 2200326, 2198952, 2197597, 2196258,
+    2194936, 2193629, 2270226, 2264013, 2258001, 2252247, 2246835, 2241870, 2237448, 2233619, 2230366, 2227610, 2225245, 2223165,
+    2221285, 2219543, 2217898, 2216321, 2214797, 2213312, 2211860, 2210435, 2209033, 2207652, 2206290, 2204946, 2203617, 2282924,
+    2276678, 2270615, 2264782, 2259253, 2254118, 2249470, 2245372, 2241834, 2238811, 2236216, 2233952, 2231930, 2230079, 2228351,
+    2226710, 2225133, 2223604, 2222115, 2220657, 2219226, 2217818, 2216431, 2215063, 2213713, 2295670, 2289394, 2283288, 2277389,
+    2271760, 2266482, 2261640, 2257303, 2253502, 2250215, 2247383, 2244919, 2242737, 2240763, 2238938, 2237222, 2235586, 2234009,
+    2232478, 2230984, 2229522, 2228085, 2226672, 2225280, 2223907, 2308457, 2302157, 2296011, 2290056, 2284343, 2278943, 2273935,
+    2269388, 2265345, 2261807, 2258736, 2256061, 2253705, 2251591, 2249657, 2247855, 2246151, 2244518, 2242942, 2241410, 2239913,
+    2238446, 2237005, 2235588, 2234191, 2321282, 2314959, 2308780, 2302775, 2296990, 2291486, 2286334, 2281603, 2277342, 2273569,
+    2270262, 2267372, 2264829, 2262561, 2260505, 2258605, 2256823, 2255129, 2253501, 2251926, 2250393, 2248893, 2247424, 2245979,
+    2244557, 2334138, 2327795, 2321587, 2315539, 2309691, 2304097, 2298822, 2293930, 2289474, 2285481, 2281949, 2278842, 2276105,
+    2273672, 2271479, 2269470, 2267601, 2265836, 2264152, 2262529, 2260955, 2259421, 2257920, 2256447, 2255000, 2347023, 2340663,
+    2334429, 2328343, 2322440, 2316769, 2311386, 2306352, 2301721, 2297528, 2293781, 2290462, 2287527, 2284920, 2282580, 2280450,
+    2278482, 2276639, 2274890, 2273214, 2271596, 2270023, 2268489, 2266987, 2265512
+];
+/// QD_P, corner mos_tt: PMOS drain charge, NEGATED, Vsg = ix * 62.5 mV, Vsd = iy * 62.5 mV, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_tt.cir, block qd_p).
+pub const TT_QD_P : [625]u32 = [
+    2097152, 2102473, 2107698, 2112834, 2117885, 2122854, 2127747, 2132565, 2137311, 2141989, 2146601, 2151150, 2155637, 2160065,
+    2164437, 2168754, 2173019, 2177233, 2181399, 2185518, 2189592, 2193622, 2197611, 2201560, 2205470, 2094093, 2099444, 2104701,
+    2109869, 2114953, 2119957, 2124884, 2129737, 2134518, 2139231, 2143877, 2148459, 2152979, 2157439, 2161842, 2166190, 2170485,
+    2174728, 2178921, 2183068, 2187168, 2191224, 2195237, 2199209, 2203142, 2091000, 2096383, 2101671, 2106870, 2111987, 2117024,
+    2121985, 2126873, 2131689, 2136436, 2141116, 2145732, 2150286, 2154779, 2159215, 2163594, 2167919, 2172192, 2176415, 2180589,
+    2184716, 2188798, 2192837, 2196834, 2200790, 2087850, 2093282, 2098599, 2103829, 2108978, 2114048, 2119042, 2123963, 2128814,
+    2133596, 2138311, 2142961, 2147549, 2152076, 2156544, 2160956, 2165312, 2169616, 2173868, 2178071, 2182226, 2186336, 2190401,
+    2194423, 2198404, 2084466, 2090067, 2095425, 2100685, 2105864, 2110966, 2115993, 2120947, 2125832, 2130648, 2135398, 2140083,
+    2144705, 2149267, 2153768, 2158213, 2162602, 2166936, 2171219, 2175452, 2179636, 2183773, 2187865, 2191913, 2195919, 2079900,
+    2086316, 2091784, 2097082, 2102293, 2107426, 2112486, 2117473, 2122392, 2127242, 2132026, 2136746, 2141403, 2145998, 2150534,
+    2155012, 2159433, 2163800, 2168114, 2172377, 2176590, 2180755, 2184875, 2188950, 2192981, 2071925, 2080644, 2086587, 2091974,
+    2097229, 2102399, 2107493, 2112516, 2117469, 2122355, 2127175, 2131930, 2136622, 2141252, 2145823, 2150335, 2154790, 2159189,
+    2163535, 2167829, 2172072, 2176267, 2180415, 2184517, 2188575, 2060048, 2071579, 2078990, 2084686, 2090040, 2095267, 2100407,
+    2105471, 2110464, 2115389, 2120247, 2125041, 2129772, 2134440, 2139048, 2143597, 2148088, 2152524, 2156904, 2161232, 2165508,
+    2169735, 2173913, 2178045, 2182132, 2045790, 2059264, 2068986, 2075612, 2081234, 2086575, 2091785, 2096904, 2101945, 2106915,
+    2111817, 2116653, 2121426, 2126135, 2130784, 2135373, 2139903, 2144377, 2148795, 2153159, 2157471, 2161732, 2165944, 2170108,
+    2174227, 2030254, 2044859, 2056665, 2064982, 2071280, 2076880, 2082219, 2087420, 2092525, 2097549, 2102500, 2107383, 2112201,
+    2116955, 2121646, 2126277, 2130850, 2135364, 2139822, 2144225, 2148574, 2152872, 2157120, 2161319, 2165470, 2013984, 2029263,
+    2042515, 2052725, 2060246, 2066395, 2071989, 2077331, 2082530, 2087625, 2092636, 2097573, 2102441, 2107242, 2111980, 2116656,
+    2121272, 2125829, 2130328, 2134772, 2139161, 2143498, 2147783, 2152018, 2156205, 1997258, 2012965, 2027156, 2038951, 2048008,
+    2055090, 2061158, 2066753, 2072102, 2077301, 2082392, 2087394, 2092320, 2097175, 2101963, 2106687, 2111349, 2115951, 2120494,
+    2124981, 2129412, 2133788, 2138113, 2142386, 2146609, 1980231, 1996223, 2011035, 2023998, 2034542, 2042845, 2049666, 2055683,
+    2061280, 2066638, 2071840, 2076929, 2081926, 2086843, 2091688, 2096465, 2101177, 2105827, 2110417, 2114949, 2119423, 2123842,
+    2128207, 2132520, 2136782, 1962999, 1979189, 1994425, 2008216, 2020000, 2029596, 2037404, 2044054, 2050037, 2055637, 2061002,
+    2066210, 2071299, 2076293, 2081204, 2086041, 2090809, 2095512, 2100152, 2104731, 2109252, 2113716, 2118124, 2122479, 2126781,
+    1945624, 1961955, 1977490, 1991869, 2004615, 2015397, 2024305, 2031776, 2038306, 2044261, 2049864, 2055237, 2060451, 2065543,
+    2070536, 2075444, 2080275, 2085036, 2089731, 2094362, 2098932, 2103443, 2107897, 2112296, 2116641, 1928150, 1944582, 1960333,
+    1975137, 1988610, 2000393, 2010377, 2018783, 2026012, 2032452, 2038386, 2043991, 2049372, 2054592, 2059688, 2064681, 2069588,
+    2074415, 2079171, 2083859, 2088483, 2093045, 2097548, 2101994, 2106385, 1910607, 1927113, 1943023, 1958140, 1972158, 1984748,
+    1995703, 2005067, 2013099, 2020147, 2026517, 2032432, 2038038, 2043427, 2048654, 2053754, 2058750, 2063656, 2068481, 2073233,
+    2077916, 2082534, 2087090, 2091586, 2096025, 1893020, 1909578, 1925607, 1940957, 1955389, 1968615, 1980400, 1990670, 1999549,
+    2007297, 2014204, 2020516, 2026416, 2032023, 2037418, 2042652, 2047757, 2052756, 2057662, 2062487, 2067236, 2071916, 2076529,
+    2081080, 2085570, 1875405, 1892000, 1908119, 1923646, 1938392, 1952118, 1964591, 1975674, 1985383, 1993880, 2001405, 2008199,
+    2014464, 2020348, 2025956, 2031357, 2036597, 2041708, 2046710, 2051619, 2056444, 2061192, 2065869, 2070479, 2075025, 1857778,
+    1874397, 1890583, 1906246, 1921235, 1935352, 1948386, 1960172, 1970654, 1979905, 1988098, 1995443, 2002145, 2008369, 2014240,
+    2019847, 2025253, 2030500, 2035616, 2040623, 2045534, 2050360, 2055108, 2059783, 2064390, 1840150, 1856783, 1873019, 1888785,
+    1903964, 1918390, 1931877, 1944258, 1955433, 1965408, 1974284, 1982229, 1989428, 1996051, 2002239, 2008098, 2013705, 2019115,
+    2024367, 2029490, 2034501, 2039415, 2044243, 2048991, 2053666, 1822529, 1839168, 1855440, 1871288, 1886614, 1901284, 1915138,
+    1928015, 1939797, 1950440, 1959986, 1968552, 1976292, 1983369, 1989925, 1996082, 2001929, 2007534, 2012948, 2018206, 2023334,
+    2028350, 2033268, 2038098, 2042847, 1804923, 1821562, 1837860, 1853769, 1869212, 1884076, 1898226, 1911516, 1923819, 1935063,
+    1945244, 1954428, 1962734, 1970303, 1977275, 1983774, 1989902, 1995737, 2001342, 2006759, 2012022, 2017155, 2022176, 2027098,
+    2031930, 1787338, 1803972, 1820286, 1836242, 1851777, 1866797, 1881188, 1894817, 1907564, 1919339, 1930104, 1939885, 1948762,
+    1956849, 1964272, 1971154, 1977602, 1983703, 1989529, 1995131, 2000551, 2005819, 2010957, 2015983, 2020908, 1769780, 1786404,
+    1802728, 1818719, 1834326, 1849471, 1864056, 1877967, 1891089, 1903325, 1914618, 1924960, 1934396, 1943010, 1950910, 1958208,
+    1965011, 1971413, 1977491, 1983306, 1988906, 1994328, 1999600, 2004744, 2009774
+];
+/// QG_P, corner mos_tt: PMOS gate charge, NEGATED, same axes, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_tt.cir, block qg_p).
+pub const TT_QG_P : [625]u32 = [
+    2097152, 2094124, 2091131, 2088171, 2085245, 2082354, 2079497, 2076675, 2073887, 2071134, 2068415, 2065730, 2063078, 2060458,
+    2057869, 2055311, 2052783, 2050284, 2047813, 2045368, 2042950, 2040558, 2038189, 2035845, 2033523, 2108913, 2105855, 2102830,
+    2099837, 2096878, 2093952, 2091061, 2088204, 2085381, 2082594, 2079841, 2077122, 2074436, 2071784, 2069164, 2066575, 2064017,
+    2061488, 2058989, 2056518, 2054073, 2051655, 2049262, 2046894, 2044549, 2120595, 2117504, 2114448, 2111424, 2108432, 2105473,
+    2102547, 2099655, 2096798, 2093976, 2091188, 2088435, 2085715, 2083030, 2080377, 2077757, 2075168, 2072610, 2070082, 2067582,
+    2065110, 2062666, 2060247, 2057854, 2055486, 2132252, 2129106, 2126020, 2122965, 2119941, 2116949, 2113990, 2111064, 2108172,
+    2105315, 2102492, 2099705, 2096951, 2094232, 2091546, 2088894, 2086273, 2083684, 2081126, 2078597, 2076098, 2073626, 2071181,
+    2068763, 2066370, 2144196, 2140840, 2137709, 2134624, 2131570, 2128546, 2125554, 2122595, 2119669, 2116777, 2113920, 2111098,
+    2108310, 2105557, 2102838, 2100152, 2097500, 2094880, 2092291, 2089733, 2087204, 2084705, 2082233, 2079789, 2077371, 2158071,
+    2153678, 2150415, 2147290, 2144203, 2141149, 2138125, 2135133, 2132175, 2129250, 2126359, 2123503, 2120682, 2117895, 2115143,
+    2112425, 2109741, 2107090, 2104471, 2101884, 2099327, 2096800, 2094302, 2091833, 2089390, 2177753, 2170446, 2166605, 2163378,
+    2160245, 2157154, 2154096, 2151071, 2148079, 2145120, 2142196, 2139306, 2136451, 2133631, 2130847, 2128097, 2125381, 2122699,
+    2120051, 2117434, 2114850, 2112296, 2109772, 2107277, 2104810, 2204137, 2193274, 2187654, 2184064, 2180821, 2177668, 2174564,
+    2171496, 2168464, 2165467, 2162504, 2159577, 2156685, 2153828, 2151007, 2148221, 2145470, 2142754, 2140072, 2137423, 2134807,
+    2132223, 2129670, 2127146, 2124652, 2234619, 2221263, 2212801, 2208110, 2204557, 2201277, 2198096, 2194971, 2191889, 2188845,
+    2185839, 2182868, 2179934, 2177036, 2174174, 2171348, 2168558, 2165804, 2163084, 2160400, 2157748, 2155130, 2152543, 2149988,
+    2147463, 2267294, 2252465, 2241377, 2234649, 2230304, 2226728, 2223402, 2220186, 2217036, 2213936, 2210877, 2207858, 2204877,
+    2201934, 2199028, 2196160, 2193328, 2190532, 2187773, 2185049, 2182359, 2179704, 2177082, 2174492, 2171933, 2301221, 2285506,
+    2272551, 2263487, 2257688, 2253473, 2249856, 2246481, 2243228, 2240050, 2236927, 2233851, 2230818, 2227826, 2224872, 2221958,
+    2219081, 2216242, 2213440, 2210674, 2207945, 2205250, 2202590, 2199963, 2197369, 2335918, 2319642, 2305448, 2294380, 2286715,
+    2281397, 2277231, 2273569, 2270146, 2266853, 2263644, 2260498, 2257403, 2254354, 2251348, 2248384, 2245459, 2242574, 2239726,
+    2236917, 2234145, 2231409, 2228709, 2226044, 2223412, 2371114, 2354465, 2339442, 2326866, 2317354, 2310568, 2305517, 2301367,
+    2297661, 2294190, 2290857, 2287617, 2284445, 2281329, 2278263, 2275243, 2272266, 2269330, 2266435, 2263580, 2260763, 2257984,
+    2255241, 2252535, 2249864, 2406643, 2389737, 2374144, 2360480, 2349398, 2341027, 2334796, 2329905, 2325755, 2322007, 2318491,
+    2315118, 2311845, 2308646, 2305508, 2302424, 2299388, 2296398, 2293451, 2290545, 2287680, 2284855, 2282069, 2279319, 2276607,
+    2442395, 2425311, 2409314, 2394866, 2382543, 2372692, 2365127, 2359260, 2354470, 2350313, 2346527, 2342967, 2339555, 2336249,
+    2333023, 2329862, 2326759, 2323706, 2320701, 2317742, 2314825, 2311951, 2309117, 2306322, 2303566, 2478294, 2461085, 2444797,
+    2429778, 2416505, 2405383, 2396489, 2389494, 2383877, 2379155, 2374988, 2371165, 2367564, 2364115, 2360776, 2357521, 2354337,
+    2351213, 2348144, 2345124, 2342151, 2339223, 2336338, 2333495, 2330692, 2514287, 2496990, 2480489, 2465048, 2451051, 2438892,
+    2428780, 2420614, 2414033, 2408594, 2403919, 2399739, 2395882, 2392243, 2388758, 2385386, 2382104, 2378895, 2375751, 2372663,
+    2369627, 2366641, 2363700, 2360805, 2357951, 2550331, 2532976, 2516317, 2500559, 2486011, 2473021, 2461856, 2452565, 2444955,
+    2438680, 2433373, 2428731, 2424538, 2420650, 2416975, 2413454, 2410051, 2406741, 2403508, 2400342, 2397236, 2394184, 2391183,
+    2388230, 2385323, 2586395, 2569004, 2552228, 2536231, 2521259, 2507612, 2495562, 2485250, 2476618, 2469436, 2463394, 2458186,
+    2453569, 2449363, 2445446, 2441737, 2438183, 2434749, 2431411, 2428154, 2424967, 2421842, 2418773, 2415757, 2412791, 2622456,
+    2605045, 2588184, 2572003, 2556703, 2542536, 2529758, 2518552, 2508957, 2500853, 2494006, 2488145, 2483015, 2478416, 2474196,
+    2470253, 2466512, 2462927, 2459463, 2456098, 2452817, 2449608, 2446463, 2443377, 2440346, 2658493, 2641075, 2624153, 2607833,
+    2592277, 2577698, 2564325, 2552354, 2541886, 2532887, 2525209, 2518629, 2512911, 2507844, 2503258, 2499026, 2495057, 2491287,
+    2487672, 2484179, 2480788, 2477482, 2474251, 2471087, 2467984, 2694490, 2677077, 2660113, 2643686, 2627930, 2613023, 2599166,
+    2586552, 2575311, 2565478, 2556976, 2549643, 2543278, 2537676, 2532659, 2528084, 2523840, 2519848, 2516050, 2512406, 2508886,
+    2505468, 2502139, 2498887, 2495704, 2730436, 2713035, 2696045, 2679538, 2663625, 2648455, 2634205, 2621052, 2609142, 2598551,
+    2589261, 2581172, 2574123, 2567934, 2562428, 2557453, 2552885, 2548629, 2544616, 2540792, 2537121, 2533574, 2530132, 2526779,
+    2523506, 2766319, 2748938, 2731934, 2715367, 2699331, 2683951, 2669380, 2655780, 2643296, 2632030, 2622008, 2613181, 2605438,
+    2598626, 2592583, 2587157, 2582216, 2577654, 2573388, 2569354, 2565506, 2561809, 2558237, 2554770, 2551395, 2802131, 2784776,
+    2767767, 2751158, 2735026, 2719478, 2704645, 2690673, 2677701, 2665842, 2655153, 2645629, 2637199, 2629747, 2623134, 2617214,
+    2611855, 2606943, 2602386, 2598109, 2594056, 2590186, 2586464, 2582867, 2579377
+];
+/// QD_N, corner mos_ss: NMOS drain charge, Vgs = ix * 62.5 mV, Vds = iy * 62.5 mV, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_ss.cir, block qd_n).
+pub const SS_QD_N : [625]u32 = [
+    2097152, 2099315, 2101447, 2103549, 2105624, 2107674, 2109698, 2111700, 2113679, 2115637, 2117575, 2119493, 2121392, 2123272,
+    2125135, 2126981, 2128810, 2130623, 2132421, 2134203, 2135971, 2137725, 2139466, 2141193, 2142907, 2095906, 2098076, 2100215,
+    2102325, 2104408, 2106465, 2108498, 2110508, 2112495, 2114462, 2116408, 2118334, 2120241, 2122130, 2124001, 2125855, 2127692,
+    2129513, 2131318, 2133108, 2134884, 2136645, 2138393, 2140127, 2141848, 2094651, 2096829, 2098975, 2101092, 2103183, 2105247,
+    2107288, 2109306, 2111302, 2113277, 2115231, 2117166, 2119081, 2120978, 2122858, 2124720, 2126565, 2128394, 2130207, 2132005,
+    2133788, 2135557, 2137312, 2139053, 2140781, 2093380, 2095572, 2097725, 2099849, 2101947, 2104019, 2106068, 2108094, 2110098,
+    2112081, 2114043, 2115986, 2117910, 2119816, 2121703, 2123574, 2125427, 2127264, 2129085, 2130891, 2132682, 2134458, 2136220,
+    2137969, 2139704, 2092037, 2094280, 2096444, 2098576, 2100680, 2102760, 2104817, 2106851, 2108863, 2110854, 2112825, 2114776,
+    2116709, 2118623, 2120519, 2122397, 2124259, 2126105, 2127934, 2129748, 2131546, 2133331, 2135100, 2136856, 2138599, 2090314,
+    2092815, 2095015, 2097157, 2099271, 2101359, 2103423, 2105466, 2107486, 2109486, 2111466, 2113426, 2115368, 2117291, 2119196,
+    2121083, 2122954, 2124808, 2126646, 2128468, 2130276, 2132068, 2133846, 2135610, 2137360, 2087518, 2090733, 2093095, 2095267,
+    2097395, 2099496, 2101572, 2103625, 2105656, 2107667, 2109657, 2111628, 2113580, 2115514, 2117429, 2119327, 2121208, 2123073,
+    2124921, 2126753, 2128570, 2130371, 2132159, 2133932, 2135691, 2083536, 2087581, 2090424, 2092706, 2094870, 2096991, 2099084,
+    2101152, 2103198, 2105223, 2107227, 2109211, 2111176, 2113123, 2115051, 2116962, 2118855, 2120731, 2122591, 2124435, 2126263,
+    2128076, 2129874, 2131658, 2133427, 2078831, 2083438, 2086982, 2089578, 2091837, 2094000, 2096120, 2098209, 2100274, 2102315,
+    2104336, 2106336, 2108316, 2110277, 2112219, 2114144, 2116051, 2117941, 2119814, 2121671, 2123511, 2125337, 2127147, 2128943,
+    2130724, 2073721, 2078663, 2082810, 2085931, 2088416, 2090669, 2092836, 2094956, 2097045, 2099108, 2101147, 2103165, 2105162,
+    2107139, 2109098, 2111037, 2112959, 2114863, 2116751, 2118621, 2120476, 2122314, 2124137, 2125946, 2127739, 2068361, 2073514,
+    2078079, 2081754, 2084615, 2087046, 2089298, 2091469, 2093592, 2095682, 2097745, 2099783, 2101799, 2103794, 2105769, 2107725,
+    2109663, 2111583, 2113485, 2115370, 2117239, 2119091, 2120928, 2122749, 2124555, 2062832, 2068128, 2072972, 2077104, 2080413,
+    2083122, 2085519, 2087770, 2089945, 2092071, 2094163, 2096225, 2098263, 2100278, 2102272, 2104245, 2106200, 2108136, 2110054,
+    2111954, 2113838, 2115705, 2117555, 2119390, 2121210, 2057180, 2062578, 2067617, 2072090, 2075824, 2078878, 2081488, 2083862,
+    2086111, 2088287, 2090416, 2092508, 2094571, 2096609, 2098623, 2100616, 2102588, 2104541, 2106476, 2108392, 2110291, 2112173,
+    2114038, 2115888, 2117721, 2051432, 2056909, 2062088, 2066812, 2070903, 2074317, 2077195, 2079737, 2082089, 2084334, 2086510,
+    2088639, 2090733, 2092796, 2094833, 2096847, 2098839, 2100810, 2102762, 2104695, 2106610, 2108507, 2110388, 2112251, 2114099,
+    2045607, 2051146, 2056432, 2061342, 2065720, 2069468, 2072640, 2075391, 2077880, 2080213, 2082451, 2084625, 2086754, 2088847,
+    2090910, 2092946, 2094959, 2096950, 2098920, 2100871, 2102803, 2104716, 2106612, 2108491, 2110353, 2039721, 2045310, 2050678,
+    2055731, 2060334, 2064373, 2067837, 2070826, 2073482, 2075928, 2078243, 2080472, 2082643, 2084770, 2086862, 2088923, 2090959,
+    2092970, 2094960, 2096929, 2098878, 2100808, 2102720, 2104615, 2106492, 2033781, 2039413, 2044847, 2050012, 2054793, 2059077,
+    2062814, 2066050, 2068898, 2071481, 2073890, 2076187, 2078407, 2080573, 2082697, 2084786, 2086846, 2088880, 2090889, 2092877,
+    2094845, 2096792, 2098721, 2100631, 2102524, 2027798, 2033464, 2038954, 2044207, 2049130, 2053617, 2057599, 2061079, 2064135,
+    2066874, 2069396, 2071774, 2074054, 2076265, 2078424, 2080544, 2082630, 2084687, 2086718, 2088726, 2090712, 2092678, 2094624,
+    2096551, 2098459, 2021776, 2027473, 2033008, 2038334, 2043372, 2048028, 2052223, 2055933, 2059201, 2062114, 2064766, 2067239,
+    2069588, 2071851, 2074052, 2076205, 2078319, 2080402, 2082456, 2084485, 2086490, 2088474, 2090437, 2092381, 2094306, 2015721,
+    2021444, 2027018, 2032404, 2037537, 2042332, 2046711, 2050633, 2054110, 2057206, 2060003, 2062585, 2065015, 2067339, 2069587,
+    2071777, 2073923, 2076032, 2078110, 2080160, 2082186, 2084189, 2086170, 2088131, 2090073, 2009637, 2015382, 2020990, 2026427,
+    2031639, 2036551, 2041088, 2045199, 2048876, 2052157, 2055112, 2057816, 2060340, 2062733, 2065034, 2067267, 2069447, 2071585,
+    2073688, 2075761, 2077808, 2079830, 2081830, 2083809, 2085767, 2003528, 2009293, 2014930, 2020410, 2025688, 2030698, 2035371,
+    2039651, 2043514, 2046979, 2050096, 2052935, 2055563, 2058037, 2060399, 2062679, 2064897, 2067067, 2069198, 2071295, 2073363,
+    2075406, 2077424, 2079421, 2081396, 1997396, 2003179, 2008841, 2014359, 2019692, 2024786, 2029576, 2034005, 2038039, 2041680,
+    2044963, 2047944, 2050688, 2053252, 2055684, 2058018, 2060279, 2062484, 2064644, 2066767, 2068858, 2070921, 2072959, 2074974,
+    2076967, 1991245, 1997043, 2002727, 2008277, 2013659, 2018824, 2023715, 2028274, 2032464, 2036273, 2039719, 2042848, 2045716,
+    2048381, 2050891, 2053287, 2055596, 2057839, 2060031, 2062182, 2064297, 2066382, 2068440, 2070473, 2072484, 1985076, 1990887,
+    1996591, 2002170, 2007594, 2012820, 2017798, 2022472, 2026802, 2030767, 2034372, 2037651, 2040650, 2043423, 2046021, 2048486,
+    2050850, 2053137, 2055364, 2057544, 2059685, 2061792, 2063871, 2065924, 2067953
+];
+/// QG_N, corner mos_ss: NMOS gate charge, same axes, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_ss.cir, block qg_n).
+pub const SS_QG_N : [625]u32 = [
+    2097152, 2095912, 2094680, 2093456, 2092239, 2091030, 2089830, 2088637, 2087453, 2086277, 2085110, 2083951, 2082800, 2081657,
+    2080523, 2079397, 2078279, 2077168, 2076066, 2074971, 2073883, 2072803, 2071729, 2070663, 2069604, 2100932, 2099685, 2098445,
+    2097213, 2095989, 2094772, 2093563, 2092363, 2091170, 2089986, 2088810, 2087643, 2086484, 2085333, 2084190, 2083056, 2081930,
+    2080812, 2079702, 2078599, 2077504, 2076416, 2075336, 2074263, 2073196, 2104688, 2103434, 2102187, 2100948, 2099716, 2098491,
+    2097275, 2096066, 2094865, 2093673, 2092488, 2091313, 2090145, 2088986, 2087835, 2086693, 2085559, 2084432, 2083314, 2082204,
+    2081101, 2080006, 2078919, 2077838, 2076765, 2108440, 2107170, 2105916, 2104670, 2103430, 2102198, 2100974, 2099757, 2098548,
+    2097347, 2096155, 2094971, 2093795, 2092627, 2091468, 2090317, 2089175, 2088041, 2086915, 2085796, 2084686, 2083583, 2082488,
+    2081401, 2080320, 2112291, 2110954, 2109687, 2108433, 2107187, 2105947, 2104715, 2103490, 2102273, 2101064, 2099863, 2098671,
+    2097486, 2096310, 2095143, 2093983, 2092833, 2091690, 2090556, 2089429, 2088311, 2087201, 2086098, 2085003, 2083915, 2116796,
+    2115115, 2113804, 2112538, 2111282, 2110034, 2108793, 2107559, 2106333, 2105115, 2103905, 2102704, 2101510, 2100325, 2099148,
+    2097980, 2096820, 2095669, 2094525, 2093390, 2092264, 2091145, 2090034, 2088931, 2087835, 2123209, 2120586, 2119067, 2117764,
+    2116491, 2115229, 2113976, 2112730, 2111493, 2110263, 2109042, 2107828, 2106624, 2105427, 2104239, 2103060, 2101889, 2100727,
+    2099573, 2098428, 2097291, 2096162, 2095042, 2093929, 2092824, 2131748, 2128026, 2125895, 2124456, 2123141, 2121855, 2120582,
+    2119320, 2118066, 2116821, 2115584, 2114356, 2113137, 2111927, 2110725, 2109533, 2108349, 2107173, 2106007, 2104849, 2103700,
+    2102559, 2101427, 2100303, 2099186, 2141589, 2137109, 2134073, 2132242, 2130811, 2129475, 2128172, 2126885, 2125610, 2124346,
+    2123092, 2121847, 2120611, 2119385, 2118167, 2116959, 2115760, 2114570, 2113390, 2112218, 2111055, 2109901, 2108755, 2107618,
+    2106490, 2152159, 2147220, 2143384, 2140887, 2139177, 2137732, 2136373, 2135050, 2133748, 2132460, 2131184, 2129919, 2128665,
+    2127421, 2126186, 2124961, 2123746, 2122540, 2121344, 2120157, 2118980, 2117811, 2116652, 2115502, 2114360, 2163180, 2157948,
+    2153544, 2150327, 2148144, 2146481, 2145019, 2143637, 2142294, 2140975, 2139674, 2138386, 2137110, 2135846, 2134593, 2133350,
+    2132118, 2130895, 2129683, 2128480, 2127287, 2126104, 2124930, 2123765, 2122609, 2174506, 2169075, 2164282, 2160457, 2157700,
+    2155692, 2154051, 2152573, 2151170, 2149809, 2148474, 2147158, 2145858, 2144572, 2143298, 2142036, 2140785, 2139545, 2138315,
+    2137096, 2135887, 2134688, 2133499, 2132319, 2131149, 2186054, 2180479, 2175412, 2171123, 2167811, 2165364, 2163460, 2161833,
+    2160341, 2158921, 2157543, 2156193, 2154864, 2153553, 2152256, 2150973, 2149702, 2148443, 2147195, 2145959, 2144733, 2143518,
+    2142313, 2141118, 2139933, 2197774, 2192088, 2186822, 2182185, 2178397, 2175486, 2173246, 2171411, 2169794, 2168292, 2166858,
+    2165466, 2164103, 2162762, 2161440, 2160133, 2158841, 2157562, 2156295, 2155041, 2153797, 2152565, 2151344, 2150133, 2148933,
+    2209629, 2203857, 2198439, 2193540, 2189364, 2186014, 2183400, 2181303, 2179519, 2177911, 2176404, 2174958, 2173554, 2172180,
+    2170829, 2169497, 2168182, 2166881, 2165594, 2164321, 2163059, 2161810, 2160571, 2159345, 2158129, 2221596, 2215753, 2210216,
+    2205115, 2200629, 2196892, 2193899, 2191503, 2189511, 2187766, 2186166, 2184656, 2183203, 2181790, 2180406, 2179046, 2177706,
+    2176382, 2175074, 2173781, 2172500, 2171233, 2169977, 2168733, 2167501, 2233656, 2227755, 2222121, 2216862, 2212127, 2208060,
+    2204708, 2201995, 2199762, 2197849, 2196135, 2194545, 2193034, 2191576, 2190155, 2188764, 2187397, 2186049, 2184718, 2183404,
+    2182104, 2180817, 2179544, 2178283, 2177034, 2245794, 2239844, 2234132, 2228744, 2223809, 2219464, 2215787, 2212758, 2210262,
+    2208153, 2206301, 2204614, 2203033, 2201523, 2200061, 2198636, 2197239, 2195865, 2194511, 2193174, 2191854, 2190548, 2189256,
+    2187978, 2186712, 2257999, 2252007, 2246229, 2240737, 2235638, 2231060, 2227096, 2223764, 2220996, 2218669, 2216657, 2214855,
+    2213191, 2211618, 2210109, 2208646, 2207217, 2205815, 2204436, 2203077, 2201735, 2200409, 2199098, 2197802, 2196519, 2270262,
+    2264234, 2258401, 2252820, 2247587, 2242813, 2238598, 2234986, 2231947, 2229388, 2227194, 2225258, 2223496, 2221852, 2220289,
+    2218782, 2217318, 2215886, 2214481, 2213097, 2211734, 2210388, 2209058, 2207742, 2206441, 2282575, 2276516, 2270635, 2264981,
+    2259635, 2254697, 2250264, 2246397, 2243095, 2240298, 2237907, 2235819, 2233944, 2232216, 2230590, 2229034, 2227531, 2226066,
+    2224633, 2223225, 2221838, 2220471, 2219121, 2217787, 2216468, 2294932, 2288846, 2282924, 2277207, 2271765, 2266687, 2262066,
+    2257970, 2254422, 2251387, 2248788, 2246532, 2244528, 2242704, 2241005, 2239394, 2237846, 2236346, 2234882, 2233447, 2232037,
+    2230648, 2229277, 2227924, 2226587, 2307327, 2301218, 2295260, 2289489, 2283965, 2278768, 2273983, 2269684, 2265906, 2262639,
+    2259829, 2257393, 2255246, 2253311, 2251528, 2249853, 2248256, 2246716, 2245219, 2243756, 2242320, 2240909, 2239517, 2238144,
+    2236788, 2319756, 2313626, 2307636, 2301818, 2296224, 2290925, 2285999, 2281518, 2277530, 2274042, 2271019, 2268395, 2266092,
+    2264034, 2262156, 2260408, 2258754, 2257169, 2255636, 2254142, 2252680, 2251245, 2249832, 2248439, 2247064, 2332214, 2326067,
+    2320049, 2314190, 2308535, 2303147, 2298097, 2293457, 2289278, 2285581, 2282349, 2279533, 2277065, 2274871, 2272886, 2271055,
+    2269337, 2267701, 2266127, 2264600, 2263110, 2261649, 2260214, 2258800, 2257406
+];
+/// QD_P, corner mos_ss: PMOS drain charge, NEGATED, Vsg = ix * 62.5 mV, Vsd = iy * 62.5 mV, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_ss.cir, block qd_p).
+pub const SS_QD_P : [625]u32 = [
+    2097152, 2102577, 2107904, 2113141, 2118291, 2123359, 2128350, 2133267, 2138112, 2142888, 2147599, 2152246, 2156833, 2161361,
+    2165832, 2170250, 2174615, 2178929, 2183196, 2187415, 2191590, 2195721, 2199810, 2203860, 2207870, 2094175, 2099628, 2104984,
+    2110249, 2115430, 2120530, 2125552, 2130500, 2135377, 2140185, 2144927, 2149605, 2154222, 2158780, 2163280, 2167726, 2172119,
+    2176461, 2180753, 2184998, 2189197, 2193353, 2197465, 2201537, 2205569, 2091172, 2096651, 2102034, 2107328, 2112538, 2117668,
+    2122721, 2127701, 2132609, 2137449, 2142223, 2146932, 2151580, 2156168, 2160699, 2165174, 2169595, 2173965, 2178284, 2182555,
+    2186780, 2190960, 2195096, 2199191, 2203245, 2088137, 2093645, 2099055, 2104376, 2109615, 2114774, 2119858, 2124868, 2129808,
+    2134680, 2139485, 2144226, 2148905, 2153524, 2158085, 2162590, 2167041, 2171439, 2175786, 2180084, 2184335, 2188540, 2192701,
+    2196820, 2200897, 2085030, 2090595, 2096032, 2101380, 2106646, 2111834, 2116947, 2121988, 2126959, 2131862, 2136698, 2141471,
+    2146182, 2150833, 2155424, 2159960, 2164440, 2168867, 2173243, 2177569, 2181847, 2186078, 2190265, 2194408, 2198509, 2081575,
+    2087387, 2092867, 2098242, 2103535, 2108751, 2113893, 2118963, 2123964, 2128898, 2133766, 2138571, 2143313, 2147995, 2152618,
+    2157184, 2161695, 2166152, 2170556, 2174911, 2179216, 2183475, 2187687, 2191856, 2195982, 2076455, 2083402, 2089034, 2094449,
+    2099772, 2105016, 2110186, 2115286, 2120317, 2125281, 2130181, 2135016, 2139790, 2144504, 2149158, 2153755, 2158296, 2162783,
+    2167217, 2171600, 2175934, 2180219, 2184458, 2188652, 2192803, 2067534, 2077060, 2083349, 2088879, 2094250, 2099530, 2104733,
+    2109865, 2114928, 2119925, 2124857, 2129725, 2134532, 2139278, 2143965, 2148594, 2153167, 2157685, 2162150, 2166563, 2170925,
+    2175239, 2179505, 2183726, 2187902, 2055153, 2067260, 2075285, 2081231, 2086722, 2092065, 2097315, 2102487, 2107588, 2112621,
+    2117589, 2122494, 2127336, 2132118, 2136840, 2141504, 2146111, 2150663, 2155161, 2159607, 2164001, 2168346, 2172642, 2176892,
+    2181096, 2040833, 2054574, 2064893, 2071961, 2077793, 2083270, 2088596, 2093824, 2098972, 2104048, 2109057, 2114001, 2118882,
+    2123702, 2128462, 2133164, 2137808, 2142397, 2146930, 2151411, 2155839, 2160217, 2164546, 2168827, 2173062, 2025466, 2040150,
+    2052335, 2061176, 2067808, 2073597, 2079071, 2084388, 2089601, 2094730, 2099786, 2104774, 2109698, 2114559, 2119359, 2124100,
+    2128784, 2133410, 2137981, 2142498, 2146963, 2151376, 2155739, 2160054, 2164321, 2009484, 2024740, 2038174, 2048811, 2056752,
+    2063172, 2068941, 2074416, 2079730, 2084932, 2090048, 2095088, 2100059, 2104965, 2109809, 2114592, 2119316, 2123982, 2128593,
+    2133148, 2137650, 2142100, 2146499, 2150848, 2155150, 1993114, 2008736, 2022983, 2035049, 2044505, 2051926, 2058225, 2063983,
+    2069462, 2074775, 2079972, 2085078, 2090105, 2095062, 2099953, 2104782, 2109549, 2114258, 2118910, 2123505, 2128046, 2132535,
+    2136971, 2141356, 2145693, 1976484, 1992355, 2007141, 2020243, 2031091, 2039741, 2046842, 2053063, 2058814, 2064297, 2069612,
+    2074806, 2079905, 2084923, 2089868, 2094747, 2099562, 2104316, 2109011, 2113649, 2118232, 2122760, 2127235, 2131659, 2136032,
+    1959674, 1975720, 1990878, 2004712, 2016691, 2026583, 2034689, 2041578, 2047743, 2053488, 2058976, 2064294, 2069487, 2074582,
+    2079593, 2084529, 2089397, 2094201, 2098944, 2103627, 2108253, 2112823, 2117340, 2121803, 2126215, 1942740, 1958910, 1974333,
+    1988689, 2001535, 2012537, 2021716, 2029442, 2036181, 2042303, 2048042, 2053535, 2058856, 2064050, 2069143, 2074149, 2079078,
+    2083938, 2088733, 2093465, 2098138, 2102754, 2107313, 2111819, 2116271, 1925717, 1941978, 1957594, 1972328, 1985829, 1997751,
+    2007956, 2016605, 2024056, 2030681, 2036768, 2042502, 2047999, 2053324, 2058520, 2063612, 2068615, 2073540, 2078393, 2083180,
+    2087904, 2092568, 2097173, 2101723, 2106219, 1908635, 1924961, 1940720, 1955735, 1969727, 1982387, 1993500, 2003072, 2011321,
+    2018564, 2025101, 2031158, 2036887, 2042386, 2047716, 2052915, 2058007, 2063008, 2067929, 2072778, 2077558, 2082275, 2086931,
+    2091529, 2096070, 1891512, 1907886, 1923753, 1938978, 1953344, 1966583, 1978466, 1988898, 1997969, 2005910, 2012991, 2019455,
+    2025486, 2031210, 2036712, 2042046, 2047248, 2052341, 2057342, 2062260, 2067105, 2071881, 2076592, 2081241, 2085832, 1874367,
+    1890773, 1906721, 1922107, 1936758, 1950453, 1962969, 1974164, 1984031, 1992705, 2000402, 2007353, 2013756, 2019764, 2025483,
+    2030986, 2036324, 2041529, 2046625, 2051626, 2056543, 2061385, 2066157, 2070863, 2075507, 1857210, 1873638, 1889648, 1905157,
+    1920030, 1934081, 1947113, 1958963, 1969563, 1978965, 1987317, 1994818, 2001662, 2008014, 2014000, 2019714, 2025218, 2030560,
+    2035768, 2040866, 2045868, 2050785, 2055625, 2060394, 2065097, 1840052, 1856493, 1872549, 1888154, 1903199, 1917533, 1930983,
+    1943385, 1954637, 1964729, 1973745, 1981836, 1989175, 1995928, 2002235, 2008202, 2013910, 2019415, 2024760, 2029972, 2035072,
+    2040076, 2044993, 2049832, 2054600, 1822902, 1839349, 1855438, 1871117, 1886298, 1900857, 1914643, 1927506, 1939326, 1950051,
+    1959711, 1968407, 1976280, 1983483, 1990159, 1996426, 2002375, 2008077, 2013583, 2018930, 2024145, 2029249, 2034254, 2039172,
+    2044012, 1805765, 1822213, 1838325, 1854062, 1869350, 1884088, 1898148, 1911392, 1923699, 1934991, 1945256, 1954549, 1962975,
+    1970664, 1977753, 1984361, 1990592, 1996524, 2002221, 2007726, 2013076, 2018294, 2023400, 2028408, 2033328, 1788649, 1805091,
+    1821220, 1837000, 1852373, 1867255, 1881537, 1895097, 1907817, 1919608, 1930428, 1940293, 1949271, 1957468, 1965001, 1971990,
+    1978540, 1984738, 1990654, 1996344, 2001849, 2007201, 2012422, 2017531, 2022541
+];
+/// QG_P, corner mos_ss: PMOS gate charge, NEGATED, same axes, at index ix * 25 + iy, units of (Q / 500 fF) in 2^-24 V plus 2^21 (deck psp103_q_ss.cir, block qg_p).
+pub const SS_QG_P : [625]u32 = [
+    2097152, 2094203, 2091286, 2088400, 2085544, 2082719, 2079926, 2077165, 2074435, 2071737, 2069070, 2066434, 2063829, 2061254,
+    2058707, 2056190, 2053700, 2051237, 2048800, 2046388, 2044002, 2041639, 2039299, 2036982, 2034687, 2108679, 2105703, 2102757,
+    2099841, 2096955, 2094099, 2091274, 2088481, 2085720, 2082990, 2080292, 2077625, 2074989, 2072383, 2069808, 2067261, 2064743,
+    2062253, 2059790, 2057353, 2054941, 2052555, 2050192, 2047852, 2045535, 2120117, 2117114, 2114141, 2111196, 2108280, 2105394,
+    2102538, 2099713, 2096920, 2094158, 2091428, 2088730, 2086063, 2083427, 2080821, 2078245, 2075699, 2073181, 2070690, 2068227,
+    2065790, 2063378, 2060991, 2058628, 2056288, 2131486, 2128452, 2125452, 2122480, 2119535, 2116620, 2113733, 2110877, 2108053,
+    2105259, 2102497, 2099767, 2097069, 2094402, 2091765, 2089160, 2086584, 2084037, 2081519, 2079028, 2076565, 2074128, 2071716,
+    2069329, 2066965, 2142864, 2139764, 2136736, 2133737, 2130765, 2127820, 2124904, 2122018, 2119162, 2116337, 2113544, 2110782,
+    2108052, 2105353, 2102686, 2100050, 2097444, 2094868, 2092321, 2089803, 2087312, 2084849, 2082411, 2080000, 2077612, 2154734,
+    2151323, 2148247, 2145220, 2142221, 2139249, 2136305, 2133389, 2130503, 2127647, 2124822, 2122029, 2119268, 2116538, 2113840,
+    2111173, 2108537, 2105931, 2103356, 2100809, 2098292, 2095802, 2093338, 2090902, 2088490, 2169381, 2164526, 2161264, 2158194,
+    2155166, 2152166, 2149194, 2146250, 2143335, 2140450, 2137595, 2134772, 2131980, 2129220, 2126492, 2123795, 2121130, 2118496,
+    2115892, 2113319, 2110774, 2108259, 2105770, 2103309, 2100875, 2190525, 2182413, 2178354, 2175152, 2172071, 2169035, 2166030,
+    2163055, 2160109, 2157193, 2154307, 2151453, 2148630, 2145840, 2143081, 2140354, 2137659, 2134996, 2132364, 2129762, 2127191,
+    2124649, 2122136, 2119650, 2117192, 2217615, 2206232, 2200065, 2196373, 2193156, 2190051, 2186998, 2183981, 2180997, 2178045,
+    2175123, 2172233, 2169375, 2166549, 2163756, 2160995, 2158267, 2155571, 2152906, 2150273, 2147671, 2145099, 2142557, 2140043,
+    2137558, 2248041, 2234551, 2225548, 2220522, 2216910, 2213654, 2210516, 2207440, 2204406, 2201409, 2198446, 2195515, 2192618,
+    2189753, 2186922, 2184123, 2181358, 2178625, 2175925, 2173257, 2170621, 2168016, 2165442, 2162897, 2160381, 2280261, 2265538,
+    2254167, 2246993, 2242441, 2238825, 2235521, 2232346, 2229242, 2226188, 2223174, 2220198, 2217257, 2214350, 2211477, 2208638,
+    2205833, 2203062, 2200324, 2197619, 2194946, 2192306, 2189697, 2187119, 2184570, 2313526, 2298052, 2285056, 2275653, 2269534,
+    2265182, 2261539, 2258184, 2254969, 2251835, 2248757, 2245726, 2242734, 2239781, 2236863, 2233981, 2231134, 2228322, 2225543,
+    2222799, 2220088, 2217411, 2214765, 2212152, 2209569, 2347445, 2331488, 2317416, 2306194, 2298226, 2292686, 2288425, 2284748,
+    2281345, 2278088, 2274921, 2271818, 2268767, 2265759, 2262791, 2259862, 2256970, 2254114, 2251293, 2248508, 2245757, 2243040,
+    2240356, 2237706, 2235087, 2381792, 2365510, 2350715, 2338149, 2328442, 2321417, 2316210, 2311995, 2308281, 2304831, 2301533,
+    2298333, 2295204, 2292130, 2289105, 2286122, 2283180, 2280277, 2277411, 2274581, 2271788, 2269030, 2266306, 2263617, 2260961,
+    2416427, 2399921, 2384624, 2371091, 2359944, 2351388, 2344978, 2339977, 2335783, 2332034, 2328540, 2325202, 2321968, 2318811,
+    2315715, 2312670, 2309671, 2306716, 2303800, 2300924, 2298085, 2295283, 2292517, 2289787, 2287091, 2451256, 2434592, 2418937,
+    2404707, 2392436, 2382490, 2374769, 2368771, 2363904, 2359720, 2355938, 2352402, 2349025, 2345758, 2342572, 2339452, 2336388,
+    2333372, 2330402, 2327473, 2324586, 2321737, 2318926, 2316152, 2313414, 2486213, 2469438, 2453523, 2438782, 2425653, 2414534,
+    2405542, 2398426, 2392714, 2387940, 2383757, 2379945, 2376370, 2372956, 2369656, 2366442, 2363297, 2360212, 2357179, 2354193,
+    2351251, 2348352, 2345492, 2342672, 2339890, 2521251, 2504399, 2488292, 2473171, 2459389, 2447317, 2437183, 2428932, 2422260,
+    2416755, 2412047, 2407862, 2404021, 2400410, 2396960, 2393627, 2390385, 2387217, 2384110, 2381058, 2378056, 2375101, 2372189,
+    2369318, 2366488, 2556333, 2539429, 2523179, 2507772, 2493491, 2480661, 2469546, 2460225, 2452548, 2446207, 2440857, 2436197,
+    2432009, 2428141, 2424496, 2421012, 2417648, 2414378, 2411185, 2408057, 2404986, 2401968, 2398998, 2396073, 2393191, 2591432,
+    2574494, 2558139, 2542514, 2527847, 2514417, 2502485, 2492202, 2483543, 2476313, 2470227, 2464994, 2460371, 2456177, 2452285,
+    2448609, 2445092, 2441698, 2438401, 2435183, 2432033, 2428944, 2425908, 2422923, 2419984, 2626525, 2609570, 2593137, 2577345,
+    2562379, 2548472, 2535868, 2524750, 2515176, 2507053, 2500175, 2494288, 2489148, 2484553, 2480353, 2476437, 2472732, 2469185,
+    2465762, 2462437, 2459195, 2456024, 2452915, 2449863, 2446862, 2661594, 2644634, 2628145, 2612225, 2597025, 2582740, 2569588,
+    2557758, 2547360, 2538380, 2530693, 2524097, 2518369, 2513303, 2508731, 2504524, 2500587, 2496854, 2493278, 2489827, 2486476,
+    2483211, 2480018, 2476890, 2473820, 2696627, 2679671, 2663143, 2647125, 2631740, 2617154, 2603555, 2591127, 2580005, 2570232,
+    2561750, 2554420, 2548052, 2542453, 2537448, 2532894, 2528679, 2524723, 2520965, 2517362, 2513883, 2510508, 2507218, 2504004,
+    2500857, 2731611, 2714666, 2698114, 2682022, 2666492, 2651663, 2637701, 2624774, 2613025, 2602535, 2593301, 2585237, 2578201,
+    2572021, 2566529, 2561574, 2557034, 2552813, 2548838, 2545056, 2541427, 2537923, 2534523, 2531211, 2527975, 2766536, 2749610,
+    2733045, 2716897, 2701253, 2686230, 2671970, 2658627, 2646342, 2635216, 2625286, 2616514, 2608803, 2602013, 2595989, 2590586,
+    2585674, 2581145, 2576917, 2572925, 2569120, 2565467, 2561939, 2558515, 2555180
+];
+
 ; --- the accessors and the tests ------------------------------------------------------------------------
 
 /// The fold of the first n entries of t after the state s. Order matters.
@@ -1237,3 +1690,99 @@ test the_grid_and_the_units_of_the_numbers_are_consistent {
     const d : f64 = amps - 0.000402305;
     assert d * d < 0.000000001 * 0.000000001;
 }
+
+test the_charge_tables_are_the_numbers_ngspice_printed {
+    // A corner's charge tables are qd_n, qg_n, qd_p, qg_p, in that order.
+    var s : u64 = fold_arr(&TT_QD_N, 625, 0);
+    s = fold_arr(&TT_QG_N, 625, s);
+    s = fold_arr(&TT_QD_P, 625, s);
+    s = fold_arr(&TT_QG_P, 625, s);
+    assert_eq(s, PIN_TT_Q);
+    s = fold_arr(&SS_QD_N, 625, 0);
+    s = fold_arr(&SS_QG_N, 625, s);
+    s = fold_arr(&SS_QD_P, 625, s);
+    s = fold_arr(&SS_QG_P, 625, s);
+    assert_eq(s, PIN_SS_Q);
+    // A pin sees one changed number.
+    var t : [625]u32 = [0] ** 625;
+    var i : usize = 0;
+    while (i < 625) {
+        t[i] = TT_QD_N[i];
+        i = i + 1;
+    }
+    t[300] = t[300] + 1;
+    assert fold_arr(&t, 625, 0) != fold_arr(&TT_QD_N, 625, 0);
+}
+
+test the_charge_tables_are_relative_to_the_origin_and_fit_the_reader {
+    // Variables, not constants: a constant assert is folded by the compiler and counts as no assert at all.
+    var off : u32 = Q_OFFSET;
+    var cref : u32 = C_REF_FF;
+    assert_eq(cref, 500);
+    assert_eq(off, 2097152);
+    // Entry 0 is the grid origin (Vgs = Vds = 0), the charge relative to which every other entry is taken: exactly the offset.
+    assert_eq(TT_QD_N[0], off);
+    assert_eq(TT_QG_N[0], off);
+    assert_eq(TT_QD_P[0], off);
+    assert_eq(TT_QG_P[0], off);
+    assert_eq(SS_QD_N[0], off);
+    assert_eq(SS_QG_N[0], off);
+    assert_eq(SS_QD_P[0], off);
+    assert_eq(SS_QG_P[0], off);
+    // No entry is below 10^6 or above 4 * 10^6: the offset is large against the most negative charge (a few fC) and small against the
+    // 2^23 the solver's u64 interpolation allows.
+    var j : u32 = 0;
+    var lo : u32 = off;
+    var hi : u32 = off;
+    while (j < 625) {
+        const a : u32 = TT_QD_N[j as usize];
+        const b : u32 = TT_QG_N[j as usize];
+        const c : u32 = TT_QD_P[j as usize];
+        const d : u32 = TT_QG_P[j as usize];
+        const e : u32 = SS_QD_N[j as usize];
+        const f : u32 = SS_QG_N[j as usize];
+        const g : u32 = SS_QD_P[j as usize];
+        const h : u32 = SS_QG_P[j as usize];
+        assert a > 1000000;
+        assert b > 1000000;
+        assert c > 1000000;
+        assert d > 1000000;
+        assert e > 1000000;
+        assert f > 1000000;
+        assert g > 1000000;
+        assert h > 1000000;
+        assert a < 4000000;
+        assert b < 4000000;
+        assert c < 4000000;
+        assert d < 4000000;
+        assert e < 4000000;
+        assert f < 4000000;
+        assert g < 4000000;
+        assert h < 4000000;
+        if (a < lo) { lo = a; }
+        if (d > hi) { hi = d; }
+        j = j + 1;
+    }
+    // The most negative entry is the tt NMOS drain charge at (1.5 V, 0): 1977380 - 2097152 = -119772 units = -3.6 fC; the largest is
+    // the tt PMOS gate charge (negated) at (1.5 V, 1.5 V): 2802131 - 2097152 = 704979 units = 21 fC.
+    assert_eq(lo, 1977380);
+    assert_eq(hi, 2802131);
+}
+
+test the_500_fF_waveform_is_the_numbers_ngspice_printed {
+    assert_eq(fold_arr(&INV500_BE, 401, 0), PIN_TRAN500);
+    // The same circuit at rest as the 10 pF run (11 units, 0.66 uV), before step 0 and after step 399: the output goes down and comes back.
+    assert_eq(INV500_BE[0], INV10_BE[0]);
+    assert_eq(INV500_BE[400], INV10_BE[400]);
+    // After step 0 the output has left the rail by 689206 units (41 mV): the PMOS charges 500 fF, and the gate edge couples charge in.
+    assert_eq(INV500_BE[1], 689206);
+    // A pin sees one changed number.
+    var t : [401]u32 = [0] ** 401;
+    var i : usize = 0;
+    while (i < 401) {
+        t[i] = INV500_BE[i];
+        i = i + 1;
+    }
+    t[200] = t[200] + 1;
+    assert fold_arr(&t, 401, 0) != fold_arr(&INV500_BE, 401, 0);
+}
diff --git a/specs/analog/mos_table.t27 b/specs/analog/mos_table.t27
index ec64300001..4667678abe 100644
--- a/specs/analog/mos_table.t27
+++ b/specs/analog/mos_table.t27
@@ -4,6 +4,13 @@
 ; the first cell with a device that is not a resistor. M3b (specs/analog/mos_psp103.t27, #10516) runs this same solver on
 ; tables made from the IHP SG13G2 PSP103 model, one NMOS and one PMOS table per corner; for that the solver takes the NMOS
 ; table and the PMOS table as two arguments (M3a passes the same table twice). Everything hand-written is this file.
+; M3c (#10777) generalizes the step once more, with the drain-charge tables of the two devices: the `*_q` functions (gq_pos, gq_neg, gq_deriv,
+; g_resid_q, newton_iter_q, inv_step_q, board_step_q, run_fix_q, and newton_iter_qf, run_tabf_q on the f64 side) take two more tables qn, qp
+; and the gate of the previous step A0, and the current-only functions of M3a and M3b are those with the zero tables NOQ, so nothing is copied and
+; every pinned number of this file is still a number of the charge path. A charge table is read and sloped like a current table (id_of, gds_cell),
+; nondecreasing along Vds. The tests of the generalization are the last three: zero charge is the old path bit for bit; a plane of charge (no
+; current) gives the closed form of the gate-edge kick and the load, for an NMOS-side and a PMOS-side charge; the mirror holds with the four tables
+; swapped. mos_psp103.t27 runs it on the PSP103 tables.
 ; -
 ; THE CIRCUIT. The reference is mirekez/cpphdl, analog/examples/cmos_inverter: Vdd = 1.2 V, an NMOS from Y to ground and a
 ; PMOS from Y to Vdd, both gates on A, A = 0 V and then 1.2 V from 2 us, Y loaded by 10 pF, dt = 10 ns, stop at 4 us (400
@@ -248,9 +255,38 @@ pub fn mag_diff(p: u32, n: u32) -> u32 {
     return n - p;
 }
 
-/// |g(y)|, the residual of the backward Euler equation at y, in units. tn is the NMOS table, tp the PMOS table.
+/// A charge table of zeros: a device that stores no charge. Passed for qn and qp it turns the charge-aware functions below
+/// (`*_q`) into the current-only ones of M3a and M3b, with every sum equal to what it was (a zero entry adds 0, a zero slope adds 0).
+pub const NOQ : [625]u32 = [0] ** 625;
+
+/// The positive terms of the charge-aware g at y: y + In(A, y) + qn(A, y) + qp(Vdd - A0, Vdd - y0). qn is the NMOS drain-charge
+/// table and qp the PMOS drain-charge table, both in the units of the current tables (the charge over C, in units of 2^-24 V) and
+/// both nondecreasing along the Vds axis (qp is stored negated, so that it grows with Vsd as qn grows with Vds). A0 is the gate of
+/// the previous step. The node's charge is C y + Qn - Qp, and the backward Euler equation is
+/// g(y) = (y - y0) + (qn(A, y) - qn(A0, y0)) - (qp(Vdd - A, Vdd - y) - qp(Vdd - A0, Vdd - y0)) + In - Ip = 0.
+pub fn gq_pos(tn: []const u32, qn: []const u32, qp: []const u32, y: u32, y0: u32, a: u32, a0: u32) -> u32 {
+    return g_pos(tn, y, a) + id_of(qn, a, y) + id_of(qp, VDD - a0, VDD - y0);
+}
+
+/// The negative terms of the charge-aware g at y: y0 + Ip(Vdd - A, Vdd - y) + qn(A0, y0) + qp(Vdd - A, Vdd - y).
+pub fn gq_neg(tp: []const u32, qn: []const u32, qp: []const u32, y: u32, y0: u32, a: u32, a0: u32) -> u32 {
+    return g_neg(tp, y, y0, a) + id_of(qn, a0, y0) + id_of(qp, VDD - a, VDD - y);
+}
+
+/// g'(y) in units of 2^20: one cell for the node's own C, plus the Vds slopes of the four tables, each read off its own table
+/// (the value and its derivative are one object, FR-002).
+pub fn gq_deriv(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y: u32, a: u32) -> u64 {
+    return (CELL as u64) + (gds_cell(tn, a, y) as u64) + (gds_cell(tp, VDD - a, VDD - y) as u64) + (gds_cell(qn, a, y) as u64) + (gds_cell(qp, VDD - a, VDD - y) as u64);
+}
+
+/// |g(y)| of the charge-aware equation, in units.
+pub fn g_resid_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y: u32, y0: u32, a: u32, a0: u32) -> u32 {
+    return mag_diff(gq_pos(tn, qn, qp, y, y0, a, a0), gq_neg(tp, qn, qp, y, y0, a, a0));
+}
+
+/// |g(y)|, the residual of the backward Euler equation at y, in units. tn is the NMOS table, tp the PMOS table. No charge.
 pub fn g_resid(tn: []const u32, tp: []const u32, y: u32, y0: u32, a: u32) -> u32 {
-    return mag_diff(g_pos(tn, y, a), g_neg(tp, y, y0, a));
+    return g_resid_q(tn, tp, &NOQ, &NOQ, y, y0, a, a);
 }
 
 /// round(g * 2^20 / d): the Newton move for a residual of magnitude g and a derivative of d / 2^20. d >= 2^20.
@@ -259,11 +295,13 @@ pub fn newton_move(g: u32, d: u64) -> u32 {
 }
 
 /// One Newton iteration of backward Euler at the output node: from y toward the root of g, clamped to [0, Vdd]. The NMOS
-/// reads tn at (A, Y), the PMOS reads tp at (Vdd - A, Vdd - Y); M3a passes one table twice, M3b (mos_psp103.t27) two.
-pub fn newton_iter(tn: []const u32, tp: []const u32, y: u32, y0: u32, a: u32) -> u32 {
-    const p : u32 = g_pos(tn, y, a);
-    const n : u32 = g_neg(tp, y, y0, a);
-    const d : u64 = (CELL as u64) + (gds_cell(tn, a, y) as u64) + (gds_cell(tp, VDD - a, VDD - y) as u64);
+/// reads tn at (A, Y) and the PMOS reads tp at (Vdd - A, Vdd - Y); the drain charges read qn at (A, Y) and qp at (Vdd - A, Vdd - Y),
+/// and the previous step's at (A0, y0) and (Vdd - A0, Vdd - y0). M3a passes one table twice, M3b (mos_psp103.t27) two, M3c adds
+/// the charge tables.
+pub fn newton_iter_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y: u32, y0: u32, a: u32, a0: u32) -> u32 {
+    const p : u32 = gq_pos(tn, qn, qp, y, y0, a, a0);
+    const n : u32 = gq_neg(tp, qn, qp, y, y0, a, a0);
+    const d : u64 = gq_deriv(tn, tp, qn, qp, y, a);
     if (p >= n) {
         const down : u32 = newton_move(p - n, d);
         if (down >= y) { return 0; }
@@ -274,29 +312,44 @@ pub fn newton_iter(tn: []const u32, tp: []const u32, y: u32, y0: u32, a: u32) ->
     return y + up;
 }
 
-/// One backward Euler step from y0 with the gate at a: k Newton iterations.
-pub fn inv_step(tn: []const u32, tp: []const u32, y0: u32, a: u32, k: u32) -> u32 {
+/// One Newton iteration without charge: `newton_iter_q` with the zero charge tables.
+pub fn newton_iter(tn: []const u32, tp: []const u32, y: u32, y0: u32, a: u32) -> u32 {
+    return newton_iter_q(tn, tp, &NOQ, &NOQ, y, y0, a, a);
+}
+
+/// One backward Euler step from y0 with the gate at a, coming from the gate a0: k Newton iterations.
+pub fn inv_step_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y0: u32, a: u32, a0: u32, k: u32) -> u32 {
     var y : u32 = y0;
     var i : u32 = 0;
     while (i < k) {
-        y = newton_iter(tn, tp, y, y0, a);
+        y = newton_iter_q(tn, tp, qn, qp, y, y0, a, a0);
         i = i + 1;
     }
     return y;
 }
 
-/// The step the board runs: K_NEWTON iterations. The trip count is a constant, so the loop unrolls (`t27b verilog` takes
-/// a `while` only with a static bound); inv_step above takes k as an argument so the tests can run K = 1 and K = 2.
-pub fn board_step(tn: []const u32, tp: []const u32, y0: u32, a: u32) -> u32 {
+/// One backward Euler step without charge: `inv_step_q` with the zero charge tables.
+pub fn inv_step(tn: []const u32, tp: []const u32, y0: u32, a: u32, k: u32) -> u32 {
+    return inv_step_q(tn, tp, &NOQ, &NOQ, y0, a, a, k);
+}
+
+/// The step the board runs, with charge: K_NEWTON iterations. The trip count is a constant, so the loop unrolls (`t27b verilog`
+/// takes a `while` only with a static bound); inv_step_q above takes k as an argument so the tests can run K = 1 and K = 2.
+pub fn board_step_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y0: u32, a: u32, a0: u32) -> u32 {
     var y : u32 = y0;
     var i : u32 = 0;
     while (i < K_NEWTON) {
-        y = newton_iter(tn, tp, y, y0, a);
+        y = newton_iter_q(tn, tp, qn, qp, y, y0, a, a0);
         i = i + 1;
     }
     return y;
 }
 
+/// The step the board runs without charge: `board_step_q` with the zero charge tables.
+pub fn board_step(tn: []const u32, tp: []const u32, y0: u32, a: u32) -> u32 {
+    return board_step_q(tn, tp, &NOQ, &NOQ, y0, a, a);
+}
+
 ; --- the f64 side: references that only the tests reach -------------------------------------------------
 ; A module-level f64 fn is carried to Verilog only when the hardware reaches it; nothing below is reached by `step`.
 
@@ -412,18 +465,48 @@ pub fn monotone(t: []const u32) -> bool {
     return ok;
 }
 
-/// The fixed-point run: 400 steps of `inv_step` with k iterations from Y = y_init, the gate mirrored (Vdd - A) when
-/// mirror is set. out gets every Y, res the last residual of every step. Returns the fold of the 400 samples.
-pub fn run_fix(tn: []const u32, tp: []const u32, k: u32, y_init: u32, mirror: bool, out: []u32, res: []u32) -> u64 {
+/// Is the table nondecreasing along one axis over all 600 neighbouring pairs of that axis: ax 0 along the first index (Vgs), ax 1
+/// along the second (Vds)? A drain charge table is nondecreasing along Vds (a capacitance is positive), a gate charge table along Vgs.
+pub fn monotone_along(t: []const u32, ax: u32) -> bool {
+    var ok : bool = true;
+    var ix : u32 = 0;
+    while (ix < NODES) {
+        var iy : u32 = 0;
+        while (iy < NODES) {
+            const here : u32 = t[(ix * NODES + iy) as usize];
+            if (ax == 0 and ix + 1 < NODES) {
+                if (t[((ix + 1) * NODES + iy) as usize] < here) { ok = false; }
+            }
+            if (ax == 1 and iy + 1 < NODES) {
+                if (t[(ix * NODES + iy + 1) as usize] < here) { ok = false; }
+            }
+            iy = iy + 1;
+        }
+        ix = ix + 1;
+    }
+    return ok;
+}
+
+/// The gate of step n, mirrored (Vdd - A) when mirror is set.
+pub fn gate_in(n: u32, mirror: bool) -> u32 {
+    if (mirror) { return VDD - gate_at(n); }
+    return gate_at(n);
+}
+
+/// The fixed-point run: 400 steps of `inv_step_q` with k iterations from Y = y_init, the gate mirrored (Vdd - A) when
+/// mirror is set, the gate before step 0 at a_rest. out gets every Y, res the last residual of every step. Returns the fold of
+/// the 400 samples.
+pub fn run_fix_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, k: u32, y_init: u32, a_rest: u32, mirror: bool, out: []u32, res: []u32) -> u64 {
     var y : u32 = y_init;
     var s : u64 = 0;
     var n : u32 = 0;
     while (n < RUN_STEPS) {
-        var a : u32 = gate_at(n);
-        if (mirror) { a = VDD - a; }
+        const a : u32 = gate_in(n, mirror);
+        var a0 : u32 = a_rest;
+        if (n > 0) { a0 = gate_in(n - 1, mirror); }
         const y0 : u32 = y;
-        y = inv_step(tn, tp, y0, a, k);
-        res[n as usize] = g_resid(tn, tp, y, y0, a);
+        y = inv_step_q(tn, tp, qn, qp, y0, a, a0, k);
+        res[n as usize] = g_resid_q(tn, tp, qn, qp, y, y0, a, a0);
         out[n as usize] = y;
         s = fold(s, y);
         n = n + 1;
@@ -431,6 +514,11 @@ pub fn run_fix(tn: []const u32, tp: []const u32, k: u32, y_init: u32, mirror: bo
     return s;
 }
 
+/// The fixed-point run without charge: `run_fix_q` with the zero charge tables (the gate before step 0 does not matter then).
+pub fn run_fix(tn: []const u32, tp: []const u32, k: u32, y_init: u32, mirror: bool, out: []u32, res: []u32) -> u64 {
+    return run_fix_q(tn, tp, &NOQ, &NOQ, k, y_init, gate_in(0, mirror), mirror, out, res);
+}
+
 /// The largest residual of a run, in units.
 pub fn max_of(r: []const u32) -> u32 {
     var m : u32 = 0;
@@ -474,27 +562,35 @@ pub fn gdsf(t: []const u32, vg: f64, vd: f64) -> f64 {
     return (d0 * (CELL_F - fx) + d1 * fx) / CELL_F;
 }
 
-/// newton_iter in real arithmetic: the same equation, the same slopes, the same clamp, no rounding.
-pub fn newton_iter_f(tn: []const u32, tp: []const u32, y: f64, y0: f64, a: f64) -> f64 {
-    const p : f64 = y + idf(tn, a, y);
-    const n : f64 = y0 + idf(tp, VDD_F - a, VDD_F - y);
-    const d : f64 = CELL_F + gdsf(tn, a, y) + gdsf(tp, VDD_F - a, VDD_F - y);
+/// newton_iter_q in real arithmetic: the same equation, the same slopes, the same clamp, no rounding.
+pub fn newton_iter_qf(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, y: f64, y0: f64, a: f64, a0: f64) -> f64 {
+    const p : f64 = y + idf(tn, a, y) + idf(qn, a, y) + idf(qp, VDD_F - a0, VDD_F - y0);
+    const n : f64 = y0 + idf(tp, VDD_F - a, VDD_F - y) + idf(qn, a0, y0) + idf(qp, VDD_F - a, VDD_F - y);
+    const d : f64 = CELL_F + gdsf(tn, a, y) + gdsf(tp, VDD_F - a, VDD_F - y) + gdsf(qn, a, y) + gdsf(qp, VDD_F - a, VDD_F - y);
     var z : f64 = y - (p - n) * CELL_F / d;
     if (z < 0.0) { z = 0.0; }
     if (z > VDD_F) { z = VDD_F; }
     return z;
 }
 
-/// The f64 model of the table run: 400 steps of k real Newton iterations from y_init, out in units.
-pub fn run_tabf(tn: []const u32, tp: []const u32, k: u32, y_init: f64, out: []f64) {
+/// newton_iter in real arithmetic, without charge.
+pub fn newton_iter_f(tn: []const u32, tp: []const u32, y: f64, y0: f64, a: f64) -> f64 {
+    return newton_iter_qf(tn, tp, &NOQ, &NOQ, y, y0, a, a);
+}
+
+/// The f64 model of the table run with charge: 400 steps of k real Newton iterations from y_init, the gate before step 0 at
+/// a_rest, out in units.
+pub fn run_tabf_q(tn: []const u32, tp: []const u32, qn: []const u32, qp: []const u32, k: u32, y_init: f64, a_rest: u32, out: []f64) {
     var y : f64 = y_init;
     var n : u32 = 0;
     while (n < RUN_STEPS) {
         const a : f64 = (gate_at(n)) as f64;
+        var a0 : f64 = a_rest as f64;
+        if (n > 0) { a0 = (gate_at(n - 1)) as f64; }
         const y0 : f64 = y;
         var i : u32 = 0;
         while (i < k) {
-            y = newton_iter_f(tn, tp, y, y0, a);
+            y = newton_iter_qf(tn, tp, qn, qp, y, y0, a, a0);
             i = i + 1;
         }
         out[n as usize] = y;
@@ -502,6 +598,11 @@ pub fn run_tabf(tn: []const u32, tp: []const u32, k: u32, y_init: f64, out: []f6
     }
 }
 
+/// The f64 model of the table run without charge: `run_tabf_q` with the zero charge tables.
+pub fn run_tabf(tn: []const u32, tp: []const u32, k: u32, y_init: f64, out: []f64) {
+    run_tabf_q(tn, tp, &NOQ, &NOQ, k, y_init, gate_at(0), out);
+}
+
 /// The analytic reference: backward Euler of the inverter with cpphdl's law, Newton to convergence (ANA_ITERS). f scales
 /// dt (1.0 is 10 ns, 0.5 is 5 ns); the gate steps at step edge. out gets Y in volts after each of nsteps steps.
 pub fn run_ana(f: f64, nsteps: u32, edge: u32, out: []f64) {
@@ -917,3 +1018,112 @@ test a_wrong_reading_of_the_table_is_caught {
     assert fold(fold(0, 1), 2) != fold(fold(0, 2), 1);
     assert_eq(fold(0x8000000000000000, 0), 1);
 }
+
+test the_charge_path_with_zero_charge_tables_is_the_current_only_path_bit_for_bit {
+    // M3c (#10777). The current-only functions are the charge-aware ones with the zero tables: not an equal run but the same
+    // function, so every pinned M3a number above (EXP_S, the first steps) is also a number of the charge path. Spot checks on the
+    // way, over several gates and states, K = 1 .. 3.
+    var t : [625]u32 = [0] ** 625;
+    build_table(&t, VTH_V);
+    var checked : u32 = 0;
+    var y : u32 = 0;
+    while (y < VDD) {
+        var a : u32 = 0;
+        while (a <= VDD) {
+            assert_eq(newton_iter(&t, &t, y, VDD - y, a), newton_iter_q(&t, &t, &NOQ, &NOQ, y, VDD - y, a, VDD - a));
+            assert_eq(inv_step(&t, &t, y, a, 2), inv_step_q(&t, &t, &NOQ, &NOQ, y, a, 0, 2));
+            assert_eq(g_resid(&t, &t, y, 777, a), g_resid_q(&t, &t, &NOQ, &NOQ, y, 777, a, VDD));
+            checked = checked + 1;
+            a = a + 6 * CELL;
+        }
+        y = y + 5 * CELL + 12345;
+    }
+    assert checked > 20;
+    assert_eq(board_step(&t, &t, 4000000, VDD), board_step_q(&t, &t, &NOQ, &NOQ, 4000000, VDD, 0));
+    assert_eq(NOQ[0], 0);
+    assert_eq(NOQ[624], 0);
+}
+
+test a_synthetic_plane_of_charge_gives_the_closed_form_feedthrough_and_load {
+    // No device current (zero tables), a drain charge that is a plane: qn(ix, iy) = 4000000 - 100000 ix + 80000 iy per node, that is a
+    // gate coupling of -K1/2^20 = -0.0954 and a drain capacitance of K2/2^20 = 0.0763 of the node's C. The backward Euler equation
+    // (y - y0)(1 + K2') = K1' (A - A0) is linear, so the gate falling from Vdd to 0 kicks the node down by K1' Vdd / (1 + K2')
+    // = 100000 * 24 / 1.0763 units, and a Newton step on a linear equation lands on it (within the rounding of the reads).
+    var qn : [625]u32 = [0] ** 625;
+    var ix : u32 = 0;
+    while (ix < NODES) {
+        var iy : u32 = 0;
+        while (iy < NODES) {
+            qn[(ix * NODES + iy) as usize] = 4000000 - 100000 * ix + 80000 * iy;
+            iy = iy + 1;
+        }
+        ix = ix + 1;
+    }
+    assert monotone_along(&qn, 1);
+    assert monotone_along(&qn, 0) == false;
+    const y0 : u32 = 20000000;
+    const want : f64 = (y0 as f64) - 2400000.0 / (1.0 + 80000.0 / CELL_F);
+    const got : u32 = inv_step_q(&NOQ, &NOQ, &qn, &NOQ, y0, 0, VDD, 3);
+    assert fabs64((got as f64) - want) <= 3.0;
+    assert got < y0 - 2000000;
+    // The same plane as a PMOS charge (stored negated, so it grows with Vsd as qn grows with Vds): the equation is
+    // (y - y0)(1 + P2') + P1' (A - A0) = 0 with the same numbers, the same kick.
+    const gotp : u32 = inv_step_q(&NOQ, &NOQ, &NOQ, &qn, y0, 0, VDD, 3);
+    assert fabs64((gotp as f64) - want) <= 3.0;
+    // With the gate still, no current and any charge, the node does not move: charge conservation at rest.
+    assert_eq(inv_step_q(&NOQ, &NOQ, &qn, &qn, y0, 5 * CELL, 5 * CELL, 3), y0);
+    assert_eq(inv_step_q(&NOQ, &NOQ, &qn, &qn, y0, VDD, VDD, 3), y0);
+    // The derivative is read off the same table: a plane's slope is its K2 per cell, one cell of Vds.
+    assert_eq(gds_cell(&qn, 7 * CELL + 5, 9 * CELL + 77), 80000);
+    // The f64 twin lands on the closed form too (two real Newton iterations).
+    const g1 : f64 = newton_iter_qf(&NOQ, &NOQ, &qn, &NOQ, y0 as f64, y0 as f64, 0.0, CELL_F * 24.0);
+    const g2 : f64 = newton_iter_qf(&NOQ, &NOQ, &qn, &NOQ, g1, y0 as f64, 0.0, CELL_F * 24.0);
+    assert fabs64(g2 - want) <= 1.0;
+}
+
+test the_mirrored_circuit_with_charge_gives_the_mirrored_step_with_the_four_tables_swapped {
+    // The mirror of the charge-aware equation: swap NMOS and PMOS (tn <-> tp and qn <-> qp, qp stored negated), A -> Vdd - A,
+    // A0 -> Vdd - A0, y -> Vdd - y: the step is the mirror of the step, bit for bit. Two unlike current tables and two unlike
+    // charge tables; the states include a gate edge.
+    var tn : [625]u32 = [0] ** 625;
+    var tp : [625]u32 = [0] ** 625;
+    build_table(&tn, VTH_V);
+    build_table(&tp, VTH_WRONG_V);
+    var qn : [625]u32 = [0] ** 625;
+    var qp : [625]u32 = [0] ** 625;
+    var ix : u32 = 0;
+    while (ix < NODES) {
+        var iy : u32 = 0;
+        while (iy < NODES) {
+            qn[(ix * NODES + iy) as usize] = 4000000 - 90000 * ix + 70000 * iy + 1000 * ix * iy / 5;
+            qp[(ix * NODES + iy) as usize] = 3000000 + 30000 * ix + 120000 * iy + 700 * iy * iy;
+            iy = iy + 1;
+        }
+        ix = ix + 1;
+    }
+    assert monotone_along(&qn, 1);
+    assert monotone_along(&qp, 1);
+    var checked : u32 = 0;
+    var y0 : u32 = 3000000;
+    while (y0 < VDD - 1000000) {
+        var a0 : u32 = 0;
+        while (a0 <= VDD) {
+            var a : u32 = 0;
+            while (a <= VDD) {
+                const m : u32 = inv_step_q(&tp, &tn, &qp, &qn, VDD - y0, VDD - a, VDD - a0, 3);
+                assert_eq(inv_step_q(&tn, &tp, &qn, &qp, y0, a, a0, 3), VDD - m);
+                assert_eq(g_resid_q(&tn, &tp, &qn, &qp, y0 + 4321, y0, a, a0), g_resid_q(&tp, &tn, &qp, &qn, VDD - y0 - 4321, VDD - y0, VDD - a, VDD - a0));
+                checked = checked + 1;
+                a = a + 4 * CELL;
+            }
+            a0 = a0 + 12 * CELL;
+        }
+        y0 = y0 + 7 * CELL + 999;
+    }
+    assert checked > 50;
+    // Not an identity of the code but of the equation: swapping only the current tables, or only the charge tables, is another step.
+    assert inv_step_q(&tn, &tp, &qn, &qp, 9000000, 0, VDD, 3) != inv_step_q(&tp, &tn, &qn, &qp, 9000000, 0, VDD, 3);
+    assert inv_step_q(&tn, &tp, &qn, &qp, 9000000, 0, VDD, 3) != inv_step_q(&tn, &tp, &qp, &qn, 9000000, 0, VDD, 3);
+    // A step that sees no previous gate (A0 = A) differs from one coming from Vdd: the gate edge is in the charge.
+    assert inv_step_q(&tn, &tp, &qn, &qp, 9000000, 0, VDD, 3) != inv_step_q(&tn, &tp, &qn, &qp, 9000000, 0, 0, 3);
+}