Changeset 2776


Ignore:
Timestamp:
Apr 27, 2006, 1:27:41 PM (18 years ago)
Author:
jack
Message:

Update to the parallel shallow water test case.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/parallel/test_parallel_sw.py

    r2739 r2776  
    11#!/usr/bin/env python
     2# Test a run of the sequential shallow water domain against
     3# a run of the parallel shallow water domain.
     4# WARNING: This assumes that the command to run jobs is mpirun.
     5# Tested with MPICH.
    26
    37mesh_filename = "merimbula_10785.tsh"
     
    178182        assert result == 0
    179183
     184# Because we are doing assertions outside of the TestCase class
     185# the PyUnit defined assert_ function can't be used.
     186def assert_(condition, msg="Assertion Failed"):
     187    if condition == False:
     188        raise AssertionError, msg
     189
    180190if __name__=="__main__":
    181191    if pypar.size() == 1:
     
    188198        l1norm_par, l2norm_par, linfnorm_par = parallel_test()
    189199        if pypar.rank() == 0:
    190             assert len(l1norm_seq) == len(l1norm_par)
    191             assert len(l2norm_seq) == len(l2norm_par)
    192             assert len(linfnorm_seq) == len(linfnorm_par)
    193             assert len(l1norm_seq) == len(l2norm_seq)
    194             assert len(l2norm_seq) == len(linfnorm_seq)
     200            assert_(len(l1norm_seq) == len(l1norm_par))
     201            assert_(len(l2norm_seq) == len(l2norm_par))
     202            assert_(len(linfnorm_seq) == len(linfnorm_par))
     203            assert_(len(l1norm_seq) == len(l2norm_seq))
     204            assert_(len(l2norm_seq) == len(linfnorm_seq))
    195205            # Anything smaller than tol we consider to be 0.
    196206            tol = pow(10, -1 * double_precision())
    197207            for x in range(len(l1norm_seq)):
    198208                for y in range(3):
    199                     assert (l1norm_seq[x][y] - l1norm_par[x][y]) < tol
    200                     assert (l2norm_seq[x][y] - l2norm_par[x][y]) < tol
    201                     assert (linfnorm_seq[x][y] - linfnorm_par[x][y]) < tol
     209                    assert_((l1norm_seq[x][y] - l1norm_par[x][y]) < tol)
     210                    assert_((l2norm_seq[x][y] - l2norm_par[x][y]) < tol)
     211                    assert_((linfnorm_seq[x][y] - linfnorm_par[x][y]) < tol)
     212                    if x > 0:
     213                        assert_((l1norm_seq[x][y] - l1norm_seq[x-1][y]) < tol)
     214                        assert_((l2norm_seq[x][y] - l2norm_seq[x-1][y]) < tol)
     215                        assert_((linfnorm_seq[x][y] - linfnorm_seq[x-1][y]) < tol)
     216                        assert_((l1norm_par[x][y] - l1norm_par[x-1][y]) < tol)
     217                        assert_((l2norm_par[x][y] - l2norm_par[x-1][y]) < tol)
     218                        assert_((linfnorm_par[x][y] - linfnorm_par[x-1][y]) < tol)
    202219        pypar.finalize()
     220
Note: See TracChangeset for help on using the changeset viewer.