def test_gauss_seidel_2(self):
print 'test gauss_seidel 2'
g1,g2,d1,d2 = self.braess_heterogeneous(.25, .25)
fs = gauss_seidel([g1,g2], [d1,d2], solver_2)
a = np.array([[.125,.25],[.125,.0],[.0, .25],[.125, .0],[.125, .25]])
self.check(fs, a, 1e-3)
g1,g2,d1,d2 = self.braess_heterogeneous(1., 1.)
a = np.array([[.5,.5],[.5,.5],[.0, .0],[.5, .5],[.5, .5]])
fs = gauss_seidel([g1,g2], [d1,d2], solver_2)
self.check(fs, a, 1e-3)
g1,g2,d1,d2 = self.braess_heterogeneous(.75, .75)
fs = gauss_seidel([g1,g2], [d1,d2], solver_2, max_iter=200)
a = np.array([[.375, .625],[.375, .125],[.0, .5],[.375, .125],[.375, .625]])
self.check(fs, a, 1e-3)
def test_gauss_seidel_2(self):
print 'test gauss_seidel 2'
g1, g2, d1, d2 = self.braess_heterogeneous(.25, .25)
fs = gauss_seidel([g1, g2], [d1, d2], solver_2)
a = np.array([[.125, .25], [.125, .0],
[.0, .25], [.125, .0], [.125, .25]])
self.check(fs, a, 1e-3)
g1, g2, d1, d2 = self.braess_heterogeneous(1., 1.)
a = np.array([[.5, .5], [.5, .5], [.0, .0], [.5, .5], [.5, .5]])
fs = gauss_seidel([g1, g2], [d1, d2], solver_2)
self.check(fs, a, 1e-3)
g1, g2, d1, d2 = self.braess_heterogeneous(.75, .75)
fs = gauss_seidel([g1, g2], [d1, d2], solver_2, max_iter=200)
a = np.array([[.375, .625], [.375, .125],
[.0, .5], [.375, .125], [.375, .625]])
self.check(fs, a, 1e-3)
def braess_parametric_study():
'''
parametric study of heterogeneous game on the Braess network
'''
g1, g2, d1, d2 = braess_heterogeneous(.0, 1.5)
fs = solver_2(g1, d1, display=1)
print '.0, 1.5'
np.savetxt('data/braess/test_1.csv', fs, delimiter=',')
g1, g2, d1, d2 = braess_heterogeneous(.5, 1.)
fs = gauss_seidel([g1, g2], [d1, d2], solver_2, display=1)
print '.5, 1.'
np.savetxt('data/braess/test_2.csv', fs, delimiter=',')
g1, g2, d1, d2 = braess_heterogeneous(.75, .75)
fs = gauss_seidel([g1, g2], [d1, d2], solver_2, display=1)
print '.75, .75'
np.savetxt('data/braess/test_3.csv', fs, delimiter=',')
g1, g2, d1, d2 = braess_heterogeneous(1., .5)
fs = gauss_seidel([g1, g2], [d1, d2], solver_2, display=1)
print '1., .5'
np.savetxt('data/braess/test_4.csv', fs, delimiter=',')
g1, g2, d1, d2 = braess_heterogeneous(1.5, .0)
fs = solver_2(g2, d2, display=1)
print '1.5, .0'
np.savetxt('data/braess/test_5.csv', fs, delimiter=',')
def braess_parametric_study():
'''
parametric study of heterogeneous game on the Braess network
'''
g1,g2,d1,d2 = braess_heterogeneous(.0, 1.5)
fs = solver_2(g1, d1, display=1)
print '.0, 1.5'
np.savetxt('data/braess/test_1.csv', fs, delimiter=',')
g1,g2,d1,d2 = braess_heterogeneous(.5, 1.)
fs = gauss_seidel([g1,g2], [d1,d2], solver_2, display=1)
print '.5, 1.'
np.savetxt('data/braess/test_2.csv', fs, delimiter=',')
g1,g2,d1,d2 = braess_heterogeneous(.75, .75)
fs = gauss_seidel([g1,g2], [d1,d2], solver_2, display=1)
print '.75, .75'
np.savetxt('data/braess/test_3.csv', fs, delimiter=',')
g1,g2,d1,d2 = braess_heterogeneous(1., .5)
fs = gauss_seidel([g1,g2], [d1,d2], solver_2, display=1)
print '1., .5'
np.savetxt('data/braess/test_4.csv', fs, delimiter=',')
g1,g2,d1,d2 = braess_heterogeneous(1.5, .0)
fs = solver_2(g2, d2, display=1)
print '1.5, .0'
np.savetxt('data/braess/test_5.csv', fs, delimiter=',')
