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=',')