def test_gauss_seidel_2(self):
print 'test gauss_seidel 2'
g1,g2,d1,d2 = 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 = 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 = 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_fw_heterogeneous_2(self):
print 'test fw_heterogeneous 2'
g1,g2,d1,d2 = braess_heterogeneous(.25, .25)
fs = fw_heterogeneous_2([g1,g2], [d1,d2], q=10, past=10)
a = np.array([[.125,.25],[.125,.0],[.0, .25],[.125, .0],[.125, .25]])
self.check(fs, a, 1e-2)
g1,g2,d1,d2 = braess_heterogeneous(1., 1.)
a = np.array([[.5,.5],[.5,.5],[.0, .0],[.5, .5],[.5, .5]])
fs = fw_heterogeneous_2([g1,g2], [d1,d2], q=10, past=10)
self.check(fs, a, 1e-2)
g1,g2,d1,d2 = braess_heterogeneous(.75, .75)
fs = fw_heterogeneous_2([g1,g2], [d1,d2], q=200, past=10, max_iter=200)
a = np.array([[.375, .625],[.375, .125],[.0, .5],[.375, .125],[.375, .625]])
self.check(np.sum(fs,1), np.sum(a,1), 1e-2)
def test_gauss_seidel_3(self):
print 'test gauss_seidel 3'
g1,g2,d1,d2 = braess_heterogeneous(.25, .25)
fs = gauss_seidel([g1,g2], [d1,d2], solver_3)
a = np.array([[.125,.25],[.125,.0],[.0, .25],[.125, .0],[.125, .25]])
#print fs
#print a
self.check(fs, a, 1e-3)
g1,g2,d1,d2 = 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_3)
self.check(fs, a, 1e-3)
g1,g2,d1,d2 = braess_heterogeneous(.75, .75)
fs = gauss_seidel([g1,g2], [d1,d2], solver_3, q=50)
a = np.array([[.375, .625],[.375, .125],[.0, .5],[.375, .125],[.375, .625]])
self.check(fs, a, 1e-3)
