Numerical Performance Results From The Shallow Water Equation Test Suite

Numerical Performance Results from the Shallow
Water Equation Test Suite

John B. Drake, Ed.
Updated: 3 May 1993
Introduction
The DOE Computer Hardware, Advanced Mathematics and Model Physics

(CHAMMP) program seeks to provide climate researchers with an advanced
modeling capability for the study of global change issues and is interested in the
development of new methods for the study of climate dynamics. The shallow
water equations have been used as a kernel for both oceanic and atmospheric
general circulation models and are useful in evaluating numerical methods for
weather forecasting and climate modeling. To promote development of new
methods, a set of test cases has been proposed [?] and example software and
reference solutions provided [?].
This report summarizes the performance of methods that have been applied
to the test cases. Promising schemes should be subjected to other tests appropriate to their intended application. It is hoped that the bibliography provided
herein will offer pointers to the appropriate literature for more comprehensive
studies of the strengths and weaknesses of individual methods.
Comparison of Algorithms and Computing Platforms
Table 1 gives a comparison of methods by accuracy and computational performance. In Table 1, execution time is given in seconds and represents the best
measurement of how long it takes to perform a 5 day integration on a dedicated
machine. Dedicated time is not always available and measurements of time are
often peculiar to a given installation. The accuracy reported is the normalized
l2 (h) error as requested in [?, eq. 83]. Gflops is an estimate of the number of
send correspondence to bbd@ornl.gov or Mathematical Sciences Section, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831-8083
floating point operations performed per second during the integration. Hardware performance monitors are the preferred measurement method.
Algorithm
Spectral
Spectral
Spectral
Spectral
TIG Model
A-L
Icosohedral PIC
Resol.
.
T42
T42
T213
T213
2562
72 44
10242
Machine
Accuracy
Gflops
Y-MP
Y-MP
Y-MP
Y-MP
C90
C90
Y-MP
1
6
1
6
1
1
1
1010
1010
1010
1010
2.5 104
2.5 104
7.5 104
0.162
0.567
0.215
1.210
0.087
0.351
0.103
Execution
Time (sec)
3.5
1.0
690.0
130.0
20.4
3.9
26.0
Notes
Table 1: Best CPU time - Accuracy for Test Case 2

Table 2 compares the parallel performance of methods. P is the number
of processors used in the computation and SP is the parallel speed up with P
processors over a single processor time. If TP denotes the execution time for P
S
processors, then SP = TTP1 . The parallel efficiency is given by EP = Pp .
Algorithm
Resol.
Machine
SP
EP
Spectral
Spectral
Spectral
Spectral
Spectral
T42
T213
T21
T42
T85
Y-MP
Y-MP
iPSC/860
iPSC/860
iPSC/860
6
6
64
128
128
3.5
5.4
5.6
18.4
49.6
0.58
0.90
0.08
0.14
0.39
Execution
Time (sec)
1.0
130.0
1.37
3.92
16.9
Notes
2
2
6
6
6
Table 2: Parallel Performance on Test Case 2
Notes
1. Results of STSWM [?]. Solution exactly representable in spectral expansion so accuracy not representative. The Y-MP results were calculated in
64bit arithmetic.
2. Rudy Jacobs results of multitasked STSWM reported at the Third CHAMMP
Workshop on Numerical Solution of PDEs in Spherical Geometry.
3. TIG is the twisted icosahedral grid method described in [?]. Execution
time estimated from 600 sec timesteps at 0.0284 sec/step on test case 5.
4. Arakawa-Lamb as described in [?]. Execution time estimated from 600
sec timesteps at 0.0284 sec/step on test case 5.
2
1
2
2
2
3
4
5
5. The PIC method is applied on an icosahedral grid of 10242 points. 90

timesteps were taken for the 5 day simulation. Results presented by John
Baumgardner at the Third CHAMMP Workshop on Numerical Solution
of PDEs in Spherical Geometry.
6. The Intel iPSC/860 results are 32bit arithmetic with accuracy O(105 ).
The T21 case required 90 timesteps, T42 180, and T85 360, for the
five day integration.
Literature
Seven test cases were proposed in [?]. These cases collect several tests common
in the literature but particularly follow work in [?]. A code to solve the shallow
water equations using the spectral transform method (STSWM) is described
in [?]. High resolution test case solutions using the spectral code STSWM are
given in [?].
The report [?] compares solutions using an icosahedral grid twisted to maintain grid symmetry between hemispheres. Parallel algorithms for the spectral
transform are discussed in [?, ?].

Numerical Performance Results From The Shallow Water Equation Test Suite

Uploaded by

Copyright:

Available Formats

Numerical Performance Results From The Shallow Water Equation Test Suite

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Numerical Performance Results From The Shallow Water Equation Test Suite

Uploaded by

Copyright:

Available Formats

Numerical Performance Results from the Shallow

Water Equation Test Suite

Updated: 3 May 1993

The DOE Computer Hardware, Advanced Mathematics and Model Physics

Comparison of Algorithms and Computing Platforms

Table 1: Best CPU time - Accuracy for Test Case 2

Table 2: Parallel Performance on Test Case 2

5. The PIC method is applied on an icosahedral grid of 10242 points. 90

You might also like