BT_SERIAL
NAS Parallel Benchmark BT, Serial Version, FORTRAN90

BT_SERIAL is a FORTRAN90 program which is a serial implementation of the NAS Parallel Benchmark known as BT ("block tridiagonal solver for 3D PDE").

The NAS Parallel Benchmarks are "paper and pencil" specifications; that is, the proposers have laid down what the programs must do algorithmically, but only very general rules about computer language implementation have been given. Thus, a user or vendor has enormous freedom in implementing the benchmark program in a way that is optimal for a given architecture or parallel programming scheme.

As a starting point for implementors, the proposers have supplied certain initial versions of the benchmark programs.

This directory contains material derived from one such initial benchmark code. Persons interested in the original source of the benchmark codes should NOT work from this material, but should go to the NAS Parallel Benchmark website directly!

The benchmark programs are available in SERIAL, MPI and OPEN_MP versions, and in a sequence of "classes" of increasing size: S, W, A, B, C, D and E.

Statistics for the classes S, A, and B, are available, as run on one processor of a Cray Y-MP, with a FORTRAN77 version of the program.

Class	Size	Memory (Mw)	Time (sec)	MFLOPS
S	12x12x12	0.3	7.2	34
A	64x64x64	24.0	923.0	192
B	102x102x102	96.0	3554.0	203
C	162x162x162	?	?	?

Statistics for classes S, W, A, B, and C are available, as run on an Apple PowerPC G5, using this SERIAL FORTRAN90 version of the program:

Class	Size	Memory (Mw)	Time (sec)	MFLOPS
S	12x12x12	0.3	1	230
W	24x24x24	?	33	233
A	64x64x64	24.0	735	228
B	102x102x102	96.0	3192	219
C	162x162x162	?	?	?

This directory contains the serial version of the benchmark, along with include files that allow it to run in some of the smaller sizes.

The web site for the NAS Parallel Benchmarks is http://www.nas.nasa.gov/Resources/Software/npb.html.

Languages:

BT_SERIAL is available in a FORTRAN90 version.

Related Data and Programs:

CG_SERIAL, a FORTRAN program which is a serial version of the NAS Parallel Benchmark CG (conjugate gradient solver).

EP_SERIAL, a FORTRAN90 program which runs the serial version of the NAS Parallel Benchmark known as "EP" (embarassingly parallel).

FIRE_SERIAL, a FORTRAN90 program which simulates a forest fire over a rectangular array of trees, starting at a single random location. It is intended as a starting point for the development of a parallel version.

IS_SERIAL, a C program which is a serial version of the NAS Parallel Benchmark IS (integer sort).

LINPACK_BENCH, a FORTRAN77 program which measures the time taken by LINPACK to solve a particular linear system.

MATMUL, a FORTRAN90 program which is an interactive matrix multiplication benchmark.

MDBNCH, a FORTRAN77 program which is a benchmark molecular simulation calculation.

MXM_SERIAL, a FORTRAN90 program which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.

POISSON_SERIAL, a FORTRAN90 program which computes an approximate solution to the Poisson equation in a rectangle, and is intended as the starting point for the creation of a parallel version.

PRIME_SERIAL, a FORTRAN90 program which counts the number of primes between 1 and N, intended as a starting point for the creation of a parallel version.

SEARCH_SERIAL, a FORTRAN90 program which searches the integers from A to B for a value J such that F(J) = C. this version of the program is intended as a starting point for a parallel approach.

SUBSET_SUM_SERIAL, a FORTRAN90 program which seeks solutions of the subset sum problem, in which it is desired to find a subset of a set of integers which has a given sum; this version of the program is intended as a starting point for a parallel approach.

Reference:

1. David Bailey, Eric Barszcz, John Barton, D Browning, Robert Carter, Leonardo Dagum, Rod Fatoohi, Samuel Fineberg, Paul Frederickson, Thomas Lasinski, Robert Schreiber, Horst Simon, V Venkatakrishnan, Sisira Weeratunga,
   The NAS Parallel Benchmarks,
   RNR Technical Report RNR-94-007, March 1994.

Source Code:

• bt_serial.f90, the source code.
• header.h, an include file for BT_SERIAL.
• work_lhs.h, an include file for BT_SERIAL.
• wtime.c, a small routine used for timing.
• wtime.h, an include file for WTIME.

Examples and tests:

BT_SERIAL_S runs the benchmark in serial mode with problem size S.

• nbparams_S.h, the include file, which should be renamed nbparams.h before compilation.
• bt_serial_S_output.txt, the output file.

BT_SERIAL_W runs the benchmark in serial mode with problem size W.

• nbparams_W.h, the include file, which should be renamed nbparams.h before compilation.
• bt_serial_W_output.txt, the output file.

BT_SERIAL_A runs the benchmark in serial mode with problem size A.

• nbparams_A.h, the include file, which should be renamed nbparams.h before compilation.
• bt_serial_A_output.txt, the output file.

BT_SERIAL_B runs the benchmark in serial mode with problem size B.

• nbparams_B.h, the include file, which should be renamed nbparams.h before compilation.
• bt_serial_B_output.txt, the output file.

BT_SERIAL_C runs the benchmark in serial mode with problem size C.

• nbparams_C.h, the include file, which should be renamed nbparams.h before compilation.
• bt_serial_C_output.txt, the output file.

List of Routines:

• MAIN is the main program for the BT_SERIAL benchmark.
• ADD updates the vector U.
• ADI
• ERROR_NORM computes the norm of the error.
• RHS_NORM
• EXACT_RHS computes the right hand side based on exact solution.
• EXACT_SOLUTION returns the exact solution at point (xi, eta, zeta).
• INITIALIZE initializes the field variable u.
• LHSINIT
• PRINT_RESULTS
• COMPUTE_RHS
• SET_CONSTANTS
• MATVEC_SUB subtracts bvec=bvec - ablock*avec
• MATMUL_SUB subtracts a(i,j,k) X b(i,j,k) from c(i,j,k)
• BINVCRHS
• BINVRHS
• TIMER_CLEAR
• TIMER_START
• TIMER_STOP
• TIMER_READ
• ELAPSED_TIME
• VERIFY is the verification routine
• XSOLVE performs line solves in the X direction.
• Y_SOLVE performs line solves in the Y direction.
• ZSOLVE performs line solves in the Z direction.

You can go up one level to the FORTRAN90 source codes.