Basic HTML version of Foils prepared 19 September 98

Foil 77 Jacobi iteration: data structures

From MPI Message Passing Interface Computational Science for Simulations -- Fall Semester 1998. by Geoffrey C. Fox, Nancy McCracken
CALL MPI_Type_vector( nyblock,
1, nxlocal+1, MPI_REAL, rowtype, ierr)
CALL MPI_Type_commit(
rowtype, ierr );
dx = 1.0/nx; dy = 1.0/ny; err = tol * 1.0 e6
DO j = 0, nylocal
  • DO i = 0, nxlocal
    • f(i,j) = -2*(dx*(i+coords(2)*nxblock))**2
    • + 2*dx*(i+coords(2)*nxblock)
    • u(i,j) = 0.0
    • unew(i,j) = 0.0
  • END DO
END DO

© 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 Apr 11 1999