==============================================================================
Theorem 2: canonical relations level by level (max partial-isometry defect)
==============================================================================
A = 1_4 (Cuntz O_4)                           primitive=True  noncomm=0.943  defects by level [0.0, 0.513, 0.513]  first failure: level 2  (bound 2 <= N=4)
    pair (p_11, p_22) via alpha=(1, 2): ||p_ab p_cd - sum_beta T|| = 0.0e+00; ||sum TT^*-1|| = 3.8e-15; ||sum T^*T-1|| = 3.7e-15; max|T| over non-admissible beta = 0.0e+00
A = J_4 - I_4                                 primitive=True  noncomm=0.943  defects by level [0.0, 0.513, 0.512]  first failure: level 2  (bound 2 <= N=4)
    pair (p_11, p_22) via alpha=(1, 2): ||p_ab p_cd - sum_beta T|| = 0.0e+00; ||sum TT^*-1|| = 3.8e-15; ||sum T^*T-1|| = 3.7e-15; max|T| over non-admissible beta = 2.5e-17
hub: 4 independent twins + looped hub (N=5)   primitive=True  noncomm=0.943  defects by level [0.0, 0.0, 0.513, 0.513]  first failure: level 3  (bound 3 <= N=5)
    pair (p_11, p_22) via alpha=(1, 5, 2): ||p_ab p_cd - sum_beta T|| = 0.0e+00; ||sum TT^*-1|| = 3.8e-15; ||sum T^*T-1|| = 3.7e-15; max|T| over non-admissible beta = 0.0e+00
twins at distance 3 (N=6)                     primitive=True  noncomm=0.943  defects by level [0.0, 0.0, 0.0, 0.513, 0.513]  first failure: level 4  (bound 4 <= N=6)
    pair (p_11, p_22) via alpha=(1, 5, 6, 2): ||p_ab p_cd - sum_beta T|| = 0.0e+00; ||sum TT^*-1|| = 3.8e-15; ||sum T^*T-1|| = 3.7e-15; max|T| over non-admissible beta = 0.0e+00
Example 3: A_D (QAut = dual D_inf)            primitive=True  noncomm=0.707  defects by level [0.0, 0.354, 0.375, 0.375]  first failure: level 2  (bound 3 <= N=4)
    pair (p_11, p_33) via alpha=(1, 3): ||p_ab p_cd - sum_beta T|| = 0.0e+00; ||sum TT^*-1|| = 0.0e+00; ||sum T^*T-1|| = 0.0e+00; max|T| over non-admissible beta = 0.0e+00

==============================================================================
Example 3: A_D = [[1,0,1,1],[0,1,1,1],[1,1,0,0],[1,1,0,0]]
==============================================================================
A_D^2 = [[3, 2, 1, 1], [2, 3, 1, 1], [1, 1, 2, 2], [1, 1, 2, 2]] -> primitive: True
B = J - A_D = [[0, 1, 0, 0], [1, 0, 0, 0], [0, 0, 1, 1], [0, 0, 1, 1]] (adjacency of K_2 (+) K_2-with-loops); out-degrees of B: [1, 1, 2, 2]
Aut(A_D) = [(1, 2, 3, 4), (1, 2, 4, 3), (2, 1, 3, 4), (2, 1, 4, 3)]  (order 4 )
w entries partial isometries: True ; range/source projections magic: True ; A Pr = So A: True
max partial-isometry defect of w_{alpha,beta} over admissible words of length <= 6: 0.0
||Delta(2e-1) - (2e-1)(x)(2e-1)|| (evaluated in the block model) = 0.0
non-commutativity [s,t] in the model: 2.8284271247461903
DONE s4
