J = E gamma_q^+ = 3.1664009357e-04  (grid quadrature; rigorous value 3.1663998209e-04)
check E gamma_q = -0.24999999, E gamma_3 = -0.25000001 (both -1/4); max gamma_3 = -0.052736 (< 2*0.49788-1 = -0.00424)
k = 1: best two-role value - 5/8 = 1.187458e-04 = 0.3750 J  (bound (1-4^-k) = 0.7500); f = 10; largest contribution of nu with undetermined role pattern, over all f: 5.73e-09
    nu = 0: weights of t = 0..3: 0.2000, 0.0000, 0.6000, 0.2000
    nu = 1: weights of t = 0..3: 0.0000, 1.0000, 0.0000, 0.0000
k = 2: best two-role value - 5/8 = 1.187458e-04 = 0.3750 J  (bound (1-4^-k) = 0.9375); f = 1100; largest contribution of nu with undetermined role pattern, over all f: 5.19e-06
    nu = 00: weights of t = 0..3: 0.2000, 0.0000, 0.6000, 0.2000
    nu = 01: weights of t = 0..3: 0.2000, 0.0000, 0.6000, 0.2000
    nu = 10: weights of t = 0..3: 0.0000, 1.0000, 0.0000, 0.0000
    nu = 11: weights of t = 0..3: 0.0000, 1.0000, 0.0000, 0.0000
k = 3: best two-role value - 5/8 = 1.187458e-04 = 0.3750 J  (bound (1-4^-k) = 0.9844); f = 11110000; largest contribution of nu with undetermined role pattern, over all f: 3.70e-05
k = 4: best two-role value - 5/8 = 1.187458e-04 = 0.3750 J  (bound (1-4^-k) = 0.9961); f = 1111111100000000; largest contribution of nu with undetermined role pattern, over all f: 7.28e-05
