-- Mathlib/.lake/packages/lean4/src/lean/Init/Data/BitVec/Bitblast.lean
theorem lawful_divSubtractShift (qr : DivModState w) (h : qr.Poised args) :
    DivModState.Lawful args (divSubtractShift args qr)

Proof state:
case neg.hrWidth
w : Nat
qr : DivModState w
n d : BitVec w
h : DivModState.Poised { n := n, d := d } qr
hrwn : qr.wr + qr.wn = w
hd : 0 < { n := n, d := d }.d
hrd : qr.r.toNat < { n := n, d := d }.d.toNat
hr : qr.r.toNat < 2 ^ qr.wr
hn : qr.q.toNat < 2 ^ qr.wr
hrnd : { n := n, d := d }.n.toNat >>> qr.wn = { n := n, d := d }.d.toNat * qr.q.toNat + qr.r.toNat
hwn_lt : 0 < qr.wn
this : d.toNat * (qr.q.toNat * 2) = d.toNat * qr.q.toNat * 2
rltd : ¬qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) < d
⊢ { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.r.toNat <
    2 ^
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.wr

case neg.hqWidth
w : Nat
qr : DivModState w
n d : BitVec w
h : DivModState.Poised { n := n, d := d } qr
hrwn : qr.wr + qr.wn = w
hd : 0 < { n := n, d := d }.d
hrd : qr.r.toNat < { n := n, d := d }.d.toNat
hr : qr.r.toNat < 2 ^ qr.wr
hn : qr.q.toNat < 2 ^ qr.wr
hrnd : { n := n, d := d }.n.toNat >>> qr.wn = { n := n, d := d }.d.toNat * qr.q.toNat + qr.r.toNat
hwn_lt : 0 < qr.wn
this : d.toNat * (qr.q.toNat * 2) = d.toNat * qr.q.toNat * 2
rltd : ¬qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) < d
⊢ { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.q.toNat <
    2 ^
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.wr

case neg.hdiv
w : Nat
qr : DivModState w
n d : BitVec w
h : DivModState.Poised { n := n, d := d } qr
hrwn : qr.wr + qr.wn = w
hd : 0 < { n := n, d := d }.d
hrd : qr.r.toNat < { n := n, d := d }.d.toNat
hr : qr.r.toNat < 2 ^ qr.wr
hn : qr.q.toNat < 2 ^ qr.wr
hrnd : { n := n, d := d }.n.toNat >>> qr.wn = { n := n, d := d }.d.toNat * qr.q.toNat + qr.r.toNat
hwn_lt : 0 < qr.wn
this : d.toNat * (qr.q.toNat * 2) = d.toNat * qr.q.toNat * 2
rltd : ¬qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) < d
⊢ { n := n, d := d }.n.toNat >>>
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.wn =
    { n := n, d := d }.d.toNat *
        { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
              r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.q.toNat +
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
            r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.r.toNat

Tactic:
case neg.hrWidth =>
  simp only
  have hdr' : d ≤ (qr.r.shiftConcat (n.getLsbD (qr.wn - 1))) :=
    BitVec.not_lt_iff_le.mp rltd
  have hr' : ((qr.r.shiftConcat (n.getLsbD (qr.wn - 1)))).toNat < 2 ^ (qr.wr + 1) := by
    apply toNat_shiftConcat_lt_of_lt <;> bv_omega
  rw [BitVec.toNat_sub_of_le hdr']
  omega

Result:
case neg.hqWidth
w : Nat
qr : DivModState w
n d : BitVec w
h : DivModState.Poised { n := n, d := d } qr
hrwn : qr.wr + qr.wn = w
hd : 0 < { n := n, d := d }.d
hrd : qr.r.toNat < { n := n, d := d }.d.toNat
hr : qr.r.toNat < 2 ^ qr.wr
hn : qr.q.toNat < 2 ^ qr.wr
hrnd : { n := n, d := d }.n.toNat >>> qr.wn = { n := n, d := d }.d.toNat * qr.q.toNat + qr.r.toNat
hwn_lt : 0 < qr.wn
this : d.toNat * (qr.q.toNat * 2) = d.toNat * qr.q.toNat * 2
rltd : ¬qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) < d
⊢ { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.q.toNat <
    2 ^
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.wr

case neg.hdiv
w : Nat
qr : DivModState w
n d : BitVec w
h : DivModState.Poised { n := n, d := d } qr
hrwn : qr.wr + qr.wn = w
hd : 0 < { n := n, d := d }.d
hrd : qr.r.toNat < { n := n, d := d }.d.toNat
hr : qr.r.toNat < 2 ^ qr.wr
hn : qr.q.toNat < 2 ^ qr.wr
hrnd : { n := n, d := d }.n.toNat >>> qr.wn = { n := n, d := d }.d.toNat * qr.q.toNat + qr.r.toNat
hwn_lt : 0 < qr.wn
this : d.toNat * (qr.q.toNat * 2) = d.toNat * qr.q.toNat * 2
rltd : ¬qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) < d
⊢ { n := n, d := d }.n.toNat >>>
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
          r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.wn =
    { n := n, d := d }.d.toNat *
        { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
              r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.q.toNat +
      { wn := qr.wn - 1, wr := qr.wr + 1, q := qr.q.shiftConcat true,
            r := qr.r.shiftConcat (n.getLsbD (qr.wn - 1)) - d }.r.toNat