A[3:0] ──┐  (dividend, 0-15)
         ├──► Divider ──┬──► Q[3:0]  (quotient)
B[3:0] ──┘  (divisor, 1-15)  └──► R[3:0]  (remainder)

A = Q × B + R, where R < B

P = 0  (partial remainder)
for i = 3 downto 0:
    P = (P << 1) | A[i]     // shift in next dividend bit
    diff = P - B
    if diff >= 0:           // P >= B
        Q[i] = 1
        P = diff            // keep subtraction result
    else:
        Q[i] = 0            // restore P (don't subtract)
R = P  (final remainder)

13 ÷ 3:
  Stage 3: P=0, shift in 1 → P=1, 1<3, Q[3]=0
  Stage 2: P=1, shift in 1 → P=3, 3≥3, Q[2]=1, P=0
  Stage 1: P=0, shift in 0 → P=0, 0<3, Q[1]=0
  Stage 0: P=0, shift in 1 → P=1, 1<3, Q[0]=0
  Result: Q=0100=4, R=0001=1
  Verify: 4×3+1=13 ✓

threshold-divider/
├── model.safetensors
├── create_safetensors.py
├── config.json
└── README.md