November 01, 2002
Scientific Computing Using POOMA
Jeffery D. Oldham
A very effective, MP-capable library for mathematical computation. If it works for Los Alamos, it'll work for you.
November 2002/Scientific Computing Using POOMA/Algorithm 1
// Let N(t,i,j) represent the neutron density at time t at grid position (i,j).
// Initial conditions
for all grid positions (i,j)
N(0,i,j) = 0 if not at center of grid
N(0,i,j) = 10^7 at center of grid
for each time step t,
for all grid positions (i,j) (except the outside boundary)
N(t+1,i,j) = stencil(i,j)
// Boundary conditions
N(t+1,boundary) = 0
Inspect the simulation if desired.
endfor
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
Next Page