def test_init(self):
a = [0.0, 0.0]
b = [0.0, 2.0]
c = [2.0,0.0]
d = [0.0, 4.0]
e = [2.0, 2.0]
f = [4.0,0.0]
points = [a, b, c, d, e, f]
#bac, bce, ecf, dbe, daf, dae
vertices = [ [1,0,2], [1,2,4], [4,2,5], [3,1,4]]
domain = Advection_Domain(points, vertices)
domain.check_integrity()
assert domain.quantities.has_key('stage')
assert domain.get_conserved_quantities(0, edge=1) == 0.
