def test2by2(self):
a = 2e0
S11 = 1e0
S12 = (1 + a) * exp(-a)
S = array([[S11, S12], [S12, S11]])
V, Vdag, Vinv = part_a(S)
print reduce(dot, (Vdag, S, V))
print Vinv
print inv(V)
print Vdag
print(dot(Vinv, V))
print(dot(V, Vinv))
def test2by2(self):
a = 2e0
S11 = 1e0
S12 = (1 + a) * exp(-a)
S = array([[S11, S12], [S12, S11]])
V, Vdag, Vinv = part_a(S)
print reduce(dot, (Vdag, S, V))
print Vinv
print inv(V)
print Vdag
print(dot(Vinv, V))
print(dot(V, Vinv))
def test5by5(self):
print('Test 5 by 5')
S = array([[1.06666667, 0., 0.15238095, 0., .05079365],
[0., 0.15238095, 0., 0.05079365, 0.],
[0.15238095, 0., 0.05079365, 0., 0.02308802],
[0., 0.05079365, 0., 0.02308802, 0.],
[0.05079365, 0., 0.02308802, 0., 0.01243201]])
V, Vdag, Vinv = part_a(S)
print reduce(dot, (Vdag, S, V))
print Vinv
print inv(V)
print Vdag
print(dot(Vinv, V))
print(dot(V, Vinv))
def test5by5(self):
print('Test 5 by 5')
S = array([[ 1.06666667, 0., 0.15238095, 0., .05079365],
[ 0., 0.15238095, 0., 0.05079365, 0.],
[ 0.15238095, 0., 0.05079365, 0., 0.02308802],
[ 0., 0.05079365, 0., 0.02308802, 0.],
[ 0.05079365, 0., 0.02308802, 0., 0.01243201]])
V, Vdag, Vinv = part_a(S)
print reduce(dot, (Vdag, S, V))
print Vinv
print inv(V)
print Vdag
print(dot(Vinv, V))
print(dot(V, Vinv))
