Basic HTML version of Foils prepared 15 March 1996

Foil 5 Matrix Formalism for NAS Benchmarks

From Spatial Differencing and ADI Solution of the NAS Benchmarks CPSP713 Case studies in Computational Science -- Spring Semester 1996. by Geoffrey C. Fox
1 We can write the right hand side of equation (3.26) as the matrix vector product
2 The vector Uvec has Nvec = 5 Nx Nh Nz components and RHS is a Nvec by Nvec matrix. The forcing term Hvec is also a vector with Nvec components
  • The components of vectors and rows and columns of matrices are labelled by i,j,k and CFD index index running from 1 to 5
3 Formally we can write full equations:
4 If we had used the explicit Euler method, then equations would have read the simple:
  • which only needs sparse matrix vector multiplication as its computational kernel.
in Table To:

Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Sun Feb 22 1998