Changeset 2776 for inundation/parallel/test_parallel_sw.py

Timestamp:

Apr 27, 2006, 1:27:41 PM (19 years ago)

Author:

jack

Message:

Update to the parallel shallow water test case.

File:

• inundation/parallel/test_parallel_sw.py (modified) (3 diffs)

inundation/parallel/test_parallel_sw.py

@@ -1 +1 @@
 #!/usr/bin/env python
+# Test a run of the sequential shallow water domain against
+# a run of the parallel shallow water domain.
+# WARNING: This assumes that the command to run jobs is mpirun.
+# Tested with MPICH.
 
 mesh_filename = "merimbula_10785.tsh"
@@ -178 +182 @@
         assert result == 0
 
+# Because we are doing assertions outside of the TestCase class
+# the PyUnit defined assert_ function can't be used.
+def assert_(condition, msg="Assertion Failed"):
+    if condition == False:
+        raise AssertionError, msg
+
 if __name__=="__main__":
     if pypar.size() == 1: 
@@ -188 +198 @@
         l1norm_par, l2norm_par, linfnorm_par = parallel_test()
         if pypar.rank() == 0:
-            assert len(l1norm_seq) == len(l1norm_par)
-            assert len(l2norm_seq) == len(l2norm_par)
-            assert len(linfnorm_seq) == len(linfnorm_par)
-            assert len(l1norm_seq) == len(l2norm_seq)
-            assert len(l2norm_seq) == len(linfnorm_seq)
+            assert_(len(l1norm_seq) == len(l1norm_par))
+            assert_(len(l2norm_seq) == len(l2norm_par))
+            assert_(len(linfnorm_seq) == len(linfnorm_par))
+            assert_(len(l1norm_seq) == len(l2norm_seq))
+            assert_(len(l2norm_seq) == len(linfnorm_seq))
             # Anything smaller than tol we consider to be 0.
             tol = pow(10, -1 * double_precision())
             for x in range(len(l1norm_seq)):
                 for y in range(3):
-                    assert (l1norm_seq[x][y] - l1norm_par[x][y]) < tol
-                    assert (l2norm_seq[x][y] - l2norm_par[x][y]) < tol
-                    assert (linfnorm_seq[x][y] - linfnorm_par[x][y]) < tol
+                    assert_((l1norm_seq[x][y] - l1norm_par[x][y]) < tol)
+                    assert_((l2norm_seq[x][y] - l2norm_par[x][y]) < tol)
+                    assert_((linfnorm_seq[x][y] - linfnorm_par[x][y]) < tol)
+                    if x > 0:
+                        assert_((l1norm_seq[x][y] - l1norm_seq[x-1][y]) < tol)
+                        assert_((l2norm_seq[x][y] - l2norm_seq[x-1][y]) < tol)
+                        assert_((linfnorm_seq[x][y] - linfnorm_seq[x-1][y]) < tol)
+                        assert_((l1norm_par[x][y] - l1norm_par[x-1][y]) < tol)
+                        assert_((l2norm_par[x][y] - l2norm_par[x-1][y]) < tol)
+                        assert_((linfnorm_par[x][y] - linfnorm_par[x-1][y]) < tol)
         pypar.finalize()
+