  ok   delta = 0: J_4 = <p, q> coprime, D_4 = J_4-perp, D_5 = (S_1 J_4)-perp, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 0: (f, P, Q + 7 l^6) local length 32 <= 32
  ok   delta = 0: J_4 = <p, q> coprime, D_4 = J_4-perp, D_5 = (S_1 J_4)-perp, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 0: (f, P, Q + 7 l^6) local length 32 <= 32
  ok   delta = 1: J_4 = y <y^3 + z^3, y^2 z + 2 z^3>, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 1: (f, P, Q + 7 l^6) local length 33 <= 36
  ok   delta = 1: J_4 = y <y^3 + z^3, y^2 z + 2 z^3>, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 1: (f, P, Q + 7 l^6) local length 33 <= 36
  ok   D_5 = <Y^5, Z^5>, D_4 = <Y^4, Z^4, X Y^2 Z + Y^2 Z^2>, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 34 <= 40
  ok   D_5 = <Y^5, Z^5>, D_4 = <Y^4, Z^4, X Y^2 Z + Y^2 Z^2>, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 34 <= 40
  ok   D_5 = <Y^5 + 5 X Y^4, Z^5> (non-binary), seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 35 <= 40
  ok   D_5 = <Y^5 + 5 X Y^4, Z^5> (non-binary), seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 36 <= 40
  ok   D_5 = <Y^5, Y^4 Z>, D_4 = <Y^4, Y^3 Z, Z^4 + X Y^3>, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 34 <= 40
  ok   D_5 = <Y^5, Y^4 Z>, D_4 = <Y^4, Y^3 Z, Z^4 + X Y^3>, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 3
  ok      delta = 2: (f, P, Q + 7 l^6) local length 34 <= 40
  ok   D_5 = <Y^5, Y^4 Z>, D_4 = <Y^4, Y^3 Z, Y^2 Z^2> (binary), seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = z^3 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  ok   D_5 = <Y^5, Y^4 Z>, D_4 = <Y^4, Y^3 Z, Y^2 Z^2> (binary), seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = z^3 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  ok   Gamma = V(x, yz(y+z)): D_4 = nu_4, D_5 = <Y^5 + Z^5, (Y-Z)^5>, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = y^2*z + y*z^2 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  ok   Gamma = V(x, yz(y+z)): D_4 = nu_4, D_5 = <Y^5 + Z^5, (Y-Z)^5>, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = y^2*z + y*z^2 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  ok   Gamma = V(x, z^2 y): D_4 = nu_4 = <Y^4, Y^3 Z, Z^4>, D_5 = <Y^5 + Z^5, Y^4 Z>, seed 0: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = y*z^2 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  ok   Gamma = V(x, z^2 y): D_4 = nu_4 = <Y^4, Y^3 Z, Z^4>, D_5 = <Y^5 + Z^5, Y^4 Z>, seed 1: closed, H = [1, 3, 5, 7, 3, 2]
  ok      dim J_4 = 2 + dim(X-parts of D_4) = 2
  ok      delta = 3: h = y*z^2 kills D_5, lies in (I : m^2), and h in pi(cubic parts of (I : m^2))
  accepted samples: 16 of 16; deltas seen: [0, 1, 2, 3]
  ok   all samples accepted, all deltas 0..3 represented
(A) extreme models: f = x^2, P = p, Q = q + nu l^6 (common roots of p, q of total multiplicity delta in p)
  ok   delta = 1: (x^2, y^4 + 2*y^3*z - y^2*z^2 - 2*y*z^3, y^4 + y^2*z^2 + 3*y*z^3 + 7 z^6): 36 (stable at N = 10) <= 36
  ok   delta = 2: (x^2, y^4 + y^3*z - 2*y^2*z^2, y^4 + y^3*z + 3*y^2*z^2 + 7 z^6): 40 (stable at N = 10) <= 40
  ok   delta = 2: (x^2, y^4 + 2*y^3*z - y^2*z^2 - 2*y*z^3, y^4 + 2*y^2*z^2 - 3*y*z^3 + 7 z^6): 40 (stable at N = 10) <= 40
time 17.0 s
ALL CHECKS PASSED
