(1) Q1 - Q2 + N_3/512 == 0 identically: True
(2) 16*(conic - 3/16) == ((ri-rj)^2-1)((ri+rj)^2-1): True
(3) N=400000 uniform points of [0,1/4]^3; exact locus fraction 0.3674
    Thm 1.4 with '>= 3/16' (proof version): disagreements with locus = 0
    Thm 1.4 with '<= 3/16' (printed):       disagreements with locus = 183284  (locus-only 4416, printed-only 178868)
    region1 = {Q1>=Q2} contained in locus: violations = 0
    region 2 (proof version) vs corner region {r1+r2+r3 <= 1}: disagreements = 0 (corner fraction 0.0110)
(4) GHZ d=(1/4,1/4,1/4): Q1-Q2 = -1/512 ; conic value = 1/4 > 3/16: True
    numpy Gram determinants of GHZ: [0.25, 0.25, 0.25]
ALL CHECKS DONE
