def test_flux_1_triangle0(self):
a = [0.0, 0.5]
b = [0.0, 0.0]
c = [0.5, 0.5]
points = [a, b, c]
vertices = [ [0,1,2] ]
domain = Advection_Domain(points, vertices)
domain.check_integrity()
#Populate boundary array with dirichlet conditions.
domain.neighbours = num.array([[-1,-2,-3]])
domain.quantities['stage'].boundary_values[:] = 1.0
domain.order = 1
domain.distribute_to_vertices_and_edges() #Use first order default
domain.check_integrity()
domain.compute_fluxes()
U = -domain.quantities['stage'].explicit_update
R = -0.5/domain.areas[0]
assert U==R, '%s %s' %(U, R)
