Changeset 2730

Timestamp:

Apr 19, 2006, 6:34:07 PM (18 years ago)

Author:

ole

Message:

More explanation of flux calculations in test

File:

• inundation/pyvolution/test_shallow_water.py (modified) (1 diff)

inundation/pyvolution/test_shallow_water.py

@@ -437 +437 @@
 
         #Flux across right edge of volume 1
-        normal = domain.get_normal(1,0)
+        normal = domain.get_normal(1,0) #Get normal 0 of triangle 1
+        assert allclose(normal, [1, 0])
+        
         ql = domain.get_conserved_quantities(vol_id=1, edge=0)
+        assert allclose(ql, [val1, 0, 0])
+        
         qr = domain.get_conserved_quantities(vol_id=2, edge=2)
+        assert allclose(qr, [val2, 0, 0])
+        
         flux0, max_speed = flux_function(normal, ql, qr, zl, zr)
+
+        #Flux across edge in the east direction (as per normal vector)
         assert allclose(flux0, [-15.3598804, 253.71111111, 0.])
         assert allclose(max_speed, 9.21592824046)
 
+
+        #Flux across edge in the west direction (opposite sign for xmomentum)
+        normal_opposite = domain.get_normal(2,2) #Get normal 2 of triangle 2
+        assert allclose(normal_opposite, [-1, 0])
+        flux_opposite, max_speed = flux_function([-1, 0], ql, qr, zl, zr)
+        assert allclose(flux_opposite, [-15.3598804, -253.71111111, 0.])
+        
 
         #Flux across upper edge of volume 1